Question: How To Use Tophat With Strand Specifc
0
gravatar for camelbbs
6.2 years ago by
camelbbs650
China
camelbbs650 wrote:

Hi,

While I use tophat to align rnaseq, if my data is strand-specific, do I need to add the option --library-type?

What is the default settings in tophat for that.

thanks,

Ch

rna-seq • 5.8k views
ADD COMMENTlink modified 6.2 years ago by Ido Tamir5.0k • written 6.2 years ago by camelbbs650

What do you mean by strand specific?

ADD REPLYlink written 6.2 years ago by KCC3.9k

The library was prepared by strand-specific assay

ADD REPLYlink modified 6.2 years ago • written 6.2 years ago by camelbbs650
5
gravatar for Ido Tamir
6.2 years ago by
Ido Tamir5.0k
Austria
Ido Tamir5.0k wrote:

yes library-type. Then you need to know the way the library was prepared. Currently the dUTP protocol is mostly used and for this you need to specify --library-type fr-firststrand. Default is fr-unstranded.

ADD COMMENTlink written 6.2 years ago by Ido Tamir5.0k
1

Thanks a lot !! I found some public data don't have the library preparation info, such as ERP000546. If I don't know that info, how to set the parameter. Is there big affection to mapping result?

ADD REPLYlink written 6.2 years ago by camelbbs650

The reads will be always map to the same locations irrespective of correctness of library-type (not tested). Downstream analysis will be affected: cufflinks and htseq count when the data is in the wrong orientation or no orientation is given.

Map 1. with unstranded and then you see the result (i.e. if the data is really stranded in a genome browser).

ADD REPLYlink written 6.2 years ago by Ido Tamir5.0k

Hi Ido, Thanks very much for your help. I just want to ask how to identify if the data are stranded through genome browser. I've mapped 1. with unstranded and got the bam file. I am using the IGV.

ADD REPLYlink written 6.2 years ago by camelbbs650
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 968 users visited in the last hour