Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 26.3% | 69.1 / 10.4 | 38 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (25/3); new base A (0/10); total (25/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.05e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824922‑1825196 | tatCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑Gagaaaag > 1:873845/1‑140 (MQ=255) ttCCATACTCGCCCCCGGCCCTTTTTTTTCTTTGTTATCGTTCAGGTAACCGTGGATGTTTATCGGCACCGGGGGTGTTTTACGTTAGTGGTGATCGGTT‑TTTTTTTACAACTGTTTGAGAA‑GAAAAAAAAAAAccgccg < 2:616980/140‑1 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGa > 2:240304/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGG‑TATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGAt > 2:837860/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACGGCCGATCCTGTc > 2:1062914/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTc > 2:655952/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCa > 1:902097/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCa > 1:840348/1‑140 (MQ=255) cccGGCCTTTTTTTTTCATTGTTATCGTTAGGGTACCGGGGGATGGTTAGCGGCAACGGGGGTAATTTATGTTAGTGGTGATCGGGT‑TTTTTTTACATCTTTTAGAGAA‑AAGAAAAGAAAGCCGCCGATCCTGTCCACCg < 2:784060/140‑1 (MQ=255) cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCg > 1:335137/1‑140 (MQ=255) gccccTTTTTTTTCTTTTTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGT‑TTTTTTTACAACTTTTTGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:958253/138‑1 (MQ=255) cccccTTTTTTTTCTTTGTTTTTGTTGAGTTAACGGGGGTTGGTTTTCGGAAGCTGGGGTG‑TTTTCCTTTTTGGGTATCAGGCTTTTTTTTACAATTTTTGGAAAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:902097/138‑1 (MQ=255) cGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 1:1171941/1‑140 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:997026/1‑140 (MQ=255) tCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa > 1:87556/1‑140 (MQ=255) gTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:172973/1‑140 (MQ=255) tATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa > 1:42573/1‑140 (MQ=255) tCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg > 2:515726/1‑140 (MQ=255) tCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAg > 2:526795/1‑140 (MQ=255) gaGGTACCGGTGGTGGTTTATCGGCACCGGGGGTGTTTTCGGTTGTTGGTGTTCTGTT‑TTTTTTTACAGTTTTTTAAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:282923/139‑1 (MQ=255) ggTTATCGGCAACTGGGGTGATTTACTTTGGTGGTGATCAGTT‑TTTTTTTACAATTGTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt < 2:336649/140‑1 (MQ=255) ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg > 1:6837/1‑140 (MQ=255) tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:780135/1‑140 (MQ=255) tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:400974/1‑140 (MQ=255) aTCGACACCTGGGGTGTTTTACTTTAGTGGTGAACAGGT‑TTTTTTCACAACTTTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGc < 2:335137/140‑1 (MQ=255) aCTGGGGTAATTTCCGTTAGTGGTGATCAGTT‑TTTTTTTAAAAGTTTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa < 2:158052/140‑1 (MQ=255) ctggggTGATTTACGTTGGTGGTGATCAGGT‑TTTTTTTACAACTGTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat < 1:493841/137‑1 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:1037685/1‑140 (MQ=255) tCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt > 1:158052/1‑140 (MQ=255) ggTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGattatt > 2:294831/1‑140 (MQ=255) tGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGGTGAAAGATTATTca > 2:30813/1‑140 (MQ=255) taTTTACGTTAGTGGTGATCAGTT‑TTTTTTTACAACTTTTAGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:42573/139‑1 (MQ=255) tACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt > 2:1120661/1‑140 (MQ=255) ctttAGGGGTGATCAGCT‑TTTTTTTACAACT‑TGTGAAAAGGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg < 2:6837/138‑1 (MQ=255) cGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTg > 1:1140056/1‑140 (MQ=255) ttAGTGGTGATCAGGT‑TTTTTTTACAATTGTTTGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:778984/140‑1 (MQ=255) ttAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:305005/140‑1 (MQ=255) gtggtgATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGc > 2:998828/1‑140 (MQ=255) ttttttACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc > 1:935352/1‑140 (MQ=255) tGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 1:1120661/140‑1 (MQ=255) aaaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc < 2:840348/136‑1 (MQ=255) gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 1:998828/140‑1 (MQ=255) gagaaAAGAAAACCGCCGATCCTGT‑CACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc > 1:78418/1‑140 (MQ=255) gaaaaaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACg < 1:240304/137‑1 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGc > 2:999897/1‑140 (MQ=255) aaaagaaaaCCGCCGATCCTGTGCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 1:146142/1‑140 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGC‑TTTTTTTACAACTGTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTT > NC_000913/1824922‑1825196 |
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: ‑ |