Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie - What am I missing here?

    Ok so I am a newbie here, so please bear with me. I am trying out Bowtie. And following the manual and I issue the command as it is in the manual (I just removed the --suppress 1,5,6,7 to make all the columns visible.

    ./bowtie -a -v 2 e_coli -c ATGCATCATGCGCCAT

    The first two lines of the output are as follows...

    0 - gi|110640213|ref|NC_008253.1| 148810 ATGGCGCATGATGCAT IIIIIIIIIIIIIIII 0 10:A>G,13:C>G
    0 - gi|110640213|ref|NC_008253.1| 2852852 ATGGCGCATGATGCAT IIIIIIIIIIIIIIII 0 8:T>A

    If I am reading this right, it found matches on the negative reference strand, with TWO mismatches on the first result (10:A>G,13:C>G). But when I look at the alignments closely,

    Test - ATGCATCATGCGCCAT
    Result 1 - ATGGCGCATGATGCAT

    I see mismatches at position 4,5,6,11,12,13.

    Same thing with result TWO

    Test - ATGCATCATGCGCCAT
    Result 1 - ATGGCGCATGATGCAT

    Bowtie says ONE mismatch in position 8, while I see mismatches in positions 4,5,6,11,12,13

    What AM I missing here? Please apologize for relative ignorance of the question.
    Thanks
    Quantrix

  • #2
    If the match is to the negative strand, you probably need to take the reverse complement of one of your sequences. With that being said, it looks to me like your "Test" and "Result 1" are in fact identical taking the reverse complement, so there would be no mismatches rather than two.

    Comment


    • #3
      Dear Kopi,
      You are right, in fact, I did the reverse complement for BOTH the sequences and they are exactly the same. So does this mean when Bowtie outputs the results, it DOES NOT output the actual result sequences? (And merely outputs the TEST sequence?)

      I do not know. I am more confused than ever. Can someone clarify this for me please?
      Thank you
      Q

      Comment


      • #4
        OK, I should have thought of that. What you are seeing is default Bowtie output, which according to the manual (http://bowtie-bio.sourceforge.net/ma...-bowtie-output) does not show the reference sequence, but the read sequence.

        Comment


        • #5
          Damnit!
          RTFM...R-T-F-M! Now I feel like an idiot.
          Thank you Kopi. But it is still weird though. Why even bother outputting the READ sequence when there is a match? I am sure there is an answer to it. But never mind.
          Q

          Comment


          • #6
            All of the aligners I know output the read sequence (or reverse-complemented if it is mapped to reverse strand). This is helpful for the post-processes, e.g. you can convert the output back to fastq or you can calculate the GC content of the mapped reads.

            Comment


            • #7
              Hi Nguyen,
              Thanks for the reply. It DOES make sense now. I think I got thrown off by the fact there was only ONE read in the test sequence. If there are millions of reads, it makes sense to output the reads and for other purposes such as GC content etc
              Q

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Best Practices for Single-Cell Sequencing Analysis
                by seqadmin



                While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                06-06-2024, 07:15 AM
              • seqadmin
                Latest Developments in Precision Medicine
                by seqadmin



                Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                Somatic Genomics
                “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                05-24-2024, 01:16 PM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Today, 07:24 AM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 08:58 AM
              0 responses
              11 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-12-2024, 02:20 PM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-07-2024, 06:58 AM
              0 responses
              184 views
              0 likes
              Last Post seqadmin  
              Working...
              X