Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,886 G→A 9.4% G273D (GGT→GAT)  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,8860GA 9.4% 146.7 / 5.0 54G273D (GGT→GAT) pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (3/2);  total (28/26)
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.84e-01

ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAG  >  NC_002947/1975742‑1976033
                                                                                                                                                |                                                                                                                                                   
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)
           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                                                                                                                                         >  5:247114/1‑75 (MQ=25)
                                                                                 tGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGGTTACCGATCATCAGGAAgc                                                                                                                                         <  6:247114/76‑2 (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                                                                           >  6:671/1‑122 (MQ=255)
                                                                                                 gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCAGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGt                                                                           <  5:671/122‑1 (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                                                                        <  3:160963/114‑1 (MQ=255)
                                                                                                            aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc                                                                        >  4:160963/1‑114 (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)
                                                                                                                                             tCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgcg     >  4:249930/1‑148 (MQ=255)
                                                                                                                                               ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg  >  8:132460/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAG  >  NC_002947/1975742‑1976033

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: ‑