Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,826,470 | (A)8→9 | coding (79/960 nt) | srlE → | sorbitol‑specific PTS enzyme IIBC1 component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,826,462 | 1 | . | A | 95.8% | 80.6 / ‑3.0 | 24 | coding (71/960 nt) | srlE | sorbitol‑specific PTS enzyme IIBC1 component |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base A (12/11); total (12/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 |
GAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑AAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGG > NC_000913/2826382‑2826535 | gAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCggg < 2:857529/81‑1 (MQ=255) aGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑a < 2:507342/81‑1 (MQ=255) aGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑a < 2:507341/81‑1 (MQ=255) catcatGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑aaaaaa < 2:507343/81‑1 (MQ=255) catcatGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑aaaaa > 2:507389/1‑80 (MQ=255) tcatGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑aaaaaaaa < 1:507344/81‑1 (MQ=255) aCGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑aaaaaaaaatcg < 2:507345/80‑5 (MQ=255) cGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTc < 2:507346/81‑1 (MQ=255) atatTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCtat < 2:507347/81‑1 (MQ=255) tatTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCtata > 1:507392/1‑81 (MQ=255) ttCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATc < 2:507348/81‑1 (MQ=255) ggATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCAcc > 2:507379/1‑81 (MQ=255) gAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCAccgccg < 1:507349/80‑1 (MQ=255) aaaaaGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGt < 1:507350/81‑1 (MQ=255) aaaaGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTa < 2:507351/81‑1 (MQ=255) aaaaGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTa > 2:507396/1‑81 (MQ=255) aaaGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTAc > 1:507382/1‑81 (MQ=255) ggTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGc < 1:507352/81‑1 (MQ=255) gTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCg < 2:507353/81‑1 (MQ=255) ggCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGAt > 2:507385/1‑81 (MQ=255) gCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGAtt > 2:507399/1‑81 (MQ=255) ggggCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGAttgtt > 2:507391/1‑81 (MQ=255) cggcCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGAc > 1:507401/1‑80 (MQ=255) ggcCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACaa < 1:507354/81‑1 (MQ=255) ccGCTTGAGCTGAAAGCCACGCCGGG‑AAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTg < 1:507355/81‑1 (MQ=255) aGCTGAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcaca < 1:507356/80‑1 (MQ=255) gAAAGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTa > 2:507400/1‑81 (MQ=255) aaaGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTAc > 2:507383/1‑81 (MQ=255) aaaGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTAc < 2:507357/81‑1 (MQ=255) aaGCCACGCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTAct > 2:507402/1‑81 (MQ=255) gCCGGGAAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActggctgg > 2:507386/1‑80 (MQ=255) | GAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGG‑AAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGG > NC_000913/2826382‑2826535 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |