| 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; DUF1376 domain‑containing protein YdaU |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,421,643 | 1 | . | C | 94.7% | 59.4 / ‑3.0 | 19 | coding (525/858 nt) | ydaU | Rac prophage; DUF1376 domain‑containing protein YdaU |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (13/5); total (13/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.16e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA > NC_000913/1421515‑1421760 | caaacaaaACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCcggg < 2:41104/140‑1 (MQ=255) aCAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTg < 2:116456/140‑1 (MQ=255) cAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGa > 1:304838/1‑140 (MQ=255) gtTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCctc > 1:672456/1‑140 (MQ=255) gATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTa > 1:293320/1‑140 (MQ=255) aaGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGcc > 1:789222/1‑140 (MQ=255) gtggtgTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGa > 2:732804/1‑140 (MQ=255) gatgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtct > 2:30243/1‑140 (MQ=255) atgTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAGCAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGAGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGtctc > 2:574374/1‑140 (MQ=255) cAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGAAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCtt > 2:390690/1‑140 (MQ=255) aGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTgg > 1:633962/1‑139 (MQ=255) aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTggg > 1:102142/1‑140 (MQ=255) aCTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTggg < 2:789222/140‑1 (MQ=255) tACAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGc > 2:805866/1‑140 (MQ=255) aCAGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGcc < 2:102142/140‑1 (MQ=255) aGATAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAg > 1:415390/1‑140 (MQ=255) tAAAGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCa < 2:293320/140‑1 (MQ=255) aaGATCTAAAAACAGATCCCCCCCTAAATCCCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCAtt > 1:365687/1‑140 (MQ=255) aaaaaCAGATCCCCCCCTAAATCCCCCCCCGGGGAAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGaaaa < 1:390690/140‑1 (MQ=255) | CAAACAAAACAAGAGAGAAATTTAACAGGTGTTCAAACAGATGTTGAAGTGGTGTTTGAACATGATGTCAACACAAAGGCAACTAATAAAGATACAGATAAAGATCTAAAAACAGATCCCCCCCTAAAT‑CCCCCCCGGGGGAATCGAGGTGTCAAAAAGTTTGACCCTCTGGATATTACTTTGCCGAACTGGATTTCTGTCTCGCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAA > NC_000913/1421515‑1421760 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |