Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Understanding Maq log output

    Hi,

    I have a couple of questions about my Maq map output if anyone knows the answers I'd really appreciate them:

    Code:
    [match_search]  98% processed in 20717.101 sec: 31618863 / 21294334 = 1.485
    [match_search]  99% processed in 20765.187 sec: 28388678 / 20109661 = 1.412
    [match_search] 100% processed in 20768.725 sec: 1164880 / 1048335 = 1.111
    [match_core] sorting the hits and dumping the results...
    [ma_load_reads] loading reads...
    [ma_load_reads] 2000000*2 reads loaded.
    [mapping_count_single] 7, 11, 16, 27
    [maq_indel_pe] the indel detector only works with short-insert mate-pair reads.
    [match_data2mapping] 1922258 out of 4000000 raw reads are mapped with 0 in pairs.
    -- (total, isPE, mapped, paired) = (2000000, 0, 1922258, 0)
    First, what is the "mapping_count_single" - what do these 4 numbers mean?

    Second, this file contained 2,000,000 reads (8 million lines in fastq format), so how come the logfile says x out of 4000000 raw reads are mapped?

    Many Thanks
    Anna

  • #2
    I have a similar issue. My output log reports n_reads*2 because they are paired end, but it looks like that is not the case for your data. My issue is the "x out of n raw reads are mapped" part:

    [match_data2mapping] 1976956 out of 2000000 raw reads are mapped with 1969400 in pairs.
    -- (total, isPE, mapped, paired) = (2000000, 1, 1976956, 1969400)

    What happened to the remaining 23,044 reads? It is not that they are simply unmapped, as I have a number of unmapped reads in my .sam alignment file.

    Any help much appreciated.

    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, Yesterday, 06:37 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    8 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    66 views
    0 likes
    Last Post seqadmin  
    Working...
    X