vg map, paired-end reads, multi-mapping and MAPQ
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.
• 1,149 views
•
link
written
by
cristian.groza.mcgill •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
warning[vg::giraffe]: Cannot cluster reads with a fragment distance smaller than read distance
written by cassiwatt •
-
Mapping quality in `vg giraffe`
written by Wenhai •
-
Mate reads are aligned to different chromosomes
written by wormball •
-
Variation in mapping percentage with genome
written by onkar •
-
Exclude reads with mate mapped to a different chr
written by emma.a •
-
Getting uniquely aligned properly paired reads
written by nadiabeg.comsats •
-
Dragen alignment failing/terminating
written by Floydian_slip
-
Filtering BAM file to include PE reads that did not map, even though their mate did
written by nute11a •
-
Calculating percentage of reads mapped to a reference genome
written by rujiporn.sun •
-
map methylome data with GKNO mosaik program
written by zhiguangli88 •
