compute normal-tumor coverage ratio from exome BAMs
1
0
Entering edit mode
8.1 years ago
tayebwajb ▴ 110

Could someone please suggest a quick way to compute the data ratio of uniquely mapped reads in
the normal to uniquely mapped reads in the tumor, as required by Varscan in the command below? I have over 50 exome BAMs.

(normal_unique_mapped_reads/tumor_unique_mapped_reads).

java -jar VarScan.jar copynumber normal-tumor.mpileup
output.basename -min-coverage 10 --data-ratio
[data_ratio] --min-segment-size 20
--max-segment-size 100

varscan gatk exome samtools bedtools • 1.7k views
ADD COMMENT
0
Entering edit mode
8.0 years ago
Chirag Nepal ★ 2.3k

Hi tayebwajb,

Did you manage how to extract the uniquely mapping reads from normal-tumor.mpileup or from bam files.

thanks !

ADD COMMENT

Login before adding your answer.

Traffic: 938 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6