| 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.1% | 108.7 / ‑3.1 | 35 | 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 (22/12); total (23/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
CGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTA > NC_000913/1427960‑1428106 | cgcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcg > 1:262665/1‑81 (MQ=255) gcgGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTcgc < 1:262633/81‑1 (MQ=255) ggATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCgcgc < 1:262634/81‑1 (MQ=255) gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCt < 2:262635/81‑1 (MQ=255) gATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCt > 1:262684/1‑81 (MQ=255) tCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCtgtg > 1:262668/1‑81 (MQ=255) cTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGgcgc < 1:262636/81‑1 (MQ=255) aGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCtt > 1:262673/1‑81 (MQ=255) gTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTgg > 2:262672/1‑81 (MQ=255) tttttACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGt < 2:262637/81‑1 (MQ=255) aCCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCg > 1:262680/1‑81 (MQ=255) cTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCgg > 2:262678/1‑80 (MQ=255) ttGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCa > 2:262689/1‑81 (MQ=255) tGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCAt > 1:262667/1‑81 (MQ=255) gCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATc > 1:262687/1‑81 (MQ=255) ccGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCg > 2:262663/1‑81 (MQ=255) gggACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGc < 1:262638/81‑1 (MQ=255) cGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCt > 2:262669/1‑78 (MQ=21) cGCTCGGCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTgag < 2:262639/81‑1 (MQ=18) cTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTgaga > 2:408215/1‑80 (MQ=34) ggCGATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACaaa > 1:262685/1‑81 (MQ=255) gATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTg > 1:262674/1‑81 (MQ=255) gATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTg < 2:262641/81‑1 (MQ=18) gATCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTg < 1:262640/81‑1 (MQ=18) aTCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGc > 1:262690/1‑81 (MQ=18) tCAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGct < 1:262642/81‑1 (MQ=18) cAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctc > 1:262676/1‑81 (MQ=18) cAGCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctc > 2:262693/1‑81 (MQ=18) gCCAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGctcct > 1:262697/1‑81 (MQ=18) cAGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCctct > 2:262671/1‑81 (MQ=18) aGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCtctc < 2:262643/81‑1 (MQ=18) aGTCCACATCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCtctc > 2:262675/1‑81 (MQ=18) aTCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAg < 1:262644/81‑1 (MQ=255) aTCCACCTCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAg > 2:262696/1‑81 (MQ=255) tCGGGCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTa > 1:262686/1‑81 (MQ=255) | CGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAGTCCACATCCACCTCGGCCAGCTCCTCGCGCTGTGGCGCTCCTTGGTAGCCGGCATCGGCTGAGACAAATTGCTCCTCTCCATGAAGCAGATTA > NC_000913/1427960‑1428106 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |