Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 497,020 | 0 | A | C | 68.4% | 7.0 / 17.1 | 19 | *559E (TAA→GAA) | ybaL | Kef family K(+) transporter |
Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/13); total (5/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.16e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCG > NZ_CP009273/496891‑497153 | cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAc > 1:232930/1‑131 (MQ=255) cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAc < 2:232930/131‑1 (MQ=255) cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac > 2:34991/1‑139 (MQ=255) cACGCCGCATCGGGCATTGGGCTCTCGTTGCTTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCAACCACCTTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCcag < 2:287842/138‑1 (MQ=37) aTCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACTCCCGGCACCCCCCTTTTCCCCCGCCCCCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg < 2:381046/139‑1 (MQ=37) gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa < 1:254164/136‑1 (MQ=255) tGCCTGAGGCGCCGCTGGCGCGTCTTATCATTCCTACATATTTTTCATATTTTACATCCGGCACCCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc < 1:316568/139‑1 (MQ=255) tGCCTGAGGCGCCCCTGGCGCGTCTTACCATCCCTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc < 1:254081/139‑1 (MQ=255) cTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt < 1:195361/139‑1 (MQ=38) aTGCGCCGCTGGCGCGTCTTATCATCCCTACATATTTTTCATATTTTCCATCCGGCAACCCCCTTTTCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg < 1:34991/138‑1 (MQ=255) gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc > 2:428049/1‑138 (MQ=255) gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt > 2:347906/1‑139 (MQ=255) cTTATCATGCCTACATATTTTTCATATTTTACATCCGCCAACCCCCGTTTCCCCCGTCACCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt < 2:369781/138‑1 (MQ=255) aTGCCTACATATTTTTCATATTTTACATCCGCCAACCACCGTTTCCCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc < 1:92956/138‑1 (MQ=255) ccTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg < 2:206622/138‑1 (MQ=255) cATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg < 2:436694/138‑1 (MQ=255) cGGCAACCCCCTTTTCCCCCGCCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt < 1:178615/138‑1 (MQ=255) ggCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa > 1:354948/1‑139 (MQ=255) cGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg < 2:210662/139‑1 (MQ=255) | CGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCG > NZ_CP009273/496891‑497153 |
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: ‑ |