| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,743,669 | 0 | G | A | 39.4% | 33.7 / 18.5 | 33 | E313K (GAA→AAA) | yiaK | 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent |
| Reads supporting (aligned to +/- strand): ref base G (19/1); new base A (0/13); total (19/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.44e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATT > NC_000913/3743523‑3743812 | ttGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGaaa > 2:319517/1‑150 (MQ=255)ttGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCAGGAATTTACTACCCTGCTGGCCGaaa > 1:132371/1‑150 (MQ=255) aaGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATc > 2:618150/1‑150 (MQ=255) tATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGtt > 1:84388/1‑151 (MQ=255) aTCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtg > 2:536146/1‑151 (MQ=255) aTCGACGGCCCCCCCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACAAAAATCAGCCCATTCGCTT‑CCCCGGCCATAAATTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTtg < 2:80115/151‑1 (MQ=255) gctCCCCCCCGCGAGGCCAGGCTGCAACGCTTCATGGTTTACGTTCCTAGGCCCGGGCGTGCGAACAAAAATCGGCCCATCCGTTT‑CCCCGGCAATAATTTTACTCCCCTGCTGGCCAAAATCCGCCGTAACGGCATCACTGTtgatga < 1:77059/147‑1 (MQ=255) cccACCCGCGATCCCAACTTCAAATGAATCGTGGATTCCGTTGCGTGGCCGGGGCGTGCTAACAAAATTCAGCCCTTTCGCTC‑ACCCGCCAAATATTTTACTCCCCTTCTGGCCAAAAACCGCCGTAACGGCATCACTGTtgatga < 2:632011/146‑1 (MQ=255) cACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg > 2:2548/1‑151 (MQ=255) aaGCTCAAACCCATCATGGTTTCCTTTTCTAGTGCCGGGCGTGTGGACAAAAAGCAGGCCATCCTCTC‑CCCCGGCTTTGATTTTCCTCCCCGGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGcc < 2:267531/144‑1 (MQ=255) gCTGCAACGCAGCATGGATTACGTTACGAGCGCCGGGCGTGCGGAAAAAAACCGGGCCTTTCTTTC‑CCCCGGCAATGATTTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCa < 2:204680/150‑1 (MQ=255) cggggggggCTGACAAAAATCAGCCCATCCCCTC‑CCCCGGCTATAATTTTATCCCCCTGCTGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTgg < 1:532946/142‑1 (MQ=255) cggccGTGCTGAGAAAAAACA‑GCCCTTCCCTTCCCCCGCCCATAAATTTCCTCCCCGGCCGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAttttt < 2:496946/143‑1 (MQ=255) gCGTGCTGAAAAAAACCGGCCCTTTCTCTT‑CCCCGCCAATAATTTTCCTCCCCGGTTGCCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACa < 2:592337/149‑1 (MQ=255) cTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcgcg > 1:115903/1‑151 (MQ=255) aaaaaaCCAGCCCATCCGCTT‑CCCCGGCCTGTATTTTCCTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAg < 1:202532/148‑1 (MQ=255) aaaaaCCAGGCCTTTCGTTT‑CCCCGGCCATAATTTTACTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCa < 2:311887/147‑1 (MQ=255) aaaTCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGccc > 1:534304/1‑150 (MQ=255) aaaTCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGccc > 2:440129/1‑150 (MQ=255) aaaCCAGCCCTTCCGCTC‑CCCCGCCAATAAATTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCa < 2:425057/151‑1 (MQ=255) aTCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAAt > 2:463049/1‑151 (MQ=255) aTTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCg > 1:451900/1‑151 (MQ=255) tCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGtt > 2:98534/1‑151 (MQ=255) cccGCCAATGATTTTCCTCCCCGGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt < 1:115155/145‑1 (MQ=255) ccgccaataattttACCCCCCTGTGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgcc < 1:299641/140‑1 (MQ=255) ccGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgcc > 2:283825/1‑151 (MQ=255) ccGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgcc > 2:8803/1‑151 (MQ=255) tGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTg > 1:289565/1‑151 (MQ=255) tttACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGaaa > 1:650995/1‑150 (MQ=255) ccccTGCGGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCt < 2:506780/150‑1 (MQ=255) ctgctgGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCg > 2:597520/1‑150 (MQ=255) gctgGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGAtt > 2:601940/1‑151 (MQ=255) gctgGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCCCGAtt > 1:401349/1‑151 (MQ=255) | TTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTT‑ACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCTCGATT > NC_000913/3743523‑3743812 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |