Python和Perl处理fasta文件对比_2018-09-11

从fasta中抽取需要的序列:
Python方法一:(这种方式最慢)
由于平时工作比较忙,工作中遇到的问题都没有及时记录,现在有空时间了所以整理一下一些常见的问题及技巧供大家参考:

!/usr/bin/python

import sys
def Fasta(inputfile) :
f = open(inputfile,”r”)
fastadic = {}
while True :
line = f.readline().rstrip()
if len(line) == 0 :
break
else :
pass
if line.startswith(“>”) :
name = line.split()[0]
fastadic[name] = “”
else :
fastadic[name]+=line
f.close()
return fastadic

def usage():
print (“1,fasta\n2,list”)

if len(sys.argv) != 3 :
usage()
sys.exit()

print (sys.argv[0])

fastafile = sys.argv[1]
listfile = sys.argv[2]

dic = Fasta(fastafile)
ff=open(listfile,”r”)
while 1:
line1 = ff.readline().rstrip()
if len(line1) == 0 :
break
else :
pass
array = line1.split(“\t”)
if “>”+array[0] in dic.keys() :
print (“>”+array[0]+”\n”+dic[“>”+array[0]])
else :
print (“Error:”+line1)
ff.close()
Python方法二:(次慢)

!/usr/bin/python

import sys
def Fasta(inputfile) :
f = open(inputfile,”r”)
fastadic = {}
while True :
line = f.readline().rstrip()
if len(line) == 0 :
break
else :
pass
if line.startswith(“>”) :
name = line.split()[0]
fastadic[name] = []
else :
fastadic[name].append(line)
f.close()
return fastadic

def usage():
print (“1,fasta\n2,list”)

if len(sys.argv) != 3 :
usage()
sys.exit()

print (sys.argv[0])

fastafile = sys.argv[1]
listfile = sys.argv[2]

dic = Fasta(fastafile)
ff=open(listfile,”r”)
while 1:
line1 = ff.readline().rstrip()
if len(line1) == 0 :
break
else :
pass
array = line1.split(“\t”)
if “>”+array[0] in dic.keys() :
print (“>”+array[0]+”\n”+””.join(dic[“>”+array[0]]))
else :
print (“Error:”+line1)
ff.close()
Perl 方法三(最快)

!/usr/bin/perl -w

die “perl 《Python和Perl处理fasta文件对比_2018-09-11》ARGV[0] ;
while (<I>) {
chomp;
if (/^>/) {
@name = split/\s+/,《Python和Perl处理fasta文件对比_2018-09-11》fasta{《Python和Perl处理fasta文件对比_2018-09-11》fasta{《Python和Perl处理fasta文件对比_2018-09-11》_ ;
}
}
close I;

open T,《Python和Perl处理fasta文件对比_2018-09-11》_ ;
if (exists 《Python和Perl处理fasta文件对比_2018-09-11》tax[0]}) {
print “>”.《Python和Perl处理fasta文件对比_2018-09-11》fasta{“>”.$tax[0]}, “\n” ;
}
}

close T;

    原文作者:凤舞琦天
    原文地址: https://www.jianshu.com/p/009ea69015f7
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞