Hisat2 mapping stringency
1
0
Entering edit mode
4 months ago
MA ▴ 10

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 • 452 views
ADD COMMENT
0
Entering edit mode
4 months ago
rfran010 ★ 1.3k

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..

ADD COMMENT
1
Entering edit mode

Thanks a lot! The answer is really helpful!

ADD REPLY

Login before adding your answer.

Traffic: 2251 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