Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • sampie
    Member
    • Jun 2014
    • 10

    Bowtie2 output with paired end reads

    Hi,

    I am trying to use bowtie2 to align paired end Illumina read files.

    I am using following command:

    ./bowtie2 -a --threads 8 --n-ceil 42 --very-sensitive-local --quiet --no-sq --no-hd --np 0 --ma 1 --mp 0,0 --rdg 1,1 --rfg 1,1 --score-min L,0,0.95 -x rdp/bt2_index/rdp.db -q -1 sR1.fastq -2 sR2.fastq > bt2out.txt

    The index is made from RDP II database and currently sR1.fastq and sR2.fastq contain only one read:

    sR1.fastq:
    @1_M01363:14:000000000-A5K8J:1:1101:8617:1211 1:N:0:16
    TACGGAGGATCCGAGCGTTATCCGGA...

    sR2.fastq:
    @2_M01363:14:000000000-A5K8J:1:1101:8617:1211 2:N:0:16
    TTCACCGTTGCCGGCGTACTCCCCAGGTGGA

    The bt2out.txt file shows mostly results with sR1.fastq read header.

    1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
    1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
    1_M01363:14:000000000-A5K8J:1:1101:8617:1211....
    .
    .
    .

    I am not sure how to read bowtie2 results. How can I obtain information how each matching read aligned as a pair?

    Thanks
  • westerman
    Rick Westerman
    • Jun 2008
    • 1104

    #2
    Bowtie2 should be outputting BAM/SAM files. Read up on that format. If you have further questions after reading about SAM files -- especially the bitwise flag field -- then ask a question.

    Comment

    • GenoMax
      Senior Member
      • Feb 2008
      • 7142

      #3
      It appears that sampie omitted "-S filename_to_hold_alignments.sam" option on the command line to write the alignments to a file.

      Comment

      • westerman
        Rick Westerman
        • Jun 2008
        • 1104

        #4
        Originally posted by GenoMax View Post
        It appears that sampie omitted "-S filename_to_hold_alignments.sam" option on the command line to write the alignments to a file.
        While good practice the '-S' is optional. As per the docs:


        -S <hit>
        File to write SAM alignments to. By default, alignments are written to the "standard out" or "stdout" filehandle (i.e. the console).
        So the stdout file of bt2out.txt should be a SAM file.

        Comment

        • jwfoley
          Senior Member
          • Jun 2009
          • 183

          #5
          As long as we're talking best practices, I would definitely pipe the output through samtools and save a BAM rather than a SAM file. There's no reason a SAM file should ever exist on a hard drive.

          Comment

          • sampie
            Member
            • Jun 2014
            • 10

            #6
            Thanks for replies.

            I took a look at SAM format specification. It was more clear to me than bowtie2 manual.

            I am now testing bits 0x1, 0x2 and 0x8 from flags and if they are on I am determining that current match is R2 and previous match is R1 (of the same pair).

            Hopefully it works correctly now.

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Pathogen Surveillance with Advanced Genomic Tools
              by seqadmin




              The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
              03-24-2025, 11:48 AM
            • seqadmin
              New Genomics Tools and Methods Shared at AGBT 2025
              by seqadmin


              This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

              The Headliner
              The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
              03-03-2025, 01:39 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-20-2025, 05:03 AM
            0 responses
            42 views
            0 reactions
            Last Post seqadmin  
            Started by seqadmin, 03-19-2025, 07:27 AM
            0 responses
            51 views
            0 reactions
            Last Post seqadmin  
            Started by seqadmin, 03-18-2025, 12:50 PM
            0 responses
            38 views
            0 reactions
            Last Post seqadmin  
            Started by seqadmin, 03-03-2025, 01:15 PM
            0 responses
            194 views
            0 reactions
            Last Post seqadmin  
            Working...