What is the fold 80 base penalty of Mitochondrial Data?
0
0
Entering edit mode
23 months ago
NikhilP ▴ 20

Hello, I have sequenced whole genome mitochondrial samples. I calculated the fold 80 base penalty using CollectHsMetrics(GATK). But I am getting above 100, which I think is incorrect. The code which I used for calculating is the below code

gatk CollectHsMetrics -I /mnt/NGS2/MitoRes/Batch_20/10255168_26/BAM_Files/10255168_26_MtMkDupSorted.bam \
-O /mnt/NGS2/MitoRes/Batch_20/10255168_26/BAM_Files/10255168_26_MtMkDupSorted_metrics.txt \
-R /mnt/NGS2/Database/Homo_sapiens_assembly38.chrM.fasta \
-BI /mnt/NGS2/Database/Twist_MitoPanel_chrM_all_hg38_target.interval_list \
-TI /mnt/NGS2/Database/Twist_MitoPanel_chrM_all_hg38_target.interval_list

For this particular sample, I am getting Fold-80 as 40. So then I used the parameter -covMax 2000 and the score went down to 4.

So I wanted to know the correct value for -covMax while calculating fold 80 for mitochondrial Data. And what is the optimum result we should get for mitochondrial data?

Thank you.

NGS mitochondiral • 497 views
ADD COMMENT

Login before adding your answer.

Traffic: 2478 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