Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,960,176 G→A R540W (CGG→TGG)  ptsP ← phosphoenolpyruvate‑‑protein phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,960,1760GA100.0% 140.3 / NA 35R540W (CGG→TGG) ptsPphosphoenolpyruvate‑‑protein phosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/20);  total (15/20)

ACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGAC  >  NZ_CP009273/2960042‑2960311
                                                                                                                                      |                                                                                                                                       
aCCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCacac                                                                                                                                      <  2:576284/138‑1 (MQ=255)
 ccGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTACCGGCATTAGCACGCAGCATCGCCCACACCt                                                                                                                                    <  1:506915/139‑1 (MQ=255)
    atcatcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGa                                                                                                                                 >  1:236416/1‑139 (MQ=255)
       atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCa                                                                                                                              <  2:236416/139‑1 (MQ=255)
               gACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGAtc                                                                                                                      >  2:610204/1‑139 (MQ=255)
                aCTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGAtct                                                                                                                     >  2:621126/1‑139 (MQ=255)
                          cGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATc                                                                                                           <  2:97137/139‑1 (MQ=255)
                          cGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGACCACACCTGGATCAAGAAGATCTCCGGCTGAt                                                                                                            <  1:555242/138‑1 (MQ=255)
                                gTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGt                                                                                                     <  2:86433/139‑1 (MQ=255)
                                           cggcgTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATcc                                                                                          >  1:371906/1‑139 (MQ=255)
                                                tGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGc                                                                                     >  1:178238/1‑139 (MQ=255)
                                                  cTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCa                                                                                   <  1:610204/139‑1 (MQ=255)
                                                    tcatcaACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAAc                                                                                 >  1:108565/1‑139 (MQ=255)
                                                               aTCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGa                                                                      >  1:228738/1‑139 (MQ=255)
                                                                         gtgACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCg                                                             <  1:580902/138‑1 (MQ=255)
                                                                         gtgACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCg                                                             <  1:163419/138‑1 (MQ=255)
                                                                              cATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATtct                                                                  >  1:311121/1‑128 (MQ=255)
                                                                              cATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATtct                                                                  <  2:311121/128‑1 (MQ=255)
                                                                               aTCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCg                                                      <  2:637517/139‑1 (MQ=255)
                                                                                tCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACAcc                                                                                                                                     <  1:116240/59‑1 (MQ=255)
                                                                                tCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACAcc                                                                                                                                     >  2:116240/1‑59 (MQ=255)
                                                                                 cGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCg                                                                                               <  2:30710/96‑1 (MQ=255)
                                                                                 cGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCg                                                                                               >  1:30710/1‑96 (MQ=255)
                                                                                 cGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCgg                                                     >  2:403215/1‑138 (MQ=255)
                                                                                          aaTATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGc                                           <  1:376822/139‑1 (MQ=255)
                                                                                          aaTATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGc                                           <  1:289267/139‑1 (MQ=255)
                                                                                                cAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgctg                                      <  2:103523/138‑1 (MQ=255)
                                                                                                  ggTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAgctgct                                    >  2:60216/1‑138 (MQ=255)
                                                                                                  ggTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCtt                                   <  2:529652/139‑1 (MQ=255)
                                                                                                    ttGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGc                                           >  2:336945/1‑129 (MQ=255)
                                                                                                    ttGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGc                                           <  1:336945/129‑1 (MQ=255)
                                                                                                    ttGCCCGTAGCGGCATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTAt                                 <  2:108565/139‑1 (MQ=255)
                                                                                                                 cATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATc                    >  2:600086/1‑139 (MQ=255)
                                                                                                                 cATTAGCACGCAGCATCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATc                    >  1:209265/1‑139 (MQ=255)
                                                                                                                           cagcaTCGCCCACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGc          <  1:485518/139‑1 (MQ=255)
                                                                                                                                   cccACACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGAc  <  1:621126/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
ACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGAC  >  NZ_CP009273/2960042‑2960311

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