Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,967,533 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,5330AG97.3% 95.0 / NA 37L301P (CTG→CCG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/16);  minor base C (0/1);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

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

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: