5、提取snp indel 位点

标签：perl   print   ane   提取   length   lan   final   gre   list   

le final.snp.list | perl -lane ‘{$a+=1;print "$a\t$F[0]\t$F[1]\t$F[1]"}‘ | less >snp_site

le final.indel.vcf |grep -v ‘^#‘ | less -S|perl -lane ‘{$a+=1;$b=$F[1]+length($F[3]);print "$a\t$F[0]\t$F[1]\t$b"}‘ | less -S >indel_site

5、提取snp indel 位点

标签：perl   print   ane   提取   length   lan   final   gre   list   

原文地址：http://www.cnblogs.com/renping/p/7356586.html