Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,422,159 A→G N151S (AAC→AGC)  ulaC → L‑ascorbate specific PTS enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,1590AG98.2% 199.1 / ‑5.0 56N151S (AAC→AGC) ulaCL‑ascorbate specific PTS enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/33);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACC  >  NC_000913/4422021‑4422293
                                                                                                                                          |                                                                                                                                      
tGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTA‑CGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGcg                                                                                                                                      >  2:1007481/1‑140 (MQ=255)
  gcggcgGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAgcgc                                                                                                                                     <  1:262821/140‑1 (MQ=255)
   cggcggTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGcgcg                                                                                                                                    <  1:441318/140‑1 (MQ=255)
    ggcggTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCgg                                                                                                                                   <  1:817488/140‑1 (MQ=255)
    ggcggTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCgg                                                                                                                                   <  1:494561/140‑1 (MQ=255)
        gTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGc                                                                                                                               >  1:455472/1‑140 (MQ=255)
              gCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTaaga                                                                                                                         <  1:394286/140‑1 (MQ=255)
                        cccAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGaa                                                                                                               <  2:372204/139‑1 (MQ=255)
                        aCCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGaa                                                                                                               <  2:598105/140‑1 (MQ=255)
                         ccAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAAc                                                                                                              >  1:132725/1‑140 (MQ=255)
                              aaGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTc                                                                                                         <  2:136324/140‑1 (MQ=255)
                                  gggCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTa                                                                                                     >  2:185078/1‑140 (MQ=255)
                                   ggCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTAc                                                                                                    <  1:778240/140‑1 (MQ=255)
                                    gCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTAcc                                                                                                   >  2:184033/1‑140 (MQ=255)
                                    gCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTAcc                                                                                                   >  2:1073762/1‑140 (MQ=255)
                                     catcatGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCg                                                                                                  <  1:277974/140‑1 (MQ=255)
                                        catGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATg                                                                                               >  1:73937/1‑140 (MQ=255)
                                          tGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGtt                                                                                             >  1:301876/1‑140 (MQ=255)
                                               aTCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGTTCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAg                                                                                        >  2:681778/1‑140 (MQ=255)
                                               aTCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAg                                                                                        <  1:285269/140‑1 (MQ=255)
                                                    cAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTcgcg                                                                                   >  2:781388/1‑140 (MQ=255)
                                                     aaCCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCgcgc                                                                                  >  1:96011/1‑140 (MQ=255)
                                                       ccTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTg                                                                                <  2:501374/140‑1 (MQ=255)
                                                       ccTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTg                                                                                <  1:185078/140‑1 (MQ=255)
                                                       ccTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTg                                                                                <  2:829473/140‑1 (MQ=255)
                                                         tGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGa                                                                              <  1:395389/140‑1 (MQ=255)
                                                             tgaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAc                                                                          <  2:132725/140‑1 (MQ=255)
                                                             tgaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAc                                                                          <  1:3731/140‑1 (MQ=255)
                                                             tgaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAc                                                                          <  1:1073762/140‑1 (MQ=255)
                                                              gaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  2:301876/140‑1 (MQ=255)
                                                              gaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  1:781388/140‑1 (MQ=255)
                                                              gaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  1:512254/140‑1 (MQ=255)
                                                              gaagatgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  1:213297/140‑1 (MQ=255)
                                                                  atgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  2:621826/136‑1 (MQ=255)
                                                                  atgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         >  1:621826/1‑136 (MQ=255)
                                                                  atgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         <  2:530385/136‑1 (MQ=255)
                                                                  atgaagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAAcc                                                                         >  1:530385/1‑136 (MQ=255)
                                                                     aagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTAt                                                                  >  1:788653/1‑140 (MQ=255)
                                                                     aagaGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTAt                                                                  >  2:32438/1‑140 (MQ=255)
                                                                             tttGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAgcg                                                          <  2:73937/140‑1 (MQ=255)
                                                                              ttGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAgcgc                                                         <  2:788653/140‑1 (MQ=255)
                                                                                     tttACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGaa                                                  >  2:330604/1‑140 (MQ=255)
                                                                                         cgcgCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCa                                              <  1:184033/140‑1 (MQ=255)
                                                                                                gCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCt                                       <  2:502348/140‑1 (MQ=255)
                                                                                                              gAAGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGC‑AGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGt                        <  1:1007481/140‑1 (MQ=255)
                                                                                                               aaGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGt                        <  1:316351/140‑1 (MQ=255)
                                                                                                                aGTACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTc                       <  1:616926/140‑1 (MQ=255)
                                                                                                                  tACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGa                     <  2:455472/140‑1 (MQ=255)
                                                                                                                  tACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGa                     <  2:526767/140‑1 (MQ=255)
                                                                                                                  tACTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGa                     <  1:316466/140‑1 (MQ=255)
                                                                                                                    cTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACa                   >  2:835888/1‑140 (MQ=255)
                                                                                                                    cTGGATCTCATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACa                   >  2:502427/1‑140 (MQ=255)
                                                                                                                            cATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGaa           <  1:32438/140‑1 (MQ=255)
                                                                                                                            cATTGACCGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGaa           >  1:768439/1‑140 (MQ=255)
                                                                                                                                   cGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCa    >  1:433982/1‑140 (MQ=255)
                                                                                                                                   cGCACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCa    >  1:603729/1‑140 (MQ=255)
                                                                                                                                     cACCAGCGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCAcc  >  2:260218/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACGCGCCTGCCGTACCGAGCAGGAAGTACTGGATCTCATTGACCGCACCAACGCGGCAGCTTAAGAAGGAATTGAACATGTCATTACCGATGTTGCAAGTCGCGCTGGACAACCAGACTATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACC  >  NC_000913/4422021‑4422293

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