Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 256,444 C→T 10.2% A122T (GCC→ACC)  U712_01195 ← Putative PTS system glucosamine‑specific EIICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881256,4440CT10.2% 110.6 / 3.3 39A122T (GCC→ACC) U712_01195Putative PTS system glucosamine‑specific EIICBA component
Reads supporting (aligned to +/- strand):  ref base C (12/23);  new base T (2/2);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

GCCAAACGAATGAGAAAATGACGCCGATCACCAATGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCGAGGGTGACGGTTAAAATTAAGTATCCGATGACAG  >  CP006881/256309‑256578
                                                                                                                                       |                                                                                                                                      
gCCAAAGCACTGACAAATTGACGCCGATCACCAATGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCgag                                                                                                                                     <  1:260183/139‑1 (MQ=255)
ccaaaaCGATTCACAAAATCCCGCCGATCACCAATGACCTGACGGAAGTAATGATCGGGACAACGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCgag                                                                                                                                     <  2:309370/136‑1 (MQ=255)
                 aTGACGCCGATCACCAATGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATtc                                                                                                                    >  6:33026/1‑139 (MQ=255)
                        cGATCACCAATGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGaa                                                                                                             <  5:287233/139‑1 (MQ=255)
                             acaaaTGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTTATGAAAATCTCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAAccc                                                                                                        <  1:342761/136‑1 (MQ=255)
                                aaTGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCa                                                                                                     <  5:33026/139‑1 (MQ=255)
                                       tGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATatcaat                                                                                              >  7:259634/1‑139 (MQ=255)
                                          cGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAg                                                                                           <  5:307456/139‑1 (MQ=255)
                                             aaGTAATGATCGGGACAAAGCGGTTTCCACATAAAACACCACGGCCCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTg                                                                                        <  6:32653/139‑1 (MQ=255)
                                                  aTGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCgg                                                                                   <  7:38641/139‑1 (MQ=255)
                                                  aTGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCgg                                                                                   <  4:195550/139‑1 (MQ=255)
                                                  aTGATCGGGACAAAGCGTGTTCCAGAGAAAAATCCACGGACCGGATGCAGCTCGATTGATGAAACTCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCgg                                                                                   <  5:70200/139‑1 (MQ=255)
                                                   tGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGc                                                                                  >  7:40839/1‑139 (MQ=255)
                                                     aTCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCt                                                                                 <  7:109138/138‑1 (MQ=255)
                                                         ggACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCtt                                                                             >  7:4475/1‑138 (MQ=255)
                                                               cccGTTTCCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATCGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      <  3:1039/137‑1 (MQ=255)
                                                               aGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      <  7:256750/139‑1 (MQ=255)
                                                               aGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACg                                                                       >  8:256750/1‑138 (MQ=255)
                                                                    tttCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      >  8:372440/1‑134 (MQ=255)
                                                                    tttCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      <  7:372440/134‑1 (MQ=255)
                                                                    tttCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCt                                                                 >  3:378257/1‑139 (MQ=255)
                                                                         agagAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAgg                                                            <  8:40839/139‑1 (MQ=255)
                                                                         acacaaaaaaCGAAGGACCGGATGCAGCTCGAATGATGAAAATCGAATATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAgg                                                            <  1:159643/135‑1 (MQ=255)
                                                                                  tCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAAt                                                   <  3:191939/139‑1 (MQ=255)
                                                                                      aGGACCGGATGCAGCTCGATTGAAGAAAATGGCTTATATAAATACGCGGCGAGAAGCCCGATAAGGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTc                                               <  5:301663/139‑1 (MQ=255)
                                                                                           cGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCata                                          <  8:4475/139‑1 (MQ=255)
                                                                                            ggATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCgg                                                                                   >  8:65389/1‑97 (MQ=255)
                                                                                            ggATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCgg                                                                                   <  7:65389/97‑1 (MQ=255)
                                                                                             gATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTc                                               <  6:102200/132‑1 (MQ=255)
                                                                                             gATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTc                                               >  5:102200/1‑132 (MQ=255)
                                                                                                 cAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCg                                    <  4:378257/139‑1 (MQ=255)
                                                                                                    cTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGt                                                      >  7:198511/1‑118 (MQ=255)
                                                                                                    cTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGt                                                      <  8:198511/118‑1 (MQ=255)
                                                                                                            atgaAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCGAGGGTGACGGt                         >  4:110203/1‑139 (MQ=255)
                                                                                                                         atataAATAGGCGGTGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCat                                           >  7:234877/1‑108 (MQ=255)
                                                                                                                         atataAATAGGCGGTGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCat                                           >  7:229321/1‑108 (MQ=255)
                                                                                                                         atataAATAGGCGGTGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCat                                           <  8:229321/108‑1 (MQ=255)
                                                                                                                         atataAATAGGCGGTGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCat                                           <  8:234877/108‑1 (MQ=255)
                                                                                                                                 aGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCGAGGGTGACGGTTAAAATTAAGTATCCGATGa     >  2:246902/1‑138 (MQ=255)
                                                                                                                                   gcggcgAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCGAGGGTGACGGTTAAAATTAAGTATCCGATGACAg  >  4:113670/1‑139 (MQ=255)
                                                                                                                                    cggcgAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      <  6:272946/70‑1 (MQ=255)
                                                                                                                                    cggcgAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACgg                                                                      >  5:272946/1‑70 (MQ=255)
                                                                                                                                     ggcgAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATatcaatc                                                                                             <  6:189553/46‑1 (MQ=255)
                                                                                                                                     ggcgAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATatcaatc                                                                                             >  5:189553/1‑46 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GCCAAACGAATGAGAAAATGACGCCGATCACCAATGAACTGACGGAAGTAATGATCGGGACAAAGCGTTTTCCAGAGAAAAATCCAAGGACCGGATGCAGCTCGATTGATGAAAATCGCTTATATAAATAGGCGGCGAGAAGCCCGATAATGATTCCTCCGAAAACCCCCATATCAATCAGGTGCTCGGCTCCTTCATACGGAGGCTGAAGGCCGAGTAATTTTCCCATATTGTCGAGGGTGACGGTTAAAATTAAGTATCCGATGACAG  >  CP006881/256309‑256578

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