Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • mbk0asis
    Member
    • Dec 2011
    • 41

    #1

    BISMARK paired-end result sam file contains on read1 info

    Hello, I ran bismark with paired-end reads.
    According to the manual, the result SAM file should had contained both read 1 and 2 sequences, mapping qual, and etc. But I didn't get them both.
    I used TrmGalore for adapter removal with the paired-end option to make sure both mates were paired.

    Command I used :
    Code:
    bismark --non_directional --bowtie2 -p 8 -N 1 -L 10 ../mm10/ -1 MEF_R1_val_1.fq -2 MEF_R2_val_2.fq
    Here's some examples from my result,

    Code:
    HISEQ:284:HAK1MADXX:1:1101:1451:2148_1:N:0:CGATGT       99      chr7    48002502        23      101M    =       48002628        227     AAATGTGGTATATTTATATAATGGAGTATTATTTAGTTATTAAAAAGAGTGTATTTATGAAATTTTTAGTTAAATGGATGGATTTGGAGGGTATTATTTTG   C@CFFFEFFHHHHJJIJIJIJIJJJJJHJJJJJJJIJIJIIIIIJHIEGBDEGHIJGJJJIGJJJJJJIJJIHIIJJHHGGGHIJJHAEFF5@@D>CEEEA   NM:i:19 MD:Z:10C2C2C1C9C2C1C2C11A15C0C3C0C11C0C7C2C2C0C2        XM:Z:..........h..h..h.h.........h..h.x..h...........................hh...hh...........hx.......h..h..hx..      XR:Z:CT XG:Z:CT
    HISEQ:284:HAK1MADXX:1:1101:1451:2148_1:N:0:CGATGT       147     chr7    48002628        23      101M    =       48002502        -227    ATTTATATAATATGTATTTATTGATAAGTGTTTATTAGTTTAAAATTTAGGATATTTAAGATATAAGATATAATTTATTAAATATATGAAATTTAAGAAGA   FDCHHHHHHHIIGGIGDGFGGIHDIHJJJHIJJJJJIIIIJJHIIEJJJJJJIGGIJJIJJJJJIHFIJJJJIHHC,IJIGJJIJJJJHHHHDDFFFFCCC   NM:i:23 MD:Z:1C1C1C1C8C1C1C9C7C0C0C5C7C0C0C13C5G0C4C1C6C1C5A1   XM:Z:.h.h.h.h........h.h.x.........h.......hhh.....h.......hhh.............h......h....h.h......h.h.......      XR:Z:GA XG:Z:CT
    Thank you.
  • frozenlyse
    Senior Member
    • Sep 2008
    • 135

    #2
    The first row is read 1, and the second row is read 2 - you can tell from the SAM flag (second column)

    99 == 0x01 + 0x02 + 0x20 + 0x40 which translates to "template having multiple segments in sequencing" + "each segment properly aligned according to the aligner" + "SEQ of the next segment in the template being reversed" + "the first segment in the template"

    147 == 0x01 + 0x02 + 0x10 ("SEQ being reverse complemented") + 0x80 ("the last segment in the template").

    Comment

    • mbk0asis
      Member
      • Dec 2011
      • 41

      #3
      Thank you for your answer!
      I think 'frozenlyse' is right.
      I tried SAM flag interpreter (?) and confirmed they were indeed properly paired reads.

      '99' means;
      read paired
      read mapped in proper pair
      mate reverse strand
      first in pair

      '147' means;
      read paired
      read mapped in proper pair
      read reverse strand
      second in pair

      By the way, what does the 'read reverse strand' means?

      Comment

      • dpryan
        Devon Ryan
        • Jul 2011
        • 3478

        #4
        "Read reverse strand" means that the reverse complement of the read aligned, rather than the read itself. In other words, it mapped to the "-" strand rather than the "+" strand (leaving out that things aren't that simple in bisulfite sequencing...the fragment that gave rise to those example reads originated from the original + strand...if you want to know more about this then have a look at the trim_galore manual where Felix went over things quite nicely).

        Comment

        Latest Articles

        Collapse

        • 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
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-22-2026, 11:13 AM
        0 responses
        60 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 10:48 AM
        0 responses
        42 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        75 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        74 views
        0 reactions
        Last Post SEQadmin2  
        Working...