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

Missing parameter required by StringTie if using STAR as aligner #857

Closed
xianggenti opened this issue Aug 4, 2022 · 1 comment
Closed
Labels
bug Something isn't working
Milestone

Comments

@xianggenti
Copy link

xianggenti commented Aug 4, 2022

Description of the bug

Ref: https://github.com/gpertea/stringtie
Missing XS tag from STAR.

bam file Header:
@pg ID:STAR PN:STAR VN:STAR_2.6.1d CL:STAR --runThreadN 12 --runRNGseed 0 --genomeDir STARIndex --readFilesIn EngTreg2_1_1_val_1.fq.gz EngTreg2_1_2_val_2.fq.gz --readFilesCommand zcat --outFileNamePrefix EngTreg2_1. --outSAMtype BAM Unsorted --outSAMattributes NH HI AS NM MD --outSAMattrRGline ID:EngTreg2_1 SM:EngTreg2_1 --outFilterMultimapNmax 20 --alignSJDBoverhangMin 1 --sjdbGTFfile genes.gtf --quantMode TranscriptomeSAM --quantTranscriptomeBan Singleend --twopassMode Basic

Quote from Stringtie:
"Any SAM record with a spliced alignment (i.e. having a read alignment across at least one junction) should have the XS tag to indicate the transcription strand, i.e. the genomic strand from which the RNA that produced this read originated. TopHat and HISAT2 alignments already include this tag, but if you use a different read mapper you should check that this tag is also included for spliced alignment records. STAR aligner should be run with the option --outSAMstrandField intronMotif in order to generate this tag."

Command used and terminal output

No response

Relevant files

No response

System information

No response

@xianggenti xianggenti added the bug Something isn't working label Aug 4, 2022
@xianggenti xianggenti changed the title Use STAR as alignor parameters for Stringtie Missing parameter for Stringtie if use STAR as alignor, Aug 4, 2022
@xianggenti xianggenti changed the title Missing parameter for Stringtie if use STAR as alignor, Missing parameter for Stringtie if use STAR as alignor Aug 4, 2022
@drpatelh drpatelh added this to the 3.9 milestone Sep 25, 2022
@drpatelh drpatelh changed the title Missing parameter for Stringtie if use STAR as alignor Missing parameter required by StringTie if using STAR as aligner Sep 25, 2022
drpatelh added a commit to drpatelh/nf-core-rnaseq that referenced this issue Sep 27, 2022
@drpatelh
Copy link
Member

Fixed in 3716c62

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants