Background
The identification of genetic variants has great significance in genetic research. To call variants using next-generation sequencing data, current methods rely primarily on mapped reads produced by a separate read aligner without taking into account existing genetic variants[1]. Thus, these methods usually require a large number of reads (high coverage) to be able to detect variants accurately[2]. Moreover, the separation of read alignment and variant calling results in a workflow is complex and involves many separate steps and different tools[3].