Quote:
Originally Posted by michy
Hello,
I've used MAQ to call SNPs in Illumina data, it worked fine apart from the read depth. When I looked at the alignment using IGV the read depth is often alot more. Im sure theres a reason, can anyone tell me why?.
Thank you
|
Look at the default filtering: mapping quality etc. These may be reducing coverage (filtering reads) during SNP calling. This is a good thing, since poorly mapped reads or high-error reads may introduce artifacts.