| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,867,455 | G→A | Q33Y (CAG→TAT) | 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,455 | 0 | G | A | 95.5% | 50.1 / NA | 22 | Q33Y (CAG→TAT) | rpoS | RNA polymerase sigma factor RpoS |
| Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (12/9); minor base C (1/0); total (13/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01 | |||||||||||
AACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGG > NC_000913/2867227‑2867648 | aaCGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTc < 1:68007/250‑1 (MQ=255) tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTc < 1:154328/236‑1 (MQ=255) ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTGGGTTCATATGCTACTAAGGc < 1:82542/223‑1 (MQ=255) natcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATACCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAttt > 1:72926/2‑250 (MQ=255) ggCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg > 1:213209/1‑224 (MQ=255) naCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACAAAGGTAAAGCTGCGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAACTCGTTATCACTGGGTTCATATTCTACTAAgg > 1:52458/2‑179 (MQ=255) gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCATGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGt > 1:250886/1‑251 (MQ=255) cgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTg > 1:41909/1‑251 (MQ=255) aaaCTTCTTCTTCGGACGTTAACAGTGGTGAATAACCAATCTCACCAACGGAAAGCTGAGCCGCGTCCAACAAAAGCTGTCTGGCTCCCTCCGATACCAGCTCCCCTTAGGCCAAATCGTTATCACTGGGGTCATATTCTACTAAGGCCTTTTCGTCacac > 1:40376/1‑158 (MQ=255) ccAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAAc > 1:108222/1‑119 (MQ=255) ncAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGCTCCTATTCTACTAACGACTTTTCGTCCAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCACAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGc > 1:32245/2‑250 (MQ=255) cAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCCAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc > 1:278868/1‑138 (MQ=255) nacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGACAAATCGTTTTCACTGGGTTCATATTCTCCTAAGGACCTTTCGTCAAACACCTCAACTCCGTTCTCATCAAATTCCGCATCTTAATTTACATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGCTCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCAGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGt > 1:31454/2‑251 (MQ=255) cGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTTGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATaaag < 1:268859/241‑3 (MQ=255) gtgtgGCTCCCTGCGATAACAGTTCATCTTTGGCCCAATCGTTATCACTGGGTTCATCTTCTACTACCGCCATTTCGTCAACAACCTCAACTCCGTTCTCATCAAATTCCGCACCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGCTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACgg > 1:116835/1‑251 (MQ=255) cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCtttt < 1:9916/69‑1 (MQ=255) cTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCATATTCTACTAAGGCCtttt > 2:9916/1‑69 (MQ=255) ncNGTNCCTCTNCGGCNAAATCGTTATCACTGNGTTCATATTCTACTAAGGCNNTNGCGGCNAAAAc < 2:108222/66‑1 (MQ=255) cAAATCGTTATCGCTGGGTTCATATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAAGTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCGGACTCATAAGGTGGCTCCTACCCGTTATCCCTTGACGGAACATTCAAGCAAAAGCCTGGGTCCGCCGATTTATc < 1:38330/187‑1 (MQ=255) nnnnnncANNNNNTTCNTATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATATAAATCATGAACTTTCAGCGTATTCTGACTCATAATGTGGCTCCTACCCGTGATCCCTTGACGGAACa < 2:177444/144‑1 (MQ=255) tCACTGGGTTCATATTCTACAAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCAGCTTCATTTAAATCATGAACTTTCAGCGGATTCTGTCTCATAAGGTGGCTCTTACCCGTTATCCCTTTACGGAACATTCAAGCAAAAGCCTGGTTccgcc < 1:198327/169‑1 (MQ=255) tCACTGGGNTNNNANNCTACTAAGGCCTTTTCGTCCAAAACCTCAACACCGTTCTCATCANATTACGAATCTGCATNTANATCANGATNTTTCAGCGTCTTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgca > 2:198327/1‑168 (MQ=255) | AACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGG > NC_000913/2867227‑2867648 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |