Mappability and read length

Abstract

Power-law distributions are the main functional form for the distribution of repeat size and repeat copy number in the human genome. When the genome is broken into fragments for sequencing, the limited size of fragments and reads may prevent an unique alignment of repeat sequences to the reference sequence. Repeats in the human genome can be as long as 104 bases, or 105 − 106 bases when allowing for mismatches between repeat units. Sequence reads from these regions are therefore unmappable when the read length is in the range of 103 bases. With a read length of 1000 bases, slightly more than 1% of the assembled genome, and slightly less than 1% of the 1 kb reads, are unmappable, excluding the unassembled portion of the human genome (8% in GRCh37/hg19). The slow decay (long tail) of the power-law function implies a diminishing return in converting unmappable regions/reads to become mappable with the increase of the read length, with the understanding that increasing read length will always move toward the direction of 100% mappability.

Exome sequencing of 457 autism families recruited online provides evidence for novel ASD genes

“…We first filtered our variant call set for rare heterozygous coding variants (MAF<=1x10 -4 across all populations represented in gnomAD and ExAC databases). To account for regions in the reference genome that are more challenging to resolve, we removed variant sites found in regions of non-unique mappability (score<1; 300bp), likely segmental duplication (score>0.95), and known low-complexity 82 . We then excluded sites located in MUC and HLA genes and imposed a maximum variant read depth threshold of 500.…”

From Sequence Data to Returnable Results: Ethical Issues in Variant Calling and Interpretation

“…This can be due to a variety of reasons, including ''pseudogenes'' (areas of the DNA that look like a particular gene, but are slightly different), repetitive sequences, and regions of duplications or deletions, all of which naturally occur in all individuals. At read lengths typically employed in NGS, 5-10% of the genome cannot be confidently mapped due to these factors [reviewed in Li and Freudenberg (2014)]. …”

PAFFT: A new homology search algorithm for third-generation sequencers

“…Repetitive elements in the genome, because of the presence of similar or identical sequences, sometimes cause mapping error of the short reads [6][7][8]. The chloroplast genome of higher plants contains the inverted repeat region (IR) [9].…”