How to detect False positive variant in repeat region ?
1
0
Entering edit mode
6.6 years ago
istdasklar • 0

My variant caller returns variants like the following one in repeat region. ( PGM / Amplicon ) .
How can I know if there are true of false ? I mean caused by sequencing error

IGV IGV2

bam vcf calling • 2.4k views
ADD COMMENT
2
Entering edit mode
6.6 years ago

Hey, didn't you have another account? - your profile photo looks familiar.

The only way to confirm a variant is through the use of an ancillary method, like Sanger sequencing. NGS always struggles to correctly call variants in repeat regions whose length approaches the average read length that you're using. Why? - in part, it is due to the issue of mis-alignment in these regions. Even prior to in silico alignment, homopolymers like AAAAAAA, GGGGGGGG, etc., can be difficult to faithfully sequenced during the sequence run itself.

To guard against errors in repeat regions, you can do some basic QC thresholds:

  • Prior to alignment, trim bases at read ends whose average base qualities fall below 30
  • Prior to alignment, eliminate short reads
  • Prior to variant calling, eliminate reads with MAPQ<40, 50, or 60
  • Require that variants are called at minimum of 18 read depth
  • Require that variants have 'high' genotype qualities (at least 30)
  • Only look at variants that pass a threshold for strand bias (given by PV4 tag)
ADD COMMENT
0
Entering edit mode

To guard against errors in repeat regions, you can do some basic QC thresholds:

  • Prior to alignment, trim bases at read ends whose average base qualities fall below 30
  • Prior to alignment, eliminate short reads
  • Prior to variant calling, eliminate reads with MAPQ<40, 50, or 60
  • Require that variants are called at minimum of 18 read depth
  • Require that variants have 'high' genotype qualities (at least 30)
  • Only look at variants that pass a threshold for strand bias (given by PV4 tag)
  • Take a variant caller that do denovo assembly like freebayes or GATK's HaplotypeCaller
ADD REPLY
0
Entering edit mode

Thanks for your reply. I will check all the things !

ps : Yes, I logged to a wrong account to post my question .. sorry !

ADD REPLY

Login before adding your answer.

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