Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,531 A→G L301P (CTG→CCG)  ptsP ← PEP‑protein phosphotransferase enzyme I, GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,5310AG97.6% 105.3 / NA 41L301P (CTG→CCG) ptsPPEP‑protein phosphotransferase enzyme I, GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/17);  minor base C (0/1);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT  >  NC_000913/2967392‑2967663
                                                                                                                                           |                                                                                                                                    
gACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCATGATGATTCAGCAATCGCTGACCCGgcg                                                                                                                                   >  2:627275/1‑143 (MQ=255)
   aacaCCGACTAAGCGATCCTGGGGCAGCCCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCAc                                                                                                                                <  2:302260/143‑1 (MQ=255)
         gACTGAGCGAGCCTGGGGCAGCTCAGCAGGCGGTGCCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAAt                                                                                                                          <  1:623012/143‑1 (MQ=255)
           cTAAGCCATCCTGGGGCATCTCAACAAGCGGCGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCg                                                                                                                        <  1:186942/143‑1 (MQ=255)
               gCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAg                                                                                                                    >  1:411645/1‑143 (MQ=255)
                   tCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGCTCGGCCCTTGATTAGCGCCATCAAGCTGAAACAGCAATCGCTGACCCGGCGCCCCGAACTCGCCAGCCCg                                                                                                                >  1:260499/1‑143 (MQ=255)
                      tGGGGCAGCTCAGCAAGCGTTGTCGTTGACAGTTCATCTGCTACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGCTGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGctc                                                                                                             <  2:474377/143‑1 (MQ=255)
                         ggcgggtccgcgaGCGTTGTGGCTGACCGTTCGTCTGCCACCAGTATGAAACGTTCCTGCCAGGCGTTCGGCCCTTGACTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCttt                                                                                                          <  1:166643/131‑1 (MQ=255)
                              cTCAGCAAGCGTTGTCGCTAACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTGCGGCCCCTGATTAGCGTCATCAAGATGAAACAGAAAACGCTGACCCGGCGCACGCAAATCGCCAGCGCGCGCTTCGAGAt                                                                                                     >  1:553546/1‑143 (MQ=255)
                                  gCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGCGAtgctt                                                                                                 >  2:692502/1‑139 (MQ=255)
                                    aaGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTAGGCCCTAGATTAGCGTCATCAAGCTGAAACAGCAATCACTGAATCGACGCACGTAAATCGCCCGCTCGCTCCTTGAGATactcct                                                                                               >  1:315295/1‑138 (MQ=255)
                                               cggccAGTTCATCCGCCACCAGAACGACACCTTCCGGCCAGGCGTCCGGCCCTTGATTAGCCTCATCAAGATGAAACAGCAACCGCCGCCCCCGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGcgcg                                                                                    <  2:251698/139‑1 (MQ=255)
                                                tGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCgg                                                                                   >  2:16522/1‑143 (MQ=255)
                                                 gACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATACAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGGCCGCTCTGTGAGATAGTGATCTCttcggtcggc                                                                                  >  1:248633/1‑135 (MQ=255)
                                                       tCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCCAACTg                                                                            >  2:691294/1‑143 (MQ=255)
                                                             gCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGCGATAGTTATCACTTAGCGCGGCAAACTGTTCGGc                                                                      <  2:137057/143‑1 (MQ=255)
                                                             gCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTCACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGCTAGGTATCACTCAGCGCGGCCAACTGTCCGGc                                                                      >  1:562517/1‑143 (MQ=255)
                                                              ccaccaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGGTAGTGATCACTTAGCGGGGCAAACTGTGCGGCa                                                                     >  1:103173/1‑143 (MQ=255)
                                                               caccaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCTAACTGTTCGGCaa                                                                    >  2:672628/1‑143 (MQ=255)
                                                               caccaGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCACACTGTTCGGCac                                                                    >  2:299678/1‑142 (MQ=255)
                                                                accaGAATGAAACGTTCCGGCCAGGCGTTCGGACCTTAATTAGCGTCATCCAGATGAAACAGCAATCCCTGACCCGGCGCACGTCACTCGCCAGCCCGCTCTTTCAGATACTTACCACCTAGCGCGCCCACCTGTTCCGCaca                                                                   >  1:57715/1‑141 (MQ=255)
                                                                      aTGAAAGGTTCTGGCCAGCCGTGCGGCCCTTGATTAGCGTCAGCAAGATGAAACAGCAGTCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTc                                                             <  2:256831/143‑1 (MQ=255)
                                                                        gAAACGTGCCGGCAAGGCGTTCGGCCCTTGATTAGCGTCATCAAGAGGAAACAGCAATCGCTGTCCCGGCGCCCGTAAATCGCCAGCCCCCTCTTTGAGATAGTTATCACTTCGCGCGGCAAACTGTTCCGCAAATTTTTCaa                                                           <  2:390774/143‑1 (MQ=255)
                                                                               tCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCCCGTAAATAGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAAGTTTTCAATGACCGt                                                    >  2:57692/1‑143 (MQ=255)
                                                                                  ggCCAGGCGTTGGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGtttt                                                 <  2:139606/143‑1 (MQ=255)
                                                                                    ccAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTa                                               >  1:293915/1‑143 (MQ=255)
                                                                                    ccAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTa                                               <  2:21936/143‑1 (MQ=255)
                                                                                     cAGGCGTCCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAAGCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCCAACTGTCAGTCAAATTTTTCAATGACCGTTTTTAc                                              <  1:627275/143‑1 (MQ=255)
                                                                                        gCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGCTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTATTCAATACCCGTTTTTACCGc                                           >  2:405473/1‑143 (MQ=255)
                                                                                          gTTCGGCCCTTGATTAGTGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAATTCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCGAATGTTTCGATGACCGgttgtccggctc                                         >  2:234063/1‑134 (MQ=255)
                                                                                                  cTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGCGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGcc                                 >  2:559455/1‑143 (MQ=255)
                                                                                                   ttGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCa                                >  2:312246/1‑143 (MQ=255)
                                                                                                        gaGCGTCATTGAGATGAAACAGCAATCGCTGACCCGGCGCGCGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAg                           <  1:157148/142‑1 (MQ=255)
                                                                                                              catcaACATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTa                     <  2:411645/143‑1 (MQ=255)
                                                                                                                tcgaggtGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATc                   <  1:568867/137‑1 (MQ=255)
                                                                                                                   aGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCACTGACCGTTTTTACCGCCCCCTCTGCCACCGAGCCTTTATCAAc                >  1:282514/1‑143 (MQ=255)
                                                                                                                   aGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAAc                <  1:16522/143‑1 (MQ=255)
                                                                                                                   aGATGAAACAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAAc                >  1:657572/1‑143 (MQ=255)
                                                                                                                         aaCAGCAATCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACTGTTTTTACCGCCCCCTCTGCCACCGAGCCTTGATCacgcgcggc          >  2:228141/1‑135 (MQ=255)
                                                                                                                            acccaTCGCGGCCCCGGCGCACGTAAATTGCCAGCCCGCTCTTTGAGATACTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCaaa       <  2:200173/139‑1 (MQ=255)
                                                                                                                                 tCGCTGACCCGGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAAtt  <  2:293915/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
GACAACACCGACTAAGCGATCCTGGGGCAGCTCAGCAAGCGTTGTCGCTGACAGTTCATCTGCCACCAGAATGAAACGTTCCGGCCAGGCGTTCGGCCCTTGATTAGCGTCATCAAGATGAAACAGCAATCGCTGACCCAGCGCACGTAAATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATT  >  NC_000913/2967392‑2967663

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