Structural Variation Discovery: Delly's Output
I have a question about Delly's MAPQ value. Before calling SVs on my data, I filtered all reads with MAPQ 20 or less. In Delly's output, however, many of the SV predictions have MAPQ values less than 20. So my question is how does Delly calculate MAPQ score?
Thanks.
1 answer
For SV discovery, delly uses MAPQ and filters reads with MAPQ < q where q is by default 1.
For SV genotyping, delly hashes the mapping qualities of both pairs and uses the minimum of these 2 mapping qualities as the paired-end quality. For bwa alignments it uses the AS tag (alignment score) instead of MAPQ because that worked better in genotyping benchmarks.
Log in to answer this question.
More posts like this
-
How to connect DELLY output to ShatterSeek input
written by Xiaofan •
-
featureCounts generates low assigned rate
written by yuabrahamliu •
-
BLAST bit-score/score to E-value comparison for HSPs
written by cklinger •
-
The interpretation of SV output file
written by nina •
-
SNP calls all seem to have quality of zero
written by novice
-
Qual score in VCF file
written by swatipuraanik •
-
Pipeline of structural variation calling using multiple tools for multiple samples
written by michealsmith
-
What is your opinion on the network score of IPA (Ingenuity Pathway Analysis) analysis?
written by itsunoda •
-
Puzzle on size of Bam file
written by devinliao0918 •
-
BWA MEM -T option. What does really do and how?
written by cpcantalapiedra
