Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,579	T→C	T299A (ACC→GCC)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,579	0	T	C	98.2%	197.2 / NA	57	T299A (ACC→GCC)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/25);  minor base A (0/1);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCA  >  NC_000913/2446447‑2446713                                                                                                                                     |                                                                                                                                       gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTg                                                                                                                                 >  1:342772/1‑140 (MQ=255)          aaaCGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGt                                                                                                                        >  1:637238/1‑140 (MQ=255)               aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                   >  1:1007645/1‑140 (MQ=255)               aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                   >  1:280546/1‑140 (MQ=255)                 cGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgcg                                                                                                                 >  2:63264/1‑140 (MQ=255)                 cGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgcg                                                                                                                 >  1:246161/1‑140 (MQ=255)                  gCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgcgc                                                                                                                <  2:246161/140‑1 (MQ=255)                     aGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCat                                                                                                             <  1:1011119/140‑1 (MQ=255)                        aTCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATg                                                                                                          >  1:572715/1‑140 (MQ=255)                          cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATggg                                                                                                        >  2:1515/1‑140 (MQ=255)                            cTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCa                                                                                                      <  1:1079091/140‑1 (MQ=255)                            cTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCa                                                                                                      >  1:397768/1‑140 (MQ=255)                                      cGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgct                                                                                            <  1:1515/140‑1 (MQ=255)                                      cGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGAAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgct                                                                                            <  2:354974/140‑1 (MQ=255)                                         cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                         <  1:267765/140‑1 (MQ=255)                                         cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                         >  1:441864/1‑140 (MQ=255)                                         cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                         <  2:572715/140‑1 (MQ=255)                                          aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGccc                                                                                        <  2:1007645/140‑1 (MQ=255)                                          aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGccc                                                                                        <  1:927652/140‑1 (MQ=255)                                              cgcgGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgctg                                                                                    >  1:200522/1‑140 (MQ=255)                                                 ggATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGctgctg                                                                                 >  1:175608/1‑140 (MQ=255)                                                   aTCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGAt                                                                               >  1:437446/1‑140 (MQ=255)                                                             aGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCgaga                                                                     >  1:47599/1‑140 (MQ=255)                                                               ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                   >  2:1113655/1‑140 (MQ=255)                                                               ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                   >  1:1005290/1‑140 (MQ=255)                                                                gATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATg                                                                  >  2:1109038/1‑140 (MQ=255)                                                                 aTCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgc                                                                 >  1:1110057/1‑140 (MQ=255)                                                                  tCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgcc                                                                <  1:696305/140‑1 (MQ=255)                                                                  tCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgcc                                                                <  1:1109038/140‑1 (MQ=255)                                                                         ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCa                                                         >  2:480066/1‑140 (MQ=255)                                                                          gCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCATTGATTTGCTGCCCGCTGCCGATACCGCCGAGAATGCCGCCCGCAt                                                        <  1:838692/140‑1 (MQ=255)                                                                           ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                       >  2:552009/1‑140 (MQ=255)                                                                           ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                       >  2:171782/1‑140 (MQ=255)                                                                            cTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATgg                                                      <  1:1113655/140‑1 (MQ=255)                                                                                ggTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgcccg                                                           >  2:828460/1‑131 (MQ=255)                                                                                ggTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgcccg                                                           <  1:828460/131‑1 (MQ=255)                                                                                  tGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGctc                                                <  1:225720/140‑1 (MQ=255)                                                                                   gATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGctct                                               <  2:441864/140‑1 (MQ=255)                                                                                     tcatcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTgg                                             >  1:443314/1‑140 (MQ=255)                                                                                      catcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGa                                            <  1:171782/140‑1 (MQ=255)                                                                                       atcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaa                                           <  2:280546/140‑1 (MQ=255)                                                                                        tctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaaa                                          <  2:175608/140‑1 (MQ=255)                                                                                         ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                         <  2:47599/140‑1 (MQ=255)                                                                                         ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                         <  2:200522/140‑1 (MQ=255)                                                                                           cAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                       >  2:727057/1‑140 (MQ=255)                                                                                           cAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                       >  1:754654/1‑140 (MQ=255)                                                                                           cAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                       >  1:786034/1‑140 (MQ=255)                                                                                              cttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCt                                    <  1:63264/140‑1 (MQ=255)                                                                                                    cGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTggcg                              >  2:929115/1‑138 (MQ=255)                                                                                                    cGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTg                              >  1:276610/1‑140 (MQ=255)                                                                                                            cGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATc                      <  1:552009/140‑1 (MQ=255)                                                                                                                            ccGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcg      <  2:397768/140‑1 (MQ=255)                                                                                                                              ggCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgcg    <  2:1005290/140‑1 (MQ=255)                                                                                                                               gCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGTCgcgc   >  2:899451/1‑140 (MQ=255)                                                                                                                               gCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCgcgc   >  1:526188/1‑140 (MQ=255)                                                                                                                               gCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCgcgc   <  2:637238/140‑1 (MQ=255)                                                                                                                                cACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCa  >  1:542155/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                       GATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCA  >  NC_000913/2446447‑2446713

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