tophat- paired end reads combined in one file
Hi All,
Just a quick question re tophat- How does tophat handles the paired-end data when both the reads (left and right) are in one file. I mapped this data (still running) using tophat without separating the paired-end reads in two different files. Will top handle this or I should separate the files?
Thanks,RT
1 answer
You should separate the FASTQ reads into two files for Tophat. It will not recognize them as paired-end reads if they are together in the same file.
From the manual:
##Using TopHat
Usage: `tophat [options]* <genome_index_base> <reads1_1[,...,readsN_1]> [reads1_2,...readsN_2]`
When running TopHat with paired reads it is **critical** that the `*_1` files an the `*_2` files appear in separate comma-delimited lists, and that the order of the files in the two lists is the same.
• 0 views
•
link
updated
by
Ram
•
written
by
Ashutosh Pandey
•
written
by
Ashutosh Pandey
Log in to answer this question.
More posts like this
-
How bedtools genomecov handles reads without both ends mapped in paired sequencing
written by rimelof •
-
Separate unaligned reads in fastq format using bowtie
written by MAPK
-
reads diffreence in fastqc and Tophat input
written by zz105 •
-
Count right and left mate pairs from Sam files
written by PAn •
-
Tophat vs Tophat-fusion concordant alignment rate
written by IP
-
Tophat, low alignement left read
written by claire.morandin •
-
How to calculate FPKM from the reads mapped to two genomes
written by GR
-
Mapping paired end and single end reads together
written by GR
-
need to separate read1 and read2 to run tophat2
written by Ann
-
tophat output files containing the reads which mapped uniquely as a pair
written by trakhtenberg
