Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • duhaimj
    Junior Member
    • Sep 2009
    • 2

    #1

    Tophat2: bam file has no quality sequences for a lot of reads

    • With the version 2 of tophat I get an '*' as quality sequences in the bam file for more that 70% of the reads.

      Ex:
      Read 1
      GCCCTCGCGGTGCTCTCCGGGTCTGTGCTGAGGAGAACGCAACTCCGCCG DDDDDDDFFFGHHHJJIJIJIIIJJIGGHIGJJJJJJHGHHHFFFFFCC@
      Read 2
      GGCGTGTGTTGC *
      Read 3
      CCGGCCCTCGCTCCAGCAGCTGGACCCCTACCTGCCGTCTGCTGC *

      Read 4
      CAACTGCTGGCCTGTGCCAGGGTGCAAGCTGAGCACTGGAGTGGAGTTTTC C@CFFFFFHGHHHIGHJJIIJICEHIGIIIIGHJGIJJJBGDFGGHHIIIH

      The bam file for this same sample with tophat 1.4 is ok.

      Has someone encountered this problem? THank you
  • Jon_Keats
    Senior Member
    • Mar 2010
    • 279

    #2
    Why are your read lengths different. Looks like the two that are the same size have proper quality strings

    Comment

    • pzumbo
      Member
      • Mar 2009
      • 11

      #3
      I encountered the same problem.
      Problem was alleviated by switching back to 1.4.1.

      Comment

      • duhaimj
        Junior Member
        • Sep 2009
        • 2

        #4
        My reads were filtered for quality. But the same problems happens without filtering. And version 1.4 corrects the problem as mentionned by pzumbo.

        Comment

        • caddymob
          Member
          • Apr 2009
          • 36

          #5
          I'm having the same issue with tophat2. I just posted over here: http://seqanswers.com/forums/showpos...66&postcount=4

          Comment

          • QTLdum
            Junior Member
            • Feb 2011
            • 5

            #6
            I had the same problem when running tophat2 with bowtie2. However, I did not see this issue after running tophat2 with the --bowtie1 option. Thus, it might be problem with bowtie2.

            Comment

            Latest Articles

            Collapse

            • SEQadmin2
              16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by SEQadmin2




              NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
              Today, 06:46 AM
            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              10-14-2026, 09:55 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 10-28-2026, 09:57 AM
            0 responses
            22 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-28-2026, 09:51 AM
            0 responses
            14 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-22-2026, 11:13 AM
            0 responses
            67 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-22-2026, 10:48 AM
            0 responses
            47 views
            0 reactions
            Last Post SEQadmin2  
            Working...