Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Olalla
    Member
    • Aug 2014
    • 11

    #16
    Thanks a lot to all for your suggestions. Finally I got the program running, by deleting the bloody white spaces in the read names and then setting -NW:mrnl=80 in MIRA, after that all worked perfectly

    Comment

    • maubp
      Peter (Biopython etc)
      • Jul 2009
      • 1544

      #17
      If you had several lines starting GIMXFMA02G21Y1 but the descriptions were different, something strange has happened. I suspect you still have a problem in your data, and all you have done by editing the read names is avoid MIRA's sanity test without fixing the real problem.

      Comment

      • maubp
        Peter (Biopython etc)
        • Jul 2009
        • 1544

        #18
        Originally posted by Olalla View Post
        So, the read names in my fastq files appear as follows:
        GIMXFMA02G21Y1 length=60 xy=2788_0299 region=2 run=R_2010_06_09_09_17_36_ that is, long names with spaces. So I think that what MIRA is doing is juts taking the first 14 characters as the read name (e.g. GIMXFMA02G21Y1),
        Yes, that is the identifier - everything after the space is an optional description in FASTQ.
        Originally posted by Olalla View Post
        which in fact are repeated among many of the files that I do have (I have found common lines when comparing files including only these names in the lines). However, when I search for common lines in files including all information in headers (like GIMXFMA02G21Y1 length=60 xy=2788_0299 region=2 run=R_2010_06_09_09_17_36_), the output of the search is that there are no common lines between any of the files, so I think that I should first eliminate spaces, maybe replacing them by ":".
        Something is wrong, GIMXFMA02G21Y1 should be unique. What does this give?
        Code:
        grep GIMXFMA02G21Y1 your_file.fastq

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM
        • SEQadmin2
          New Genomics Technologies Take Aim at Long-Standing Limits
          by SEQadmin2


          Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

          We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
          ...
          09-28-2026, 10:25 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-22-2026, 11:13 AM
        0 responses
        31 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 10:48 AM
        0 responses
        22 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-15-2026, 10:52 AM
        0 responses
        39 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-14-2026, 05:39 AM
        0 responses
        45 views
        0 reactions
        Last Post SEQadmin2  
        Working...