Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SAM file after Bowtie is messed up

    Hi,



    I have chip-seq data from E. coli (51 bp). I mapped my fastq file to my reference genome (custom build) using Bowtie in Galaxy. In the SAM output file, some rows have the sequence in the quality score columns, and the quality scores in the OPT column. Some rows are fine.



    Anyone would know what is causing it and how to fix that ?



    Thanks

  • #2
    Are you sure that's actually the case? It's incredibly more likely that you're just miscounting the columns.

    Comment


    • #3
      See step 14

      Comment


      • #4
        Have you set the file format to "fastqsanger" for your original data files (I can't tell from the history you shared). Here is how you would do it: https://wiki.galaxyproject.org/Suppo...ognize_dataset Then you should not have to groom your data. If this is recent data correct so it should already be in sanger fastq format.

        It appears that part of illumina fastq header (1:N:0:18) is missing from the reads that appear to have an alignment (at least that is what it looks like in the web page).

        Comment


        • #5
          Hi,

          thank you for looking at my data.

          I have tried without grooming, just changing data type (my reads are illumina 1.9 encoding) and I have the exact same result.

          The illumina fastq header (1:N:0:18) is present for all reads in the fastq file.

          I have tried galaxy GVL instance and galaxy main. Same results.

          I don't have this problem when I use BWA mapping. But it's better to use Bowtie for E. coli reads since BWA looks for intron so better used for eukaryotes is that right ?

          Comment


          • #6
            I don't have this problem when I use BWA mapping. But it's better to use Bowtie for E. coli reads since BWA looks for intron so better used for eukaryotes is that right ?
            No, BWA, like Bowtie does not take into account the introns.
            Only splice-junction aware aligners, like TopHat and STAR do, in which case you have to provide them with the genome annotation indicating the location of the junctions.
            TopHat actually delegates the alignment to Bowtie1 or 2, and only handles the splicing.

            In the link to the Galaxy instance that you posted, you are using a version of Bowtie that dates back to 2010, version 0.12.7. It's not clear from your post if you've already tried this, but the first troubleshooting step I would take would be to upgrade to a more modern version of Bowtie. There is a long list of bugs that have been fixed in Bowtie since 2010.

            Comment


            • #7
              Thanks for that ! That's really helpful.

              I didn't check which version of Bowtie I was using thinking that the Galaxy main instance would display the most up to date version. I will have a look at that.

              Thanks a lot.

              Comment

              Latest Articles

              Collapse

              • 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
              • seqadmin
                The Impact of AI in Genomic Medicine
                by seqadmin



                Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                02-26-2024, 02:07 PM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 03-14-2024, 06:13 AM
              0 responses
              33 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-08-2024, 08:03 AM
              0 responses
              72 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-07-2024, 08:13 AM
              0 responses
              81 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-06-2024, 09:51 AM
              0 responses
              68 views
              0 likes
              Last Post seqadmin  
              Working...
              X