Biostar Test Site

This is site is used for testing only. Visit: https://www.biostars.org to ask a question.

vg map, paired-end reads, multi-mapping and MAPQ
0
0
Entering edit mode
6 months ago

Hi,

I have a question about how vg treats reads that are multi-mapping to different chromosomes but are then anchored to the correct location by their mate read. Normally, a multi-mapping read is assigned a MAPQ of 0. However, if vg manages to rescue the pair, will it adjust the MAPQ to be above 0 for both fragments? I know this is the case with bwa mem. I just want to make sure that the following command will not discard multi-mapping reads that are anchored to the correct location by their mate: vg filter aligned.gam -t 40 -r 0.90 -fu -m 1 -q 10 -D 999 -x index.xg > filtered.gam

Thank you in advance.

vg • 165 views
ADD COMMENT

Login before adding your answer.

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