temp1

1.从bam文件中提取位点的测序深度

samtools depth -b intervals.bed 1.bam 

 

intervals.bed文件包括需要提取的位点

chr1 39717581

chr1 39898741

chr1 40928207

chr1 91405097

chr1 99225626

chr1 150812033

chr1 155448184

chr1 160098527

 

2.提取所有位点的ref_depth和var_depth

 

bam-readcount -b 15 -f hg19.fa 1.bam -l intervals.bed