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
          Latest Developments in Precision Medicine
          by seqadmin



          Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

          Somatic Genomics
          “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
          05-24-2024, 01:16 PM
        • 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

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 05-24-2024, 07:15 AM
        0 responses
        16 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 10:28 AM
        0 responses
        18 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-23-2024, 07:35 AM
        0 responses
        22 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 05-22-2024, 02:06 PM
        0 responses
        11 views
        0 likes
        Last Post seqadmin  
        Working...
        X