SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   MAQ and SNP calling (http://seqanswers.com/forums/showthread.php?t=3333)

michy 12-09-2009 08:41 AM

MAQ and SNP calling
 
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

nilshomer 12-09-2009 04:48 PM

Quote:

Originally Posted by michy (Post 11307)
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.

bioenvisage 12-10-2009 02:21 AM

can any one tell me the details regarding the alignment output of maq

BertieWooster 12-11-2009 09:23 AM

Quote:

Originally Posted by bioenvisage (Post 11348)
can any one tell me the details regarding the alignment output of maq



Well, the maq map output itself is binary. You can run the pileup command to give you the read bases and the quality at each position. Check the maq documentation and that should elaborate further what the columns in the output are.

BertieW


All times are GMT -8. The time now is 11:46 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.