Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • sam flag 97 and 145

    Hello,

    After alignment, I got reads in the sam format with flag 97 and 145.
    flag 97:
    the read is paired in sequencing;
    strand of the mate reverse;
    the read is the first read in a pair;

    flag 145:
    the read is paired in sequencing;
    strand of the query;
    the read is the second read in a pair;

    These reads are not mapped in a proper pair but both query and mate are mapped. how to explain this? could you give me an example?

    Thanks
    Yanju

  • #2
    They might be mapping to different contigs? What are the full SAM lines for these two reads?

    Comment


    • #3
      The reads are looking like the following. I only have one reference sequence.

      ENST000004134650_104_988_445_609_200# 97 ENSG00000141510 13261 37 36M = 14182 957 CGGTCAACCGTTTTGTAGAACAACTCCCGTCCCCTC 22222
      2222222222222222222222222222222 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36

      ENST000004134650_104_988_445_609_200# 145 ENSG00000141510 14182 37 36M = 13261 -957 GCGGCCACATCCTCGACGACCACGTCCCCGGTGCCC 22222
      2222222222222222222222222222222 XT:A:U NM:i:0 SM:i:37 AM:i:37 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:36

      Comment


      • #4
        anybody help?

        Comment


        • #5
          First read, 97 == 1 + 32 + 64 == 0x1 + 0x20 + 0x40

          0x0001 the read is paired in sequencing, no matter whether it is mapped in a pair
          0x0020 strand of the mate (0 for forward; 1 for reverse strand)
          0x0040 the read is the first read in a pair

          Second read, 145 == 1 + 16 + 128 == 0x1 + 0x10 + 0x80

          0x0001 the read is paired in sequencing, no matter whether it is mapped in a pair
          0x0010 strand of the query (0 for forward; 1 for reverse strand)
          0x0080 the read is the second read in a pair

          So first read (flag 97) is on the forward strand (since it does not have 0x10 set), second read (flag 145) is on the reverse strand (since it does have 0x10 set).
          Last edited by maubp; 05-14-2010, 01:31 PM. Reason: corrected forward/reverse as discussed in next two posts

          Comment


          • #6
            I think it is the other way round:
            the first read (97) is on the forward strand, since the strand of the mate is 1 (reverse);
            the second read is on the reverse strand.

            and in both cases (97,145)
            0x0002 the read is mapped in a proper pair =0
            0x0004 the query sequence itself is unmapped=0 (0 means it is mapped ???)
            0x0008 the mate is unmapped=0 (0 means it is mapped ???)
            This means both of them are mapped? but not proper?
            Last edited by hollandorange; 05-14-2010, 01:18 PM.

            Comment


            • #7
              Originally posted by hollandorange View Post
              I think it is the other way round:
              the first read (97) is on the forward strand, since the strand of the mate is 1 (reverse);
              the second read is on the reverse strand.
              Yes - I had it right the first time, but flipped in back in a hurry before leaving the computer

              Comment


              • #8
                I had a detail look at the samtool manu. I found, for these reads, the inferred insert size is very large like 957. but in my setting, the insertsize_mean is 200 and std is 20.

                is it the reason that they are mapped, but their mapped position is too far away and then samtool report not proper mapped?

                Comment


                • #9
                  Originally posted by hollandorange View Post
                  I had a detail look at the samtool manu. I found, for these reads, the inferred insert size is very large like 957. but in my setting, the insertsize_mean is 200 and std is 20.

                  is it the reason that they are mapped, but their mapped position is too far away and then samtool report not proper mapped?
                  The flags are set by the aligner, not samtools, so checkout the aligner's documentation. Check out the following to help you digest individual tags:


                  Also, you can use the "-X" option in "samtools view" to convert the flag field to a string format.

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Strategies for Sequencing Challenging Samples
                    by seqadmin


                    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                    03-22-2024, 06:39 AM
                  • 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

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Yesterday, 06:37 PM
                  0 responses
                  10 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 06:07 PM
                  0 responses
                  9 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-22-2024, 10:03 AM
                  0 responses
                  51 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-21-2024, 07:32 AM
                  0 responses
                  67 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X