Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 453,912 A→G E344G (GAA→GGA)  clpX → ATP‑dependent protease ATP‑binding subunit ClpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273453,9120AG96.9% 109.7 / ‑3.2 32E344G (GAA→GGA) clpXATP‑dependent protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (15/16);  total (15/17)
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.06e-01

TCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC  >  NZ_CP009273/453799‑454029
                                                                                                                 |                                                                                                                     
tcgtcTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGAc                                                                                              >  1:4357/1‑139 (MQ=255)
  gtcTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACg                                                                                             <  2:401914/138‑1 (MQ=255)
          gTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTg                                                                                     >  2:105510/1‑138 (MQ=255)
             gTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATg                                                                                 >  2:333126/1‑139 (MQ=255)
                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             >  2:37695/1‑139 (MQ=255)
                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             >  1:246555/1‑139 (MQ=255)
                 cAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTAt                                                                             <  1:41552/139‑1 (MQ=255)
                        gaatgaaCTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTaa                                                                       <  1:105510/138‑1 (MQ=255)
                                        gaagCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTa                                                      <  1:37695/139‑1 (MQ=255)
                                            ctctGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAAc                                                  >  2:51463/1‑139 (MQ=255)
                                              ctGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCg                                                >  2:283404/1‑139 (MQ=255)
                                                     aGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGccc                                          <  2:257982/138‑1 (MQ=255)
                                                      gATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGt                                        >  2:238106/1‑139 (MQ=255)
                                                           ccAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCt                                   >  1:322253/2‑139 (MQ=255)
                                                                gagCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGtt                              <  2:424488/139‑1 (MQ=255)
                                                                 agCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTc                             >  1:87075/1‑139 (MQ=255)
                                                                   ccGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGgcg                                                                                        <  1:378301/78‑1 (MQ=255)
                                                                   ccGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGgcg                                                                                        >  2:378301/1‑78 (MQ=255)
                                                                               cTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCg               <  1:238106/139‑1 (MQ=255)
                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  2:246555/138‑1 (MQ=255)
                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:51463/138‑1 (MQ=255)
                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:333126/138‑1 (MQ=255)
                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:398248/138‑1 (MQ=255)
                                                                                  aCCAAGCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  2:382140/138‑1 (MQ=255)
                                                                                       gCAGTATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTc       >  2:52517/1‑139 (MQ=255)
                                                                                           tATCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  2:87075/139‑1 (MQ=255)
                                                                                           tATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  1:314926/139‑1 (MQ=255)
                                                                                            aTCAGGCGCTGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGGGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAc  >  1:115228/1‑139 (MQ=255)
                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             >  2:181540/1‑119 (MQ=255)
                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCa             <  1:181540/119‑1 (MQ=255)
                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   >  1:434576/1‑129 (MQ=255)
                                                                                                     tGTTTAATCTGGGAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATa   <  2:434576/129‑1 (MQ=255)
                                                                                                                 |                                                                                                                     
TCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGCAGTATCAGGCGCTGTTTAATCTGGAAGGCGTGGATCTGGAATTCCGTGACGAGGCGCTGGATGCTATCGCTAAGAAAGCGATGGCGCGTAAAACCGGTGCCCGTGGCCTGCGTTCCATCGTAGAAGCCGCACTGCTCGATAC  >  NZ_CP009273/453799‑454029

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