| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 580,851 | G→A | K6K (AAG→AAA) | nohD → | DLP12 prophage; DNA packaging protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 580,851 | 0 | G | A | 97.9% | 145.4 / ‑4.2 | 47 | K6K (AAG→AAA) | nohD | DLP12 prophage; DNA packaging protein |
| Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (26/20); total (26/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.67e-01 | |||||||||||
AGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCG > NC_000913/580757‑580930 | agGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCTGTGGAATGTGCAATGGAAGTCAACAAAAAACAgct > 2:127523/1‑100 (MQ=255)agGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgct > 1:546664/1‑100 (MQ=255) gAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgctgg > 1:539404/1‑100 (MQ=255) aCAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACa > 1:842051/1‑100 (MQ=255) cAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt < 2:1154537/100‑1 (MQ=255) gCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTc > 1:693022/1‑99 (MQ=255) gCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCg < 1:467936/100‑1 (MQ=255) gCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCg > 1:337526/1‑100 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCgg > 1:1479398/1‑100 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCgg < 1:1381847/100‑1 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCgg < 1:1339412/100‑1 (MQ=255) tGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGGTGACATTTTCg > 2:1496250/1‑98 (MQ=255) aaaGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGc < 1:1421198/100‑1 (MQ=255) aaaGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGc < 2:1013254/100‑1 (MQ=255) aGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGa > 2:1218458/1‑100 (MQ=255) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAg < 1:341292/100‑1 (MQ=255) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAg < 1:598184/100‑1 (MQ=255) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAg < 2:834818/100‑1 (MQ=255) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAAACAAAAAAAAACAGCTGGCTGACATTTTCGGTGCGAg < 2:430479/100‑1 (MQ=17) cGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGt > 2:1083289/1‑100 (MQ=255) gAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTa > 1:311900/1‑100 (MQ=255) gAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTa > 2:525771/1‑100 (MQ=255) tttttGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCg > 2:725697/1‑100 (MQ=255) ttttGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATAAAAAAAAAAAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGt < 2:1113040/100‑1 (MQ=14) ggCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCa < 1:31332/100‑1 (MQ=255) ctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAg < 2:1126269/100‑1 (MQ=255) gTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGaa > 1:1215528/1‑100 (MQ=255) tCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCa < 2:1215521/100‑1 (MQ=255) cTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg > 1:1068141/1‑100 (MQ=255) cTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg > 1:1354003/1‑100 (MQ=255) cTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGACAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg < 2:546661/100‑1 (MQ=255) tctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAg > 1:1225404/1‑100 (MQ=255) tctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAg < 2:427262/100‑1 (MQ=255) ctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAgg > 1:534886/1‑100 (MQ=255) ttttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAAt < 2:595167/100‑1 (MQ=255) gTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGccc < 1:440611/100‑1 (MQ=255) gTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGccc > 2:19296/1‑100 (MQ=255) cGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGtt > 2:1418711/1‑100 (MQ=255) gTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTCTCAGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc > 2:1101553/1‑100 (MQ=255) tGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCt < 1:725700/100‑1 (MQ=255) ggAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTg > 1:394325/1‑100 (MQ=255) aaTGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGTAATGCCCGTTCTGCg > 1:1546325/1‑100 (MQ=255) tgtgCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAg < 2:1546315/100‑1 (MQ=18) gtgCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAgg < 1:680639/100‑1 (MQ=17) tgCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc > 1:786694/1‑100 (MQ=14) tgCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc > 1:1218633/1‑100 (MQ=14) gCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTCCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCg > 2:845320/1‑100 (MQ=2) | AGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCG > NC_000913/580757‑580930 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |