Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,835	G→A	G213G (GGC→GGT)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,835	0	G	A	96.2%	250.1 / ‑2.3	79	G213G (GGC→GGT)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (36/40);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTC  >  NC_000913/2446740‑2446930                                                                                                |                                                                                                aTTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGCCAAAGACCGGCTCGCCAAGTCCGGCAGGAAcaccac                                                                                             >  1:1266005/1‑100 (MQ=255)   ttCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTg                                                                                           <  2:1830896/100‑1 (MQ=255)    tCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTgg                                                                                          <  2:1029988/100‑1 (MQ=255)      cACGCCTTTCACCGCGTTGATGCTCATTAGCGCATGGGCGATGTCGGCATCCAGGCGGTCAAAGACAGGCTCGCCCAGTCCGGCGGGGAAGGCCCTGGca                                                                                        >  1:651536/1‑100 (MQ=255)      cACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGca                                                                                        >  1:2552732/1‑100 (MQ=255)        cGCCTTTCACCGCGTTGATGCTCATCAGCGCATGAGCGATGTCAGAAACCAAGCGGTCAAAGACCGACTCGCCAAGTCCGGCCGTACCACCACTGGcaac                                                                                      <  2:245703/100‑1 (MQ=255)         gCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGcaaca                                                                                     <  1:1199489/100‑1 (MQ=255)           cTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCaacaac                                                                                   <  1:596657/100‑1 (MQ=255)               cACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACggag                                                                               >  2:59134/1‑98 (MQ=255)               cACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTg                                                                               <  1:125444/100‑1 (MQ=255)               cACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTg                                                                               >  2:1480920/1‑100 (MQ=255)               cACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCCGCAGGAACACCACTGGCAACAACGGTg                                                                               <  2:1546521/100‑1 (MQ=255)               cACCGCGTTGATGCTCATCAGCGCAGGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTg                                                                               <  2:1933916/100‑1 (MQ=255)               cACCGCGTAGATGCTCCTCAGCGCATGGGCGCTGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTg                                                                               <  1:1446308/100‑1 (MQ=255)                 ccGCGTTGATTCTCATCAGCGCATGGGCGATGTCAGCAGCCAGGCGGTCAAGGACCGGCTCGCCAGGTCCGGCAGGAACACCACTGGGAACAACtgtggc                                                                             >  2:1423394/1‑98 (MQ=255)                 ccGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGAc                                                                             >  2:1208950/1‑100 (MQ=255)                   gcgTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACtt                                                                           >  1:856367/1‑100 (MQ=255)                         aTGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAgcgc                                                                     >  1:1798718/1‑100 (MQ=255)                          tGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGcc                                                                    <  1:1799547/100‑1 (MQ=255)                          tGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGcc                                                                    <  2:1592589/100‑1 (MQ=255)                           gCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCg                                                                   <  1:1232731/100‑1 (MQ=255)                            cTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGa                                                                  >  1:1319600/1‑100 (MQ=255)                            cTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGa                                                                  <  2:1466050/100‑1 (MQ=255)                             tcatcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGAt                                                                 <  2:308085/100‑1 (MQ=255)                             tcatcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGAt                                                                 <  1:2021990/100‑1 (MQ=255)                              catcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATg                                                                >  1:1168258/1‑100 (MQ=255)                              catcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATg                                                                >  2:33836/1‑100 (MQ=255)                              catcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATg                                                                >  2:34472/1‑100 (MQ=255)                               atcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATgg                                                               >  2:721378/1‑100 (MQ=255)                               atcaGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATgg                                                               >  1:2628998/1‑100 (MQ=255)                                 caGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAg                                                             <  2:88414/100‑1 (MQ=255)                                   gcgcATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTc                                                           >  1:1698606/1‑100 (MQ=255)                                     gcATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGc                                                         <  1:2016980/100‑1 (MQ=255)                                      cATGGGCGGTGTCAGCATCCAGGCGGTCAAAGAACGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGcc                                                        <  1:1654740/100‑1 (MQ=255)                                         gggCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCctc                                                     >  1:1826858/1‑100 (MQ=255)                                         gggCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCctc                                                     >  1:2152720/1‑100 (MQ=255)                                          ggCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCctct                                                    >  1:1123228/1‑100 (MQ=255)                                             gATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCtttt                                                 <  1:922316/100‑1 (MQ=255)                                              aTGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCttttt                                                <  1:630417/100‑1 (MQ=255)                                              aTGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCttttt                                                <  2:1927706/100‑1 (MQ=255)                                               tGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCtttttt                                               >  1:1543164/1‑100 (MQ=255)                                                gTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTc                                              <  1:2397252/100‑1 (MQ=255)                                                gTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTc                                              <  1:2396522/100‑1 (MQ=255)                                                 tCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCa                                             >  2:2094692/1‑100 (MQ=255)                                                      aTCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGcgca                                        >  2:721902/1‑100 (MQ=255)                                                      aTCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGcgca                                        <  2:1854843/100‑1 (MQ=255)                                                      aTCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGcgca                                        <  1:125420/100‑1 (MQ=255)                                                       tCCAGGCGCTCAAAGACCTTCTTTTTTAGTTGGGCAGGAGCACCACGTGCAATAACTGTGACGGGAGCTCCGATGGAGTCGCCCTCTTGTTTCAGcgcac                                       <  2:2589847/100‑1 (MQ=255)                                                         cAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCGCTCTTTTTTCAGcgcacgc                                     >  1:254683/1‑100 (MQ=255)                                                            gCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATc                                  <  1:973435/100‑1 (MQ=255)                                                               gTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAAc                               >  1:2169321/1‑100 (MQ=255)                                                                tCAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGCCGCCCTCTTTTTTCAGCGCACGCATCAACt                              >  2:850886/1‑100 (MQ=255)                                                                 cAAAGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTc                             >  1:1635429/1‑100 (MQ=255)                                                                  aaaGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCg                            >  1:1247563/1‑100 (MQ=255)                                                                   aaGACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGt                           >  1:1520888/1‑100 (MQ=255)                                                                     gACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCTTCTTTTTTCAGCGCAGGCATCAACTCGTCt                         >  2:495763/1‑100 (MQ=255)                                                                     gACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCt                         <  1:2418068/100‑1 (MQ=255)                                                                     gACCGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCt                         >  2:2472779/1‑100 (MQ=255)                                                                      aCCGGCTCGCCAAGGCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTGTTCAGCGCACGCATCAACTCGTCTa                        <  2:925999/100‑1 (MQ=255)                                                                       ccGGCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTaa                       >  1:724327/1‑100 (MQ=255)                                                                         ggCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAAcg                     >  1:1403136/1‑100 (MQ=255)                                                                         ggCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAAcg                     >  1:75704/1‑100 (MQ=255)                                                                         ggCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAAcg                     <  1:7236/100‑1 (MQ=255)                                                                         ggCTCGCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAAcg                     <  2:2540550/100‑1 (MQ=255)                                                                              gCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCg                >  2:2127616/1‑100 (MQ=255)                                                                              gCCAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCg                <  1:640161/100‑1 (MQ=255)                                                                                cAAGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAt              <  2:932419/100‑1 (MQ=255)                                                                                 aaGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAtt             <  2:1817390/100‑1 (MQ=255)                                                                                 aaGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAtn             <  1:2067303/100‑2 (MQ=255)                                                                                  aGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAttt            >  1:1783272/1‑100 (MQ=255)                                                                                  aGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAttt            <  1:34372/100‑1 (MQ=255)                                                                                  aGTCCGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGAttt            <  1:35020/100‑1 (MQ=255)                                                                                    tCCGGCAGGAACACCGCTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTg          <  2:1907272/100‑1 (MQ=255)                                                                                     ccGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGt         >  1:623965/1‑100 (MQ=255)                                                                                     ccGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGGCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGACGATTTTGt         >  1:612761/1‑100 (MQ=255)                                                                                      cGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTTTTTTTTTTCATCAACTCGTCTAACGCGTCGATTTTGTc        >  2:743957/1‑100 (MQ=255)                                                                                      cGGCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTc        <  2:2300638/100‑1 (MQ=255)                                                                                        gCAGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCgg      <  1:1622589/100‑1 (MQ=255)                                                                                          aGGAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCgggg    >  2:1510721/1‑100 (MQ=255)                                                                                            gAACACCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTc  <  1:1025850/100‑1 (MQ=255)                                                                                                |                                                                                                ATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACGCCACTGGCAACAACGGTGACTTTAGCGCCGATGGAGTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTC  >  NC_000913/2446740‑2446930

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