Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 281,662 | Δ1 bp | coding (322/1155 nt) | yagA ← | CP4‑6 prophage, putative DNA‑binding transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 281,660 | 0 | G | . | 98.0% | 194.0 / ‑2.7 | 49 | coding (324/1155 nt) | yagA | CP4‑6 prophage, putative DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (21/27); minor base C (0/1); total (21/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGA > NC_000913/281525‑281779 | aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca < 1:1464319/139‑3 (MQ=255) gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtgt > 2:360348/1‑139 (MQ=255) gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGgtgt > 1:1395774/1‑139 (MQ=255) cTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCG‑CCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGccc > 2:829891/1‑139 (MQ=255) tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGcccc < 2:715110/139‑1 (MQ=255) cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCa < 1:550226/139‑1 (MQ=255) gcggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCCGGGCATCACGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAgc > 2:1474552/3‑139 (MQ=255) ggCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc < 2:1158697/139‑1 (MQ=255) gCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCt < 2:199086/139‑1 (MQ=255) tCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCa < 2:1000649/139‑1 (MQ=255) cccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATg > 1:1244526/1‑139 (MQ=255) gTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACg < 1:360348/139‑1 (MQ=255) aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg < 1:459169/111‑1 (MQ=255) aTGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCggg > 2:459169/1‑111 (MQ=255) ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATa < 2:820475/139‑1 (MQ=255) cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCTCAGCGTTCATGACGGTCATGGGCCATACg < 2:871856/132‑1 (MQ=255) cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg > 1:1089704/1‑132 (MQ=255) cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg < 2:1089704/132‑1 (MQ=255) cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACg > 1:871856/1‑132 (MQ=255) ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa < 1:821637/110‑1 (MQ=255) ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa > 2:821637/1‑110 (MQ=255) gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg < 1:829891/139‑1 (MQ=255) gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg < 2:915950/139‑1 (MQ=255) ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt > 1:1494913/1‑81 (MQ=255) ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCtt < 2:1494913/81‑1 (MQ=255) tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg > 1:1241969/1‑108 (MQ=255) tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg < 2:1241969/108‑1 (MQ=255) gggCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg < 2:758504/139‑1 (MQ=255) gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca < 1:1497177/114‑1 (MQ=255) gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCGTACgcagca > 2:1497177/1‑114 (MQ=255) gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc > 2:1082070/1‑138 (MQ=255) gCCATCAGGTTATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc < 1:1082070/138‑1 (MQ=255) ttATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa < 1:1474552/139‑1 (MQ=255) tATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc > 2:58930/1‑71 (MQ=255) tATGGACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc < 1:58930/71‑1 (MQ=255) ggACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc > 2:894413/1‑106 (MQ=255) ggACGGTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGcc < 1:894413/106‑1 (MQ=255) gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg < 2:179035/126‑1 (MQ=255) gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg > 1:179035/1‑126 (MQ=255) gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa > 2:1206972/1‑130 (MQ=255) gTGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa < 1:1206972/130‑1 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg > 1:859363/1‑83 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg < 2:859363/83‑1 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg > 2:923992/1‑82 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg > 1:59511/1‑118 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg < 1:1276991/118‑1 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg < 2:59511/118‑1 (MQ=255) tGCTGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg > 2:1276991/1‑118 (MQ=255) tGCTGAA‑GCCGGCATGGTGTGCCCCTGGTCTTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATgg < 1:923992/82‑1 (MQ=255) tGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa < 2:1194712/126‑1 (MQ=255) tGAAGGC‑GGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa > 1:1194712/1‑126 (MQ=255) | AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGA > NC_000913/281525‑281779 |
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: ‑ |