Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 645,916 | 0 | T | G | 52.2% | 26.4 / 34.7 | 23 | T11P (ACT→CCT) | citE | citrate (pro‑3S)‑lyase subunit beta |
Reads supporting (aligned to +/- strand): ref base T (2/9); new base G (12/0); total (14/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.73e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATT > NZ_CP009273/645780‑646050 | cgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGttt < 2:169606/139‑1 (MQ=255) aCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTTGGTTTTGCCGCGag > 2:70486/1‑138 (MQ=255) gCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTTTTGCAGCGAAGc > 1:365663/1‑139 (MQ=255) tGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTCGTTGTTGCAGCGAAGCGGa > 1:353229/1‑138 (MQ=255) tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTGTTGCAGCGTAACGGAAACCATTGGc > 2:216772/1‑138 (MQ=255) tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTGTAGCGAAGAGGGAAACAtttgg > 2:269087/1‑135 (MQ=255) tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGGAGCGGAAATCATTGGc > 2:353356/1‑138 (MQ=255) ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGGGGTTTTCGTTGTTGCCGCGGAAAGGAAAACATTGGCAAt > 2:131244/1‑139 (MQ=255) ttCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTGTTGCAGCGAAGCGGGAAACATTGGCCAt > 2:186644/1‑139 (MQ=255) cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAt < 2:359743/137‑1 (MQ=255) cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAAt < 1:353356/137‑1 (MQ=255) gTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGTTTTACGTTGTTTCCGCGAAGCGGAAATCATTGGCAATCCTccc > 1:318009/1‑139 (MQ=255) aCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGc < 1:85794/139‑1 (MQ=255) gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGgcg < 2:269123/138‑1 (MQ=255) gCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa < 1:269087/138‑1 (MQ=255) cAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATg < 2:370209/139‑1 (MQ=255) aGGCACAAACAACATGCTGCGGCGGGTGCGCGGTTTTTGTTGTTGCCGCGAAGCGGAAATCCTTGGCAATCCTCCCACGGGCGA‑GGCGGGGTACCGCTGGGGCGTGGCCGCAGGGCTTCCCGTCGGGGACGGaaaaag > 2:32569/1‑136 (MQ=255) aGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTCGTTGTTGCAGCGAAGCGGGAAACATTTGCAATCCTCCCCTTGCCGAGGCCGGGA‑ACCCCTGGCGCGGGCCCGCCGGGGTTCCCCTCGGGCACCTAAAATg > 2:387880/1‑138 (MQ=255) acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTc < 2:318009/139‑1 (MQ=255) aacaacATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAg < 1:70486/139‑1 (MQ=255) gcggcgGGTGCGAGGTTTACGTTGTTGCAGCGAAGCGGAAATAATTGGCAATCCCCCCCTGGCAGA‑GCCGGGGGACCGCTGGCGCGTGCCAGACGGGGTTCCCGTCGTGCCCGGAAAATGGAGGCCAGGGCGCCctttt > 2:341951/1‑137 (MQ=255) cggGGGCGAGTTTTTCGTTGTTGCAGCGAAGCGGAAATAATTGGCAATCCTCCCCTGGCAGA‑GGCGGGGTACCGCTGGCGCGGGCCAGCAGGGCTTCCAGTCGTGCCCGGAAAATGCAGTCCAGGGCGCCTTTGTCATc > 2:335097/1‑139 (MQ=255) cggGGGCGAGGTTTTCGTTGTTGCAACGGAGAGGAAAACATTGGGAATCCTCCCAAGGCAGA‑GCCGGGGGACCCCCGGCGCGGGCCAGCAGGGCTTCCAGTCGGGCGCGTAAAATGCAGTCCAGTGCGCCTTTGGCATc > 2:14631/1‑139 (MQ=255) cGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACAtt < 1:387880/139‑1 (MQ=255) | CGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGA‑GCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATT > NZ_CP009273/645780‑646050 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |