Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,348,801 | 0 | G | A | 55.6% | 6.0 / 22.2 | 18 | E240K (GAA→AAA) | gltB | glutamate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base G (7/1); new base A (0/10); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCC > NZ_CP009273/3348731‑3348890 | gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGTGCATAACGGTGAAATCAACACCATCACCgg > 2:172408/1‑90 (MQ=255) ccGTCCCGCGCTGCCCGCGGGCCCACCCTTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgc < 2:513698/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac < 2:473374/63‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac > 1:473374/1‑63 (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) | GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCC > NZ_CP009273/3348731‑3348890 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |