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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,447,002	0	T	C	100.0%	57.0 / NA	17	T158A (ACC→GCC)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/7);  total (10/7)

GTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATCGCGCAGA  >  NC_000913/2446871‑2447126                                                                                                                                    |                                                                                                                             gTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGc                                                                                                                      <  2:1062913/140‑1 (MQ=255)                    gcacgcaTCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCaaa                                                                                                   >  2:1065763/1‑140 (MQ=255)                                    cTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGatat                                                                                   >  2:1010500/1‑140 (MQ=255)                                          cgTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTg                                                                             <  2:343834/140‑1 (MQ=255)                                          cgTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTg                                                                             <  2:244625/140‑1 (MQ=255)                                                                             ttttGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTc                                          <  1:877685/140‑1 (MQ=255)                                                                                         tGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGCCAGGCGGCCACGGATTTCATGCCCAAATTTCTCGGCGAGATATTTTTTGGAAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGgaaga                              <  2:980671/140‑1 (MQ=255)                                                                                                            ccAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATaa           >  1:269758/1‑140 (MQ=255)                                                                                                            ccAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATaa           >  2:288754/1‑140 (MQ=255)                                                                                                            ccAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATaa           >  1:989082/1‑140 (MQ=255)                                                                                                            ccAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATaa           >  1:834072/1‑140 (MQ=255)                                                                                                            ccAGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCACGGCGGTTTCGCGGGCGGAAGAAACTCCACCGCCGCGACaa           >  1:658718/1‑140 (MQ=255)                                                                                                              aGCGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATc         >  1:344774/1‑140 (MQ=255)                                                                                                                cGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATcgc       >  1:773500/1‑140 (MQ=255)                                                                                                                cGGAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGATTCGCCGGCGGAAAAACAGCCACCGCCTCCATAATccc       >  2:525764/1‑138 (MQ=255)                                                                                                                 ggAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATcgcg      <  2:773500/140‑1 (MQ=255)                                                                                                                 ggAATGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATcgcg      <  2:989082/140‑1 (MQ=255)                                                                                                                     tGTCGCCCATCTGGGCCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGGCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATCGCGCAGa  <  1:1065763/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                             GTCGCCCTCTTTTTTCAGCGCACGCATCAACTCGTCTAACGCGTCGATTTTGTCGGGGTCCGGGCAAAAAAACGGATTTTGCTCGACCTGCGACCAGTCTTTGATATCCAGCGGAATGTCGCCCATCTGGGTCAGGCAGCCACGGATTTCAATACCAAATTTCTCGGCGAGATATTTTTTGGCAATAGCTCCTGCCGCCACGCGCATGGCGGTTTCGCGGGCGGAAGAACGTCCACCGCCGCGATAATCGCGCAGA  >  NC_000913/2446871‑2447126

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