gff文件从gene的pos位点在基因组里面抽取序列

大佬教程收集整理的这篇文章主要介绍了gff文件从gene的pos位点在基因组里面抽取序列，大佬教程大佬觉得挺不错的，现在分享给大家，也给大家做个参考。

open FASTA,"$ARGV[0]" or die "can not open GFF file,$!";
open GFF,"$ARGV[1]" or die "can not open FASTA file,$!";
 while (<FASTA>)
{
chomp;
if (/>/){$key=$_;$key =~ s/>//g;}
else {$hash{$key}.=$_}
}#把基因序列文件扫描进去hash表
 while (<GFF>)
{
chomp;
@F=split;
next unless $F[2] eq "gene";
$out=substr($hash{$F[0]},$F[3]-1,$F[4]-1);
print ">$F[1]\n$out\n";
}

大佬总结

以上是大佬教程为你收集整理的gff文件从gene的pos位点在基因组里面抽取序列全部内容，希望文章能够帮你解决gff文件从gene的pos位点在基因组里面抽取序列所遇到的程序开发问题。

如果觉得大佬教程网站内容还不错，欢迎将大佬教程推荐给程序员好友。

本图文内容来源于网友网络收集整理提供，作为学习参考使用，版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ：384754419，请注明来意。