SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Tophat paired end alignment help: comma or blank to separate two pairs? (http://seqanswers.com/forums/showthread.php?t=17268)

mediator 01-31-2012 09:25 AM

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)

jbrwn 01-31-2012 09:40 AM

a space between file_1 and file_2.

mediator 01-31-2012 09:41 AM

@jbrwn
Thanks! Hope that improves my alignment rate.

jbrwn 01-31-2012 09:44 AM

Quote:

Originally Posted by mediator (Post 63377)
@jbrwn
Thanks! Hope that improves my alignment rate.

i typically include --mate-std-dev <int> as well, though i'm not sure how much difference it makes. this information should be available from the bioanalyzer output or someone-from-the-lab output.


All times are GMT -8. The time now is 09:47 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.