Python根据ID从fasta文件提取序列
根据ID从fasta文件提取序列,为什么结果只有最后一个ID的结果,之前的被覆盖了?以前用这个提取过序列,没有出现问题,这次不知道为什么,只有一个ID的结果
#!/usr/bin/python3
import sys
args=sys.argv
fr=open(args[1], 'r')
fw=open('./out.fasta', 'w')
dict={}
for line in fr:
if line.startswith('>'):
name1=line.split(">")[1]
name=name1.split(":")[0]
dict[name]=''
else:
dict[name]+=line.replace('\n','')
fr.close()
ls=open(args[2], 'r')
for lin in ls:
for name in dict.keys():
if name ==lin:
fw.write(">"+ name)
fw.write('\n')
fw.write(dict[name])
fw.write('\n')
fw.close()
...全文
相关推荐
FASTA序列提取方法[代码]
本文介绍了两种从FASTA序列文件中提取指定序列的方法。第一种方法是使用seqkit命令行工具,通过seqkit grep命令从FASTA文件中提取指定ID的序列,并详细说明了安装和使用步骤。第二种方法是使用Python脚本,通过读取输入FASTA文件,检查目标ID列表,并将匹配的序列写入输出文件。两种方法均适用于处理包含大量序列的FASTA文件,能够高效地提取所需序列。文章还提供了详细的代码示例和使用说明,帮助读者快速掌握这些技术。
python学习——通过命令行参数根据fasta文件中染色体id提取染色体序列
提取fasta文件genome_test.fa中第14号染色体的序列,其内容如下: >chr1 ATATATATAT >chr2 ATATATATATCGCGCGCGCG >chr3 ATATATATATCGCGCGCGCGATATATATAT >chr4 ATATATATATCGCGCGCGCGATATATATATCGCGCGCGCG >chr...
extract_fasta_序列提取_
给定物种名,批量分子矩阵中提取给定的物种名称对应的序列,并生成新矩阵用于系统发育分析。
cbbl:从NCBI核苷酸数据库下载所有cbbl基因并以fasta格式输出的脚本
cbbl 从NCBI核苷酸数据库下载所有cbbl基因的脚本,并以fasta格式输出。 只需更改搜索词和输出即可将其通用化。
perl and python.zip_DNA_freedom178_序列提取_提取DNA
从DNA序列文件里按照序列号提取文件,可根据实际情况更改


其他技术讨论专区
