Question: Preventing alignment to new Splice Junctions with STAR
0
gravatar for MarVi
9 days ago by
MarVi20
MarVi20 wrote:

Hello to all,

I hope everyone is well! I want to ask if anyone knows if there are any parameters in STAR that prevent the alignment of reads to unannotated splice junctions? I want to speed the alignment task by limiting alignments to splice junctions that are not known. I already included --sjdbScore to 0, --alignSJoverhangMin to 10. But, it still takes some time and I still get splices that for the moment don't interest me.

Thanks in advance for your comments!

ADD COMMENTlink modified 8 days ago by benformatics2.0k • written 9 days ago by MarVi20
1
gravatar for benformatics
8 days ago by
benformatics2.0k
ETH Zurich
benformatics2.0k wrote:

As per the manual:

The parameter: --alignSJoverhangMin sets the:

minimum overhang for unannotated junctions

Whereas, the parameter: --alignSJDBoverhangMin sets the:

minimum overhang for annotated junctions

Note the presence of SJDB vs. SJ in the latter, this refers to your GTF database and/or 2-pass splice junctions, etc... Essentially, what you consider to be annotated.

In other words, for your specific case you should change --alignSJoverhangMin <N> so that N is a number greater than your read length.

You should no longer get any spliced reads to unannotated sites with overhangs shorter than N (which means no reads should splice to unannotated junctions).

Assuming standard RNA-seq you may be fine just setting it to something like 1000.

ADD COMMENTlink modified 8 days ago • written 8 days ago by benformatics2.0k

You can read more here: https://groups.google.com/u/1/g/rna-star/c/J6qH9JCysZw/m/VmgoGKSE9qQJ

If you want to modify the actual outputs in the splice junction files.

ADD REPLYlink written 8 days ago by benformatics2.0k

Thanks so much! You were of much help!

ADD REPLYlink written 8 days ago by MarVi20
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: 1561 users visited in the last hour