| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 580,851 | G→A | 40.4% | 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 | 40.4% | 49.9 / 70.3 | 69 | K6K (AAG→AAA) | nohD | DLP12 prophage; DNA packaging protein |
| Reads supporting (aligned to +/- strand): ref base G (15/26); new base A (17/11); total (32/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.49e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-01 | |||||||||||
AAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGG > NC_000913/580752‑580940 | aagaaagGGATCTAGAGGTGCTGAAAGCGAGCGTTTTGGCGGCTGTCGTTTCCTGTCTCTGTGTTTGGCCGTGGAATGTGCAATGGAAGTCAACAAAAag > 2:198006/1‑100 (MQ=14)aagaaagGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAag < 1:1955348/100‑1 (MQ=255) AGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc > 1:265105‑M1/1‑100 (MQ=255) aaggaaacgacagaggccaaaaagctcgctttcagcacctgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc > 1:655507‑M2/55‑100 (MQ=255) AAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc > 2:1173486‑M1/1‑100 (MQ=255) AAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc < 1:380887‑M1/100‑1 (MQ=255) ggaaacgacagaggcccaaaagctcgttttcagcacctgtcgtttcctttctCTGTTTTTGGCCGTGGAATGAACAATGGAAGTCAACAAAAAGCAgctg > 1:1487793‑M2/53‑100 (MQ=255) ggaaacgacagaggccaaaaagctcgctttcagcacctgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctg > 2:2334414‑M2/53‑100 (MQ=255) gaaacgacagaggccaaaaagctcgctttcagcacctgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGTAATGGAAGTCAACAAAAAACAgctgg > 2:1492695‑M2/52‑100 (MQ=255) AAACGACAGGTGCTGAAAGCGAGCTTTTTGGTCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTAAACAAAAAACAgctggc > 2:1009104‑M1/1‑100 (MQ=255) AAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgctggc > 1:1240739‑M1/1‑100 (MQ=255) AAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTCTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgctggc > 1:142467‑M1/1‑100 (MQ=255) GACAGGTGCTGAAAGGGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAgctg > 1:2677‑M1/1‑94 (MQ=255) GACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGAc > 2:2376320‑M1/1‑100 (MQ=255) CAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACAt < 1:144672‑M1/100‑1 (MQ=255) gaggccaaaaagctcgctttcagcacctgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAttt > 1:1925106‑M2/43‑100 (MQ=255) GGTGCTGAAAGCTAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAttt > 2:2039785‑M1/1‑100 (MQ=255) GGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAttt > 2:1699112‑M1/1‑100 (MQ=255) GGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGATGACAt < 1:614405‑M1/98‑1 (MQ=255) gTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt < 1:727988/100‑1 (MQ=255) tGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGAAAGACAACAAAAAACAGCTGGCTGACATTTTc < 2:560248/100‑1 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg < 1:748369/100‑1 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCgg > 1:1451571/1‑100 (MQ=255) cTGAAAGCGAGCTTTTTGGCCTCTGTAGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGTCTGACATTTTCgg > 1:1453468/1‑100 (MQ=255) tGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt < 2:832819/100‑1 (MQ=255) gAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt > 1:671891/1‑99 (MQ=255) gAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTg > 2:272433/1‑100 (MQ=255) gcaaGCTTTATGGCCTCTGTCGTTTCCTGTTTCTGTTTTTGTCCGTGGATTGTGCAATAGAAAAAAAAAAAAAACAGCTGGCTGACATTTTCGGTGCGAg < 2:99728/97‑1 (MQ=17) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAg < 2:2199441/100‑1 (MQ=255) gCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAg > 2:213243/1‑100 (MQ=255) cGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGt < 1:1803798/100‑1 (MQ=255) aGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCTGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTAt > 2:1556077/1‑100 (MQ=255) cTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATcc > 1:1035925/1‑100 (MQ=255) tttGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTa < 1:1264482/100‑1 (MQ=255) tGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTAcc < 1:1112446/100‑1 (MQ=255) gCCTCTGTCGTTTCCTTTCTCTGATTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAt > 2:870403/1‑100 (MQ=255) ccTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCAtt < 1:649727/100‑1 (MQ=255) ctctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTc > 2:1533030/1‑100 (MQ=255) tctctcGTTTCCTTTCTCTGTTTTTGTCCGTGTTATGTGCAATAAAAAAAAAAAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTc < 2:1557891/99‑1 (MQ=14) tctGTCGTTTCTTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCATCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCa < 2:2342734/100‑1 (MQ=255) tctGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCa < 1:1283536/100‑1 (MQ=255) ctgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCa > 1:392855‑M2/17‑99 (MQ=255) tgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGa < 1:907335‑M2/85‑1 (MQ=255) tgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGa < 2:489898‑M2/85‑1 (MQ=255) tgtcgtttcctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGa > 2:757860‑M2/16‑100 (MQ=255) ntttcctttctCTGTTTTTGTCAGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACt > 2:487846‑M2/12‑99 (MQ=255) ctttctCTGTTTTTGTTTGTTTAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAg > 2:1406660‑M2/7‑99 (MQ=255) ctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg < 1:573900‑M2/94‑1 (MQ=255) ctttctCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAg > 2:2168645‑M2/7‑99 (MQ=255) tttatCTGTTTTTGTTCGTGGAATGTGCAATGGCAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAgg < 1:1304611‑M2/94‑1 (MQ=255) tCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAg < 1:2123058‑M2/99‑1 (MQ=255) ctGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAgg < 2:1232053/100‑1 (MQ=40) tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg > 1:912464/1‑100 (MQ=255) tGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAggg > 1:1661394/1‑100 (MQ=255) ttttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAAt < 1:1928904/100‑1 (MQ=40) tttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATg < 2:1617806/100‑1 (MQ=255) ttGTCCGTGGAATGTGCAATGGAAGTCAATAACAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATg < 2:1133852/99‑1 (MQ=255) ttGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGc < 1:1533048/100‑1 (MQ=255) tGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGcc < 1:1937420/100‑1 (MQ=255) ccGTGGAATGTGCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCg > 1:363854/1‑99 (MQ=255) cGTGGAATGTGCAATGGAAATCAACAAAAAGCAGCTGGCTGACAGTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGtt < 1:219993/100‑1 (MQ=255) gTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc < 1:1883990/100‑1 (MQ=35) tGGAATGTGCAATGGAAGTCAACAAATAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTc < 2:856879/99‑1 (MQ=255) tGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCt < 1:1895037/100‑1 (MQ=35) tGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCt > 2:241069/1‑100 (MQ=38) aaTGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCg < 1:241072/100‑1 (MQ=35) aaTGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCg < 2:1826094/100‑1 (MQ=35) tgtgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAg < 2:2256522/100‑1 (MQ=33) gtgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAgg < 2:1925084/100‑1 (MQ=32) tgCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc < 1:1012645/100‑1 (MQ=31) gCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGc < 1:1850082/99‑1 (MQ=11) gCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCg > 1:940667/1‑100 (MQ=11) gCAATGGAAGTCAACAAAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCg > 1:2019601/1‑100 (MQ=11) gtgaacaaAAAGCAGCTGGCTGACATTGTCGGTGCGAGTATCCGTAACAATCAGAACGGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAgg < 2:244754/97‑1 (MQ=1) caccaaAAAACAGCTGGCTGACATTTTCGGTGCGAGTATCAGTACCATTCAGAACTGGCAGGAACGGGGAATGCCCGTTATGCGAGGCGGCGACAAggg > 2:2286711/4‑99 (MQ=1) | AAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGTACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGG > NC_000913/580752‑580940 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |