| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 497,013 | 0 | A | C | 52.0% | 18.2 / 30.7 | 25 | intergenic (+128/+5) | gsk/ybaL | inosine/guanosine kinase/Kef family K(+) transporter |
| Reads supporting (aligned to +/- strand): ref base A (7/5); new base C (0/13); total (7/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.65e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAA > NZ_CP009273/496878‑497147 | aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCg > 1:231997/1‑139 (MQ=255) 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 > 1:321957/1‑139 (MQ=255) ttttCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGata > 2:206918/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) | ATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTT‑ACCCCGTCACCACCTCA‑CCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAA > NZ_CP009273/496878‑497147 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |