Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,920,205 | (T)7→8 | intergenic (+62/‑18) | yebT → / → rsmF | MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,920,198 | 1 | . | T | 98.0% | 162.0 / ‑3.7 | 51 | intergenic (+55/‑25) | yebT/rsmF | MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base T (24/26); total (24/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01 |
AAAGCCCAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑TTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCC > NC_000913/1920060‑1920329 | aaaGCCCAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTg < 2:1398923/139‑1 (MQ=255) ccAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttt < 1:60967/139‑1 (MQ=255) aaGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc > 1:1342545/1‑136 (MQ=255) cAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGt > 1:573183/1‑139 (MQ=255) acTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTAc > 1:1193693/1‑139 (MQ=255) cTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCt < 2:573183/139‑1 (MQ=255) gTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCg < 1:1150611/139‑1 (MQ=255) gCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTgg > 1:512509/1‑139 (MQ=255) gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTg < 2:463059/134‑1 (MQ=255) gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTg > 1:463059/1‑134 (MQ=255) tGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGTGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacaca < 2:544014/139‑1 (MQ=255) tGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacaca > 2:191847/1‑139 (MQ=255) taaCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCTCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacacac < 1:1467052/138‑1 (MQ=255) gTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc < 2:381175/99‑4 (MQ=255) gTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc > 1:381175/1‑96 (MQ=255) tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacacac > 2:1181199/1‑127 (MQ=255) tGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacacac < 1:1181199/127‑1 (MQ=255) ggCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGc < 1:166774/117‑1 (MQ=255) ggCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGc > 2:166774/1‑117 (MQ=255) tGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg > 1:72072/1‑139 (MQ=255) gAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc < 1:191847/139‑1 (MQ=255) tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCtt < 1:826932/131‑1 (MQ=255) tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCtt > 2:826932/1‑131 (MQ=255) tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc > 2:1314093/1‑78 (MQ=255) tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑ttttttttgc < 1:1314093/81‑4 (MQ=255) cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACaa > 1:526771/1‑139 (MQ=255) cTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTg > 2:286108/1‑95 (MQ=255) cTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTg < 1:286108/95‑1 (MQ=255) cTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAAt > 2:1363090/1‑139 (MQ=255) cccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcg < 2:1193693/139‑1 (MQ=255) ccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcgc < 2:526771/139‑1 (MQ=255) gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacaca > 2:1090688/1‑92 (MQ=255) gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAacaca < 1:1090688/92‑1 (MQ=255) aCTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCtt < 2:678405/139‑1 (MQ=255) tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACa > 2:421724/1‑68 (MQ=255) tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACa < 1:421724/68‑1 (MQ=255) gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGa < 1:407892/139‑1 (MQ=255) tCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACg < 2:72072/139‑1 (MQ=255) tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCg > 1:185518/1‑64 (MQ=255) tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCg < 2:185518/64‑1 (MQ=255) tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAATGGTACATGCTCg < 2:1315234/64‑1 (MQ=255) tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAATGGTACATGCTCg > 1:1315234/1‑64 (MQ=255) ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAg < 1:1297722/49‑1 (MQ=39) ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAg > 2:1297722/1‑49 (MQ=39) ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGAcaca < 1:593642/105‑1 (MQ=255) ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGAcaca > 2:593642/1‑105 (MQ=255) cTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtg < 2:512509/139‑1 (MQ=255) gcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtgatga > 1:775948/1‑139 (MQ=255) gcCGGAGCGTTTATGCTAAACTGCGCGCCTG‑TTTTTTTCGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtgatga < 2:1342545/139‑1 (MQ=255) ccGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTgatgat > 2:1400699/1‑139 (MQ=255) aCTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTcc < 2:1264707/62‑1 (MQ=255) aCTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTcc > 1:1264707/1‑62 (MQ=255) cTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCt > 2:332721/1‑35 (MQ=255) cTGCGCGCCTGTTTTTTTT‑GCCAGTGGTACATGCt < 1:332721/35‑1 (MQ=255) cgcgcCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGt < 2:187334/50‑1 (MQ=255) cgcgcCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGt > 1:187334/1‑50 (MQ=255) cgcgcCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc < 1:1363090/139‑1 (MQ=255) gcgcCTGTTTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGccc > 2:573531/1‑139 (MQ=255) | AAAGCCCAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑TTTTTTT‑GCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCC > NC_000913/1920060‑1920329 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |