Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,130,531 A→G I540T (ATC→ACC)  ptsP ← phosphoenolpyruvate‑‑protein phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,130,5310AG100.0% 145.5 / NA 39I540T (ATC→ACC) ptsPphosphoenolpyruvate‑‑protein phosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/17);  total (22/17)

ACCTTAGCGTTATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCAT  >  NZ_CP009273/4130405‑4130649
                                                                                                                              |                                                                                                                      
aCCTTAGCGTTATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTgg                                                                                                             <  2:279380/138‑1 (MQ=255)
aCCTTAGCGTTATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGa                                                                                                            >  1:5008/1‑139 (MQ=255)
            tCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCtcatca                                                                                                 >  2:352634/1‑138 (MQ=255)
            tCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCtcatca                                                                                                 >  1:308677/1‑138 (MQ=255)
               cgATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACg                                                                                             <  2:151132/139‑1 (MQ=255)
                             cagATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTGCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTa                                                                               <  1:101660/139‑1 (MQ=255)
                                      cGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACg                                                                             <  2:477533/132‑1 (MQ=255)
                                      cGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACg                                                                             >  1:477533/1‑132 (MQ=255)
                                      cGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGc                                                                      >  1:117286/1‑139 (MQ=255)
                                       gTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCt                                                                     >  1:206474/1‑139 (MQ=255)
                                           ggTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGttt                                                                 >  1:567636/1‑139 (MQ=255)
                                            gTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTg                                                                >  2:144509/1‑139 (MQ=255)
                                                      ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGc                                                        >  1:131119/1‑137 (MQ=255)
                                                      ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGc                                                        >  1:214777/1‑137 (MQ=255)
                                                      ccAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGc                                                        >  2:314899/1‑137 (MQ=255)
                                                                aGAAATCAATCTCTTCTCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTtttcttt                                            <  2:117286/139‑1 (MQ=255)
                                                                aGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTg                                                                       >  1:484875/1‑112 (MQ=255)
                                                                aGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTg                                                                       <  2:484875/112‑1 (MQ=255)
                                                                 gAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTtttctttc                                           <  2:415966/139‑1 (MQ=255)
                                                                  aaaTCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCa                                          <  1:314899/139‑1 (MQ=255)
                                                                         tctcTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGaa                                                                                          <  1:457673/84‑1 (MQ=255)
                                                                         tctcTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGaa                                                                                          >  2:457673/1‑84 (MQ=255)
                                                                              tCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGAtct                              >  1:8228/1‑139 (MQ=255)
                                                                                    caTTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCa                        <  2:567636/139‑1 (MQ=255)
                                                                                      ttGATCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCAtt                      <  1:238756/139‑1 (MQ=255)
                                                                                         aTCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgag                   >  2:98114/1‑139 (MQ=255)
                                                                                          tCGATGATGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgag                   >  1:19774/1‑138 (MQ=255)
                                                                                            gatgatGAACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgagga                 <  1:144509/138‑1 (MQ=255)
                                                                                              tgatgaACATCACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTgag                   >  2:86472/1‑134 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCtcat                                                                                                   <  1:152732/44‑1 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCtcat                                                                                                   >  2:152732/1‑44 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTtttct                                              >  2:321/1‑97 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTtttct                                              <  1:321/97‑1 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAt     >  1:551842/1‑138 (MQ=255)
                                                                                                        cACCGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGAt     >  1:527791/1‑138 (MQ=255)
                                                                                                          ccGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCAt  <  2:308677/139‑1 (MQ=255)
                                                                                                          ccGACGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCAt  <  2:131119/139‑1 (MQ=255)
                                                                                                             aCGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCg                                                         >  1:549234/1‑81 (MQ=255)
                                                                                                             aCGGCACTTCCAGCATGGTGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCg                                                         <  2:549234/81‑1 (MQ=255)
                                                                                                                              |                                                                                                                      
ACCTTAGCGTTATCGCGATCCACCGCCAGCAGATACTGCGTCAGGTCATTACTACCAATACTAAAGAAATCAATCTCTTCGCAGCATTGATCGATGATGAACATCACCGACGGCACTTCCAGCATGATGCCGAGCTGGATTTTCTCATCAAACGGAATGTGTTCGTTACGTAGTTGCTGTTTGGCTTCCGCCAGTTTTTCTTTCACCCATAAGATCTCTTCCATTGAGGAGATCATCGGGATCAT  >  NZ_CP009273/4130405‑4130649

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