Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • length of uniquely mapped seq in Bismark

    Dear all

    Bismark reports the number of allignments with a unique best hit, is there a way to know what was the Length of those unique sequences ? for instance I have an RRBS library which after trimming has seq length of 20 - 140 bp, this gave me about 40% unique mapping in Bismark, I would like to know the length of those fragments which mapped uniquely.

    Best regards

  • #2
    Hi Hedi,

    Every alignment that is reported by Bismark are unique hits (ambiguosly aligning reads are discarded, and not-aligning reads - well - don't align at all). So you can simply look at the aligned read length distribution using your program of choice. It would only be 2 clicks in SeqMonk for example. As a general rule though, longer reads tend to align better.

    40% mapping efficiency doesn't sound great to be honest. Which organisma are you working with, and did you perform appropriate adapter/quality trimming etc?

    Comment


    • #3
      thank you Felix

      so i just imported the BAM files generated by Bismark into Seqmonk, without any quantification i can click on read length distribution histogram on data sets. but trimmed files should have 20 - 140 bp length, why still i can see seq with for example 400 bp in BAM files?

      we are working with bull data, and yes we did pay attention to evrything when it comes to trimming and quality assessment. but still this is what we have unfortunatelly.

      Appreciate your halp

      Comment


      • #4
        That sounds like you are working with paired-end data. In that case the Read 1 and Read 2 will be ends of the fragment you are sequencing, and they may well be further apart. In this case it doesn't really matter how long each read is, but it is the combination of the mappability of R1 and R2 that counts.

        Comment


        • #5
          thank you

          we actually used two different kits for RRBS library preparation, I figured out that there is lots of short reads under 70 bp in the first method and using another method the peak of read length was about 130 bp which we would like to have, as you mentioned longer reads tend to align better. but surprisingly mapping was almost the same, that's why i was curious to look at the length of those aligned reads in both libraries.

          Best

          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
          13 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