Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,029,596 | A→G | V125A (GTG→GCG) | rlmI ← | 23S rRNA m(5)C1962 methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,029,596 | 0 | A | G | 97.4% | 132.1 / ‑4.8 | 39 | V125A (GTG→GCG) | rlmI | 23S rRNA m(5)C1962 methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (18/20); total (19/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.87e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01 |
CTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGA > NC_000913/1029466‑1029732 | cTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGATCATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAAt > 2:754384/1‑140 (MQ=255) ttCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCg > 2:363441/1‑140 (MQ=255) cttcttTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGaaa < 2:89281/140‑1 (MQ=255) ttttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgat < 1:919886/140‑1 (MQ=255) tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 2:1001333/140‑1 (MQ=255) tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 2:472825/140‑1 (MQ=255) tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 2:1062030/140‑1 (MQ=255) tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 2:1033911/140‑1 (MQ=255) tttttACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 1:408201/140‑1 (MQ=255) gTACCGCGACGTCGCTGCGATCGTAAATCGATCATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTa < 1:754384/140‑1 (MQ=255) tACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGTCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatcg > 2:263028/1‑133 (MQ=255) tACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGTCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatcg < 1:263028/133‑1 (MQ=255) ccGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATa < 2:1142323/140‑1 (MQ=255) cgcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc > 1:1359580/1‑129 (MQ=255) cgcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAAcgatc < 2:1359580/129‑1 (MQ=255) cgcgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATAc > 2:1263058/1‑140 (MQ=255) cgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGaaa < 1:861850/122‑1 (MQ=255) cgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGaaa > 2:861850/1‑122 (MQ=255) cgACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATAccc > 1:1143281/1‑140 (MQ=255) cTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGcc < 1:1123748/140‑1 (MQ=255) cGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAAc < 2:284794/140‑1 (MQ=255) tACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATa < 2:419818/140‑1 (MQ=255) cAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAAc < 1:17440/134‑1 (MQ=255) cAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAAc > 2:17440/1‑134 (MQ=255) cGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGt > 1:62643/1‑140 (MQ=255) ttGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGa > 2:1306984/1‑140 (MQ=255) ggCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCAt < 1:774729/140‑1 (MQ=255) tGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCgcca > 1:789023/1‑140 (MQ=255) tGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCgcca > 2:799667/1‑140 (MQ=255) gCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCgccag < 1:656591/140‑1 (MQ=255) ccccAGCACTCAGCAGTTGCAGCGCCGGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTtgt > 1:89303/1‑140 (MQ=255) ccccAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTtgt > 1:391649/1‑140 (MQ=255) ccccAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTtgt > 2:842082/1‑140 (MQ=255) ccAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACg > 1:801649/1‑128 (MQ=255) ccAGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACg < 2:801649/128‑1 (MQ=255) aGCACTCAGCAGTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCtt < 1:1306984/140‑1 (MQ=255) cagcagTTGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCa > 2:381014/1‑140 (MQ=255) ttGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAAcggc > 1:558296/1‑140 (MQ=255) ttGCAGCGCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAAcggc > 1:1240289/1‑140 (MQ=255) gcagcgCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAAcggcgg > 1:87976/1‑140 (MQ=255) agcaCCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAACCGATAGCTGTCGAGGCCATCTTTTGGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGa < 2:1143281/140‑1 (MQ=255) | CTCCATTCCTTCTTTTTTACGTACCGCGACGTCGCTGCGATCGTAAATCGAACATTCCGGGTACAGCGTTTGCAGGGCACTAATTAATGCCGCGCGCTGATATTCTGCCCCAGCACTCAGCAGTTGCAGCACCAGAAAATTACCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGA > NC_000913/1029466‑1029732 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |