| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,825,061 | 0 | G | A | 32.1% | 38.0 / 11.0 | 28 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
| Reads supporting (aligned to +/- strand): ref base G (16/3); new base A (0/9); total (16/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.19e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT > NC_000913/1824922‑1825195 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGagaaaag > 1:1090816/1‑140 (MQ=255)tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGagaaaag > 1:491270/1‑140 (MQ=255) tCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGagaaaagaa > 1:20689/1‑140 (MQ=255) tATCCCATCCTCGCCCCCGCCCTTTTTTTTTCTTTTTTATCGTTCAGGTACCGGTGGATGTTTAGCGGCAACTGGGGTTATTTACGTTAGTGGTGATCAGTTTTTTTTTAC‑AACTTTTTGAGAAGAAAAAAAAAAAccgc < 2:381223/140‑1 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 2:74803/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 2:600203/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:340888/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:759376/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:739866/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGa > 1:645487/1‑140 (MQ=255) aCTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTg > 1:599669/1‑140 (MQ=255) ccGCCCCCGGCCCTTTTTTTTCTTTGTTATCGTTCAGTTAACGGTGGATGTTTATGGGCAACTGGGGTGTTTTACTTTAGGGGTAAGCAGTTTTTTTTTCC‑ACTTTTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTc < 2:1090816/139‑1 (MQ=255) cTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:887603/1‑140 (MQ=255) cccGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:1220528/1‑140 (MQ=255) ttttttCTTTGTTTTCGTTCGGTTAACGGTGATTGGTTATCGGAACCGGGGGTGATTTCCGTTGGTGTTGATCA‑GCTTTTTTTACAAACTGTTTAAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:16930/138‑1 (MQ=255) ttttttCATTGTTATCTTTCAGGTAACGGTGGTTGGTTATCGGCAACGGGGGTGATTTACGTTGGGGGTGATCAGGTTTTTTTTAC‑AACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCa < 2:759376/138‑1 (MQ=255) cAGGTAACCGTGGATGTTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAAAAAAGAAAACCGCCGATCATGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 1:32485/140‑1 (MQ=255) agggTACCGGTGGGTGGTTAGGGGAACCTGGGGTTTTTTAGGTTGGTGGTGATGAGTTTTTTTTAAA‑AGTTGTTAGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:287653/138‑1 (MQ=255) cGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGTTTTTTTTTAC‑AACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTa < 1:1136593/140‑1 (MQ=255) aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttt > 2:348567/1‑140 (MQ=255) tCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCt > 1:1116522/1‑140 (MQ=255) cGGCACCTGGGGTATTTTACGTTAGTGGTGACCAGTTTTTTTTTAC‑AACTGTTGGAGAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTg < 2:599669/140‑1 (MQ=255) taTTTACGTTAGTGGTGATCAGGTTTTTTTAAC‑AATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:491270/139‑1 (MQ=255) tACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt > 1:600226/1‑140 (MQ=255) ttgaggggtgATCAGTTTTTTTTTAC‑AATTGTTGGAAAAAAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGc < 2:1187808/134‑1 (MQ=255) ttAGTGGTGATCGGGTTTTTTTTAC‑AATTTTTTAAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:593097/140‑1 (MQ=255) ttAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa > 1:47538/1‑140 (MQ=255) ggtggtgACCAGGTTTTTTTTAC‑AACTTTTTAAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAgg < 1:365067/139‑1 (MQ=255) tggtgATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 1:593097/1‑140 (MQ=255) gCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:662371/1‑140 (MQ=255) tGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt > 2:522631/1‑140 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:1028866/1‑140 (MQ=255) gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCTGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 1:1220528/140‑1 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCt > 1:1090720/1‑140 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTAC‑AACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCT > NC_000913/1824922‑1825195 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |