Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,866,592 | +C | coding (960/993 nt) | rpoS ← | RNA polymerase, sigma S (sigma 38) factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,866,588 | 1 | . | C | 98.6% | 249.2 / ‑3.3 | 71 | coding (964/993 nt) | rpoS | RNA polymerase, sigma S (sigma 38) factor |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (30/40); total (30/41) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 |
TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA > NC_000913/2866491‑2866682 | tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc < 2:182065/100‑1 (MQ=255) ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cccc > 1:273089/1‑100 (MQ=255) gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑ccccc > 2:851767/1‑99 (MQ=255) ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGc < 1:127415/100‑1 (MQ=255) ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGc > 2:733713/1‑100 (MQ=255) ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGc < 1:249794/100‑1 (MQ=255) tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCa < 2:1203798/100‑1 (MQ=255) tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAgg > 2:10268/1‑100 (MQ=255) tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAgg < 2:1395173/100‑1 (MQ=255) cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGAt < 1:1211825/100‑1 (MQ=255) cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGAt > 2:590074/1‑100 (MQ=255) acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGCtt > 1:1800498/1‑100 (MQ=255) caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGAttt < 2:1359683/100‑1 (MQ=255) gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGCTATTCAGCCCCC‑TGCGTTTGCAGGATTTc > 1:430584/1‑99 (MQ=255) gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTcg > 2:1015193/1‑100 (MQ=255) aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTcg < 2:1496530/99‑1 (MQ=255) ggCCAGCCTCGCTTGAGACCGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCaa > 2:1645870/1‑100 (MQ=255) ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAAc < 2:1692025/100‑1 (MQ=255) cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACg > 2:749583/1‑100 (MQ=255) aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACgg > 2:939598/1‑100 (MQ=255) gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGgc > 2:143916/1‑100 (MQ=255) cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGgcgc > 1:1773705/1‑100 (MQ=255) tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCa < 1:164485/100‑1 (MQ=255) tCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCa < 1:624975/100‑1 (MQ=255) cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAg < 1:407851/100‑1 (MQ=255) cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAg < 2:980544/100‑1 (MQ=255) cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc < 1:1812687/100‑1 (MQ=255) ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc < 2:546889/100‑1 (MQ=255) gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt > 1:1185380/1‑100 (MQ=255) aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGCAACAGCGCTTCGATATTCAG‑CCCCTTGCGTTTTCAGTATTTCGCGCAAACGGCCCAGGCCTTCaa < 1:1686584/100‑1 (MQ=255) cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc > 1:1353389/1‑100 (MQ=255) cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc > 1:945531/1‑100 (MQ=255) tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc < 2:833453/100‑1 (MQ=255) gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg < 2:1606885/100‑1 (MQ=255) ccTTTCTGACAGATGCTTACTTATTGGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTTGCGTTAATGGCGCAGGCCTTCAACCTg < 1:1038538/99‑1 (MQ=255) cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa < 1:953863/100‑1 (MQ=255) cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa < 1:1336498/100‑1 (MQ=255) cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGGTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt > 2:1202860/1‑99 (MQ=255) tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg < 2:2438/100‑1 (MQ=255) aCAGATGCTTACTTACTCGCTGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg < 2:75547/100‑1 (MQ=255) cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa < 1:631972/100‑1 (MQ=255) gCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg < 1:1411509/100‑1 (MQ=255) cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt > 1:1392959/1‑100 (MQ=255) cttactCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc < 1:733709/100‑1 (MQ=255) cttactCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc > 2:318379/1‑100 (MQ=255) ttactCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg < 2:947751/100‑1 (MQ=255) tactCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg < 1:1015183/100‑1 (MQ=255) actCGCGGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg < 1:939592/100‑1 (MQ=255) cgcgGAACAGCGCTTCGATGTTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa > 2:1709003/1‑100 (MQ=255) cgcgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa > 2:170182/1‑100 (MQ=255) cgcgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa > 1:1709270/1‑100 (MQ=255) gcgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGACGAACACGTTCACGGGTGa > 2:155412/1‑99 (MQ=255) gcgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCGTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg < 1:750265/100‑1 (MQ=255) cgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg < 2:32036/100‑1 (MQ=255) cgGAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg > 1:1651006/1‑100 (MQ=255) ggAACAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc < 1:1788362/100‑1 (MQ=255) cAGCGCTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt > 2:95676/1‑100 (MQ=255) gcgcTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt > 1:706361/1‑100 (MQ=255) gcgcTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt > 1:126391/1‑100 (MQ=255) cgcTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTc < 2:1767404/100‑1 (MQ=255) gcTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa < 2:83639/100‑1 (MQ=255) gcTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa < 1:590070/100‑1 (MQ=255) cTTCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc > 2:674858/1‑100 (MQ=255) tCGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCTCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa < 2:1563728/100‑1 (MQ=255) cGATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc < 2:1627045/100‑1 (MQ=255) gATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa > 2:355835/1‑98 (MQ=255) gATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc > 1:1787497/1‑100 (MQ=255) gATATTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc < 2:791027/100‑1 (MQ=255) atatTCAGCCCCC‑TGCGTCTGCAGCAATTCGCGCACACGGCGCAGGCCTTCAACCAGAATTTGGCCCACACGTTCACGGGTGAGGGCAATTTCa > 1:1653567/1‑94 (MQ=255) atatACAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTAAGGCCAATTTCAc < 1:441308/95‑1 (MQ=255) tatTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa < 2:1687314/100‑1 (MQ=255) atTCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc < 2:664667/100‑1 (MQ=255) ttCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa > 1:1178300/1‑100 (MQ=255) ttCAGCCCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa < 2:1800511/100‑1 (MQ=255) | TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CCCC‑TGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA > NC_000913/2866491‑2866682 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |