This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Detecting Human Contamination in Low Coverage NGS sample

Hi,

I have created contaminated human sequence R1, R2 fq.gz files by adding 10% random sequences from another human sample. When I align to the reference genome, the average coverage is around 1.0x.

I have tried to detect the contamination with verifybamid and also GATK GetPileupSummaries/CalcualteContamination but the the contamination is not detected.

Has anyone been successful detecting contamination on low coverage .bam files with coverage around 1.0x?

Thanks.

verifybamid calculatecontamination
• 1,271 views • link
written by igorm • igorm

1 answer

SOLVED: The .bam had to be duplication-marked. When using the .bam produced by Picard's MarkDuplicates with verifybamid, I get the correct contamination estimate.

• 0 views • link
written by igorm • igorm

Log in to answer this question.

More posts like this