学习
实践
活动
专区
工具
TVP
写文章

算法(三)列举所有k-mer组合

关键词:k-mer; recursive; trick; 什么是k-mer? 比如,“ATGC”的所有1-mer是:’A’, ‘T’, ‘G’, ‘C’。共4^1=4种组合。 而“ATGC”的所有2-mer是 “AA”, “AT”, “AG”, “AC” “TA”, “TT”, “TG”,“TC” “GA”, “GT”, “GG”,“GC” “CA”, “CT”, “CG”, 那么如何打印出所有的k-mer组合呢?如果是2-mer,我们可以用两个for循环来列出所有组合,如果是3-mer,可以用三个for循环。但是如果是10-mer呢?岂不是要10个for循环? 从而会生成不同的k-mer。 细细研读这段代码后,可以发现这种方法只适用于字符串长度为2的指数的情况。 最后 我们再给出列举“ABCDEFGH”的所有k-mer组合的代码: ? 如果任何问题欢迎交流!

51420
  • 广告
    关闭

    2023新春采购节

    领8888元新春采购礼包,抢爆款2核2G云服务器95元/年起,个人开发者加享折上折

  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    k-mer分析:你的基因组有没有被污染?

    其中Kmergenie常用于预测de novo组装中最优组装k-mer大小,根据reads分割k-mers并绘制k-mer深度分布曲线。 现在所有的测序reads均产生测序k-mers,由于测序深度较高,k-mers出现的频次也即k-mer深度较大,在去除错误率影响的前提下,可以认为其中完全不同的k-mer数目即是genomic k-mer 数目,使完全不同的k-mer的数目最大的k值就是最佳k-mer大小。 hash的二进制文件结果给出,需要统计出k-mer总数(Total),特异的k-mer数目(Distinct),只出现过一次的k-mer数(Unique),频数最高的k-mer数目(Max_count) 等信息,以stats命令来运行,如下所示: jellyfish stats -o mer_counts_stats.txt mer_counts.jf 对k-mer的计数结果有个直观的认识,则需要统计出现了

    68830

    关于k-mer与基因组(组装)的那些事

    在这个过程中,我们经常会遇到k-mer这个名词,然而这个抽象的名词是什么意思呢?它又有什么用呢?接下来,就随着小编一起去探究这k-mer背后的含义吧! k-mer是什么? 通过将reads切割成以k为单位的k-mer,由于测序错误具有随机性,这些由于测序错误生成的k-mer绝大多数都是原测序物种中不存在的k-mer,因此都只出现了1次,要是将这些k-mer去掉,那么就会较大的可能除去测序错误 我们用k-mer做什么? 在了解了k-mer是什么以及通过去掉低频率的k-mer能够使得组装结果更加准确以后,k-mer就没有别的用途了吗?当然不是! 下图是在k-mer=15、17、19时分别作的k-mer深度分布曲线。 说了那么多使用k-mer分析的优点,好像忘了一个重要的点:k-mer怎么好像只有奇数呢? 是的,k-mer只能是奇数,就是为了防止通过k-mer组装时,正反链混淆。

    7K63

    基于深度学习微表情识别的调查

    早期的MER方法主要是基于传统的外观和几何特征。最近,随着深度学习(DL)在各个领域的成功,神经网络在MER中受到越来越多的关注。 由于上述ME字符,基于DL的MER变得具有挑战性。对于数据,人们提出了各种DL方法来解决ME问题并提高MER性能。 在这项调查中,我们对深度微表情识别(MER)进行了全面的回顾,包括数据集、深度MER管道和最具影响力的方法的基准标记。这项调查为该领域定义了一个新的分类法,包含了基于DL的MER的所有方面。 此外,我们总结了剩余的挑战和设计强大的深度MER系统的潜在方向。据我们所知,这是第一次对深度MER方法的调查,这个调查可以作为未来MER研究的参考点。 as a reference point for future MER research.

    35700

    miRNA 靶向预测软件targetscan

    01 Targetscan靶向预测思想 TargetScan 基于序列互补原则,找到比对到靶 3'UTR 的保守性 8 mer、7 mer 或 6 mer 位点(seed match 序列),进一步根据热力学稳定性筛选得到 seed 序列配对主要考虑三种类型:7 mer-1a(miRNA 的第 2-7nt 与靶基因互补配对, 而且 UTR 上与 miRNA 1nt 互补配对的位置是 A);7 mer-m8 (miRNA 2 -8nt 与靶基因完全配对);8 mer (miRNA 2-8nt 与靶基因完 全配对,而且 UTR 上与miRNA 1nt 互补配对的位置是 A)。 主要包括如下几部分: Site Type 8 mer > 7 mer-m8 > 7 mer-1a; 3' pairing contribution:除了与 miRNA seed 区域配对,与 miRNA12 其中标题各列的含义如下: Gene ID:基于 ID Species ID:物种 ID Mirbase ID:miRbase 中 miRNA 的 ID Site Type:配对类型(8mer、7 mer-m8

    4.2K20

    AAAI | 联合建模医学命名实体识别和标准化的神经多任务学习框架

    使层次化任务(MER和MAN)在保持任务间相互支持的同时,转化为并行多任务模式成为可能。 在本文中,作者将MER和MEN看做两个并行的任务。MER和MEN采用相同的输入但具有不同的输出。 因此,可以将MEN视为具有与MER相同的输入的序列标记任务。 和MEN任务的反馈进行的改进,证明了反馈策略都能提高这两个任务的性能,其中对MER效果尤为明显;最后一部分展示了将Bi-LSTM、多任务学习(MTL)和MER,MEN的反馈策略联合起来的模型效果,最终证明了文章提出的模型效果是最优的 最后作者还分析了普通模型和本文模型的边界不一致误差,实验结果表明MTL可以显著的缓解MER和MEN边界不一致问题,从而提高模型性能。 为了更先进、更智能地利用两者之间的关系,文章提出了一种新的具有两种显式反馈策略的深层神经多任务学习框架来联合建模MER和MEN。

    46660

    扫码关注腾讯云开发者

    领取腾讯云代金券