Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Anthony.287
    Member
    • Dec 2010
    • 95

    RLMIDs vs GSMIDs

    Hello all!
    I'm finally starting to get the hang of sfffile, and am finding some things that I can't explain. For instance, one particular region (a pool of 4 RLMID tagged rapid libraries) on the most recent PTP had 43,949 reads. 43,792 of those got sorted in sff's based on the RLMIDs. Out of curiosity, I sorted using the GSMIDs, and ~15,000 were sorted then. What could be happening?

    [PMGF@PMGF-454 ~]$ sfffile -s GSMIDs /data/R_2012_11_29_16_12_43_FLX05100626_Administrator_Shotgun/sff/HYPJBHH12_Shotgun.sff
    Reading the input SFF file(s)...
    Generating the split SFF file(s)...
    MID1: 3 reads written into the SFF file.
    MID2: 0 reads found.
    MID3: 13 reads written into the SFF file.
    MID4: 10898 reads written into the SFF file.
    MID5: 4660 reads written into the SFF file.
    MID6: 0 reads found.
    MID7: 0 reads found.
    MID8: 0 reads found.
    MID9: 1 reads written into the SFF file.
    MID10: 0 reads found.
    MID11: 0 reads found.
    MID12: 0 reads found.
    MID13: 0 reads found.
    MID14: 0 reads found.
    [PMGF@PMGF-454 ~]$ sfffile -s RLMIDs /data/R_2012_11_29_16_12_43_FLX05100626_Administrator_Shotgun/sff/HYPJBHH12_Shotgun.sff
    Reading the input SFF file(s)...
    Generating the split SFF file(s)...
    RL1: 24344 reads written into the SFF file.
    RL2: 11022 reads written into the SFF file.
    RL3: 3692 reads written into the SFF file.
    RL4: 4734 reads written into the SFF file.
    RL5: 0 reads found.
    RL6: 1 reads written into the SFF file.
    RL7: 0 reads found.
    RL8: 0 reads found.
    RL9: 0 reads found.
    RL10: 0 reads found.
    RL11: 0 reads found.
    RL12: 0 reads found.


    Can the RLMIDs be misinterpreted as GSMIDs, and vice-versa? Any ideas?

    Thank you!
  • flxlex
    Moderator
    • Nov 2008
    • 412

    #2
    Yes, presumably because by default, up to two mismatches are allowed between the MID and the read. If you try allow only for one error, this should not happen. Either copy the Installation_path/454/config/MIDConfig.parse file to a new location and edit it:
    Code:
    mid = "MID4", "AGCACTGTAG", 2;
    becomes
    Code:
    mid = "MID4", "AGCACTGTAG", 0;
    Then run
    Code:
    sfffile -s GSMIDs -mcf modified_MIDConfig.parse /data/R_2012_11_29_16_12_43_FLX05100626_Administrator_Shotgun/sff/HYPJBHH12_Shotgun.sff
    Or to indicate specific MIDs and allow no errors (adding '/0') run
    Code:
    sfffile -s GSMIDs:mid4/0,mid5/0@/data/R_2012_11_29_16_12_43_FLX05100626_Administrator_Shotgun/sff/HYPJBHH12_Shotgun.sff
    sfffile -s RLMIDs:mid1/0,mid2/0,mid3/0,mid4/0@/data/R_2012_11_29_16_12_43_FLX05100626_Administrator_Shotgun/sff/HYPJBHH12_Shotgun.sff

    Comment

    Latest Articles

    Collapse

    • seqadmin
      New Genomics Tools and Methods Shared at AGBT 2025
      by seqadmin


      This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

      The Headliner
      The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
      03-03-2025, 01:39 PM
    • seqadmin
      Investigating the Gut Microbiome Through Diet and Spatial Biology
      by seqadmin




      The human gut contains trillions of microorganisms that impact digestion, immune functions, and overall health1. Despite major breakthroughs, we’re only beginning to understand the full extent of the microbiome’s influence on health and disease. Advances in next-generation sequencing and spatial biology have opened new windows into this complex environment, yet many questions remain. This article highlights two recent studies exploring how diet influences microbial...
      02-24-2025, 06:31 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-20-2025, 05:03 AM
    0 responses
    17 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-19-2025, 07:27 AM
    0 responses
    18 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-18-2025, 12:50 PM
    0 responses
    19 views
    0 reactions
    Last Post seqadmin  
    Started by seqadmin, 03-03-2025, 01:15 PM
    0 responses
    185 views
    0 reactions
    Last Post seqadmin  
    Working...