| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,857,146 | 0 | C | G | 48.4% | ‑0.1 / 56.5 | 31 | L18L (CTC→CTG) ‡ | mutS | methyl‑directed mismatch repair protein |
| Reads supporting (aligned to +/- strand): ref base C (8/8); new base G (11/4); total (19/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.73e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTG > NC_000913/2857062‑2857245 | cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGAtttt > 2:115801/1‑139 (MQ=255) tatCAGTGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGc > 1:735658/1‑94 (MQ=255) tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAAc < 1:115801/139‑1 (MQ=255) cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa < 2:293832/134‑1 (MQ=255) cAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa > 1:293832/1‑134 (MQ=255) aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGtt > 1:509376/1‑139 (MQ=255) gggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGttt > 2:263802/1‑139 (MQ=255) ccGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 2:113759/128‑1 (MQ=255) ccGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa > 1:113759/1‑128 (MQ=255) ccGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATga > 2:192531/1‑139 (MQ=255) gACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg > 1:749563/1‑123 (MQ=255) gACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATg < 2:749563/123‑1 (MQ=255) cATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGacgacg < 1:192531/139‑1 (MQ=255) aaCCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 1:215278/120‑1 (MQ=255) aaCCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa > 2:215278/1‑120 (MQ=255) aTGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgcg > 1:830840/1‑139 (MQ=255) tGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAAcgcgc > 2:427051/1‑139 (MQ=255) gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCg > 2:353059/1‑72 (MQ=255) gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑AGTC‑‑GCTGAAAGCCCAGCATCCCg < 1:353059/72‑1 (MQ=255) gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGAT‑‑‑AGTCGCTGAAAGCCCAGCATCCCg < 1:891988/72‑1 (MQ=255) gTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCg > 2:891988/1‑72 (MQ=255) cAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCa > 2:124489/1‑139 (MQ=255) tAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCctg < 2:865131/74‑1 (MQ=255) tAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCctg > 1:865131/1‑74 (MQ=255) aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa > 2:659599/1‑100 (MQ=255) aaaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGa < 1:659599/100‑1 (MQ=255) aaaTTTCGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAActgctg > 2:775051/1‑139 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGTGATAG‑TC‑‑GCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 2:614649/94‑1 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGT‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa > 1:614649/1‑94 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGC‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa > 2:167989/1‑94 (MQ=255) cGACGCCCATACGCCCATGATGCAGCAGC‑GA‑T‑AGTCGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGaa < 1:167989/94‑1 (MQ=255) | CATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGT‑AT‑C‑TCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGTCGCAACTGCTG > NC_000913/2857062‑2857245 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |