Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,348,816 | 0 | A | C | 59.1% | 17.0 / 25.1 | 22 | T245P (ACC→CCC) | gltB | glutamate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/13); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC > NZ_CP009273/3348730‑3348884 | aGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCg > 1:39135/1‑90 (MQ=255) gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCgg > 2:497257/1‑90 (MQ=255) cacCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgc > 1:337796/1‑90 (MQ=255) ccGTCCCGCGCTGGCCGCGGGCCCACCCGTCCCCCTTTCGGGCGCAAACGGGTAAAACCACCCCCATCCCCGGTAACCGCCAATGGgcgc < 2:455078/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 1:569633/1‑90 (MQ=255) cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc > 1:455078/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:492460/1‑90 (MQ=255) cgcgggccgcgggcGCAACCGTCCCGCTATCGGGCGCATAACGGTAAAATCAACCCCATCCCGGGTAACCGCCAATGGgcgcgtgcgcgt < 1:544713/90‑1 (MQ=255) cgcgggccgcggccGCACCCGTCCCGTTTTCGGGCGCATAAGGGTAAAATCACCCCCACCCCCGGTAACCGCCAATGGgcgcgtgcgcgt < 2:238123/90‑1 (MQ=255) tgggccgcgggcGCACCCTTCCCGCTGTTTGGGGAATAGGGGTAAAATCAACCCCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTAc < 1:470237/88‑1 (MQ=255) cgggccgcgggcCCACCCTTCCCGCTTTCGGGCGCAAAAGGGTAAAATCACCCCCCTCCCCGGTAACCGCCAATGGGCGCGTGCGCGTAc < 2:163577/88‑1 (MQ=255) cgcgggcCCACCCGTCCCGCTTTCGGGCGAATAAGGGTAAAATAAACCCCATCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTtata < 2:243058/90‑1 (MQ=255) ctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaaa > 2:544713/1‑90 (MQ=255) cGTTCCGCTATGGGGCGCATAGGGGTAAAATAACCACCTCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 2:309244/90‑1 (MQ=255) cGTCCCGCTATCGGGCGCATAGGGGAAAAACAAACCCCCCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 2:588988/90‑1 (MQ=255) ttCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCg > 1:29705/1‑90 (MQ=255) cGCTTTCGGGCGAAAAACGGTAAAAACACCCCCACCCCCGTTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:47396/90‑1 (MQ=255) cGCTGTCTGGCGCATAAGGGTAAATTCAACCCCCTCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:41643/90‑1 (MQ=255) cGCTGTCGGCCGCAAAACGGAAAAATCACCCCCCTCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 1:497257/90‑1 (MQ=255) cTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCt > 1:431296/1‑90 (MQ=255) cgcATAACGGAAAAACCACCCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc < 1:269485/90‑1 (MQ=255) aagggTAAAACCAACACCACCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:569633/87‑1 (MQ=255) | AGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC > NZ_CP009273/3348730‑3348884 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |