Question: STAR Aligner minimum read-length
0
gravatar for informatics bot
3.5 years ago by
United States
informatics bot570 wrote:

Is there a parameter in STAR for setting the minimum length of reads to be mapped?

In the "Log.final.out" file, there is a category, "% of reads unmapped: too short," how do I change the "too short" variable?

rna-seq star • 5.8k views
ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by informatics bot570

Hi,

Maybe this question can help you ? Do you know the length of those unmapped reads ?

ADD REPLYlink written 3.5 years ago by Thibault D.510

That question is a decent start, but it doesn't specify what I'm trying to do... I'm comparing percent of reads mapped between STAR and GSNAP. GSNAP automatically throws out reads less that 17 bp. I was wondering if there is a parameter similar in STAR.

ADD REPLYlink written 3.5 years ago by informatics bot570
1
gravatar for informatics bot
3.5 years ago by
United States
informatics bot570 wrote:

I was able to answer my own question:

The following parameters reduced my "Too short" percent from 10.8% to ~4%. Most of the reads went to multi-mapped or "other," but at least it worked. 

--seedSearchStartLmax 30  # increase overall mapping sensitivity

--outFilterScoreMinOverLread 0 --outFilterMatchNminOverLread 0  --outFilterMatchNmin 50 # allow output of short alignments

There apparently isn't a minimum read length threshold for STAR. 

ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by informatics bot570
3

STAR version 2.4.2a provides this option:

    --outFilterMatchNmin 20

In this case, it requires a minimum 20 bp match.

ADD REPLYlink written 3.5 years ago by robm9119150
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: 845 users visited in the last hour