Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,348,810 | 0 | A | C | 54.0% | 43.1 / 51.7 | 37 | T243P (ACC→CCC) | gltB | glutamate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base A (16/1); new base C (0/20); total (16/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGC > NZ_CP009273/3348722‑3348889 | gTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac > 2:77452/1‑90 (MQ=255) caGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCAcc < 2:442483/90‑1 (MQ=255) gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCgg > 1:824740/1‑90 (MQ=255) gcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTa > 1:83145/1‑90 (MQ=255) aaCACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATgg > 2:732365/1‑90 (MQ=255) ccGTCCCGCGCGGGCCGCTGGCGCACCCGTCCCGCTATCTGGCGCATAACGGTGAAATCAACCCCATCCCGGGTAACCGCCAATGGgcgc < 2:544559/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:831516/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:811885/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 1:451395/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:864060/1‑90 (MQ=255) cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc > 2:1076389/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:591047/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:544559/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:1066875/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 2:968110/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 2:433460/1‑90 (MQ=255) ggctggccgcgggcGCACCCGTTCCGCTATCGGGCGCATAAGGGTAAAATCAACCCCACCCCGGGTAACCGCCAATGGgcgcgtgcgcgt < 1:378028/89‑1 (MQ=255) gtggcggcgggcGCACCCTTCCCCCTATCGGGCCAAAAAGGGAAAAACCACCCCCCCCCCCGGAAACCGCCAATGGGCGCGTGCGCGTAc < 1:732365/89‑1 (MQ=255) cgggccgctggcGCACCCGTCCCGCTATCTGGCGCAAAACGGTGAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTAc < 2:83145/88‑1 (MQ=255) ggccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACtt > 2:409128/1‑90 (MQ=255) cgggcGCACCCGTCCCGCTATCGGGCGCAAAAGGGAAAAAACAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaaa < 1:811885/88‑1 (MQ=255) gcgcACCCGTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTc < 2:42996/90‑1 (MQ=255) gcAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCa > 2:439203/1‑90 (MQ=255) cGTTCCGGTTTCTGGCGCAAAACGGTAAAAACAACCCCATCTCCGATAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 2:201268/90‑1 (MQ=255) ccGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCt > 1:965790/1‑90 (MQ=255) cGCTTTCGGGCGCATAGCGGTAAAATCACCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:80131/90‑1 (MQ=255) cGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:68892/90‑1 (MQ=255) cGCTATCTGGCGAATAAGGGTAAAATCAACCCCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:987376/90‑1 (MQ=255) cGCTATCGGGCGCATAACGGTGAAATCAACCCCATCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:722741/90‑1 (MQ=255) tcgggCGCATAGGCGTAAAAACACCCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATccc < 1:1101253/87‑1 (MQ=255) cgcATAAGGGTAAAACCACCCCCCCCCCGGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc < 2:1173923/90‑1 (MQ=255) cgaaTAACGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc < 2:369989/87‑1 (MQ=255) aagggTAAAATCACCCCCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:954715/87‑1 (MQ=255) aagggTAAAATCACCCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:326189/87‑1 (MQ=255) aaCGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 1:409128/90‑1 (MQ=255) aaCGGTAAAACCACCCCCCCCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:15506/90‑1 (MQ=255) aaCGGAAAAACCAACCCCACCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 1:439203/90‑1 (MQ=255) aaaaaCCAACCCCACCCCGGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGAcgc < 2:607982/88‑1 (MQ=255) | GTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGC > NZ_CP009273/3348722‑3348889 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |