foreach my $i (@region){
$all++;
@Rgene = split "\t",$i;
system "samtools depth -r $Rgene[1] $bam > $all.txt";
my $file = "$all.txt";
# print $file,"\n";
open(IN2,"$file");
while(<IN2>){
chomp;
$sum++;
my @line = split "\t";
转载于:https://www.cnblogs.com/blueicely/p/3385945.html