| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,745,019 | 0 | C | G | 38.7% | 38.5 / 10.1 | 31 | G72G (GGC→GGG) | yiaN | 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein |
| Reads supporting (aligned to +/- strand): ref base C (6/13); new base G (9/3); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.90e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTAT > NC_000913/3744783‑3745240 | gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCA‑TGAATgtgggggg > 1:1085326/1‑239 (MQ=255) cGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACA‑TTAATgggggggggctt > 1:1744986/1‑234 (MQ=255) aaGGAGCTTAATCGTGGCGGTGGTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTCCCTATTGCCTGGGCACTGTTGTTGTGGGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCa < 2:248307/238‑1 (MQ=255) gCTTAATCTGGGCTGTGCTGATTTTTCTGTGCTGTCTGTTGGGTGGTATTGCTAGCTGTTTGCCTAGTGCCTGGGCACTGGTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACg < 2:384825/238‑1 (MQ=255) gCTTAATCATGGCGGTGCTGATTTGTCTGTGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACg < 2:1897833/238‑1 (MQ=255) cgttAGTCATAGCTGTGCTTGTTTTGCTGGGTTGTCTGTGGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACg < 2:1439534/236‑1 (MQ=255) gtctgttttgttttttGTGAGCTGTCTTTTGGTTGGTATGGTGTGCGGTTGGCCTATTGCCGGGGCACTTTTGTTGGGCGGGGCGGCGTTAATGTTCTGGCTGGACATTTTTGATGTCCAGTTCATGGCACAACCGCTGGTGGACGGAGCCGATAGTTTCTCCCTGCGGGCGTTTCCGTTCTTTGGGCGGGCGGGGGAATTCT‑TGAATGGGGGGGTTTTGCCAAAACG‑CATTTTTTAcc < 2:1369677/228‑1 (MQ=255) tgtgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCTTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGACCTGc < 1:259331/239‑1 (MQ=255) tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATAC‑TGAATGCGGGGGGCCTGTTAAAAAG‑CCTTTTTGACCTGccc > 1:1369677/1‑238 (MQ=255) tgCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACA‑TGAATTCGGGGGGCTTGTCAAAACA‑CAATGTTGACCTccc > 1:502280/1‑238 (MQ=255) tCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATTC‑TGAAAgcgggc > 1:756345/1‑199 (MQ=255) tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAAACA‑TGAATGCGGGCGGCTTGGCAAAACG‑CATTGTTGGACTGCCGGTGAAAATGGGGGGGGA‑Ta > 1:1086258/1‑239 (MQ=255) ggtggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATAA‑TGAAAGAGGGGGGCTTGTCCAAACG‑CATTGTTTACCTGCCGCTGAAACAGGGGGGGGA‑AAAACCAg > 1:1618789/1‑238 (MQ=255) tATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGGGGCTTGGCAAAAAG‑CAATGTTTACCTGCCGCTGAAAATGGGGGGGGA‑TAAACCAGGCgggg > 1:129221/1‑238 (MQ=255) cTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCt > 1:644491/1‑177 (MQ=255) cTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCt < 2:644491/177‑1 (MQ=255) ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAATCA‑TGAATGCGGGCGGCCTGTCAAAACG‑CCTTGTTGAGCTGCCGGTGAAACTGGGGGGGGA‑TAAACCAGGGgggctgtgccgcgc > 1:1567220/1‑234 (MQ=255) ggTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGGGGCTTGTCAAAACG‑CATTGTTGAACTGCCGATGAAACTGGGGGGGCA‑TAAACCAGGCgggctgggggtag > 1:2030279/1‑233 (MQ=255) ccTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCaaa < 1:1605045/173‑1 (MQ=255) ccTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCt < 1:1439399/159‑1 (MQ=255) gggCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGGAAATA‑TTAATTCGGGGGGCTTGTCAAAAAG‑CCttgtt > 2:1605045/1‑173 (MQ=255) tgttgttgtGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTTTGCTGGCGGGGGGAATCATTGAATCGGGGGGGCTT‑TCCAAACGCCTTTTTTTACCGCCCGTTGAAA‑TGGGGGGGCACAAAACCACGGGGGC‑GGGGTACGACGGCGCGCTGGCAggcctggta > 2:2002126/1‑231 (MQ=255) gCGGGGCGGCGTTAATGTTCTGGCTGGACTTGTTTGTTGTCCAGTTCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAACCG‑CATTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGc < 1:1872328/239‑1 (MQ=255) tgtggcgctgtCATGGCAACACCGCTGGTGAGGGGACGCGATCACTGCTCCCGGGCGGCCTTGCTGTTCTGGTGGCTGGCGGGGGAGATCA‑TGAGTGGGGGCGGGTGGTAAAACCG‑CTGTGTTCACCTGCGGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGTCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGTATACCGCCGCCGt < 2:1744986/229‑1 (MQ=255) aaaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGGGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGAACTGCCGATGAAACTGGGGGGGCA‑TAAAc > 1:787376/1‑135 (MQ=255) tGACCGGAGCCGATAGCTTCCCCCTGCTGGCGTTTCCGTTTTTTGTGCGGGCGGGTGAAATCA‑TGAATGCGGGGGGCTTGTAAAAACG‑CATTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGAt < 1:1018728/166‑1 (MQ=255) ggCCGATGGCTTCCCCCGGCGGGCGATCCCTTTTTTTGTGGTGGGGGGTAAAATCA‑TGAATGCGGGGGGTTTGTCAAAACG‑CATTGTTGCCCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGcc < 1:2002126/237‑1 (MQ=255) ctacccccGGCTGGGGTTGCGGTTCGTTGTGCTGGCGGGTGAAATCA‑GGAATGCGGGCGGCTTGTCAAACAG‑AAGTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTg < 2:1594327/134‑1 (MQ=255) ctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATgcgggcgg < 1:413626/57‑1 (MQ=255) ctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCC‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGACCTGCCGGTGAAACTGGTGGGGCA‑TAAACCCGGCGGGCTGGGCTACGTCGGCGGGCTGGCCGCCATGa > 1:1594327/1‑139 (MQ=255) ttCGGTTTTTTGTGCGGGGGGGGAAAATCA‑TGAATGCGGGCGGTTTGTAAAAACG‑CTTTGTTGCCCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGc < 2:1070508/127‑1 (MQ=255) tCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACg < 2:1445566/54‑1 (MQ=255) tCCGTTCTTTGTGCTGGCGGGTGAAAACA‑TGAATGCGGGGGGCTTGTCAAAACg > 1:1445566/1‑54 (MQ=255) tgtgCTGGCGGGTGAAATCA‑TGAATGCGGGGGGCTTGTCAAAACG‑CATTGTTGACCTGCCGATGAAACTGGGGGGGCA‑TAAACCAGGCGGGGTGGGCGACGTCGGGGGGCCGGGAGCTATGAttttt > 1:1070508/1‑124 (MQ=255) tCA‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCTTTGCGGGTACCGCCGCCGGTGCCGCTCTAATGGTGCCGATGATGCGCCGCGCCCACTATCCCGGCAACCCGGCGGGGGGGCTGATTGCCTCTTGCCGCAtttt > 2:1759758/1‑233 (MQ=255) cgggcggCTTGTCAAAACG‑CATTGTTAACATGCCGATGAAACTGGTGGGGCA‑TAAAACAGGCGGGCTGGGCTACGTAGGCCTGCCGGCAGCTATg > 1:1400465/1‑95 (MQ=255) | GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCA‑TGAATGCGGGCGGCTTGTCAAAACG‑CATTGTTGACCTGCCGATGAAACTGGTGGGGCA‑TAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTAT > NC_000913/3744783‑3745240 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |