Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Paired read names / SAM qname format

    Hi everyone,
    I am trying to figure out what possible paired read name formats are in use and good to support for open source alignment programs. The most common paired read name format i encounter is something like the following,
    >READ_A/1
    AAAAACCCCCCC
    >READ_A/2
    TTTTTGGGGGAA
    Is this the only format used by different sequencing platforms? Are there other formats such as,
    >READ_B:1
    AAAAACCCCCCC
    >READ_B:2
    TTTTTGGGGGAA
    Also with respect to SAM format, does anyone know what the proper query name for these reads should be?
    Most aligners that I have tried seem to agree that for paired reads that end in "/1" or "/2", these last two characters are truncated to get the SAM query name.
    For example the first read pair would be reported as,
    READ_A
    But the second read pair would have two query names,
    READ_B:1 and READ_B:2 in the SAM file
    I was thinking that maybe the safest approach would be to find the largest common prefix and truncate there to get the pair read name. But I am not sure if this will always work, for example, if the read names are,
    >READ_C_Z12411/1
    AAAAACCCCCCC
    >READ_C_Z12516/2
    TTTTTGGGGGAA
    The pair read names based on largest common prefix would be,
    READ_C_Z12
    Where as most aligners would report two separate names as,
    READ_C_Z12411 and READ_Z12516
    Thanks for the help!

    Misko

  • #2
    As you've observed, Sam QNAME isn't what traditional assemblers and viewers refer to as a read name. Instead it's really the template name - in your example qname is READ_A or READ_B. I think the automatic detection of /1 and /2 to trim down a read-name to a qname is there simply because /1 and /2 was the MAQ convention. Ideally no one should rely on this convention now.

    A common prefix is prone to error too. When dealing with fastq there's a huge variety of options:

    - 1 fastq line per template, with an implicit assumption we chop in half to get both reads (the original illumina format IIRC)

    - 1 fastq file with fwd and rev reads alternating, so they always come in pairs.

    - 2 fastq files with each end stored in its own file.

    - 1 fastq file for single-ended data. This obviously causes confusion with the first two conventions.

    Every aligner seems to want something marginally different, and fastq frankly is hopeless as a format for embedding such meta data so it's tricky for tools to work out which read layout you use without either explicit command line options or trying to "be clever". In my experience that nearly always ends up shooting you in the foot sooner or later.

    As for SAM. I'd like to see the true read suffix (or entire name if inappropriate) having a standard auxillary key:value tag to go from template names back to the original read names. Adding that and the ability to have more than 2 reads per qname and we gain the ability to use SAM for mixed assemblies with capillary finishing reads in it.

    James

    Comment


    • #3
      Thanks for the answer, this helps clarify

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Quality Control Essentials for Next-Generation Sequencing Workflows
        by seqadmin




        Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

        Nucleic Acid Quality Control
        Preparing for NGS starts with isolating the...
        02-10-2025, 01:58 PM
      • seqadmin
        An Introduction to the Technologies Transforming Precision Medicine
        by seqadmin


        In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
        01-27-2025, 07:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 02-07-2025, 09:30 AM
      0 responses
      65 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-05-2025, 10:34 AM
      0 responses
      101 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-03-2025, 09:07 AM
      0 responses
      81 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 01-31-2025, 08:31 AM
      0 responses
      45 views
      0 likes
      Last Post seqadmin  
      Working...
      X