Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,497,341 (C)5→4 coding (624/2496 nt) ptsP ← phosphoenolpyruvate‑‑protein phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,497,3370C.92.5% 154.7 / 7.2 41coding (628/2496 nt)ptsPphosphoenolpyruvate‑‑protein phosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base . (21/16);  total (23/18)
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.84e-01

CGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTT  >  NZ_CP009273/2497200‑2497466
                                                                                                                                         |                                                                                                                                  
cgGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCgccgcc                                                                                                                                   >  2:96516/1‑139 (MQ=255)
cgGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCgccgcc                                                                                                                                   <  1:1270377/139‑1 (MQ=255)
cgGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCgccgcc                                                                                                                                   >  2:781553/1‑139 (MQ=255)
 ggCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGccc                                                                                                                                  >  1:1264264/1‑139 (MQ=255)
 ggCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGccc                                                                                                                                  >  2:425568/1‑139 (MQ=255)
   cTTCAGTCCCGGCCATGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGATAACTCGTTTGCGCAGATAATCGCTGGCTGACGCAGAAAGTTTGGCGCAAACCTGCTCCATATTCCTGATGAGCGCCGcccca                                                                                                                                <  2:663729/139‑2 (MQ=255)
      cAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGTTGATCGCCGCCCCCAGc                                                                                                                             >  2:271866/1‑139 (MQ=255)
             ggCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     <  1:27193/139‑1 (MQ=255)
             ggCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     <  2:531885/139‑1 (MQ=255)
             ggCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     >  1:531885/1‑139 (MQ=255)
             ggCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     >  2:969110/1‑139 (MQ=255)
                  ggTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGtt                                                                                                                <  2:415464/139‑1 (MQ=255)
                  ggTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGtt                                                                                                                <  1:951253/139‑1 (MQ=255)
                   gTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTc                                                                                                               >  2:194103/1‑139 (MQ=255)
                             gcaACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGccc                                                                                                                                  <  2:1280526/111‑1 (MQ=255)
                             gcaACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGccc                                                                                                                                  >  1:1280526/1‑111 (MQ=255)
                             gcaACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     >  2:347771/1‑123 (MQ=255)
                             gcaACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATg                                                                                                                     <  1:347771/123‑1 (MQ=255)
                                aCTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGa                                                                                                 <  1:194103/139‑1 (MQ=255)
                                aCTGGTCGCTGATATCGCGAATGTCGATAGCACGTTCGCGCAGAGAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGGTTATGATCGCGGCCACAGCCCCTGATGCTGTTCTGTCATCAGGGGa                                                                                                 <  2:674418/139‑1 (MQ=255)
                                    ttCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACgg                                                                                              <  1:999075/138‑1 (MQ=255)
                                    ttCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGa                                                                                             >  1:761440/1‑139 (MQ=255)
                                                 cgAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCaaa                                                                                >  1:443941/1‑139 (MQ=255)
                                                       tCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcat                                                                           <  1:425568/138‑1 (MQ=255)
                                                                  ttCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAg                                                                                                     >  1:901182/1‑101 (MQ=255)
                                                                  ttCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAg                                                                                                     <  2:901182/101‑1 (MQ=255)
                                                                   tCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGa                                                                                                 <  1:766996/104‑1 (MQ=255)
                                                                   tCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGa                                                                                                 >  2:766996/1‑104 (MQ=255)
                                                                    cgcgCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGc                                                             >  2:449764/1‑139 (MQ=255)
                                                                                     gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCgccgcc                                                                                                                                   <  1:705567/54‑1 (MQ=255)
                                                                                     gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCgccgcc                                                                                                                                   >  2:705567/1‑54 (MQ=255)
                                                                                     gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGtt                                                                                                                <  1:1203411/72‑1 (MQ=255)
                                                                                     gcagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGtt                                                                                                                >  2:1203411/1‑72 (MQ=255)
                                                                                      cagaggcagaAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAgg                                           <  2:761440/139‑1 (MQ=255)
                                                                                                         gcAAACCTGCTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGt                        <  1:393564/139‑1 (MQ=255)
                                                                                                                  cTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtcacgctca                >  2:612957/1‑139 (MQ=255)
                                                                                                                  cTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcat                                                                           >  1:970570/1‑79 (MQ=255)
                                                                                                                  cTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATtcat                                                                           <  2:970570/79‑1 (MQ=255)
                                                                                                                  cTCCATATTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtcacgctcac               >  2:930779/1‑139 (MQ=255)
                                                                                                                       tatTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAAt          >  1:198880/1‑139 (MQ=255)
                                                                                                                       tatTGCTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGAAAGTTTTGATTTCGCCGTCACGCTCACGCAAt          >  1:1160578/1‑139 (MQ=255)
                                                                                                                            cTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTg        <  1:1251901/138‑1 (MQ=255)
                                                                                                                            cTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTg     >  1:854354/1‑139 (MQ=255)
                                                                                                                            cTGATGATCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTg     >  1:616820/1‑139 (MQ=255)
                                                                                                                               atgatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtc                       >  1:80797/1‑118 (MQ=255)
                                                                                                                               atgatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGtc                       <  2:80797/118‑1 (MQ=255)
                                                                                                                               atgatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGtt  >  1:530921/1‑139 (MQ=255)
                                                                                                                                 gatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGc       >  2:907057/1‑132 (MQ=255)
                                                                                                                                 gatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGc       <  1:907057/132‑1 (MQ=255)
                                                                                                                                 gatCGCCGCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGtt  >  2:831771/1‑137 (MQ=255)
                                                                                                                                 gatCGCCGCCCCAGCCCCTGATGCTGTTCTGCCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCGCACGCAATTGCTGGtt  >  2:546518/1‑137 (MQ=255)
                                                                                                                                    cgccgcCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                                                    >  1:873457/1‑84 (MQ=255)
                                                                                                                                    cgccgcCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATgg                                                    <  2:873457/84‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CGGCTTCAGTTCCGGCCAGGTGATATGCAGCAACTGTTCGCTGATATCGCGAATGTCGCTAACACGTTCGCGCAGATAATCGCTGGCAGAGGCAGAAAGTTTGGCGCAAACCTGCTCCATATTGCTGATGATCGCCGCCCCCAGCCCCTGATGCTGTTCTGTCATCAGGCGACGGATATTGCCTGCAAATTCATCATCCTGAATCAGCGACAAATGGGCGCTGAGGATAGTTTTGCTTTCGCCGTCACGCTCACGCAATTGCTGGTT  >  NZ_CP009273/2497200‑2497466

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