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
          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
          ...
          Yesterday, 10:25 AM
        • SEQadmin2
          How Immunogenomics Decodes Immunity’s Genetic Blueprint
          by SEQadmin2




          The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

          This convergence of genetics, immunology, and computation...
          09-01-2026, 05:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Today, 09:51 AM
        0 responses
        9 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 09-25-2026, 09:06 AM
        0 responses
        32 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 09-23-2026, 11:05 AM
        0 responses
        27 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 09-18-2026, 11:37 AM
        1 response
        47 views
        0 reactions
        Last Post pekgio
        by pekgio
         
        Working...