| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 578,568 | G→A | *154* (TAG→TAA) | rzpD → | DLP12 prophage, putative murein endopeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 578,568 | 0 | G | A | 97.2% | 216.9 / ‑1.9 | 72 | *154* (TAG→TAA) | rzpD | DLP12 prophage, putative murein endopeptidase |
| Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (34/36); total (35/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01 | |||||||||||
GGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAA > NC_000913/578434‑578700 | ggCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt < 1:1444206/139‑1 (MQ=255) tGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAcc < 1:1492969/139‑1 (MQ=255) tGCAACCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACAATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGa > 2:984420/1‑139 (MQ=255) aGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATgg > 1:590151/1‑139 (MQ=255) ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc < 1:984420/139‑1 (MQ=255) ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc < 1:272460/139‑1 (MQ=255) cTGGCAGACACCGCTGAACGGGATTATTGCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGTAGGAACCCAGAAGTATATTAATGGGCTGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa < 1:231885/139‑1 (MQ=39) gCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAatta > 2:910802/1‑139 (MQ=255) acacCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt < 2:739087/139‑1 (MQ=255) cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa < 2:110227/139‑1 (MQ=255) cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa > 1:110227/1‑139 (MQ=255) cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt < 2:1113268/81‑1 (MQ=255) cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCt > 1:1113268/1‑81 (MQ=255) ccTCAGAGAGAGGCTGATCACTATGCAAACACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttta > 2:543100/1‑114 (MQ=255) ccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTg < 1:543100/115‑1 (MQ=255) agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg < 1:71883/113‑1 (MQ=255) agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg > 2:71883/1‑113 (MQ=255) agagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaa < 2:590151/139‑1 (MQ=255) agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt < 2:1510434/104‑1 (MQ=255) agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTAtttt > 1:1510434/1‑104 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc < 2:1628341/117‑1 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc > 1:1628341/1‑117 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg > 1:177882/1‑125 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg < 2:177882/125‑1 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg > 2:1599163/1‑128 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg < 1:1599163/128‑1 (MQ=255) gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGc > 1:305914/1‑139 (MQ=255) tGATCACTATGCAAAAACAACTGGAAGGAGCCCAGAAGTATATTACTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAAt > 2:976485/1‑107 (MQ=255) tGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAAt < 1:976485/107‑1 (MQ=255) gATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATaaag > 2:92837/1‑63 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATg > 1:572905/1‑77 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATg < 2:572905/77‑1 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc < 1:1016245/136‑1 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc > 2:1016245/1‑136 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa < 2:1031547/139‑1 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa > 1:1198541/1‑139 (MQ=255) aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATAATAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa > 1:1031547/1‑139 (MQ=255) tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCa > 2:1330217/1‑98 (MQ=255) tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCa < 1:1330217/98‑1 (MQ=255) aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATggg < 1:1263141/69‑1 (MQ=255) aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa > 2:730359/1‑139 (MQ=255) aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGACAAAGCTGACCATATCGATggg > 2:1263141/1‑69 (MQ=255) aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATaaag < 1:324477/50‑1 (MQ=255) aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc > 1:1341944/1‑101 (MQ=255) aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc < 2:1341944/101‑1 (MQ=255) aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg > 2:1111091/1‑109 (MQ=255) aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg < 1:1111091/109‑1 (MQ=255) acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa > 1:548073/1‑139 (MQ=255) tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACAcgc > 1:1692481/1‑94 (MQ=255) tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATCGATGGGCAACTCATGCAATTATTGTGAGCAATACACAcgc < 2:1692481/94‑1 (MQ=255) tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt < 1:910802/139‑1 (MQ=255) cccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg > 2:32314/1‑46 (MQ=255) cccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg < 1:32314/46‑1 (MQ=255) cccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc < 2:211984/85‑1 (MQ=255) cccAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc > 1:211984/1‑85 (MQ=255) aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt < 2:774198/65‑1 (MQ=255) aGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt > 1:774198/1‑65 (MQ=255) aaGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATga > 2:822452/1‑138 (MQ=255) atatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGc < 2:1198541/139‑1 (MQ=255) atatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGc < 1:822452/139‑1 (MQ=255) tatTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCt < 2:305914/139‑1 (MQ=255) atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg < 1:749262/84‑1 (MQ=255) atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg > 2:749262/1‑84 (MQ=255) aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc > 2:863202/1‑56 (MQ=255) aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc < 1:863202/56‑1 (MQ=255) cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg > 2:1491344/1‑53 (MQ=255) cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg < 1:1491344/53‑1 (MQ=255) cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGt > 2:1362124/1‑139 (MQ=255) cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGt > 2:1536929/1‑139 (MQ=255) tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa > 1:1060376/1‑105 (MQ=255) tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa < 2:1060376/105‑1 (MQ=255) cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACgcgc > 2:567583/1‑62 (MQ=255) cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACgcgc < 1:567583/62‑1 (MQ=255) cAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTaa < 1:1362124/139‑1 (MQ=255) | GGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAA > NC_000913/578434‑578700 |
| 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: ‑ |