This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Weird Duplicated Lines With Gatk

Hi all,

I came across weird lines produced by GATK and I was wondering if anyone ever came across a similar problem.

I have two lines for the same position (but no indel):

chr19    9413330    .    G    T    529.48    .GT:DP:GQ:PL  0/1:111:99:559,0,2802
chr19    9413330    .    G    .    16.10    LowQual    GT:DP:GQ:PL  0/0:106:0.02:23,0,2904

I understand having multiple lines for indels but what do these lines even mean ? I mean the sample is probably heterozygous at that position.

Thanks !

gatk genotyping
• 2,326 views • link
updated by Alex Paciorkowski Alex Paciorkowski • written by Gabriel R. Gabriel R.

1 answer

It looks like your second read at chr19:9413330 is of low quality and there is no call. I would filter that read out of further analyses.

• 0 views • link
written by Alex Paciorkowski Alex Paciorkowski

It's not a read, it's the output of the genotyper.

• 0 views • link
written by Gabriel R. Gabriel R.

Same idea. There's still no call at that site. I would still toss it.

And your sample is heterozygous at position. That's what 0/1 in your first "read" means.

• 0 views • link
written by Alex Paciorkowski Alex Paciorkowski

I understand, the bottom line is that this line should not exist.

• 0 views • link
written by Gabriel R. Gabriel R.

Log in to answer this question.

More posts like this