Question: Proper Bam / Tophat Output For Strand Specific Paired-End Data?
1
gravatar for user
6.9 years ago by
user810
United States
user810 wrote:

If I have a strand-specific paired-end library and align with Tophat with an option like fr-firstrand, what is the proper to encode the strand of the pairs in the BAM, and how do mappers like Tophat handle it? With fr-firstrand we know that the first strand is the strand of the read, but by definition the second mate will be on opposite orientation. If the first read in the pair has strand + in its alignment in the BAM file, then will the second read in the pair also have strand + or can we count on it being -? thanks.

ADD COMMENTlink modified 6.9 years ago by Ido Tamir5.0k • written 6.9 years ago by user810
0
gravatar for Ido Tamir
6.9 years ago by
Ido Tamir5.0k
Austria
Ido Tamir5.0k wrote:

Tophat1 should give you:

  +/- or -/+ if you specify fr-X: --> <-- 
  +/+ or -/- if you specify ff-X: --> -->

Tophat1 also gives you singletons as output (i.e. only the first or second read map). Tophat2 can also give your discordant reads i.e. +/+ although fr (--> <--) was specified.

Important is the XS flag which tells you which strand the read comes from depending on X-firststrand or X-secondstrand.

ADD COMMENTlink written 6.9 years ago by Ido Tamir5.0k
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: 1565 users visited in the last hour