Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Maq Alignment Coordinates

    I have used maq to align a set of illumina reads to a reference assembly. Is there a way to directly get the alignment coordinates for the reference assembly?

    Thanks,
    Anamika

    PS, I had inadvertently posted the above question in a different section.

  • #2
    Hey! I'm not at a computer so i cant verify, but I think it's maq mapview...

    Comment


    • #3
      I believe ECO is right. It's "maq mapview filename.map" to display a binary .map file into a human readable format. I usually end up piping it to less so that it doesn't scroll by at unreadable speeds.

      Still, I have to wonder, ECO, If you're not at a computer, how did you post that? (-;
      The more you know, the more you know you don't know. —Aristotle

      Comment


      • #4
        Originally posted by apfejes View Post
        I believe ECO is right. It's "maq mapview filename.map" to display a binary .map file into a human readable format. I usually end up piping it to less so that it doesn't scroll by at unreadable speeds.

        Still, I have to wonder, ECO, If you're not at a computer, how did you post that? (-;
        iPhone != computer in most senses...

        Comment


        • #5
          Well, there go my theories of you being an AI or having a cybernetic implant...
          The more you know, the more you know you don't know. —Aristotle

          Comment


          • #6
            maq mapview right solution

            Thanks Eco. Once I ran it with the -b option ( the read sequence and the quality are not displayed).

            From this output, I plan to get the coordinates of those regions that did not align (roughly about ~15% of the assembly.)

            Cheers,
            Anamika.

            Comment


            • #7
              Good luck!

              Comment


              • #8
                Anamika,

                It would help if you checked the coverage of your reads on the assembly using maq.

                First assemble the reads into a consensus:

                0) merge all your reads into a single .map file
                1) maq assemble output.cns genome.bfa reads.map &> asm.log
                2) maq cns2win -w 10000 output.cns > windows.csv

                If you plot columns 6 & 7 you get an indication of read depth and coverage of your sequences against the genome.

                Hope this helps you out.

                Comment


                • #9
                  Hi Zee,

                  I use the option you suggest to get the read depth for the SNP data. However, my aim is to get those coordinates of the assembly that did not align with my input reads. I have about ~15% not aligning.

                  Thanks,
                  Anamika

                  Comment


                  • #10
                    maq

                    hello,
                    when using MAQ, after the command :
                    ./maq map 1131.map onecdna.bfa coli.bfa
                    the sreen print:
                    -- maq-0.7.1
                    [ma_load_reads] loading reads...
                    [ma_load_reads] set length of the first read as 35.
                    Segmentation fault


                    what the "Segmentation fault" mean?
                    provide my best services for you

                    Comment


                    • #11
                      Originally posted by xzk421 View Post
                      hello,
                      when using MAQ, after the command :
                      ./maq map 1131.map onecdna.bfa coli.bfa
                      the sreen print:
                      -- maq-0.7.1
                      [ma_load_reads] loading reads...
                      [ma_load_reads] set length of the first read as 35.
                      Segmentation fault


                      what the "Segmentation fault" mean?

                      It means that you have not defined you data path

                      Comment


                      • #12
                        Originally posted by zee View Post
                        Anamika,

                        It would help if you checked the coverage of your reads on the assembly using maq.

                        First assemble the reads into a consensus:

                        0) merge all your reads into a single .map file
                        1) maq assemble output.cns genome.bfa reads.map &> asm.log
                        2) maq cns2win -w 10000 output.cns > windows.csv

                        If you plot columns 6 & 7 you get an indication of read depth and coverage of your sequences against the genome.

                        Hope this helps you out.
                        The read depth in column 6 is easy enough to understand but what does the value in column 7 mean? My best guess is that it gives the average mapping quality for that region. Is that even close?

                        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, 03-27-2024, 06:37 PM
                        0 responses
                        12 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-27-2024, 06:07 PM
                        0 responses
                        11 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-22-2024, 10:03 AM
                        0 responses
                        53 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, 03-21-2024, 07:32 AM
                        0 responses
                        69 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X