Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • kesner
    Member
    • Apr 2012
    • 10

    #1

    tophat 2.0.10:WARNING: read pairing issues detected (check prep_reads.log) !

    I am aligning paired end reads that have been trimmed of adapters with cutadap using tophat 2.0.10.

    Paired read log has this:
    /PHShome/bak22/src/tophat-2.0.10.Linux_x86_64/prep_reads: /lib64/libz.so.1: no version information available (required by /PHShome/bak22/src/tophat-2.0.10.Linux_x86_64/prep_reads)
    prep_reads v2.0.9 (3067M)
    ---------------------------
    WARNING: read pairing issues detected (check prep_reads.log) !
    Pair #400691 name mismatch: D64TDFP1:1:1103:2790:5963#0/1 vs D64TDFP1:1:1103:2792:5883#0/2
    11670 out of 41376603 reads have been filtered out
    149176 out of 41376603 read mates have been filtered out
    ------
    I check both records and they are fine. However the entry prior to each, although matching 1 of the paried ends is a totally trimmed records. That is both prior to and after the '+'. They have nothing.
    --------
    I am wondering if this warning is a caused by a problem in tophat which can't handle such reads.

    Has anyone else seen this? Is it a serious problem?
    barry
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    Just use a wrapper for cutadapt, like trim_galore, that can filter out orphaned mates and prevent the read pairing issue.

    Comment

    • kesner
      Member
      • Apr 2012
      • 10

      #3
      thanks I'll try it
      barry

      Comment

      • ronaldrcutler
        Member
        • May 2016
        • 80

        #4
        I ran into the same problem while trying to multithread. I found that not multithreading seemed to fix it but was then running into it again even when not using multithreading.

        I am using a script to have tophat2 go through a folder with read pairs and align those and automatically move on to the next. However, I saw that it would start to randomly use different mate pairs from different lanes instead of the same one from the same lane.

        In order to fix this I just made a new folder and copied the reads and now everything seems to work okay, and with multithreading!

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
          by SEQadmin2




          NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
          Today, 06:46 AM
        • 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

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 10-28-2026, 09:57 AM
        0 responses
        24 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-28-2026, 09:51 AM
        0 responses
        15 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 11:13 AM
        0 responses
        68 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-22-2026, 10:48 AM
        0 responses
        48 views
        0 reactions
        Last Post SEQadmin2  
        Working...