Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #46
    No, that is the correct interpretation of the flags. 8 is also the flag for "mate is unmapped", which in theory, the first read should have, since it has the 4 flag for being unmapped itself. I'd guess that the 4 flag in the second read is an error, but it's hard to explain. One weirdness that I know bwa has when aligning is it concatenates separate reference sequences, so if a read crosses over two references, it's have a mapping position, and be flagged as unmapped (PICARD yells at you when it sees this), but that wouldn't seem to explain what you have there.

    And no, unmapped paired reads might not have an * in the rname. If a read is unmapped, and its mate is mapped, it's supposed to have the rname and position of its mate. That's so they will sort together. It's not a bug, it's a feature, really. The flag is supposed to show this, maybe the ISIZE field will too.

    Comment


    • #47
      thank you for the reply.

      In the example of my previous msg the reads have different coordinates and so essentially they are mapped. I had also tried doing a blat of the sequence on the mouse whole genome (mm9) and both of them show multiple hits.

      Just for clarity, below are the best hits for both the reads. the latter read has many hits.
      ---------------------------------------------------------------------------------------------------
      browser details read1_flag99 50 1 50 50 100.0% 9 - 115259995 115260044 50
      browser details read1_flag99 50 1 50 50 100.0% 9 - 95326385 95326434 50
      browser details read1_flag99 50 1 50 50 100.0% X + 166650106 166650155 50
      browser details read1_flag99 50 1 50 50 100.0% 2 + 181747899 181747948 50
      browser details read1_flag99 50 1 50 50 100.0% 2 + 57481930 57481979 50

      browser details read2_flag151 50 1 50 50 100.0% 9 - 95326353 95326402 50
      browser details read2_flag151 50 1 50 50 100.0% 6 - 50997472 50997521 50
      browser details read2_flag151 50 1 50 50 100.0% 6 - 4873828 4873877 50
      browser details read2_flag151 50 1 50 50 100.0% 6 - 4874144 4874193 50



      Also I would expect that if a read is mapped & the mate pair is unmapped then the flag for the read could not be 99, it may be the values 75 / 91 / 139 / 155. but then I might be wrong.

      Just to clarify "mate reverse strand" means that the mate pair is mapped in reverse strand?

      Comment


      • #48
        Flag = *?

        Hi, I am mapping 454 reads with the default bwa bwasw, and I'm mostly getting good results, but a couple have me stumped...

        I understand that the FLAG of 4, with aligned chromosomal coordinates indicates an alignment possibly spanning across multiple chromosomes.

        But, I'm also getting a few alignments like this:

        GVVGYHW01CSFZ8 * chrIII 7855436 199 139M * 0 0 GATCTGCCAAGCAACAGGCTAAGTCTCGCAATCAAGCCGTCAGAAAGTTCGCAGTCAAGGCTTAAGTGGCTTGTATTCATTGTTATCCATTCATGGACATTGTTCTGGTTCAATTTCAATGAAAATTGTTGCTTATGTT @BEEA<<552275<----/8//66B<<@BIIIIHHHIIIIIIIIIIIIHCDCHHEIEEEECCDBIEEEIIII???IIIIIIIIIHHHDDB131ABB@AA=EABACB?455;IIIIIIIIIIIIIIIIIIIIIIIHDDDI AS:i:139 XS:i:0 XF:i:3 XE:i:4 XN:i:0

        As you can see, the FLAG for this alignment is "*"... is this similar to the FLAG=4 situation described above, or something else entirely?

        Thanks...

        Comment


        • #49
          Hi,

          Is it 16 refer to read map at reverse strand?
          Can I know which bitwise flag refer to those read that map in pair?

          Thanks.

          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, 03-27-2024, 06:37 PM
          0 responses
          12 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-27-2024, 06:07 PM
          0 responses
          11 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-22-2024, 10:03 AM
          0 responses
          53 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-21-2024, 07:32 AM
          0 responses
          68 views
          0 likes
          Last Post seqadmin  
          Working...
          X