| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,348,816 | 0 | A | C | 60.9% | 21.4 / 25.8 | 23 | T245P (ACC→CCC) | gltB | glutamate synthase large subunit |
| Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/14); total (9/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.22e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACC > NZ_CP009273/3348744‑3348895 | aaCACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATgg > 1:99153/1‑90 (MQ=255) accccGTCCCGCGCGGGCCGCGGGCGCACCCGTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATggg < 2:123847/87‑1 (MQ=255) ccGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgc > 2:117080/1‑90 (MQ=255) ccGCACCGCGCTGGCCCCGGGCGCACCCTTCCCGCTTTGGGGGGAAGAGGGGTGAAACCACCCCCTCCCCCGGTATCCGCCAATGGgcgc < 1:418125/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:418125/1‑90 (MQ=255) cgcgcgGGCCGCTGGCCCACCCGTCCCGCTTTCTGGCGCATAACGGTGAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgcgtgcgc < 2:367688/90‑1 (MQ=255) cgcgCTGGCCGCTGGCGCACCCTTCCCGCTATCTGGCGCATAACGGTGAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgcgtgcgc < 2:179150/90‑1 (MQ=255) cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc > 1:234073/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:441809/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:235100/1‑90 (MQ=255) cgctgcccgctggcGCACCCGTTCCGCTATCTGGCGAATAACGGTGAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgcgtgcgcgt < 2:68517/90‑1 (MQ=255) tggccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACt > 2:93471/1‑90 (MQ=255) cccGTTCCGCTATCTGGCGCATAAGGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac < 2:86784/89‑1 (MQ=255) cccTATCTGGCGCATACGGGGAAAACCACCACCCCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAAACACCGCtt < 1:896/88‑1 (MQ=255) cGCTTTCTGGCGCATAAGGGTAAAACCAACCCCCTCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:401951/90‑1 (MQ=255) cGCTTTCTGGCGAAAAAGGGAAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 1:176331/90‑1 (MQ=255) cGCTATCTGGCGCATAACGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:176823/90‑1 (MQ=255) tcgggCCCAAAAGGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATccc < 2:345077/87‑1 (MQ=255) tCTGGCGCATAACGGTAAAATCAACACCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATccc < 2:182666/90‑1 (MQ=255) gcATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCt > 1:319575/1‑90 (MQ=255) gaCGGTGAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:308363/89‑1 (MQ=255) aaCGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:97327/90‑1 (MQ=255) cAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCAcc > 1:308363/1‑90 (MQ=255) | AACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACC > NZ_CP009273/3348744‑3348895 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |