elements from a pileup. They should all be positioned at the same locus. The alleles considered will align at this locus + 1.
reference bases. The length determines the size of alleles to consider. The first element should be the reference base at locus elements.head.locus + 1.
Map from sequenced allele -> the ReadSubsequence instances for that allele.
Map from sequenced allele -> number of reads supporting that allele.
elements from a pileup.
elements from a pileup. They should all be positioned at the same locus. The alleles considered will align at this locus + 1.
Compute likelihood P(data|mixture) of the sequenced bases (data) given the specified mixture.
Compute likelihood P(data|mixture) of the sequenced bases (data) given the specified mixture.
Reads with mapping quality 0 are ignored.
Map from sequenced allele -> variant allele fraction
log10 likelihood probability, always non-positive
All sequenced alleles that are not the ref allele, sorted by decreasing allelic depth.
Map from allele to read names supporting that allele.
The reference sequence as a string.
reference bases.
reference bases. The length determines the size of alleles to consider. The first element should be the reference base at locus elements.head.locus + 1.
Alt allele with second-most reads.
The sequenced alleles at this site as ReadSubsequence instances.
The sequenced alleles at this site as ReadSubsequence instances.
Every read that is "anchored" on either side of the reference region by a matching, non-variant base is represented here.
Alt allele with most reads.
Total depth, including reads that are NOT "anchored" by matching, non-variant bases.
Fraction of reads supporting the given allele.
Statistics over a PileupElement instances (a pileup).
This is not a case class because we do not intend for it to be serialized, and we want to get an error if we try.
We can consider sequenced alleles of any reference length. For example, if we are evaluating support for a 10-base deletion, we want the keys returned by
allelicDepths
method to be the sequenced alleles aligning to that 10-base region; in particular the reference allele would be length 10. The length of therefSequence
attribute establishes the reference length to be considered.