This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Hisat2 mapping stringency

Dear all, I am running hisat2 with default settings in genomic mapping analysis. I want to run hisat2 with more strict parameters to avoid mismatched reads. Does anyone provide me some suggestions on setting relavent hisat2 options? Thank you very much!

hisat2

1 answer

You can adjust this parameter:

Scoring:
  --mp <int>,<int>   max and min penalties for mismatch; lower qual = lower penalty <6,2>
  --sp <int>,<int>   max and min penalties for soft-clipping; lower qual = lower penalty <2,1>

I would suggest a parameter sweep where you run multiple tests with different values to see how it affects mapping.

One consideration too is that hisat2 is soft-clipping reads, so you may want consider adjusting those options as well, You can see the scoring option above. There's also an option to disable soft-clipping entirely..

Thanks a lot! The answer is really helpful!

Log in to answer this question.