Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,879 A→G 9.1% M271V (ATG→GTG) ‡ pheA → bifunctional chorismate mutase/prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,975,8790AG 9.1% 150.3 / 5.0 56M271V (ATG→GTG) ‡pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (3/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCC  >  NC_002947/1975733‑1976024
                                                                                                                                                  |                                                                                                                                                 
tttCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtgat                                                                                                                                                  >  2:248219/1‑148 (MQ=255)
   cGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACTACTCCACGCGCTTCCTgatgat                                                                                                                                               <  7:47048/148‑1 (MQ=255)
    gAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATc                                                                                                                                              <  3:251345/148‑1 (MQ=255)
         aTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc                                                                                                                                        <  3:199750/149‑1 (MQ=255)
          tGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAcc                                                                                                                                       <  2:10945/149‑1 (MQ=255)
          tGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc                                                                                                                                        >  1:47770/1‑148 (MQ=255)
             cGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAgg                                                                                                                                    <  7:288161/149‑1 (MQ=255)
                    gCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATc                                                                                                                                              <  5:183564/132‑1 (MQ=255)
                    gCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATc                                                                                                                                              >  6:183564/1‑132 (MQ=255)
                    gCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAcc                                                                                                                             >  7:100438/1‑149 (MQ=255)
                                    gTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCgacgac                                                                                                             <  1:248219/149‑1 (MQ=255)
                                       tgTGCAACTCGGCGCCTTTCGCGGGCTTTATGGCGGCCACCCTGTCCGGCCTGACCCGTCTGGCCTATAATATCCAAGACCGCCCTCACAACTCCCCGCGCTTCCTGTTTAGCGTTAACCAGCAAGTACCGCCGACCGGCGACGACaa                                                                                                           <  8:100438/147‑1 (MQ=255)
                                         tGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGAc                                                                                                        <  2:47770/149‑1 (MQ=255)
                                           gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt                                                                                                      >  4:164326/1‑149 (MQ=255)
                                           gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGAACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCAGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt                                                                                                      >  4:33845/1‑149 (MQ=255)
                                                 ggcggcGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGatc                                                                                                 <  5:241256/148‑1 (MQ=255)
                                                            ccGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATg                                                                                     <  7:8838/149‑1 (MQ=255)
                                                             cGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt                                                                                      <  6:127976/147‑1 (MQ=255)
                                                             cGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGAt                                                                                      >  5:127976/1‑147 (MQ=255)
                                                               gCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGAAACGATGatc                                                                                  <  5:198496/149‑3 (MQ=255)
                                                                  atatGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGCCGACGACAAGACCTCGATCATCGTTTCGATGAGcaa                                                                                >  5:249743/1‑148 (MQ=255)
                                                                              atcttctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCg                                                                   >  6:257674/8‑149 (MQ=255)
                                                                                cTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGtt                                                                 >  2:93993/1‑149 (MQ=255)
                                                                                   ctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTg                                                                <  5:257674/145‑1 (MQ=255)
                                                                                        ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGTAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGt                                                         <  3:33845/149‑1 (MQ=255)
                                                                                        ccTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGGAACCAGGAAGTACCGCCGACCGGCGCCGACAAGAGCTCGATCATCGTTTGGATGAGCGACAAGCCAGGTGGGTTGCATGAGt                                                         >  6:102416/1‑149 (MQ=255)
                                                                                          tGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGGTTACCGATCATCAGGAAgc                                                                                                                                <  6:247114/76‑2 (MQ=25)
                                                                                          tGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGGTTACCGATCATCAGGAAgc                                                                                                                                >  5:247114/1‑75 (MQ=25)
                                                                                            aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCt                                                     >  5:158174/1‑149 (MQ=255)
                                                                                              ccGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTgg                                                   >  8:293437/1‑149 (MQ=255)
                                                                                                 tCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGc                                                >  6:118150/1‑149 (MQ=255)
                                                                                                      ccGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc                                           >  6:134148/1‑149 (MQ=255)
                                                                                                          gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATc                                       <  6:101531/149‑1 (MQ=255)
                                                                                                          gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCAGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt                                                                  <  5:671/122‑1 (MQ=255)
                                                                                                          gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCAGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt                                                                  >  6:671/1‑122 (MQ=255)
                                                                                                           aagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATTATCGGTAACCAGGAATTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCa                                      <  5:271445/149‑1 (MQ=255)
                                                                                                              aTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa                                   <  5:134148/149‑1 (MQ=255)
                                                                                                               tCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGCCCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGATTTGCTGGTGCCGTTCTATCAGAAc                                  >  2:153933/1‑149 (MQ=255)
                                                                                                                  aaGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCa                                                                         <  1:121273/107‑1 (MQ=255)
                                                                                                                  aaGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCa                                                                         >  2:121273/1‑107 (MQ=255)
                                                                                                                     aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc                                                               >  4:160963/1‑114 (MQ=255)
                                                                                                                     aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc                                                               <  3:160963/114‑1 (MQ=255)
                                                                                                                         cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGCCCTCGATCATCGTTTCGAAGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAcc                        <  5:102416/149‑1 (MQ=255)
                                                                                                                          ccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCt                       <  3:164326/149‑1 (MQ=255)
                                                                                                                                   tCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCAt              >  8:54034/1‑149 (MQ=255)
                                                                                                                                     cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc             <  6:249743/148‑1 (MQ=255)
                                                                                                                                     cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc             <  5:118150/148‑1 (MQ=255)
                                                                                                                                     cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCg            >  7:120611/1‑149 (MQ=255)
                                                                                                                                     cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAACACCTCGATCATCGTTTCGATGAGCAACAAGCCAGCTGCGTTGc                                                               >  4:104361/1‑98 (MQ=255)
                                                                                                                                      cggtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGAGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc             >  1:62039/6‑147 (MQ=255)
                                                                                                                                        gtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc                                           >  7:235462/4‑115 (MQ=38)
                                                                                                                                        gtacTTCCTGGTTACCGATCATCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTc                                           <  8:235462/112‑1 (MQ=38)
                                                                                                                                          gcTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcc       >  7:110341/1‑149 (MQ=255)
                                                                                                                                               cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTg                                                 >  1:57448/1‑102 (MQ=255)
                                                                                                                                               cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTg                                                 <  2:57448/102‑1 (MQ=255)
                                                                                                                                               cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgccc  >  7:87984/1‑149 (MQ=255)
                                                                                                                                               cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca                           >  1:252814/1‑121 (MQ=255)
                                                                                                                                               cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca                           <  2:252814/124‑4 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCC  >  NC_002947/1975733‑1976024

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