Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,326,537 C→A S133S (TCG→TCT)  ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,326,5370CA92.6% 58.2 / ‑0.8 27S133S (TCG→TCT) ftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (14/11);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTG  >  NC_000913/3326300‑3326783
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
cGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAATTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGACGTTTTGTTCTGATCTTCCGTCAGGCTGCGCGCTTTGCTCTTACCAAAAGACATGGCACCtt                                                                                                                                                                                                                                            >  1:28240/1‑251 (MQ=255)
          gcttcgcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCAccgcc                                                                                                                                                                                                                                  <  1:84033/251‑1 (MQ=255)
               cccTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCt                                                                                                                                                                                                                             <  1:137998/250‑1 (MQ=255)
                    caatctcgtgcGCCAGCAGCGTTTTACCGGTCCCCGGAGGCCCGCCCACCAAGCCGCCTTCCGGGTTCTATCCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCAtt                                                                                                                                                                                                                        <  2:114295/241‑1 (MQ=255)
                                ccAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCCTga                                                                                                                                                                                                              >  1:2318/1‑249 (MQ=255)
                                                         aggaCCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAAc                                                                                                                                                                                   >  1:105559/1‑250 (MQ=255)
                                                             ccgaccATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGGT‑CCAGACACCAATCAGCAACAGCAt                                                                                                                                                                              >  1:177834/1‑251 (MQ=255)
                                                             ccgaccATCAAGACGACTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATAATCAACCAGTTCAACAACTTCTTCTTTTGCTTCGGCGCAGCCCGCCACGTCAGCAAAGGGCGTTTTGATCTGATCTTTCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCTTTTGTCGGATGAAGAA‑ATACCAGACCGCAATCAGCAACAGCAt                                                                                                                                                                              >  2:108269/1‑251 (MQ=255)
                                                                                        ttACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGc                                                                                                                                                   >  1:289371/1‑251 (MQ=255)
                                                                                         tACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAAACCGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCGCTTACCAAAAGACATTGCACCTTTGCCACCGCCGCCCGGCATTTGACGCCTGAAGGAAAT‑CCAGACACCCCTCAGCAACCGCATCCGGAACCAGGGGATgcaggtgggagcc                                                                                                                                                  >  2:88829/1‑239 (MQ=255)
                                                                                               cgAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAgg                                                                                                                                            <  1:251209/251‑1 (MQ=255)
                                                                                                        ggAAGCGGCTCGGCTCGCGCAGATTCTCAACCAGCTCCGCAACCTCTTCTTTTGCTTCGGCGCAGCCCGCCACTTCCGCAAAGGTCGTTTTGATCTTATCTTCCGTCCGCATGCGCGCTTTGCTCTTACCCAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAATGAGAT‑CCAGACACCAATCAGCCACAGCCTCGAGAAACAGGAGATGGAGATCGACGCCAGCCGACTTGGttct                                                                                                                                   >  2:97973/1‑251 (MQ=255)
                                                                                                        ggAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGttct                                                                                                                                   <  1:177662/251‑1 (MQ=255)
                                                                                                          aaGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAACAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTCCCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACCGCCTCGGGAACCAGGAGATGAAGCTAGAAGCCAGCCGGCTTGGtttttc                                                                                                                                 >  2:8199/1‑251 (MQ=255)
                                                                                                                   ggctcgCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTc                                                                                                                        <  1:88829/251‑1 (MQ=255)
                                                                                                                                   cAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACg                                                                                                        >  1:81133/1‑251 (MQ=255)
                                                                                                                                    aaCCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTCCCAAAAGACATGGCACCttt                                                                                                                                                                                                                                           >  2:255552/1‑120 (MQ=255)
                                                                                                                                    aaCCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCtgt                                                                                                                                                                                                                                           <  1:255552/120‑3 (MQ=255)
                                                                                                                                      ccAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGACTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTc                                                                                                     <  1:280015/251‑1 (MQ=255)
                                                                                                                                         gTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGGTCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGACAGCAGGCTTGGTTCTTCAGGAGGTTAACCGACAACCTTGACGATCTTg                                                                                                  >  1:210221/1‑251 (MQ=255)
                                                                                                                                                          tttGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTCCCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGCGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCATGCCGTTCTTGGTCAACAGGTTATcc                                                                                   >  2:101974/1‑249 (MQ=255)
                                                                                                                                                           ttGCGTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGAGCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAg                                                                                 <  1:3266/250‑1 (MQ=255)
                                                                                                                                                                tcgtcgCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAAt                                                                             <  1:8199/249‑1 (MQ=255)
                                                                                                                                                                                        aaGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGTAGAT‑CCAGACACCAATCAGCAACAGCATCTGGAACCAGGAGATGAAGATACAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGCCAACCTTTACGTTCTTGGTCAACAGGTTATCCGATAATTTCGGATCCTGAACCGCAATGTAAGt                                                   >  1:125661/1‑251 (MQ=255)
                                                                                                                                                                                             cGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAAGACATGGCACCTTGGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATATAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGtat                                              <  2:241462/251‑1 (MQ=255)
                                                                                                                                                                                                                        gcgctcattGCCCTTTCCAAAAGACACAGCACCTTTGCCACCGCCGCCCTTCATATGGCGCATGAAGAAGAT‑CCAGACACCAAGCAGCAACAGCATCTGGAACCAGGAGATGAATATATAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTGCTTGGTCAACAGGTTTTCCACTAAGTTCGGATCCTGAACCGTAATGTAAGTGTTATCTCGGTCTCTATCTTTCGTgggaccgt                   <  2:210221/244‑7 (MQ=255)
                                                                                                                                                                                                                                         cAAAAGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTg  >  1:63128/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       
CGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGAT‑CCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTG  >  NC_000913/3326300‑3326783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑