Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • annavilborg
    Junior Member
    • Feb 2013
    • 8

    Input Paired-ended data to Tophat

    Hi,

    This should really be an easy thing, but I am getting conflicting advice...so my question is: What is the correct way to input paired-ended (75bp, Illumina Hi-Seq) strand specific data to Tophat? Should I enter the separate file lists for R1 and R2 fastq files? Or should I use a merged file (that is also available for download from our sequencing facility) and input the same file twice?

    The first option seems more logical to me, but I have heard and read that if one inputs the same file twice, Tophat should be able to figure out that this is paired-ended data, split it (I assume) and use that.

    Does anyone know the correct approach?

    Thank you!
  • mknut
    Member
    • Jul 2012
    • 23

    #2
    From the manual:
    Usage: tophat [options]* <index_base> <reads1_1[,...,readsN_1]> [reads1_2,...readsN_2]

    <reads1_1[,...,readsN_1]> A comma-separated list of files containing reads in FASTQ or FASTA format. When running TopHat with paired-end reads, this should be the *_1 ("left") set of files.
    <[reads1_2,...readsN_2]> A comma-separated list of files containing reads in FASTA or FASTA format. Only used when running TopHat with paired end reads, and contains the *_2 ("right") set of files. The *_2 files MUST appear in the same order as the *_1 files.
    So for example:
    Code:
    tophat -o outputDirectoryName -r 50 --library-type fr-unstranded /path/to/genome/reference /path/to/first/read/read_1.fastq /path/to/second/read/read_2.fastq
    Having a quick look at the tophat sourcecode (prep_reads.cpp), it looks like tophat processes reads sequentially, ie. the first read in read_1.fastq has to be the mate of the first read in read_2.fastq. Tophat does not seem to be doing any automatic read mate matching, at least nothing I can see. I think that if you merge the reads into one file and submit it twice as read1 and read2, you will lose synchronisation between the reads, if it works at all. Anyway, the method above is the method suggested in the manual, so I would stick to it.
    To reiterate:
    Should I enter the separate file lists for R1 and R2 fastq files?
    Yes

    Comment

    • annavilborg
      Junior Member
      • Feb 2013
      • 8

      #3
      Thank you! And thanks for explaining how tophat does this, great to get the background. I'll go ahead and run my samples as separate R1 and R2 files.

      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
      186 views
      0 reactions
      Last Post seqadmin  
      Working...