Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,900,885 | A→G | S145P (TCC→CCC) | rlmA ← | 23S rRNA (guanine(745)‑N(1))‑methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,900,885 | 0 | A | G | 80.8% | 46.1 / 9.2 | 26 | S145P (TCC→CCC) | rlmA | 23S rRNA (guanine(745)‑N(1))‑methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (4/1); new base G (14/7); total (18/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01 |
CCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAA > NZ_CP009273/1900799‑1900968 | ccAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGgaaa < 1:751282/90‑1 (MQ=255) gccgccGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACgg > 2:318454/1‑90 (MQ=255) ggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACgg < 2:338947/90‑1 (MQ=255) ggCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACgg < 1:318454/90‑1 (MQ=255) cACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCt > 2:200187/1‑90 (MQ=255) tAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGAcaca < 1:93243/90‑1 (MQ=255) aTTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACaaa > 1:116851/1‑90 (MQ=255) ctGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGa > 2:735839/1‑90 (MQ=255) ttACACGGCGTGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGc > 2:742679/1‑90 (MQ=255) gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGc > 1:265403/1‑90 (MQ=255) gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGc > 2:668585/1‑90 (MQ=255) gcgcGTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGc > 2:529093/1‑90 (MQ=255) cgcgTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCg > 2:184037/1‑90 (MQ=255) cgTAAATACGTATTATGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGtt > 1:599812/1‑90 (MQ=255) tattatGGCGTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTcgccgccgc > 1:342896/1‑90 (MQ=255) gTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc > 1:633518/1‑90 (MQ=255) gTCCATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc > 1:278600/1‑90 (MQ=255) gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc < 1:405604/90‑1 (MQ=255) gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc > 2:205622/1‑90 (MQ=255) gTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc > 1:281207/1‑90 (MQ=255) gTCCATACTGGTATCAGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGc > 1:777011/1‑90 (MQ=255) ccATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTa < 2:265403/90‑1 (MQ=255) ccATACTGGTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTa < 1:455033/90‑1 (MQ=255) ggTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGa > 2:5144/1‑90 (MQ=255) gTATCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaa < 2:599812/90‑1 (MQ=255) aTCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaaa > 2:181289/1‑89 (MQ=255) aTCGGGAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGaaa < 1:181289/89‑1 (MQ=255) | CCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAA > NZ_CP009273/1900799‑1900968 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |