Merfin: improved variant filtering, assembly evaluation and polishing via k-mer validation.

Journal Article (Journal Article)

Variant calling has been widely used for genotyping and for improving the consensus accuracy of long-read assemblies. Variant calls are commonly hard-filtered with user-defined cutoffs. However, it is impossible to define a single set of optimal cutoffs, as the calls heavily depend on the quality of the reads, the variant caller of choice and the quality of the unpolished assembly. Here, we introduce Merfin, a k-mer based variant-filtering algorithm for improved accuracy in genotyping and genome assembly polishing. Merfin evaluates each variant based on the expected k-mer multiplicity in the reads, independently of the quality of the read alignment and variant caller's internal score. Merfin increased the precision of genotyped calls in several benchmarks, improved consensus accuracy and reduced frameshift errors when applied to human and nonhuman assemblies built from Pacific Biosciences HiFi and continuous long reads or Oxford Nanopore reads, including the first complete human genome. Moreover, we introduce assembly quality and completeness metrics that account for the expected genomic copy numbers.

Full Text

Duke Authors

Cited Authors

  • Formenti, G; Rhie, A; Walenz, BP; Thibaud-Nissen, F; Shafin, K; Koren, S; Myers, EW; Jarvis, ED; Phillippy, AM

Published Date

  • June 2022

Published In

Volume / Issue

  • 19 / 6

Start / End Page

  • 696 - 704

PubMed ID

  • 35361932

Pubmed Central ID

  • PMC9745813

Electronic International Standard Serial Number (EISSN)

  • 1548-7105

Digital Object Identifier (DOI)

  • 10.1038/s41592-022-01445-y

Language

  • eng

Conference Location

  • United States