STAR Aligner minimum read-length
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?
• 17,142 views
•
link
written
by
informatics bot
1 answer
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.
• 0 views
•
link
updated
by
Ram
•
written
by
informatics bot
•
written
by
informatics bot
Log in to answer this question.
More posts like this
-
(STAR) What causes a high number of reads that are unmapped because they are too short?
written by becky.reese •
-
STAR alignment speed
written by Gilles •
-
STAR aligner can't map too short reads
written by Assa Yeroslaviz
-
STAR ALIGNMENT FOR RNA-SEQ OF PLANT SAMPLE?
written by Meeran •
-
RNA Seq analysis
written by Micro_bioinfo •
-
question about the equation of counts in STAR ReadsPerGene.out.tab
written by boymin2020 •
-
Realign Aligned.sortedByCoord.out.bam to Aligned.toTranscriptome.out.bam for RSEM input
written by mkvalluru •
-
Very low unique mapped reads - STAR alignment
written by jsl •
-
Recommended doc for interpreting STAR alignment stats
written by emblake •
-
Why Am I Getting So Many Unmapped Reads In Star, Classified As "Too Short"?
written by gaelgarcia05
Hi,
Maybe this question can help you ? Do you know the length of those unmapped reads ?
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.