Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,890 C→T 9.4% N274N (AAC→AAT) ‡ 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,8900CT 9.4% 146.6 / 4.8 53N274N (AAC→AAT) ‡pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base T (3/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGC  >  NC_002947/1975742‑1976037
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                             <  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                                                                               >  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                                                                            >  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                          <  5:118150/148‑1 (MQ=255)
                                                                                                                            cACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATc                          <  6:249743/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                                        <  2:252814/124‑4 (MQ=255)
                                                                                                                                      cTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGATGCCGTTCTGATAGAACGGCacca                                        >  1:252814/1‑121 (MQ=255)
                                                                                                                                             tCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTcgcg         >  4:249930/1‑148 (MQ=255)
                                                                                                                                               ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg      >  8:132460/1‑149 (MQ=255)
                                                                                                                                                   aCCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGTGCAGCgg   <  1:93993/148‑1 (MQ=255)
                                                                                                                                                    ccAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc  >  4:126212/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGC  >  NC_002947/1975742‑1976037

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: