美文网首页
2018-06-24(生信第二题)

2018-06-24(生信第二题)

作者: 天秤座的机器狗 | 来源:发表于2018-06-24 20:03 被阅读35次

    题目:假设有如下的fasta格式文件,想把多行序列合并成一行,怎么做?

    (多行)

    看到这个题目首先想到的是建一个空的字符串,然后把除了含“>”之外的序列合并到一起,然后打印出来,合并的方法想到的是字符串直接相加的方法。但是个有问题,就是对于多个序列就没办法操作了。所以先对单个序列试试,将序列写入到test_tmp.fa文件中。cat一下发现个问题,就是出现了空行,如下

    去除空行的方法,想到最近学习的linux命令中的sed命令,可以完成对空行的去除,如下

    对这个序列操作

    但是对于多个序列而言,这样的方法肯定不行,所以想到的是用字典的方式,将序列名那行和对应的序列构建成一个字典,然后将字典打印出来,就可以完成,如下

    这样就完成了fasta序列中多个序列多行合并成一行的问题了

    相关文章

      网友评论

          本文标题:2018-06-24(生信第二题)

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