| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 497,013 | 0 | A | C | 52.4% | 18.0 / 28.3 | 21 | intergenic (+128/+5) | gsk/ybaL | inosine/guanosine kinase/Kef family K(+) transporter |
| Reads supporting (aligned to +/- strand): ref base A (4/6); new base C (0/11); total (4/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.51e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCA > NZ_CP009273/496881‑497138 | gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGtt > 1:187186/1‑138 (MQ=255)gTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGtt < 2:187363/138‑1 (MQ=255) aGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcc < 1:382058/138‑1 (MQ=255) aGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcc < 1:572850/138‑1 (MQ=255) gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcccc < 2:355820/139‑1 (MQ=255) aGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGt > 2:591385/1‑139 (MQ=255) cccgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCTTTTACCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCcag < 1:232359/136‑1 (MQ=38) cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCcag < 1:542835/138‑1 (MQ=255) cACGCCGCATCTGGATTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCACCTTTTCCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcag < 2:328050/138‑1 (MQ=38) cACGCCGCATCGGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCACCCCCCCCGCCGGTGGCGTTTCcag < 2:559550/138‑1 (MQ=37) gggtaTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTACCATGCCTACATATTTTTCATATTTTACACCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt < 1:405794/135‑1 (MQ=38) cagtgcctctCAATGCCTAATGCGCCGCGGGCGCTTTTTATCTTGCCTACATTTTTTTCATATTTTCCTTCCGGCACCCCCCGTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa < 2:520308/133‑1 (MQ=255) gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGTATGGTACGGGCGATTTcctg > 1:520308/1‑136 (MQ=255) ccTGATGCGACGCTGGCGCTTTTTACCCTGCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc < 2:71489/139‑1 (MQ=255) gcgcCGCTGGCGCGTCTTAACATGCCAATATATTTTTCATATTATCCACCCGCCAACCCCCTTTTACCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc < 2:544482/138‑1 (MQ=255) cgcgGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa < 2:220126/139‑1 (MQ=255) gacaggTCTTATCAGCCCCACATATTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc < 2:184372/134‑1 (MQ=255) ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacctgtc > 1:184372/1‑135 (MQ=255) aCATATTTTCCATTTTTTCCACCCGCCACCCCCCGTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt < 2:12044/137‑1 (MQ=255) cATATTTTTCATATTTTACATCCGGCAACCCCCGTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg < 1:602737/138‑1 (MQ=255) cacccGGCACCCCCCGTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa < 1:549382/135‑1 (MQ=255) | GTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCA > NZ_CP009273/496881‑497138 |
| 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: ‑ |