| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 497,020 | 0 | A | C | 60.9% | 13.2 / 26.9 | 24 | *559E (TAA→GAA) | ybaL | Kef family K(+) transporter |
| Reads supporting (aligned to +/- strand): ref base A (6/3); new base C (0/14); total (6/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.32e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAA > NZ_CP009273/496892‑497147 | gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCTTTT‑ACCCCGTcac < 1:359901/138‑1 (MQ=255) agcccAGAAAGGGCTTCCACGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTCCATCCGGCACCCCCCGTTT‑CCCCCGCCCCCCCCTCC‑Cccg < 2:555525/136‑1 (MQ=25) cccgccccATCTGGCTTTTGGCTCTCGTTGCCGGATGCGCCGCTGGCCCGTCTTATCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTT‑CCCCCTTCACCCCCTCC‑CCCGCCGGTGGCGTTTCcagc < 2:63347/137‑1 (MQ=255) cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCTTTT‑CCCCCGTCACCCCCTCC‑CCCGCCGGTGGCGTTTCcag < 1:202856/138‑1 (MQ=38) aTCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAg > 1:238141/1‑139 (MQ=255) ttGGCATTTGGCTCTCGTTGCCTGATGCGCCGCTGGCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTCCATCCGGCAACCCCCGTTT‑CCCCCGTCCCCCCCTCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAg < 2:199788/136‑1 (MQ=255) ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt > 2:180725/1‑139 (MQ=255) ccTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTT‑CCCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc < 1:243928/139‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑CCCCCGTCACCACCCCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc < 1:132450/139‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGAAACCCCCGTTT‑ACCCCGTCCCCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc < 1:628787/139‑1 (MQ=255) cGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa > 2:591098/1‑139 (MQ=255) tCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt > 1:349102/1‑139 (MQ=255) cccTCCATTTTTTTCATTTTTTAC‑CCCCGCAACCCCCCTTTCCCCCCTTCCCCCCCCCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg < 2:321957/138‑1 (MQ=255) ccTACATTTTTTTCATATTTTCCCCCCGGCACCCCCCTTTT‑CCCCCGTCCCCCCCTCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg < 2:480666/138‑1 (MQ=255) ccTACATTTTTTTCATATTTTCCATCCGACCACCCCCTTTC‑CCCCCGTCCCCCCCCCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg < 2:407905/138‑1 (MQ=255) ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑CCCCCGTCCCCCCCTCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg < 2:231997/138‑1 (MQ=255) catattttacatattttaCACCCGCCACCCCCCGTTT‑CCCCCGTCACCCCCCCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg < 2:349102/139‑1 (MQ=255) cATTTTTTTCATATTTTACACCCGGCACCCCCCTTTT‑CCCCCGCCCCCCCCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg < 2:529368/139‑1 (MQ=255) aatattttccatattttccCTCCGGCAACCACCTTTT‑CCCCCGTCCCCCCCTCT‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg < 2:502094/138‑1 (MQ=255) ttttCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata > 2:206918/1‑139 (MQ=255) ttttCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata > 1:321957/1‑139 (MQ=255) ttccatttttTCCTCCCGCCACCCCCCTTT‑TCCCCCTCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata < 2:469407/132‑1 (MQ=255) atatTTTACACCCGGCACCCCCCTTTT‑CCCCCGTCCCCCCCCCC‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg < 1:180725/138‑1 (MQ=255) aaCCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATaa < 1:206918/139‑1 (MQ=255) | GTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAA > NZ_CP009273/496892‑497147 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |