Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • MAQ SNP calling

    > Hi,
    >
    > From my understanding MAQ is tested for 20-50x coverage regions for SNP
    > calling. However, we routinely have projects where depth of coverage
    > from Solexa reads is ~5000x.
    >
    > I see quite a few 'SNPs' that are not called by MAQ.
    > These are the duplicate samples, where MAQ does NOT call a SNP
    >
    > Pos Depth %a %c %g %t
    > 64 6632 0.11 92.01 0.69 7.19
    > 64 6687 0.07 90.88 0.40 8.64
    >
    > While for these duplicates it makes a SNP call!
    > 64 8594 0.05 91.91 0.54 7.49
    > 64 7889 0.13 91.15 0.49 8.23
    >
    > Any suggestions on what threshold/filter to use to make consistent
    > calling with high depth of coverage regions?
    --
    bioinfosm

  • #2
    SNP calling in high depth/coverage solexa data

    Any other options to use apart from MAQ in calling SNPs on such data?
    --
    bioinfosm

    Comment


    • #3
      Sort of on-topic....

      I'm not sure what the answer is to your question about using MAQ, but I have had the pleasure of doing SNP calls on a MAQ .map file, for which I wrote my own SNP caller. Doing your own SNP calls is actually very easy, but leaves you with the same questions about filters (though, no consistency issues, that I'm aware of.).
      The more you know, the more you know you don't know. —Aristotle

      Comment


      • #4
        Yes, MAQ SNP calling is amazingly good! I just wish it scaled to the high-depth data I have. The other alternative which I will try today is take a random subset of my data and artificially get manageable depth
        --
        bioinfosm

        Comment


        • #5
          How Maq works for small indels?

          Comment


          • #6
            Originally posted by Rao View Post
            How Maq works for small indels?
            from the PDF manual/published paper: maq for indels needs paired end reads. maps read 1, then does gapped alignment of read 2. maq indelpe will list the indels called from the read2 gapped alignment.

            try novocraft software if you want more indels from single end or paired end reads.

            Comment


            • #7
              Originally posted by dvh View Post
              from the PDF manual/published paper: maq for indels needs paired end reads. maps read 1, then does gapped alignment of read 2. maq indelpe will list the indels called from the read2 gapped alignment.

              try novocraft software if you want more indels from single end or paired end reads.
              Thanks dvh,
              How about SOAP

              Comment


              • #8
                I used SOAP to detect upto 4bp indels in the solexa data. I usually use it on the reads that MAQ is not able to align to reference.
                --
                bioinfosm

                Comment


                • #9
                  Another case of MAQ SNP calling that I would like to discuss here:

                  cns.final.snp
                  PKD1 37187 T C 255 255 1.12 63 62
                  PKD1 37189 T C 255 255 1.56 63 62

                  These are so close to each other, and look very likely to be False Positives. I was of the idea that MAQ throws such away, but not so..
                  --
                  bioinfosm

                  Comment


                  • #10
                    Missed SNP call; obvious from pileup output

                    I see this obvious SNP from the pileup file with 382 reads mapping the position of which 99.74% say its a C, as opposed to the reference T.

                    Still, the position is not reported in cns.snp file! As I understand, the maq cns2snp command on consensus.cns does not do any filtering and reports all SNPs. This ought to have been reported..

                    Anything I am missing here?
                    --
                    bioinfosm

                    Comment


                    • #11
                      MAQ SNP calling doubts

                      Originally posted by bioinfosm View Post
                      Another case of MAQ SNP calling that I would like to discuss here:

                      cns.final.snp
                      PKD1 37187 T C 255 255 1.12 63 62
                      PKD1 37189 T C 255 255 1.56 63 62

                      These are so close to each other, and look very likely to be False Positives. I was of the idea that MAQ throws such away, but not so..
                      Hi Bioinfoism,

                      This is probably 4 years late, but have you resolve the SNP calling issue with MAQ? I am actually facing the same scenario here where neighbouring SNPs were called. Any ideas on what is causing this? Also, since neighbouring SNPs were called, can the alignment be trusted?

                      Rookie question. But thanks in advance

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        Strategies for Sequencing Challenging Samples
                        by seqadmin


                        Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                        03-22-2024, 06:39 AM
                      • seqadmin
                        Techniques and Challenges in Conservation Genomics
                        by seqadmin



                        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                        Avian Conservation
                        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                        03-08-2024, 10:41 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, Yesterday, 06:37 PM
                      0 responses
                      8 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, Yesterday, 06:07 PM
                      0 responses
                      8 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-22-2024, 10:03 AM
                      0 responses
                      49 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 03-21-2024, 07:32 AM
                      0 responses
                      66 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X