| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,428,030 | C→G | A233P (GCC→CCC) | insH‑5 ← | Rac prophage; IS5 transposase and trans‑activator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,428,030 | 0 | C | G | 97.0% | 95.2 / ‑3.3 | 33 | A233P (GCC→CCC) | insH‑5 | Rac prophage; IS5 transposase and trans‑activator |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base G (15/17); total (16/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
CTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCC > NC_000913/1427956‑1428093 | cTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctc > 2:272015/1‑81 (MQ=255) gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc > 2:272026/1‑81 (MQ=255) gcgcGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcct < 2:271988/80‑1 (MQ=255) cgcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc < 1:271989/80‑1 (MQ=255) cgcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGctcctc < 2:271990/80‑1 (MQ=255) cgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgcg > 1:272040/1‑81 (MQ=255) gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCt < 1:271991/81‑1 (MQ=255) gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCt > 1:272022/1‑81 (MQ=255) gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCt > 1:272018/1‑81 (MQ=255) ttCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCtgt < 1:271992/81‑1 (MQ=255) cTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGgcg < 1:271993/80‑1 (MQ=255) cAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCt < 1:271994/79‑1 (MQ=255) ggTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTg > 2:272027/1‑81 (MQ=255) tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTgg > 1:272038/1‑80 (MQ=255) ttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTgg > 1:272016/1‑79 (MQ=255) ttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGc > 2:272029/1‑81 (MQ=255) aCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCg > 1:272028/1‑81 (MQ=255) cTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCgg > 1:272036/1‑80 (MQ=255) ttGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGAGGCGCTCCTTGGTAGCCGGCa > 1:272019/1‑81 (MQ=255) aCGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTga < 2:271995/81‑1 (MQ=25) cGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTgag < 1:271996/81‑1 (MQ=18) gCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTgaga > 1:272039/1‑81 (MQ=18) ggCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACaaa > 1:272045/1‑81 (MQ=255) gCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTTGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAAt > 2:405376/1‑81 (MQ=18) cGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAAtt > 2:272030/1‑81 (MQ=255) gATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTg < 1:271997/81‑1 (MQ=18) aGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctcc < 2:271999/81‑1 (MQ=18) aGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctcc < 2:272000/81‑1 (MQ=18) aGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctcc < 1:272001/81‑1 (MQ=18) aGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctc < 1:271998/80‑1 (MQ=21) ccAATCAACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctcctc < 1:272002/81‑1 (MQ=11) aGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCtctc < 1:272003/81‑1 (MQ=18) gTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTcc < 2:272004/81‑1 (MQ=18) | CTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCC > NC_000913/1427956‑1428093 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |