Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 56.5% | 47.3 / 61.7 | 46 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (4/16); new base G (26/0); total (30/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA > NZ_CP009273/822886‑823044 | gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAg > 2:450366/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCgggaa > 1:77407/1‑87 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg > 1:762082/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCCGAAg > 1:377527/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCC‑AAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGt > 2:1095496/1‑90 (MQ=255) ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAgg > 1:525936/1‑89 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc < 2:1227113/90‑1 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTtc > 2:355753/1‑90 (MQ=255) tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGtcgtcgg > 2:1147859/1‑89 (MQ=255) ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaa < 1:1229219/90‑1 (MQ=255) ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGaa > 1:406568/1‑90 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTTGTCGAGAAACaaaaaatc > 2:836815/1‑86 (MQ=255) ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:831310/90‑1 (MQ=255) ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg < 2:137253/90‑1 (MQ=255) cATGCTGTTGATGTGCGGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAGTCGGGTTGTCGGGAAAAAGAAAGGCCGGGTCATg > 1:1091452/1‑90 (MQ=255) cATGCTGTTGATGTGCAGCC‑AAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc < 1:1095496/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg > 2:245632/1‑39 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg < 1:245632/39‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGGGGCGGGTCGTCGGGAAACAGAAAGTCCGGGTCATGc > 1:972731/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGc > 2:1360982/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGTTTGTCGGGAAAAAAAAAGTCCGGGTCCTGc > 1:147849/1‑90 (MQ=255) gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCAt > 2:618580/1‑90 (MQ=255) cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATc > 1:549406/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 1:280952/76‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc > 2:280952/1‑76 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:1390988/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:472778/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAAGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAg > 2:1291418/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCCACGCCCGAAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTTATGCATCAACg > 1:118354/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAACg > 2:1140148/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTTCTGCATCAGCg > 2:1117332/1‑90 (MQ=255) tGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCg > 2:60335/1‑90 (MQ=255) gATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGGAGGCGGGTTGTCGGGGAAAAAAAAGTGCGGGTCATGCATCAGCGa > 1:1198033/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1079133/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:1198033/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:76206/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:549755/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 1:1147859/90‑1 (MQ=255) tgCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGGACAg > 2:606726/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATTCATCAGCGAACAgg > 1:466654/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAgg > 2:601354/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAAGCCGGGTCATGCATCAGCGAACAAgcc > 1:1220696/1‑90 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAAAAgggc > 2:727782/1‑88 (MQ=255) cAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGCGGACAGgccagg > 1:1056985/1‑89 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:814233/90‑1 (MQ=255) cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc < 2:95134/90‑1 (MQ=255) tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa < 2:466654/90‑1 (MQ=255) | GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA > NZ_CP009273/822886‑823044 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |