Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,845,829 | 0 | T | G | 34.6% | 15.7 / 22.1 | 27 | S26R (AGT→AGG) | pfeS‑II | histidine kinase |
Reads supporting (aligned to +/- strand): ref base T (9/8); new base G (5/4); total (15/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCA > NC_002947/1845686‑1845978 | gAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCCACCGCCATTCGCTGTTCTCGCAGCTGTCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATAcc > 2:398550/1‑150 (MQ=255) aCGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTg < 1:742847/150‑1 (MQ=255) cGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTga > 1:216918/1‑149 (MQ=255) aTCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCga < 2:216918/150‑2 (MQ=255) gCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGcaca > 2:777093/1‑150 (MQ=255) cACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt < 2:956741/150‑1 (MQ=255) ccGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCgatc < 1:210840/145‑4 (MQ=255) ccGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCgatc > 2:210840/1‑142 (MQ=255) ttACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCtt < 2:494865/150‑1 (MQ=255) cTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGaa > 1:477234/1‑129 (MQ=255) cTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGaa < 2:477234/129‑1 (MQ=255) tcaaagacAGCGAGGCCGAGGGAACCACGGTCTACCGCCATTTTCTGTGTTTGAAGCTGGCCATCTTGCTGGTGGGCTTATGACTGCTGATGAACGACCTCAGCGATACCCGGTGCCGGCACAGGGAAACCCAGAATGCCTTCCTCTCgg > 2:735728/8‑150 (MQ=255) cAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGa < 2:873205/150‑1 (MQ=255) ccGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGTCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgc > 1:472710/1‑150 (MQ=255) cGCCATTCGCTGTTCTAGAAGCTGGCCATCTCGCGGGTGGGCTTCTGCCTGCTGATGATCAGCCGCAGGTGTACCTGGGGCCGTCCCATGGGAACCCAGAATGCCTACCTCTCGGAGCCGACACGCGAGACCTTGCGGGGTTATgccgc‑c > 2:349343/1‑150 (MQ=255) gCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGccgc‑cg < 1:951804/150‑1 (MQ=255) gCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGa > 1:678339/1‑150 (MQ=255) tGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCATAAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAAc > 2:252994/1‑150 (MQ=255) aTCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTggc > 2:798527/1‑150 (MQ=255) tctcGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTggcg > 1:585145/1‑150 (MQ=255) cTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCga < 1:21194/51‑2 (MQ=17) cTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCga > 2:21194/1‑50 (MQ=17) ttCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTg < 1:798527/150‑1 (MQ=255) gatgatCGGCCTCAGTTATACCTGGGGCCGGAACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCAGGGGCGTTCCAGC‑CACCGACCACGCGCGAAGCTTAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtgggt > 1:235180/1‑150 (MQ=255) ggCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtg < 2:1015382/140‑1 (MQ=255) ggCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtg > 1:1015382/1‑140 (MQ=255) ccTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCGCTGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAt < 1:63653/150‑1 (MQ=255) ttATACCTGGGGCCGGCACATGGAAACCCAGAATGc < 1:212107/36‑1 (MQ=255) ttATACCTGGGGCCGGCACATGGAAACCCAGAATGc > 2:212107/1‑36 (MQ=255) ttATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGACTGGACCAGTGGGTGGCGGCCATGCATCa > 1:473464/1‑150 (MQ=255) | GAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGC‑CGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCA > NC_002947/1845686‑1845978 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |