Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SAM Format - SEQ field '='

    Hi,

    I was just looking through the SAM format spec again and I've come across something that confuses me.

    The text for the 'SEQ' field is:
    "query SEQuence; “=” for a match to the reference; n/N/. for ambiguity; cases are not maintained"

    In practice, I've never seen a "=" in the SEQ field; is this supposed to be an optional way of using the format, i.e. using a combination of "=" and "N" instead of the actual query sequence?

    Thanks,
    Bio

  • #2
    You're right, that is somehow confusing.
    By chance, I just read on the samtools manual page that samtools fillmd using option -e can "convert the read base to = if it is identical to the aligned reference base. Indel caller does not support the = bases at the moment."

    I guess reporting the actual query sequence (along with the CIGAR string) makes much more sense in terms of computation because you don't need to look up the bases in the reference. And by using N, you'd lose all SNP information.

    Comment


    • #3
      Thanks epigen, I guess it's an optional way of using the format then.

      Comment


      • #4
        Using = in the SEQ field is useful for compressing the data, see:
        In some respects the SAM/BAM specification is quite loose, in that there is more than one way to represent a given piece of information. We...


        (I know it is an old post, but the new similar thread feature of the forum brought it to my attention)

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Recent Advances in Sequencing Analysis Tools
          by seqadmin


          The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
          05-06-2024, 07:48 AM
        • seqadmin
          Essential Discoveries and Tools in Epitranscriptomics
          by seqadmin




          The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
          04-22-2024, 07:01 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 05-14-2024, 07:03 AM
        0 responses
        19 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-10-2024, 06:35 AM
        0 responses
        44 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-09-2024, 02:46 PM
        0 responses
        54 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-07-2024, 06:57 AM
        0 responses
        42 views
        0 likes
        Last Post seqadmin  
        Working...
        X