SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat v1.1 with GTF files hyjkim Bioinformatics 7 12-17-2012 07:11 AM
Error joing segment hits in tophat Wallysb01 Bioinformatics 1 03-18-2011 07:40 AM
whether tophat can find interchromosomal hits when junctions are provided... nturay Bioinformatics 0 02-09-2011 05:12 PM
Obtain Spliced Read Names/Mapping with TopHat?? mendl7 Bioinformatics 7 09-20-2010 10:53 AM
TopHat: ambiguous hits are only reported for one of the two paired end reads GoneSouth RNA Sequencing 0 01-19-2010 09:34 AM

Reply
 
Thread Tools
Old 12-08-2011, 08:07 PM   #1
liuxq
Member
 
Location: beijing, china

Join Date: Jun 2010
Posts: 36
Default why --GTF decreases spliced hits for tophat

I am a tophat user. I am now using tophat v1.3.2. My data is 2x75bp human rna-seq data.

I find that when I use --GTF option against ensembl gene set (release 64), the number of spliced hits (with N in cigar field) are less than that without using --GTF option. The number of spliced hits is as follows.
75bp_no_gtf 75bp_gtf
R1 3736037 3375768
R2 4322737 3921965
R3 4142597 3873519

My command is listed below.

tophat --mate-inner-dist 30 --bowtie-n --GTF ensemble_64.gtf --library-type fr-unstranded -o Sample_BG-R-1_after_filtering bowtie_index/hg19/hg19 r-1_R1.fq r1_R2.fq

my question is :
1. why --GTF option will decrease the spliced hits number?
2. why the number of non-spliced reads will also change when using --GTF option?

Thanks very much and look forward to your response.
liuxq is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:58 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO