Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,421,650 | (C)7→8 | coding (532/858 nt) | ydaU → | Rac prophage, conserved protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,421,643 | 1 | . | C | 95.3% | 156.5 / 0.2 | 43 | coding (525/858 nt) | ydaU | Rac prophage, conserved protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (25/16); total (26/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01 |
ACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA‑CCGATTCGAACGGAG > NC_000913/1421518‑1421775 | acaaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCcggggg > 1:632846/1‑138 (MQ=255) aaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGAT‑‑CCCCCTAAATCCCCCCCCGGGGGAATCg < 2:422272/139‑1 (MQ=255) aacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATc > 1:1062066/1‑139 (MQ=255) aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑ccccccccgg > 1:15514/1‑113 (MQ=255) aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCcgg < 2:15514/116‑1 (MQ=255) cAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg < 2:482978/139‑1 (MQ=255) ggTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAc > 2:873686/1‑139 (MQ=255) gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:335327/139‑1 (MQ=255) ttCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctc > 1:301785/1‑139 (MQ=255) gATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGttt < 1:52432/124‑1 (MQ=255) gATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGttt > 2:52432/1‑124 (MQ=255) tGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct > 1:26733/1‑128 (MQ=255) tGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:26733/128‑1 (MQ=255) tGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑cc > 1:520380/1‑90 (MQ=255) tGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑cc < 2:520380/90‑1 (MQ=255) gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg < 2:632846/139‑1 (MQ=255) gAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTg > 1:562121/1‑139 (MQ=255) gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGgt < 2:837151/100‑1 (MQ=255) gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGgt > 1:837151/1‑100 (MQ=255) gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAAc > 2:947068/1‑139 (MQ=255) gtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTg < 2:1062066/139‑1 (MQ=255) atgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTg > 2:495420/1‑139 (MQ=255) tgatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGt > 2:181045/1‑139 (MQ=255) aaaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 1:542944/95‑1 (MQ=255) aaaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct > 2:542944/1‑95 (MQ=255) aaaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCg > 1:411457/1‑139 (MQ=255) aaGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGt > 2:120404/1‑139 (MQ=255) aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg > 2:222159/1‑139 (MQ=255) aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTg > 1:873034/1‑139 (MQ=255) aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg > 1:101395/1‑139 (MQ=255) aaaGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGaa < 2:873034/139‑1 (MQ=255) aGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAAtt < 1:181045/139‑1 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCaa < 1:1009674/63‑1 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCaa > 2:1009674/1‑63 (MQ=255) aGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct > 1:212976/1‑74 (MQ=255) aGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:212976/74‑1 (MQ=255) cTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGa > 2:568754/1‑139 (MQ=255) cTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGa > 1:122957/1‑139 (MQ=255) tAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaa > 1:125795/1‑139 (MQ=255) aaaaaCAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCt < 1:863066/100‑1 (MQ=255) aaaaaCAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCt > 2:863066/1‑100 (MQ=255) gATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA‑CCGAtt > 1:804527/1‑139 (MQ=255) ccccccTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA‑CCGATTCGaa < 2:122957/139‑1 (MQ=255) ccTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCCGATTCGAACgg > 2:853371/1‑139 (MQ=255) cTAAATCCCCCCCCGGGGGAATCGGGGGGGCAAAAAGGTTGGCCCTCTGGATATTACTTTTCCGAAATGGAATTCTGTCTCGCTTTTGCGGGAGGGGGGTTAATTTCGGCCGGCATTGCGAAAA‑CCGATTCGAACgggg > 2:764403/1‑137 (MQ=255) cTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA‑CCGATTCGAACGGAg < 1:947068/139‑1 (MQ=255) | ACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA‑CCGATTCGAACGGAG > NC_000913/1421518‑1421775 |
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: ‑ |