我發現解析fasta文件的代碼。我需要計算有多少A,T,G等是每個序列中,例如: >gi|7290019|gb|AAF45486.1| (AE003417) EG:BACR37P7.1 gene product [Drosophila melanogaster]
MRMRGRRLLPIIL
在這個序列中那裏有: M - 2
R - 4
G - 1
L - 3
I - 2
P - 1
我寫了這個代碼 import sys
file = open(sys.argv[1], 'r')
string = ''
for line in file:
if line.startswith(">"):
pass
else:
string = string + line.strip()
#print (list(string))
w = i
在串聯或串在Convert1.pl使用未初始化值$Xentr4的第6行 我得到以下錯誤: Error opening no such file or directory at Convert1.pl line 6" #!/usr/bin/perl -w
# This script takes a user specified interleaved fasta input file $ARGV[