View Single Post
Old 04-25-2017, 12:12 PM   #4
Brian Bushnell
Super Moderator
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707

In a duplicated/repeat region, there is no way to know whether the read is mapped to the correct location; and thus, it's impossible to determine whether the variant call is in the correct place. Also, repeats are hard to assemble so the reference is less likely to be correct there. Typically you can just discard ambiguously-mapped reads rather than discarding variants, though.

There's a large block of Chr Y that is similar to Chr X; as an exception, variants in this region are often called on the X chromosome, since if they were discarded, much of the X/Y chromosome would be unusuable for variant detection.
Brian Bushnell is offline   Reply With Quote