| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,708,714 | C→A | 11.1% | S14I (AGC→ATC) | rseB ← | anti‑sigma factor stabilizing protein RseB |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,708,714 | 0 | C | A | 11.1% | 117.5 / 5.3 | 36 | S14I (AGC→ATC) | rseB | anti‑sigma factor stabilizing protein RseB |
| Reads supporting (aligned to +/- strand): ref base C (13/19); new base A (2/2); total (15/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01 | |||||||||||
GCCATTAAGCGTGAACGGTTCAAGTCCCGGTTCAAAATAGCTGATTTCATTGCCGCGCTGTACCACTTCCCGGCGCGGGCCATCCATTTGCAACAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGCGACGACGCTGCTCCTGTACCTGCTGCTGTTGACCATTGTTTGCGGTCGCTTCAGAAGGTACTCCCAGGCTTACCGGGCTGGCTT > NC_000913/2708443‑2708974 | gcatttAAGCGTGAACGGTTCAAGTCCCGGTTCAAAAGGGTCGATTCCATGGCCGCGCTGTACCACTTCCCGGCGCGGGCCATCCATTTGCAAAAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGc < 1:455850/268‑1 (MQ=255) aaCGGTTCAAGTCCCGGTTCAAAATAGCTGATTTCATTGCCGCGCTGTACCACTTCCCGGCGCGGGCCATCCATTTGCAACAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGACTGACAGCTCGTACTTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGCCGCGGGAGTGGCCGAGGCGTTAGCAGAGAATACCAGGCTCCCTGTCACTa > 1:289967/1‑271 (MQ=255) cGGCGCGGGCCATCCATTTGCAACAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAAtt < 1:188205/272‑1 (MQ=255) ggCGCGGGCCATCCATTTGCAACAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTc > 1:117237/1‑272 (MQ=255) aacaaTTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCTGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGATACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGc > 1:170231/1‑271 (MQ=255) tgCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAAt > 2:331587/1‑242 (MQ=255) gCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAAtt < 1:331587/242‑1 (MQ=255) aaGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTg < 1:253120/272‑1 (MQ=255) gagGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAAAAACAGGCTAc > 2:74946/1‑171 (MQ=255) agGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAAGAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAAAAACAGGCTAcc < 1:74946/171‑1 (MQ=255) ggggcgggaaggTCGTTAACGCAGGGACTCAACACCCTGTTTTTTGAGGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTg < 2:117237/261‑1 (MQ=255) ggcaattccgttAACGCAGAGATTCAAAACCCTGTTTATTGATGCTGTTGAATGACAGGTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGATACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑Ct < 2:170231/261‑1 (MQ=255) agaCTCAACTTCCTGTTTATTGATGCTTATGAATGACAGGTCGTAATTCAGTGACTGACTGGCCAAGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGGGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCg < 1:341111/272‑1 (MQ=255) acaggTTTATTTATGCTTTTGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGTTCg < 1:57760/268‑1 (MQ=255) tttATTGTTGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCCCAAACTGAAGCCTTTACAGAG‑GGGGACGGGGTTGCCGGTCGTaaaa > 2:396231/1‑270 (MQ=255) ttATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCt > 2:145771/1‑206 (MQ=255) tATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGAc > 1:135180/1‑131 (MQ=255) tATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCtt < 1:145771/206‑1 (MQ=255) tATTGATGCTGATGAATGACAGATCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGATACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTg > 1:174999/1‑188 (MQ=255) aTTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACa < 2:135180/131‑1 (MQ=255) aTTGATGCTGATGAATGACAGATCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGATACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTgg < 2:174999/188‑1 (MQ=255) atgaatgaCAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTGGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTTCGCCTGACCAACCTGAAG‑CTTTTCAGAGGGGAGCAGGACTTACACTTGGTCAACCTGCAACAGTGCa > 2:92943/1‑272 (MQ=255) aGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTAACTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTg < 1:185929/211‑1 (MQ=255) tCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGcgacga < 1:236134/272‑1 (MQ=255) gTAATTCAGTTACTGACTGGCCAGGTTCATCTGCGGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGAAGAGAATAACAGTCTACCTATCACTAATGACATGGCAAACCACAGTTGCTTCATGACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑Ct < 1:406020/211‑1 (MQ=255) gTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGtt < 1:268592/238‑1 (MQ=255) gactgactGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGACAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCg > 2:465783/1‑126 (MQ=255) actgactgGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGACAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGa < 1:465783/126‑1 (MQ=255) gctgTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACt > 1:388829/1‑173 (MQ=255) gctgTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGt > 1:273942/1‑199 (MQ=255) gctgTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGt < 2:273942/199‑1 (MQ=255) ctgTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTg < 2:388829/173‑1 (MQ=255) aataaCGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCt > 2:115707/1‑143 (MQ=255) aataaCGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGCGTGGAGTCGGCGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGAGACGACGCTGCTCCTTTACCTGCTGCTGTTGACAATTGTTTTCg > 1:189004/1‑272 (MQ=255) ataaCGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCtt < 1:115707/143‑1 (MQ=255) ggACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTAACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGGGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGCGACGACGCTGCTCCTGTACCTGCTGCTGTTGACCATTGTTTGCGGTCGCTTAAGa < 1:92943/272‑1 (MQ=255) agaATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAATTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGCGACGACGCTGCTCCTGTACCTGCTGCTGTTGACCATTGTTTGCGGTCGCTTCAGAAGGTACTCCCAGGCTTACCGGGCTGGCtt < 1:425040/272‑1 (MQ=255) | GCCATTAAGCGTGAACGGTTCAAGTCCCGGTTCAAAATAGCTGATTTCATTGCCGCGCTGTACCACTTCCCGGCGCGGGCCATCCATTTGCAACAATTGTGCAAGAGGACGGTTATCGAGGCGTGCATGTCGATAACGCAGAGACTCAACACCCTGTTTATTGATGCTGATGAATGACAGCTCGTAATTCAGTGACTGACTGGCCAGGTTCATCTGCTGTAATAACGCCCCGGACGCGGGAGTGGCCGAGGCGTTAGCAGAGAATAACAGGCTACCTGTCACTAATGACATGGCAAACCAAAGTTGCTTCATTACTGCGATTGCGTTCCTAAAGTTTGAATTCCTGGCACCTGTACAGCGGCTTGCTGGGTTTGCGCCTGCTCAAACTGAAG‑CTGTTCAGAGTGGAGTCGGCGTTGCAGTTCGTAATCCTGCAACATTGCATTAATGCGACGACGCTGCTCCTGTACCTGCTGCTGTTGACCATTGTTTGCGGTCGCTTCAGAAGGTACTCCCAGGCTTACCGGGCTGGCTT > NC_000913/2708443‑2708974 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |