title: GeneCodeq: quality score compression and improved genotyping using a Bayesian framework.
abstract: The exponential reduction in cost of genome sequencing has resulted in a rapid growth of genomic data. Most of the entropy of short read data lies not in the sequence of read bases themselves but in their Quality Scores-the confidence measurement that each base has been sequenced correctly. Lossless compression methods are now close to their theoretical limits and hence there is a need for lossy methods that further reduce the complexity of these data without impacting downstream analyses.
  
  
  
  We here propose GeneCodeq, a Bayesian method inspired by coding theory for adjusting quality scores to improve the compressibility of quality scores without adversely impacting genotyping accuracy. Our model leverages a corpus of k-mers to reduce the entropy of the quality scores and thereby the compressibility of these data (in FASTQ or SAM/BAM/CRAM files), resulting in compression ratios that significantly exceeds those of other methods. Our approach can also be combined with existing lossy compression schemes to further reduce entropy and allows the user to specify a reference panel of expected sequence variations to improve the model accuracy. In addition to extensive empirical evaluation, we also derive novel theoretical insights that explain the empirical performance and pitfalls of corpus-based quality score compression schemes in general. Finally, we show that as a positive side effect of compression, the model can lead to improved genotyping accuracy.
  
  
  
  GeneCodeq is available at: github.com/genecodeq/eval CONTACT: dan@petagene.comSupplementary information: Supplementary data are available at Bioinformatics online.
authors.affiliation: PetaGene, Ideaspace, 3 Charles Babbage Rd, Cambridge CB3 0GT, UK.
authors.forename: Daniel L
authors.initials: DL
authors.lastname: Greenfield
authors.affiliation: European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SQ, UK.
authors.forename: Oliver
authors.initials: O
authors.lastname: Stegle
authors.affiliation: PetaGene, Ideaspace, 3 Charles Babbage Rd, Cambridge CB3 0GT, UK.
authors.forename: Alban
authors.initials: A
authors.lastname: Rrustemi
mesh_terms: D001499:Bayes Theorem; D044962:Data Compression / Q000379:methods*; D005838:Genotype; D059014:High-Throughput Nucleotide Sequencing