因为fasta文件格式通常为

>name1

ATGATAGTGTCTGTAGCTGACTGT

AGTGCTGTAGATAGCTAGCTAGTC

>name2

AGTCGATCGTAGTAGCTAGCTAGC

AGTCGATGCTAGCTAGCTACGAAA

需要识别每次以'>'开头的行来区分每段编码达到处理单条序列的目的,因此需要使用正则表达式来识别每次的名字行

1.规规矩矩的写法

import re
for line in open("aa.fasta"):
    line = line.strip('\n')
    regex = re.compile('^>') 
    m = re.match(regex,line)
    if m is not None:
        print line

输出为所有开头为'>'的行

2.可以把regex那行省略

import re
for line in open("aa.fasta"):
    line = line.strip('\n')
    m = re.match('^>',line)
    if m is not None:
        print line

3.还可以把m省略

import re
for line in open("aa.fasta"):
    line = line.strip('\n')
    if re.match('^>',line) is not None:
        print line

以上三种输出结果相同,都是输出以'>'开头的名字那行。

 

相关文章:

  • 2022-12-23
  • 2021-12-18
  • 2021-08-19
  • 2021-05-21
  • 2022-01-16
  • 2021-11-27
  • 2022-12-23
  • 2021-11-26
猜你喜欢
  • 2022-02-12
  • 2022-12-23
  • 2021-12-13
  • 2022-12-23
  • 2022-12-23
  • 2022-12-23
  • 2021-09-05
相关资源
相似解决方案