| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,867,454 | T→A | Q33L (CAG→CTG) | rpoS ← | RNA polymerase sigma factor RpoS |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,867,454 | 0 | T | A | 96.7% | 91.2 / NA | 30 | Q33L (CAG→CTG) | rpoS | RNA polymerase sigma factor RpoS |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (11/18); minor base C (1/0); total (12/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.00e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 | |||||||||||
GACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTC > NC_000913/2867212‑2867696 | gACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTctact < 1:566479/251‑1 (MQ=255) gATTGCCATAACGGCTGGCAATTTTTACCACCAGACGCAAGGTACTCTCGATCATCCGGCGGCTAGAGGATACATCTCCCCGCAGTGCGCGACGCGCAAAATAAACTTCTGCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCTTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCGCTGCGGCCAAATCGTTATCACTTGGTTCCAGTGCTACTAAggtc < 2:206828/251‑3 (MQ=255) aaCGGCGGGCATATTTTACCACGAGACGCAAGTTACTCTCGATCATACGGCGGCGAGAGGCGACATCTCCACGCAGCGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTcg < 2:4400/251‑2 (MQ=255) gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaa > 2:550949/1‑246 (MQ=255) gggCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCaa < 1:550949/246‑1 (MQ=255) ttttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACt < 2:426123/250‑1 (MQ=255) tACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGt > 2:194326/1‑251 (MQ=255) caGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcat < 1:465327/251‑1 (MQ=255) aatgcaagttactcgtgagccttatgcggagAGAGGCGGCATCTCCAATCAGTGCGCGAAGCGCCAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCCATCTCACCAAAGTAAAGATTATTTGCTTCCAACACACTCTGTGTTGCTCCCTGCGATAACATTTGCTCGTTGGCCAAATCGTTAGCACTGGGTTCCAGTTCTACTAAGGCCTTGTCGTCAAAAACGACATCTCCGTTCtcatca < 2:366147/226‑1 (MQ=255) gTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg < 2:139537/251‑1 (MQ=255) cgatcctacggcggCGAAAGGCGCCATCTCCACGCAGTCCGAGACGCGCAACATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACAAATATCACCAAGGTAAAGCTGAGTAGCGACCAAAACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATAGTTATCACTGGGTTCCAGTTCTACTAAGGAATTTTCGTAAAAAACATCAACTCCGTTCTCATCAAATTACGCATCTTcg > 2:246508/9‑250 (MQ=255) cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCGCATCAAATTCCGCATCt < 1:472429/248‑1 (MQ=255) agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTc < 1:517971/251‑1 (MQ=255) aGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCAt > 1:217389/1‑251 (MQ=255) cgcctccgaccttcttcttCGGCCGTTAACAGTGGTGAATAACCAATCTCACCCAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTGCTCAGCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt < 2:192763/241‑1 (MQ=255) tcttcttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt > 1:192763/1‑239 (MQ=255) cttcGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTg > 1:448587/1‑250 (MQ=255) gtgAATAACCAATCTCACCAAGGTAAAGCTGAGACGCGTCCAACACACGCTGTTTGGCTCCCTGCGATAACAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATGCAAGCaaaa < 1:98933/251‑1 (MQ=255) tgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCATCCCACACACGCTGTGTTGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGCGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTCCCTTGACGGAACATTCAAGCaaa > 2:287308/1‑249 (MQ=255) cAATCTCACACAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGAGATACAAGATCCTATTCGGCCAAATCGTTATCACTGGGTTCCCGTTCTACTAAGGCCTTTTAGTCAAAAACCTCAACTCCGATCTAATCAAATTCCGCATCTTCATTTAAATCATGCACTTTCAGAGTATTCTGCCTCATAAGGCGGCTCCTACCCGTGATCCCTTGCCGGAACATTCAAGCAAAAGCCTGGTTc > 2:382147/1‑251 (MQ=255) cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTGATCGCTGCGGCAAATAACGCAg < 1:444598/251‑1 (MQ=255) aaCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGCTATCACTGGGTTCCAGTTCTAATAAGGCCTTTTCGTCAAAAACCTCAACTCAGTTCTCACCAAATTCCGCATCTTCATTTACATCATGAACTTTCAGCGTATTCTGACTCATCAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGCTTTATCGCTGCGGCAAATAACGCAGc > 1:351551/1‑251 (MQ=255) gCTGTGTGGCTCCCTGCGATAACAGTTCTTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATGCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGGGGCTCCTACCCGGTAGCCCCTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTa < 1:332979/251‑1 (MQ=255) ggCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAttt > 2:491234/1‑251 (MQ=255) gCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACTGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTc < 1:377409/251‑1 (MQ=255) cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGCGTTCCAGTTCTACTCAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGCCTCATAAGGTGGCTCCTACCCGTGTTCCCTTGACGGCACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTcc > 2:65429/1‑251 (MQ=255) cTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAGTTATACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGCTTTCACCtt > 2:309179/1‑251 (MQ=255) aaCAGTTCCTCTGCGGCCAAATCGTTCTCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAATCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACg < 1:538586/250‑1 (MQ=255) tCGGCCAAATCGTTATCACTGGGTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCTGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAAt < 2:223480/249‑1 (MQ=255) gggTTCCAGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTc < 2:39691/250‑1 (MQ=255) | GACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTC > NC_000913/2867212‑2867696 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |