| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,641,226 | G→A | V95V (GTG→GTA) | yjtD → | putative tRNA/rRNA methyltransferase YjtD |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,641,226 | 0 | G | A | 95.0% | 53.6 / ‑4.0 | 21 | V95V (GTG→GTA) | yjtD | putative tRNA/rRNA methyltransferase YjtD |
| Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (11/8); total (11/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01 | |||||||||||
AAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTG > NC_000913/4640982‑4641437 | aaaaTATTGGGGCAGCGGCGCGGGCAATAAAACCGATGGGTTTTAGCGTACTGCGGATTGGCGATAGTCAGGCACACCTGTAGCCAGCCACCCGCTGGGTCGCACAGTGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATAGCATTACTACGCTACTCCAGTAGAACTGGTACCGCTg < 2:350921/251‑1 (MQ=255) ggCAATGAAAACGATGGGTTTTAGCGAACTGCGGATTGTCGATAGTCAGGCACACCTTGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCATTAGAACTGGTTCCGCTGTTAGAGGAAAAATCTTCATGGa < 1:384754/251‑1 (MQ=255) aaaTGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATg < 2:158857/250‑1 (MQ=255) cAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGAc > 1:66593/1‑251 (MQ=255) aGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCCGTCGGGCGAAATATCATTACTACGCTACGCCCGTAGAACTGGTACCGCTGTTAGCGGGAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACt > 2:47813/1‑251 (MQ=255) cctgccaCCCTCTGGGTCGCACAGGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGTTCGAAATATCATTTCTACGCTACGCCAGTAGAACTGGTACCGCGGTTAGAGGAAAAATCTTCATGGAGGGGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGc < 1:236985/248‑1 (MQ=255) aCCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCACGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGACGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGTAGAGTTAGCGTTGGc > 1:306682/1‑251 (MQ=255) gggTCGCCCATGGATCGGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCTCTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGACGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGt < 1:46100/250‑1 (MQ=255) caTGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTgg > 1:313326/1‑251 (MQ=255) ggtattatatatgatatatttAAAGTTTTCCCGACACTGGCTGAATCGTTACGCGATGTCGATTTC‑ACTGTCTCCACGACGGCGCGCAGTCGGGCGAAATAGCATTACTACGCTACGCCAGTCGAACGGTTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGAGCGGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAt < 2:142672/232‑1 (MQ=255) attattTAAAATATTAAAGTTGTCCCGACACGGGTTGAATCGGTACACGATGTCG‑TTTCTACTTTCGCCACCTCTGCGCGCAGTCGGGCGAAATATCATTACTACGCTTCGCCATTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGTCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGAtggtgg < 2:304579/250‑1 (MQ=255) attattGATAATATTAAAGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCAACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGATAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGTAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGc > 1:84805/1‑251 (MQ=255) aGTTTTCCCGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCa < 1:206155/251‑1 (MQ=255) cGACACTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTgggg > 1:207890/1‑251 (MQ=255) gAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGCGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGAt > 1:18769/1‑251 (MQ=255) gAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGAt < 1:360292/251‑1 (MQ=255) gATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCt > 1:75912/1‑250 (MQ=255) gATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAAtt < 1:323728/251‑1 (MQ=255) acTGCGCGCAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATacaac > 1:363576/1‑250 (MQ=255) cAGTCGGGCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCgacg > 1:159903/1‑247 (MQ=255) gCGAAATATCATTACTACGCTACGCCAGTAGAACTGGTACCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGCTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCAGATTATCCTTCGCTCAATCTGGGGCAGGCAGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGCCGAAAAGTg > 2:90097/1‑250 (MQ=255) | AAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGATGTCGATTTC‑ACTGTCGCCACCACTGCGCGCAGTCGGGCGAAATATCATTACTACGCCACGCCAGTTGAACTGGTGCCGCTGTTAGAGGAAAAATCTTCATGGATGAGCCATGCCGCGCTGGTGTTTGGTCGCGAAGATTCCGGGTTGACTAACGAAGAGTTAGCGTTGGCTGACGTTCTTACTGGTGTGCCGATGGTGGCGGATTATCCTTCGCTCAATCTGGGGCAGGCGGTGATGGTCTATTGCTATCAATTAGCAACATTAATACAACAACCGGCGAAAAGTG > NC_000913/4640982‑4641437 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |