| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,026,694 | 0 | A | G | 52.5% | ‑4.0 / 92.3 | 61 | V63A (GTC→GCC) | drpB | cell division protein DrpB |
| Reads supporting (aligned to +/- strand): ref base A (18/11); new base G (16/16); total (34/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.41e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATG > NZ_CP009273/2026560‑2026823 | tACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAg < 2:67120/139‑1 (MQ=255) cTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGcccc > 1:240180/1‑139 (MQ=255) gCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCgg > 1:399649/1‑138 (MQ=255) ccACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTa < 2:169068/139‑1 (MQ=255) ttCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCa < 1:991284/139‑1 (MQ=255) tAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATc > 1:208307/1‑139 (MQ=255) cgaCGTGGGGGACGAGGGCGGGTCTACCCGGCAATGCCGCCTAACAATGCCACGACCGGAGCCATAAATAAAAATCCGACCAGTCCCCCCAACCACTTTCCCGCGGAAGCGCCTACATCGCTTTCACATCCCGACACgg > 2:567986/3‑139 (MQ=255) cGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACggg < 1:549757/137‑1 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAg > 1:908256/1‑78 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAg < 2:908256/78‑1 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAg < 2:280502/78‑1 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAg > 1:280502/1‑78 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCa > 1:308116/1‑139 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGt < 2:514185/102‑1 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGt > 1:514185/1‑102 (MQ=255) tGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCa > 2:1058944/1‑139 (MQ=255) tGCGGAACGAGGTCTGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCAGACCAGTGCCCCTAACTACTTTCCCGCGGTAGAGGCTAAATAGCTTTCAAATCCCGACACGGGGGCa > 2:959796/1‑139 (MQ=255) ggTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATcc > 2:72604/1‑139 (MQ=255) ggTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATcc < 1:418166/139‑1 (MQ=255) ggTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATcc > 1:51252/1‑139 (MQ=255) ggTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCt > 2:299093/1‑138 (MQ=255) ggTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCt > 2:302712/1‑138 (MQ=255) tgtACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCtcat < 1:1027447/139‑1 (MQ=255) cAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGAt > 1:1134863/1‑139 (MQ=255) tgctgcCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGa > 2:995447/1‑45 (MQ=255) tgctgcCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGa < 1:995447/45‑1 (MQ=255) tgctgcCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGt > 1:1277489/1‑138 (MQ=255) ccTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCGCTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAAcg < 2:308116/137‑1 (MQ=255) ccTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCg < 1:728754/81‑1 (MQ=255) ccTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCg > 2:728754/1‑81 (MQ=255) ccTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAAcgcg > 1:813805/1‑139 (MQ=255) aCAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCa < 2:51252/138‑1 (MQ=255) cAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCAtt > 1:1103484/1‑139 (MQ=255) ccccGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTccc > 1:303080/1‑53 (MQ=255) ccccGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTccc < 2:303080/53‑1 (MQ=255) ccccGACCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGcc > 2:507459/1‑137 (MQ=255) ccccGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGcc > 1:434337/1‑137 (MQ=255) ccccGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGcc > 1:688474/1‑137 (MQ=255) aCCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCtcatca > 2:70282/1‑109 (MQ=255) aCCAGAGCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCtcatca < 1:70282/109‑1 (MQ=255) aCCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACg > 2:211155/1‑138 (MQ=255) ccAGAGCCATAAATAAAAATCCGGCCAGTGCCTCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGa > 1:1108643/1‑138 (MQ=255) cAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGaaa < 2:852916/139‑1 (MQ=255) gagCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGaaa > 1:1127095/1‑137 (MQ=255) agCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCtca < 1:406516/101‑1 (MQ=255) agCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCtca > 2:406516/1‑101 (MQ=255) agCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAg > 2:1133229/1‑139 (MQ=255) agCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAg > 2:287840/1‑139 (MQ=255) gCCATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGc < 1:507459/139‑1 (MQ=255) ccATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGcc < 1:735331/139‑1 (MQ=255) ccATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGcc < 2:896397/139‑1 (MQ=255) ccATAAATAAAAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGcc < 1:72604/139‑1 (MQ=255) taaataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGcc < 2:340238/136‑1 (MQ=255) taaataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGcc > 1:340238/1‑136 (MQ=255) taaataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca > 1:256089/1‑139 (MQ=255) taaataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca > 1:994338/1‑139 (MQ=255) taaataaaAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca < 2:400230/139‑1 (MQ=255) taaataaaAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca < 1:1287340/139‑1 (MQ=255) taaataaaAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca < 2:688474/139‑1 (MQ=255) taaataaaAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcaca > 1:400230/1‑139 (MQ=255) aataaaAATCCGGCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCacac < 2:531480/138‑1 (MQ=255) aataaaAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACAcc < 1:174128/139‑1 (MQ=255) aaaaaTCCGGCCACTGTCCCGAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGTGGCATCGGTAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGTCGAAATAGCCACACCATg < 2:6359/139‑1 (MQ=255) aCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcac < 1:839484/124‑1 (MQ=255) aCCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCcac > 2:839484/1‑124 (MQ=255) | TACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACCAGAGCCATAAATAAAAATCCGACCAGTGCCCCTAACCACTTTCCCGCGGTAGAGCCTAAATCGCTTTCAAATCCCGACACGGGGGCATCGGGAATCCGGCTCATCACCCAGATGTAACGCGCCATTGCCCAGACGAAATAGCCACACCATG > NZ_CP009273/2026560‑2026823 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |