在perl处理fasta文件是，需要每次读入一个完整的 序列头+序列，需要利用到
local $/ = ">"
将perl默认的行分隔符切换为">"。
且从文件中第一次读取的内容为">"号，可以在<IN>进行去除，IN为定义的输入句柄。
在处理数据时：
s/>//;
my ($name, $seq) = splt/\n/,$_,2;
$seq =~ s/\n//g;
即可将一个序列的名保存到$name变量，将序列信息保留到$seq变量。
perl处理fasta文件
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 从fasta中抽取需要的序列：Python方法一：（这种方式最慢）由于平时工作比较忙，工作中遇到的问题都没有及时记...
- Perl 命令行实战1 - fasta文件的相关操作 注：本文会不断更新...... 之前的5篇简书小文已经说明了...
- 移步github 问题描述： 对宏基因组（例如肠道微生物）进行全基因组shotgun测序，对其中已知的高丰度菌进行...
