| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 803,869 | 0 | G | T | 36.4% | 42.6 / 11.0 | 33 | V123L (GTG→TTG) | ybhJ | putative hydratase YbhJ |
| Reads supporting (aligned to +/- strand): ref base G (5/16); new base T (12/0); total (17/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.26e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGC > NC_000913/803623‑804115 | aTTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgtg < 1:516833/249‑1 (MQ=255) cGCATAATACGTCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCatat < 1:78207/250‑1 (MQ=255) tCCGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTGTCTCCCCATCTTGTGGtcctccc > 2:496716/1‑250 (MQ=255) cGGAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccat‑c < 1:256855/251‑1 (MQ=255) gAAATATGGATAAAATTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGAGTTTTTGTGCCttcccccattgcgggttttcatca > 1:104132/1‑227 (MQ=255) aaaTATGGTAAAAATTAAAATTAAGTTTGAATCATTAGCCTCTCTCGATGTTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATGAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAat < 2:544806/251‑1 (MQ=255) ttAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCTTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatcc < 2:92082/231‑1 (MQ=255) ttAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGGCTCCCCATATTGCGGGcatcccggtcgggagaacacaagac > 1:92082/1‑231 (MQ=255) tGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAggcggc < 1:477579/250‑1 (MQ=255) tGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTTCCTTCACATATTGCGGGCATCCAT‑CCACATATGCGGTG‑GATGATggacggcgggg > 1:323392/1‑242 (MQ=255) tGATTCATTAGCCTCTCACGATATTACCGTTGTTGGTATTGTACAGTCCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAggcggc < 2:153836/250‑1 (MQ=255) ctcccGATATTTCCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATcc < 2:333847/245‑1 (MQ=255) atTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTATTGCGGGTATACAT‑CAATATATTCG‑TGAGATGATGGCAAGGGGCGGGAAAAAGATCCTCGGGTCAgc > 1:249741/1‑250 (MQ=255) atTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAg < 2:545983/250‑1 (MQ=255) ttACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTATTTCGGTCCTCCAA‑TAATATATGCG‑TGTGATGATGGCAGGGGGCGGGCAAATGTTCCTCGGGGCAGa > 1:365451/1‑250 (MQ=255) gTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCATATTTCGGGCATCCAT‑CAATATATGGG‑GGAGGTGATGGCAGGGGGGGGGAAAAAGATCCCCGGGGCAGACAGCCACA‑cccct > 1:62038/1‑249 (MQ=255) aCCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCAtt < 1:300013/251‑1 (MQ=255) ggTATGGACCTTTTCCCGCTGGCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGt < 2:187589/251‑1 (MQ=255) aaCGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGGGCCTCCGCATTTTGCGGTCATCCATTCAA‑ATATTCG‑TGAGATGGAGGCAGGCGGCGGGAAAATGATCCTCGGGTCAGACAACCACA‑CCCGGTACGGTGGGTTAGGGaaaaagggaggcggtg > 1:193370/1‑234 (MQ=255) cccGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGgcggcg < 2:128424/251‑1 (MQ=255) gTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTgg < 1:128453/251‑1 (MQ=255) gCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCCTATTTCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGGCAGACAGCCCCA‑CCCGTTACGGTGCATTTGGGGCAATGGCAGTCGGGGGGGGCGGCGGGGAGTTGGGaaaa > 1:440130/1‑250 (MQ=255) cAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCATATTTCGGTCCTCCAA‑CAAAATATGTG‑TGAGATGATGGCAGGGGGGGGGCAAATGATTCTCGGGTCAGAAAACCACAACCCGT‑ACGGGGCATTAAGGGCAATGGCAGGCGGCGAGGGCGGCGGGGGGTTGGGaaaaaagatgc > 1:169667/1‑244 (MQ=255) tctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTCCGCATATTGCGGTCATCCAT‑CCAAATATGCG‑TGAGATGATGGGCGGGGGGGGCCAAATGATCCCCCGGTCAGACAGCCCCA‑CCCGGTACGGGGGCTTAGGGACAATTGCAGTCGGTGGGGGGGGCGGGGGGTTTGGAAAAAAGCTGGTTAATGAAACCTggggt > 1:121349/1‑249 (MQ=255) ctGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAAAGATCCTCGGGGCAGACAGCCACA‑CCCCTTACGGGGCATTAGGGACAAAGGCAGTCGGGGGGGGGGGCGGGGAGGTGGTAAAACAGCCGCTTTATGACACCTggggt > 1:513719/1‑248 (MQ=255) cggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGGCTTCTCCTATTCCGGGCATCCCT‑CAATATATTCG‑GGAGATGATGGCACGATGTGGCCAAATGATCCTCCGGTCAGACAGCCACA‑CCCGTTACGGTGCGTTAAGGGCAAATGCAGTCGGATAGGGGGGCGGTGAGTTGGTAAAACCACTGCGTAATTACCACTGGGGTATGGGCTATCCgg > 1:397943/1‑251 (MQ=255) cACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCAAGCGTTATGGCGGGCTTTTTTTGCCTCCCCATTTTGCCGGCATCCCT‑TAATATATGCG‑TGAGATGAAGGAAGGCGGCGGCCAAATGAACCTCCGGTAAAAAAGCCACA‑CCCGTTACGGGGTATTAGGGGCAATGGCAGGCGGGGAGGGCGGGCGTGAGGTGGTAAAACAGGCGCTTAAAAAAACCCGGGGAATCGACTATCCggggg > 1:207912/1‑249 (MQ=255) gaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCt < 1:163661/251‑1 (MQ=255) gTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCgg < 1:489947/251‑1 (MQ=255) ggcggcCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCACATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGGAGGCGGCGGCAAAATGATCCTCCGGTCAAACAGCCACA‑CCCGTTTCGGGGGCTTTGGGACAATGGCAATCGGGGGGGGCGGGGGGGGGGTGGGTAAACAACTGCTTAATGAAACCCGGGGTATCGACCATCCCGGCGGGGGTGCGGGGCCTCTGCACGGGAAAACC‑GCGCCgg > 2:524743/1‑250 (MQ=255) gcggcCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTa < 2:62038/251‑1 (MQ=255) tATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTATGTGGGGCcac < 1:208933/250‑1 (MQ=255) gCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGAt > 2:287041/1‑251 (MQ=255) ttttGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGgcgc > 1:8461/1‑251 (MQ=255) | ATTCTCTCTTCGCATAATACGTCCGGAAATATGGATAAACTTAAAATTAAGTTTGATTCATTAGCCTCTCACGATATTACCTTTGTTGGTATTGTACAGACCGCTAAAGCGTCCGGTATGGAACGTTTCCCGCTGCCGTATGTGCTGACCAACTGCCATAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCAT‑CAATATATGCG‑TGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA‑CCCGTTACGGTGCATTAGGGACAATGGCAGTCGGTGAGGGCGGCGGTGAGTTGGTAAAACAGCTGCTTAATGACACCTGGGATATCGACTATCCGGGCGTGGTTGCGGTGCATCTG‑ACCGGAAAACCTGCGCCGTATGTGGGGCCACAGGATGTGGCGC > NC_000913/803623‑804115 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |