Will the --outfiltermismatchNmax in the STAR works for the single end reads?
1
0
Entering edit mode
5.9 years ago
codezy ▴ 10

Dear all,

I have some single end reads and I would like to do the mapping step with STAR. But depends on the manual of STAR the -- outFilterMismatchNmax is for setting the maximum number of mismatches per pair. So I was wondering if this --outFilterMismatchNmax will also work for the single end reads? And how would it work?

RNA-Seq alignment • 3.0k views
ADD COMMENT
0
Entering edit mode

In case you have varying read lengths , you can also use --outFilterMismatchNoverLmax:

int: alignment will be output only if its ratio of mismatches to mapped length is less than this value

ADD REPLY
0
Entering edit mode
5.9 years ago
ATpoint 86k

The manual says int: alignment will be output only if it has no more mismatches than this value., it does not seem to be an explicit paired-end option. In the worst case, take some reads, introduce a couple of mismatches to a read that perfectly matches and then see what happens.

ADD COMMENT
0
Entering edit mode

@ATpoint Thank you for your answer. I have also read this int, but in page 9 of STAR manual they said "--outFilterMismatchNmax 999 maximum number of mismatches per pair, large number switches off this lter" So I confused about that.

ADD REPLY

Login before adding your answer.

Traffic: 2447 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6