Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • samtools XS:A optional field

    Does anyone know what XS:A field in samtools output mean by any chance?

    XS:A (A is the SAM attribute type) can take the values of "+" or "-". I thought "+" or "-" might be related to how the read is mapped to the ref (forward or reverse complement), but it turned out not the case, because forward or reverse complement read can take both values.

    Thanks!

  • #2
    Originally posted by charlie_sequencing View Post
    Does anyone know what XS:A field in samtools output mean by any chance?

    XS:A (A is the SAM attribute type) can take the values of "+" or "-". I thought "+" or "-" might be related to how the read is mapped to the ref (forward or reverse complement), but it turned out not the case, because forward or reverse complement read can take both values.

    Thanks!
    The optional tags are specific to the aligner that generated them. What aligner did you use? Take a look at the manual from that aligner.

    Comment


    • #3
      There are two places in a SAM file that, on the surface, *seem* to provide equivalent information -- the second column (0 vs 16, or something) and the XS:A field.

      I have a vague memory that the first tells you about mapping relative to the reference sequence. The second (XS:A) tells you something about the direction/strand of the *read itself*.

      I wish I could find the post I read that told me that. I first came across the distinction while trying to understand the input fields required for TopHat. You might give a try searching for "TopHat" and "XS field" together. Please post back if you figure it out.

      Best of luck,

      DG

      Comment


      • #4
        Thanks so much for your reply, nilshomer and dagarfield. I found the Cufflinks manual http://cufflinks.cbcb.umd.edu/manual.html where XS field was specified. When a read is mapped to the ref across splicing junctions, the boundary of intron (GT-AG, or CT-AC) will tell you the value of XS: "+" for GT-AG and "-" for CT-AC. It tells you "which strand the RNA that produced this read came from".

        This led me to think about another question: is TopHat able to report the strand of the read if it is completely contained within an exon, i.e., without N in the CIGAR string?

        Comment


        • #5
          Early RNA-seq protocols were unstranded ie the RNA may have originated on the plus strand, but because the library was created from double stranded cDNA you obtain reads from both the plus and minus strand - in this case the samtools flag for which strand the read mapped to doesnt tell you which strand the RNA it originated from. However, if the read is across a splice junction you can work it out (as charlie_sequencing posted) hence the tophat XS flag. Reads from an unstranded library which do not cross a splice cannot be assigned a strand.

          However, lots of people are doing stranded RNA-seq these days where the reads only come from one strand which is a lot nicer.

          Comment


          • #6
            Other posts that may be of help

            Please check out these posts for a discussion of directional sequencing and XS:A tags:

            Bridged amplification & clustering followed by sequencing by synthesis. (Genome Analyzer / HiSeq / MiSeq)

            and
            Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

            Comment


            • #7
              Originally posted by frozenlyse View Post
              Early RNA-seq protocols were unstranded ie the RNA may have originated on the plus strand, but because the library was created from double stranded cDNA you obtain reads from both the plus and minus strand - in this case the samtools flag for which strand the read mapped to doesnt tell you which strand the RNA it originated from. However, if the read is across a splice junction you can work it out (as charlie_sequencing posted) hence the tophat XS flag. Reads from an unstranded library which do not cross a splice cannot be assigned a strand.

              However, lots of people are doing stranded RNA-seq these days where the reads only come from one strand which is a lot nicer.
              Hi,frozenlyse, As you said, unstranded library we don't konw reads come from which strand. Now I have a dataset of dUTP and I use --library firststrand .Tophat said,"it is assumed that only the strand generated during first strand synthesis is sequenced.",that means, the reads come from minus strand. But I found a half of reads was taged "XS:A:+"

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Genetic Variation in Immunogenetics and Antibody Diversity
                by seqadmin



                The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
                11-06-2024, 07:24 PM
              • seqadmin
                Choosing Between NGS and qPCR
                by seqadmin



                Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
                10-18-2024, 07:11 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 11-08-2024, 11:09 AM
              0 responses
              129 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-08-2024, 06:13 AM
              0 responses
              98 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 11-01-2024, 06:09 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-30-2024, 05:31 AM
              0 responses
              25 views
              0 likes
              Last Post seqadmin  
              Working...
              X