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	96.5%	95.3 / ‑6.1	29	T299A (ACC→GCC)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/8);  minor base A (1/0);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGC  >  NC_000913/2446452‑2446712                                                                                                                                |                                                                                                                                      aTTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGt                                                                                                                           >  2:343244/1‑140 (MQ=255)         gATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCa                                                                                                                   >  1:182524/1‑140 (MQ=255)           tCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgc                                                                                                                 >  1:720865/1‑140 (MQ=255)           tCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGAAATGCTGGAGGTCGGTTTCAgc                                                                                                                 >  1:514637/1‑140 (MQ=255)            cGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAgcg                                                                                                                <  2:636056/140‑1 (MQ=255)                         tCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAAt                                                                                                   >  1:207192/1‑140 (MQ=255)                                     aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGccc                                                                                       <  1:172240/140‑1 (MQ=255)                                     aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGccc                                                                                       <  1:76520/140‑1 (MQ=255)                                     aCTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGccc                                                                                       <  2:576808/140‑1 (MQ=255)                                         cgcgGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCAC‑GCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgct                                                                                           <  1:420826/131‑1 (MQ=255)                                         cgcgGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCAC‑GCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgct                                                                                           >  2:420826/1‑131 (MQ=255)                                                 ccGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAcc                                                                           >  2:647842/1‑140 (MQ=255)                                                   gACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccgc                                                                         <  2:720865/140‑1 (MQ=255)                                                                  aCGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgcccg                                                          >  2:30804/1‑140 (MQ=255)                                                                    ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCTGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCa                                                        >  1:519734/1‑140 (MQ=255)                                                                    ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCGGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCa                                                        <  2:519734/140‑1 (MQ=255)                                                                    ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCa                                                        >  2:8996/1‑140 (MQ=255)                                                                      ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                      >  1:186330/1‑140 (MQ=255)                                                                               atcatcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTg                                             >  1:393083/1‑140 (MQ=255)                                                                                  atcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaa                                          <  2:207192/140‑1 (MQ=255)                                                                                    ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                        >  1:505302/1‑140 (MQ=255)                                                                                      cAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                      >  2:264312/1‑140 (MQ=255)                                                                                            cttcGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTg                                >  2:151940/1‑140 (MQ=255)                                                                                                              tGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAttt                         <  1:709478/129‑1 (MQ=255)                                                                                                              tGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAttt                         >  2:709478/1‑129 (MQ=255)                                                                                                                   cGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCt         >  2:44841/1‑140 (MQ=255)                                                                                                                   cGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCt         >  2:711975/1‑140 (MQ=255)                                                                                                                   cGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCt         >  1:171954/1‑140 (MQ=255)                                                                                                                          gCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCgcgc  >  2:147304/1‑140 (MQ=255)                                                                                                                          gCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCgcgc  >  1:47602/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                      ATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGC  >  NC_000913/2446452‑2446712

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