Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,045 G→A 6.8% intergenic (+164/‑414) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,0450GA 6.8% 193.7 / 3.9 59intergenic (+164/‑414)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGA  >  NC_000913/374899‑375192
                                                                                                                                                   |                                                                                                                                                   
gTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGaa                                                                                                                                                   >  2:68066/1‑149 (MQ=255)
    aTTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                                                                                                                               <  2:196210/149‑1 (MQ=255)
     ttttGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt                                                                                                                                              <  2:29494/149‑1 (MQ=255)
     ttttGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCt                                                                                                                                              >  2:1019514/1‑149 (MQ=255)
              cccATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                     >  3:383512/1‑149 (MQ=255)
                 aTCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  <  1:1019514/149‑1 (MQ=255)
                  tCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                                                                                                                                 >  1:592792/1‑149 (MQ=40)
                     cAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGAt                                                                                                                              <  1:68066/149‑1 (MQ=40)
                         aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGgcg                                                                                                                          >  2:118721/1‑149 (MQ=39)
                           gACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGcgcg                                                                                                                        <  3:392508/149‑1 (MQ=39)
                            aCTGCCGGAGGCGGCGTGAACGCCTTATCCCGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCgg                                                                                                                       <  2:247062/149‑1 (MQ=255)
                             cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCggg                                                                                                                      <  4:1130949/149‑1 (MQ=40)
                                  ggATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAAttt                                                                                                                 <  2:592792/149‑1 (MQ=255)
                                    aTGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                               >  1:774231/1‑149 (MQ=40)
                                    aTGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                                               <  4:383512/149‑1 (MQ=39)
                                     tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                              <  3:630631/149‑1 (MQ=39)
                                     tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                              >  3:530387/1‑149 (MQ=39)
                                     tGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                              <  3:350407/149‑1 (MQ=40)
                                            tGAACGCCTTATTCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGa                                                                                                       <  1:633163/149‑1 (MQ=255)
                                             gAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                      <  1:118721/149‑1 (MQ=40)
                                                      aTCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGt                                                                                             <  1:361498/149‑1 (MQ=39)
                                                        ccGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTa                                                                                           <  2:774231/149‑1 (MQ=40)
                                                          ggCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                                                                                                                               <  1:951676/95‑1 (MQ=40)
                                                          ggCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                                                                                                                               >  2:951676/1‑95 (MQ=40)
                                                             cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGCCCTGATa                                                                                                     >  1:101011/1‑134 (MQ=255)
                                                             cTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGCCCTGATa                                                                                                     <  2:101011/134‑1 (MQ=255)
                                                                 ggATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                          >  2:950312/1‑93 (MQ=255)
                                                                 ggATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                          <  1:950312/93‑1 (MQ=255)
                                                                  gATGGCGTAAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                                                 >  3:344327/1‑149 (MQ=255)
                                                                   aTGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCa                                                                                <  4:267288/149‑1 (MQ=38)
                                                                      gCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                                             >  1:968303/1‑149 (MQ=37)
                                                                         tGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  >  2:630330/1‑93 (MQ=37)
                                                                         tGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                                                                                                                                  <  1:630330/93‑1 (MQ=37)
                                                                         tGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                               <  4:60102/144‑1 (MQ=35)
                                                                         tGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                                               >  3:60102/1‑144 (MQ=35)
                                                                           agaATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                              >  2:553250/1‑111 (MQ=35)
                                                                           agaATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTa                                                                                                              <  1:553250/111‑1 (MQ=35)
                                                                                    aGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACt                                                               <  2:968303/149‑1 (MQ=37)
                                                                                     ggTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTg                                                              >  2:452348/1‑149 (MQ=37)
                                                                                           aTAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGcc                                                            >  2:1117671/1‑145 (MQ=34)
                                                                                           aTAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGcc                                                            <  1:1117671/145‑1 (MQ=34)
                                                                                              aGACGCGTTAAGCGTCGCATCAGGCGTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATgc                                                      <  1:826432/149‑1 (MQ=255)
                                                                                              aGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                          >  1:464366/1‑96 (MQ=34)
                                                                                              aGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                                          <  2:464366/96‑1 (MQ=34)
                                                                                                    gTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCCTCTGCGCACGACTGCCGGATGCGGCGTGa                                               >  2:400830/1‑149 (MQ=37)
                                                                                                                       aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                      <  3:104863/76‑1 (MQ=34)
                                                                                                                       aTCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                                      >  4:104863/1‑76 (MQ=34)
                                                                                                                         cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                          <  3:236298/134‑1 (MQ=37)
                                                                                                                         cTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGcc                                          >  4:236298/1‑134 (MQ=37)
                                                                                                                               aCGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGgcg                     >  2:235554/1‑149 (MQ=40)
                                                                                                                                        ggATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAAttt            <  1:235554/149‑1 (MQ=255)
                                                                                                                                        ggATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑‑TTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                <  4:751128/122‑1 (MQ=11)
                                                                                                                                        ggATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑‑TTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                >  3:751128/1‑122 (MQ=11)
                                                                                                                                        ggATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑‑TTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                <  2:512048/122‑1 (MQ=11)
                                                                                                                                        ggATGCGGCGTAAACGCCTTATCCGGCCTAC‑‑‑‑GGCTCGG‑‑‑TTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                >  1:512048/1‑122 (MQ=11)
                                                                                                                                              ggcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGc      >  3:399021/1‑149 (MQ=255)
                                                                                                                                               gcgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGcc     >  4:1105088/1‑149 (MQ=255)
                                                                                                                                                cgTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCt    <  4:399021/149‑1 (MQ=255)
                                                                                                                                                  tGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGa  <  4:344327/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGCATTTTGTATGCCCATCCCAGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGTGAGAATTTGTAGGTCTGATAAGACGCGTT‑AGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGGTGGCGCGAGAATTTGTAGGCCTGA  >  NC_000913/374899‑375192

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