| 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 | 93.3% | 99.2 / 0.8 | 30 | coding (525/858 nt) | ydaU | Rac prophage, conserved protein |
| Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (15/13); total (15/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.83e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01 | |||||||||||
ATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGC > NC_000913/1421510‑1421776 | aTGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑ccccc < 2:630900/139‑1 (MQ=255) aaaacaaGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGaa > 2:482599/1‑139 (MQ=255) aaaTTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCaa > 1:587828/1‑139 (MQ=255) aaCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGtt < 2:577767/139‑1 (MQ=255) ggTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAc < 1:776069/139‑1 (MQ=255) gtgtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGAcc < 2:505001/139‑1 (MQ=255) gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:694256/139‑1 (MQ=255) gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:587828/139‑1 (MQ=255) gtgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAAc > 1:154542/1‑139 (MQ=255) tgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCaa > 1:319631/1‑104 (MQ=255) tgtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCaa < 2:319631/104‑1 (MQ=255) gtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg < 2:780730/111‑1 (MQ=255) gtTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg > 1:780730/1‑111 (MQ=255) cacaAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTg > 2:460615/1‑139 (MQ=255) aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTgg > 2:31116/1‑139 (MQ=255) taaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑ccccccccgg > 2:285648/1‑50 (MQ=255) taaAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑CCCCCCCcgg < 1:285648/53‑1 (MQ=255) aaaGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGaa < 1:460615/139‑1 (MQ=255) aaGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct > 1:285740/1‑81 (MQ=255) aaGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCct < 2:285740/81‑1 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACCTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGc > 2:581423/1‑139 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑ccc > 1:303990/1‑39 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑ccc < 2:303990/39‑1 (MQ=255) gATAAAGATCTAAAAACAGATCCCCCCCTAAAAC‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc < 1:539589/137‑1 (MQ=255) gATAAAGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGACTGGGTTGAATTTCGcc > 2:539589/1‑137 (MQ=255) aaaGATCTAAAAACAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCa > 1:128779/1‑139 (MQ=255) aaaaaCAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCaaa < 2:301239/50‑1 (MQ=255) aaaaaCAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCaaa > 1:301239/1‑50 (MQ=255) aaaaaCAGATCCCCCCCT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaa > 2:245566/1‑139 (MQ=255) ccccT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCg < 2:5833/67‑1 (MQ=255) ccccT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCg > 1:5833/1‑67 (MQ=255) ccccT‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACg > 1:308386/1‑139 (MQ=255) t‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGc < 2:128779/139‑1 (MQ=255) t‑AAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCAAGGCATTGCGAAAACCGATTCGAACGGAGc < 1:482599/139‑1 (MQ=255) | ATGTTCAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCT‑AAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGC > NC_000913/1421510‑1421776 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |