Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Use library information when outputting XS tag #24

Closed
jdidion opened this issue Feb 19, 2015 · 4 comments
Closed

Use library information when outputting XS tag #24

jdidion opened this issue Feb 19, 2015 · 4 comments

Comments

@jdidion
Copy link

jdidion commented Feb 19, 2015

Currently, STAR only uses intron motif information to determine the strand output in the XS tag. Other aligners (TopHat, HISAT) allow library information to be explicitly specified. When aligning reads from a stranded library (fr-firststrand), the XS tags are highly discordant between STAR and HISAT. This is not an issue when using Cufflinks for assembly (since it allows the library type to be specified), but the new StringTie assembler does not have a way to specify the library type, and relies only on the strand reported in the XS tag.

@alexdobin
Copy link
Owner

This feature is high on my TODO list, hopefully will be added in the next month or so.

@dolittle007
Copy link

dolittle007 commented Sep 20, 2016

I have exactly the same problem. I am planning to make a RNA-seq pipeline using STAR-Stringtie-edgeR.
STAR is the best choice also due to GATK RNAseq pipeline is choosing STAR.

@alexdobin
Copy link
Owner

Hi @dolittle007

please do not comment on a closed issue, since I may not see it quickly enough.
Is your data unstranded or stranded?
For the unstranded data, you can use --outSAMstrandField intronMotif option to add the XS flag to all spliced alignments.

Cheers
Alex

@dolittle007
Copy link

dolittle007 commented Sep 23, 2016

I am handling stranded library. Stringtie doesn't allow the library type to be specified, and require "XS" tag no matter the strandness of reads is. So it forces me to use Cufflinks instead.
Thank you
Ting-You Wang

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants