Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 264,717 | (A)5→6 | coding (16/477 nt) | ykfG ← | CP4‑6 prophage; RadC‑like JAB domain‑containing protein YkfG |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 264,712 | 1 | . | A | 97.4% | 133.6 / ‑3.4 | 38 | coding (21/477 nt) | ykfG | CP4‑6 prophage; RadC‑like JAB domain‑containing protein YkfG |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (16/21); total (17/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01 |
CCGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑AAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGT > NC_000913/264632‑264787 | ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt < 1:61383/81‑1 (MQ=255) ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt < 1:61382/81‑1 (MQ=255) ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt < 1:61385/81‑1 (MQ=255) ccGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt < 1:61384/81‑1 (MQ=255) cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt > 2:61427/1‑80 (MQ=255) cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGt > 2:61436/1‑80 (MQ=255) cGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑a < 1:61386/81‑1 (MQ=255) gggCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aa < 2:61387/81‑1 (MQ=255) ggCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aaa < 1:61388/81‑1 (MQ=255) tCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aaaa > 2:61447/1‑79 (MQ=255) tCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aaaa > 2:61440/1‑79 (MQ=255) tGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aaaaaagg > 2:61429/1‑77 (MQ=255) gCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑aaaaaagga > 2:61437/1‑76 (MQ=255) aGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAg > 1:61443/1‑81 (MQ=255) tGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTg < 2:61390/81‑1 (MQ=255) gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGtt > 1:61442/1‑81 (MQ=255) gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGt < 1:61391/80‑1 (MQ=255) gCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑AAAAAGGGAAGCTGttt > 1:61458/1‑81 (MQ=255) ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGttt > 1:61456/1‑80 (MQ=255) ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGttt > 2:61453/1‑80 (MQ=255) ggTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTc < 1:61392/81‑1 (MQ=255) gTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCa > 2:61454/1‑81 (MQ=255) ccAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATc > 2:61444/1‑81 (MQ=255) cAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATc < 2:61393/80‑1 (MQ=255) aGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCtt < 2:61394/81‑1 (MQ=255) gcCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCgg > 1:61432/1‑81 (MQ=255) ccGCTGAATGAGACGCCGGTCCTGTGGCGTCGTCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTc < 2:61395/81‑1 (MQ=255) ccGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTc < 1:61396/81‑1 (MQ=255) cGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCa < 1:61397/81‑1 (MQ=255) gCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAg < 1:61398/81‑1 (MQ=255) cTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGt > 2:61457/1‑81 (MQ=255) aTGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGat < 1:61399/81‑1 (MQ=255) tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata < 2:61400/81‑1 (MQ=255) tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata < 1:61401/81‑1 (MQ=255) tGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGata < 1:61402/81‑1 (MQ=255) gaCGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGataata > 2:61450/1‑81 (MQ=255) cGCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACg < 2:61403/81‑1 (MQ=255) gCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGc < 2:61405/81‑1 (MQ=255) gCCGGTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATAAGc < 1:61404/81‑1 (MQ=255) gTCCTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGaa > 1:61477/1‑81 (MQ=255) cTGTGGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATgg > 1:61472/1‑81 (MQ=255) gtgGCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCg > 2:61448/1‑81 (MQ=255) gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGt > 1:61465/1‑79 (MQ=255) gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTg < 2:61407/80‑1 (MQ=255) gCGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGa < 1:61406/81‑1 (MQ=255) cGTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAg > 1:61451/1‑81 (MQ=255) gTCATCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGc > 1:61455/1‑81 (MQ=255) aTCTCGCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCctc < 1:61408/81‑1 (MQ=255) gCCGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTg < 2:61409/78‑1 (MQ=255) gCCGGGTAAAAAAGGAAAGCTGTTTCATCTATTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGa < 2:61410/80‑1 (MQ=255) cGGGTAAAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGt > 1:61471/1‑81 (MQ=255) | CCGGGCTCATGCAGGTGGCGGTCCAGAGCCCTGAGCGCCCGCTGAATGAGACGCCGGTCCTGTGGCGTCATCTCGCCGGGT‑AAAAAGGAAAGCTGTTTCATCTGTTGCTCCTTCGGTCAGTCGATAATACGCAGAATGGCGTGAGCCTCTGGATGT > NC_000913/264632‑264787 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |