PureCN removemapq0 parameter
1
0
Entering edit mode
twtoal ▴ 10
@twtoal-15473
Last seen 13 months ago
United States

I used the option --removemapq0 when I ran Coverage.R, since it seemed reasonable that reads with 0 mapping quality should be removed.  There isn't any non-obvious reason why this is a bad idea, is there?

 

PureCN removemapq0 • 845 views
ADD COMMENT
0
Entering edit mode
@markusriester-9875
Last seen 21 months ago
United States

You'll probably get a significant drop in off-target coverage with removemapq0. There is information in those reads and on-target regions with poor mapping are ignored anyways.

I would recommend following the vignette closely and using the defaults for test runs and then work from there. 

ADD COMMENT
0
Entering edit mode

I checked a few datasets and it appears to be safe to remove those reads. I didn't see a huge impact. Let me know if it improves your results. 

ADD REPLY

Login before adding your answer.

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