Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Mapping of Paired end sequenced Amlicon?

    Hi guys,

    I am puzzled how to do mapping of my amplicons PE seq reads. Library was prepared in such a way that PE seq has been done from both sides of amplicons. Now I am slightly confused, How to map these reads over the reference amplicons using bowtie?

  • #2
    Is there a particular reason you need to use bowtie?

    Comment


    • #3
      No there is no perticular reason. I could use other tool as well.

      My question is related with the mapping. If it is single end, I can map it easily over the reference amplicon but how to map the Paired reads over the amplicons?

      Comment


      • #4
        NGS aligners are able to account for the paired-end reads. You would have no problems as long as you properly indicate input files with read pairs being mapped. What size amplicons do you have?

        Comment


        • #5
          I have reads which are around 200 bp long.

          Comment


          • #6
            See BBMap and the example for how to use paired end reads for mapping.

            Comment


            • #7
              Merge overlapping PE reads

              In continuation of the question, I have another idea of merging the overlapping reads of the PE sequenced fastq data.

              I found a very interesting page related to this: http://thegenomefactory.blogspot.com...aired-end.html

              They have suggested few tools. eg. fastq-join, FLASH etc. Have anybody tried any of such tools? So that I can have a review.


              I would appreciate if any one could suggest me for the PE sequenced data,
              1) Should I directly go for Mapping using PE option (eg. bowtie1) OR
              2) merge the PE reads R1 and R2 together to make a super-long SE read, with extra confidence of the middle bases from consensus of the overlapping sections of R1 and R2?


              Thanks all in advance.

              Comment


              • #8
                If you know that your reads are overlapping then you can merge them into a longer representation before mapping. If you find that the reads are not merging well (because insert size was bigger than you expected) then you should proceed with mapping them as paired-end.

                BBMap suite has a tool that does the merging called BBMerge. FLASH is another well regarded tool for merging reads.
                Last edited by GenoMax; 01-19-2015, 06:34 AM.

                Comment


                • #9
                  Thanks for your response.

                  I have ~350bp long DNAs fragments and after PE sequencing the length of reads are 200 (FastQC). I would expect an overlap of ~50 bases. And that is the reason, I am thinking of merging not directly mapping with PE parameter.
                  Last edited by jeni; 01-19-2015, 08:03 AM.

                  Comment


                  • #10
                    Both approaches are valid, just bear in mind that not all PE reads will merge even if they overlap, particularly if are low-complexity in the overlapping region or if the reads have high error rates in their tails. Merged reads are typically easier to work with, but unmerged reads may have lower bias if a large percentage do not merge successfully.

                    It's also possible to map paired reads, merge them based on mapping coordinates, then remap the merged reads.

                    Is your goal variant calling, quantification, or something else?

                    Comment


                    • #11
                      Basically quantification.

                      Comment


                      • #12
                        For quantification, I recommend mapping the unmerged reads as pairs, to prevent bias against low-complexity sequence.

                        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
                        10 views
                        0 likes
                        Last Post seqadmin  
                        Started by seqadmin, Yesterday, 06:07 PM
                        0 responses
                        9 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
                        67 views
                        0 likes
                        Last Post seqadmin  
                        Working...
                        X