Hello,
I am using qualimap to estimate the quality of some assemblies, and I am not sure to understand exactly what means the "general alignment error rate" given?
Thank you in advance if anyone could help me with that! Celine
1 answer
Did you mean to say alignments?
From qualimap help:
Mismatches and indels
The section reports general alignment error rate (computed as a ratio of total collected edit distance to the number of mapped bases), total number of mismatches and total number of indels (computed from the CIGAR values). Additionally fraction of the homopolymer indels among total indels is provided. Note, the error rate and mismatches metrics are based on optional fields of a SAM record (NM for edit distance, MD for mismatches). The features are not reported if these fields are missing in the SAM file.
Log in to answer this question.
-
Clarification on sequence duplication levels and Qualimap estimated duplication rate
written by keremozdel9 •
-
rRNA rate from RNA seq data- Qualimap
written by Ahmed.waraky •
-
What is -e in trim_galore
written by evcon •
-
FastTree on Mac - Segmentation Fault 11
written by celine.petitjean •
-
Question about hisat2 parameter "--no-templatelen-adjustment"
written by ChrisT •
-
Bioconductor package NOISeq: Normalization
written by KVC_bioinfo
-
Canu : how to increase the required coverage for the trimming and assembly?
written by celine.petitjean •
-
OMA orthology. Error, (in RemoveSpaces) can only assign a character
written by celine.petitjean •
-
Post alignment using RseQC
written by KVC_bioinfo
-
Base Quality Out Of Range For Specifed Quality Type
written by akshayb04 •
Hi Genomax, thank you for your answer. I have read the description of the parameter in qualimap, but I still don't understand exactly how it is calculated. It does not correspond exactly to the number of mismatches and indels reported to the number of bases in my reads, so I wanted to know exactly how it was calculated. thank you! Celine