美文网首页
Shell脚本统计DNA序列中碱基数目

Shell脚本统计DNA序列中碱基数目

作者: dming1024 | 来源:发表于2019-04-22 18:16 被阅读0次

为了不编写perl脚本统计DNA中各碱基的数目,我采用了shell脚本直接对DNA序列进行统计。我的DNA序列如下两行:
ACTGATACGTACGATCGATA
ACAGEDDADTCGATCGASTCTA。
代码很简单就一行:

 cat dna.txt | paste - - |sed  's/\t//'|grep -o .|sort|uniq -c|sort -rn
     13 A
      9 T
      8 C
      7 G
      3 D
      1 S
      1 E

相关文章

网友评论

      本文标题:Shell脚本统计DNA序列中碱基数目

      本文链接:https://www.haomeiwen.com/subject/afpugqtx.html