| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,348,816 | 0 | A | C | 55.6% | 20.6 / 24.6 | 18 | T245P (ACC→CCC) | gltB | glutamate synthase large subunit |
| Reads supporting (aligned to +/- strand): ref base A (8/0); new base C (0/10); total (8/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.29e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTC > NZ_CP009273/3348731‑3348902 | gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGTGCATAACGGTGAAATCAACACCATCACCgg > 2:172408/1‑90 (MQ=255) ccGTCCCGCGCTGCCCGCGGGCCCACCCTTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgc < 2:513698/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:422906/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:337319/1‑90 (MQ=255) cgcgggccgcgggcGCACCCGTTCCGCTATGTGGCTAATAACGGGAAAATCAACCCCACCCCCGGAAACCGCCAATGGgcgcgtgcgcgt < 2:376959/90‑1 (MQ=255) cccGTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac < 2:127472/89‑1 (MQ=255) cccGTCCCGCTATCGGGCGAATAGGGGAAAAAACACCCCCACCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac < 1:435190/89‑1 (MQ=255) cttcccGCTATCTGGCGCATAGCGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 2:252491/86‑1 (MQ=255) cGTCCCTCTATCTGCCGCAAAAGGGTAAAATAAACCCCCCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 1:436017/90‑1 (MQ=255) cGCTATCTGGCGCATAAGGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:369098/90‑1 (MQ=255) cGCTATCGGGCGCATAAGGGTAAAACCACCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:399561/90‑1 (MQ=255) tATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATc > 1:136614/1‑90 (MQ=255) aTCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATcc > 1:145130/1‑90 (MQ=255) cgcATAACGGTAAAACCACCCCCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc < 2:145130/90‑1 (MQ=255) aaCGGTAAAATAAACACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:115202/90‑1 (MQ=255) gTGAAATCAAAACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacgacg > 1:485003/1‑90 (MQ=255) gAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGAcgcc > 1:132830/1‑90 (MQ=255) aTCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTc > 1:205027/1‑90 (MQ=255) | GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTC > NZ_CP009273/3348731‑3348902 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |