Hello All,
I have been trying to figure out what would be the best parameters to use in the STAR aligner to find all possible alignments to an rna read. I realized that somehow STAR limits alignments or does not find a way to do alignments in soft mask regions. As I tried GSNAP and with this aligner I found several different positions in the genome where the read has been located. Even though, there are overlaps of alignments between STAR and GSNAP, there are some interesting positions that STAR is not able to align even though there is no mismatch or soft-clip that would be needed in order to STAR do the alignment.
I use the same parameters for both aligners, like : set the max of multimapping (very high) and set the minimum match to output the alignment. The other parameters are pretty much the basic. I have a preference for STAR as it finds a more alignments compared to GSNAP, however in some cases it doesn't retain some positions that I would be interested.
I would be glad to hear your suggestions!
Thanks in advance.
•
written
by
MarVi •
1 answer
STAR doesn't limit alignment to soft-masked regions, it more generally puts limits on multi-mapping. You can adjust these with options like --outFilterMultimapNmax and --outAnchorMultimapNmax, as suggested for tools like TETranscripts, which seem to be dealing with datasets similar to yours.
Log in to answer this question.
-
Isaac aligner V4 genome aligner questions
written by DJBill •
-
Typical percentage of multimapping reads in human rna-seq?
written by srhic •
-
How to set soft-clipped bases option in IGV batch script
written by JellyFrog •
-
Use of aligners (preferably STAR) for read-barcode matching
written by mb •
-
Can you generate a mappability mask with STAR (or another splice-aware aligner)?
written by mercerf93 •
-
Preventing alignment to new Splice Junctions with STAR
written by MarVi •
-
BWA aligning forward and reverse differently
written by dorothea.emig •
-
can I use soft-clip RNA-seq reads to summarize gene counts?
written by tujuchuanli
-
STAR unique reads
written by TEman •
-
How to pull out sequences with repeat elements from RepeatMasker output file?
written by san.san •
You may find this prior thread of interest: Hard vs soft masked references while aligning with STAR
There is an old thread in which Alex mentions (post #2) that soft-masked regions are not treated any differently by STAR.
Thanks for you reply ! So, in conclusion, STAR isn't a good aligner to align reads to repetitive regions! I didn't find any explanation of why this happens. I made some tests with Hard and Soft masked reference and actually, the result was the same. It seems that STAR simply doesn't align to repetitive zones in the genome. Then, I have either choose another aligner or use STAR + another aligner that do the job for the repetitive regions.
Well, thanks again for your reply!