...介绍一种方法。
用到软件是bedtools，具体方法如下：
Usage: bedtools getfasta [OPTIONS] -fi <fasta> -bed <bed/gff/vcf>Options: -fi Input FASTA file -bed BED/GFF/VCF file of ranges to extract from -fi -name Use the name ...
...="(chicken,((mouse,rat),(chimp,human)));";
@names=($tree=~/([a-zA-Z0-9_]+)/g);
open OUT ,">newtree.nwk" or die "can't open $name.nwk\n";
for my $name(@names){
if ( $name!~/^\d+$/){
$tree=~s/$name/$name#1/;
}
}
print OUT "$tree\n";
close(OUT);
生物信息入门...