向NCBI提交基因簇的时候需要提供sqn格式的文件，之前我在文章《NCBI上传基因簇之tbl2asn的使用》中介绍过如何使用tbl2asn生成sqn文件，遗憾的是tbl2asn官方已经不再提供软件下载了，提供的新工具为table2asn，本

使用场景假设有一个fasta格式的序列文件SRR9620252.faa，我们想要提取其中的一些序列到一个新的文件中，我们拥有这些序列的id （假设这些id存放在文件ids.txt中）。常规操作的话，可以复制id，在fasta文件中打开搜索，

向NCBI提交基因簇的时候需要提供sqn格式的文件，这个文件需要通过tbl2asn生成。 文件准备tbl2asn依赖三个文件来生成sqn文件： 文件1：fasta格式的基因组序列文件 Header处的中括号部分可以不写。 1234567

有时候需要在FASTA格式的文件中搜索短的保守序列，这个时候采用查找法比使用blast等序列比对更加人性化。但是要注意避坑，即FASTA文档中的序列一般是被打断为许多行的，如果要查找的目标序列恰好在断行处，是没有办法直接揪出它的，所以在查找

二话不说，上代码，需要安装BioPerl 123456789101112131415161718192021#!/usr/bin/perluse strict;use warnings;use Bio::SeqIO;# Author: Li

当一个基因组含有多个Contigs/Scaffolds的时候，在GenBank文件中也会产生多个LOCUS。某些软件会将GenBank文件作为输入，但仅支持一个GenBank文件中只包含一条Contig/Scaffold