The tool Calculate Unique Molecular Index Groups annotates the mapped reads with a "Unique Molecular Index group ID", that is identical for reads that are determined to belong to the same UMI. The tool can be found in the Toolbox at:
Tools | QIAseq Panel Expert Tools () | QIAseq DNA Panel Expert Tools () | Calculate Unique Molecular Index Groups ()
In the first dialog (figure 6.59), select a read mapping of reads that were previously annotated with UMI annotations.
The grouping of reads into UMI groups works as follows:
- The tool groups reads that
- start at the same position based on the end of the read to which the UMI is ligated (i.e. Read2 for paired data),
- are from the same strand, and
- have identical UMIs.
- The tool then fuzzy merges singleton groups into non-singleton groups, if the UMI of the singleton group can be made into UMI of non-singleton group by introducing a SNP, and if the non-singleton group is the biggest of such group (i.e., if two different introduced SNPs yields two different non-singleton groups, the biggest one is chosen).
- Additional merging of singletons and small groups into bigger ones can happen depending on the parameters set for the tool.
It is possible to change the following parameters (figure 6.60):
- Do not fuzzy match Groups will be merged only if they match exactly.
- Allow one mismatch Groups will be merged if they are at most one mismatch apart.
- Allow one mismatch/deletion/insertion Groups will be merged if they are at most one mismatch, deletion or insertion apart.
- Exclude ambiguously mapped reads is checked by default.
- Maximum relative size difference between merged groups will merge small groups into bigger ones if the size difference between the two groups is smaller than a certain value (set at 0.1 per default).
- Always merge singleton groups When this option is checked, a singleton UMI group is merged with a non-singleton group even if the "Maximum relative size difference between merged groups" threshold is not met.
Click Next to choose whether to Open or Save the resulting read mapping of reads which now have a "UMI group ID" annotation.
A report can also be generated. It contains:
- A summary table with the following information:
- Reads in input: Reads that were aligned to the reference
- Reads mapped multiple places (discarded): Reads that aligned to the reference in multiple places, and thus discarded
- Groups merged
- Groups not merged due to >1 candidate of equals size
- Groups not merged due to parameter thresholds
- Group size table and plots described in UMI group sizes.
Note: When the group sizes (the number of reads in UMI groups) are very large (in most cases more than 10 reads in a UMI group is not desirable), this can indicate problems, such as quality issues with the sample. It can also indicate that the sequencing depth could be reduced.