Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • problems with Bioscope output

    Hello everyone,
    I mapped my Solid data with Bioscope. In the csfasta files are around 142 million reads. This value is shown in the pairing statistics file as well. But in the bam file with the mapped and unmapped reads are together just 119 million reads. I checked that with the help of samtools. What happened to the other reads? Does Bioscope filter any reads? I couldn't find anything about that.

    Furthermore the pairing statistics show more mapped pairs than I could find in the bam file. Could anybody explain that, please?

    I would be very happy about some help or comments. Thanks

  • #2
    you do mean that you counted the mapped and unmapped files? and the numbers do not add up?
    which version of Bioscope are you using? can you post the ini file?
    http://kevin-gattaca.blogspot.com/

    Comment


    • #3
      Hello Robby,

      You mentioned that the input csfasta file contains roughly 142 million reads. As you know not all reads are mappable so you should expect some losses in the bam file. To be very specific, bam file will only have an entry for an unmapped read IF the sister read (the other half of the pair) is mapped. If both reads from a pair are not mappable, then the pair will not show up in the bam file. This should explain the discrepancy you are observing.

      Regarding the differences between pairing statistics and bam file, could you post the stats from both approaches? It is difficult to troubleshoot without something more concrete.

      Thanks

      Comment


      • #4
        Thanks for your help and for parts of the solution. As hanktu mentioned all pairs of which both reads are mappable appear in the bam file with the mapped reads as two reads. If both reads from a pair are not mappable, then the 2 reads appear in the bam file with the unmapped reads. But if one read in a pair is mappable and the other not, then one entry appears in the bam file with the mappable reads and no entry in the file with the unmappable reads. So the unmapped read of the pair doesn't appear anywhere. The problem was, that Bioscope considers pairs and I counted reads.

        But nevertheless 80-160 reads (or 40-80 pairs) are missing. In the Bioscope statistics they are counted as mapped. Do you have any ideas for that as well?

        Comment


        • #5
          What version of BioScope?

          Hi Robby,

          What version of BioScope are you running? Can you copy and paste some lines from the pairing stats file as well as the numbers you get from the bam file?

          Thanks

          Comment


          • #6
            Could this be caused by redundantly mapped reads? If you're creating the bam file from 'primary' as the config setting (I think this is the default), then there will be reads that map to many locations that are 'mappable' but will not appear in your bam files unless they can be rescued by a primary/uniquely mapping mate.

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