| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,769,122 | C→T | 13.3% | A153V (GCG→GTG) | yfjQ → | CP4‑57 prophage; uncharacterized protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,769,122 | 0 | C | T | 13.3% | 82.3 / 3.5 | 30 | A153V (GCG→GTG) | yfjQ | CP4‑57 prophage; uncharacterized protein |
| Reads supporting (aligned to +/- strand): ref base C (8/18); new base T (2/2); total (10/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTACTGAA > NC_000913/2768978‑2769267 | aCTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGGTATTGAGGGGGCGTATg > 2:31898/1‑150 (MQ=255)aCTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGGGATTGAGGGGGCGTATg > 2:91847/1‑150 (MQ=255) gATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCAGACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTg < 2:234680/150‑1 (MQ=255) gATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTg < 2:5567/150‑1 (MQ=255) gCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGTGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGa < 1:57828/150‑1 (MQ=255) gCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGa < 1:120334/150‑1 (MQ=255) gCACCAGTTCGTACCAGAGGCTGCCGGGCATGTTCAGGGTGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGa < 1:119148/150‑1 (MQ=255) cAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGGGCCACACCAGGGGGAGGTGGTGCGTCAGGGTCTTGAGGGGGCGGATGAAGTGCGGGGGAtttt > 1:67268/1‑150 (MQ=21) tACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGGGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGTGTATGAAGTGCTGGGGATTTTTGAGCgt > 1:125017/1‑150 (MQ=255) tACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGTGTATGAAGTGCTGGGGATTTTTGAGCgt > 1:125014/1‑150 (MQ=255) ggCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAgg > 2:71056/1‑150 (MQ=255) cATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGa < 1:152500/150‑1 (MQ=255) tCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGTAAAGAGGGATGccc > 2:173512/1‑149 (MQ=255) tgCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTg < 2:58640/150‑1 (MQ=255) tgCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTg < 2:58636/150‑1 (MQ=255) ggtctggtctGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCg < 2:196558/150‑1 (MQ=255) gcggcgAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTgc > 2:196025/1‑150 (MQ=255) gcgAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTgcagc > 2:166900/1‑150 (MQ=255) gcgAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTAACGCCACCAGTgcagc > 2:123285/1‑150 (MQ=255) tCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCa < 2:236802/150‑1 (MQ=255) gCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCaa < 1:182151/150‑1 (MQ=255) aGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAag < 1:71056/150‑1 (MQ=255) tGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAagcag < 2:200201/150‑1 (MQ=255) aTGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTg < 2:62034/150‑1 (MQ=255) tggtgAGTCAGGTTATTGAGGGGGTGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAgg < 2:125017/150‑1 (MQ=255) tggtgAGTCAGGTTATTGAGGGGGTGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAgg < 2:125014/150‑1 (MQ=255) atGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTa < 1:95919/149‑1 (MQ=255) gAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTACt < 2:44321/150‑1 (MQ=255) gAGGGGGCGTATGAAGTGCTGGGGATTTGTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTGTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTACt < 2:44320/150‑1 (MQ=255) ggggCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGTGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTACTGaa < 2:189923/150‑1 (MQ=255) | ACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAGGGCGGTGTGCCAGAACGGTCTGGTCTGCGGCGAGTCGTTTGGCGAGGTGCGGGTGCCACACAAGGGGGATGTGGTGAGTCAGGTTATTGAGGGGGCGTATGAAGTGCTGGGGATTTTTGAGCGTGTGGAAGAAAAGAGGGATGCCATGCAGTCGTTGCTTTTACCGCCACCTGTGCAGCAGGCACTGGCAAAAGCAGCGCTTACGTATCGCTTTGGTGAGGACCACCAGCCGGTTACTGAA > NC_000913/2768978‑2769267 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |