View Single Post
Old 01-31-2012, 09:25 AM   #1
Location: New England

Join Date: Nov 2010
Posts: 27
Default Tophat paired end alignment help: comma or blank to separate two pairs?

Hi All,
I am using the current command to run a PE alignment in Tophat:
$ tophat -p 4 -m 2 -r 200 /index/hg19 /data/raw/N1_1.fa,N1_2.fa

my question is, both N1_1.fa and N1_2.fa are the paired end reads generated from a single human tissue sample. Should I use comma or blank to separate those two? And here is the log for the alignment, is 76% alignment looks normal? Thanks!

# reads processed: 242980185
# reads with at least one reported alignment: 185782139 (76.46%)
# reads that failed to align: 56208030 (23.13%)
# reads with alignments suppressed due to -m: 990016 (0.41%)
Reported 457696000 alignments to 1 output stream(s)
mediator is offline   Reply With Quote