Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,738,236	T→C	G164G (GGT→GGC)	pheA →	fused chorismate mutase/prephenate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,738,236	0	T	C	100.0%	167.8 / NA	47	G164G (GGT→GGC)	pheA	fused chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/30);  total (17/30)

GCGGCGCGCCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGTGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTAC  >  NC_000913/2738099‑2738365                                                                                                                                          |                                                                                                                                  gcggcgCGCCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGCCTATGTCGTCGTACCGTTTGATAATACCAGCTCCGgcgc                                                                                                                                 >  2:330932/1‑140 (MQ=255)  cggcgCGCCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTCTGCCGTCCTACCGATTGAAAAAACCAGCTCCGGCGcc                                                                                                                                >  1:312722/1‑140 (MQ=255)  cggcgCGCCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGcc                                                                                                                                >  2:401377/1‑140 (MQ=255)       cgcCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATaaa                                                                                                                           >  2:414443/1‑140 (MQ=255)       cgcCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATaaa                                                                                                                           <  1:414443/140‑1 (MQ=255)        gcCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAac                                                                                                                          <  2:5204/140‑1 (MQ=255)        gcCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAac                                                                                                                          <  1:401377/140‑1 (MQ=255)        gcCAGTATGCTGCCCATCACTTCGACCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCTATAAac                                                                                                                          <  1:237823/140‑1 (MQ=255)          cAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAacga                                                                                                                        <  2:435984/140‑1 (MQ=255)            gTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAacgacg                                                                                                                      >  2:595648/1‑140 (MQ=255)                gctgcCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTa                                                                                                                  <  1:547400/140‑1 (MQ=255)                    cccGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGAt                                                                                                              >  1:681698/1‑140 (MQ=255)                       gTCACTTTGAGCAATTCATTGATAGTGGCTGTGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATctg                                                                                                           <  2:232750/140‑1 (MQ=255)                           cTTTGGGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATTCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCtgctgc                                                                                                       <  1:27086/140‑1 (MQ=255)                           cTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCtgctgc                                                                                                       <  1:16020/140‑1 (MQ=255)                            tttGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCa                                                                                                      <  2:671425/140‑1 (MQ=255)                            tttGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCa                                                                                                      <  2:52230/140‑1 (MQ=255)                              tGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAAc                                                                                                    >  2:648881/1‑140 (MQ=255)                                aGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGATCCCGCGACATAAACGACGTTAACGATCTGCTGCAACAt                                                                                                  >  2:257993/1‑140 (MQ=255)                                     ttCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAg                                                                                             <  1:595648/140‑1 (MQ=255)                                     ttCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAg                                                                                             <  1:125092/140‑1 (MQ=255)                                     ttCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCACCATAAAGGACGTTGACGATCGGCTGCAACATACCAg                                                                                             >  1:281759/1‑140 (MQ=255)                                      tCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGc                                                                                            >  2:587242/1‑140 (MQ=255)                                               gTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGAtt                                                                                   >  2:513686/1‑140 (MQ=255)                                                 ggCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGAttgt                                                                                 >  1:50446/1‑140 (MQ=255)                                                 ggCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGAttgt                                                                                 >  2:62747/1‑140 (MQ=255)                                                          aaaTTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGAt                                                                        <  1:392650/140‑1 (MQ=255)                                                          aaaTTTGCCGATATTTTTAATCAGGGGGAAACCGGCCAGGCCGACTATGCATTCGTACCGATTGAAAATTCCAGCTCCGGCGCCATAAACGACGTTTACGACCTTGTGCAACATGCCATCTTGTCGATTGTTGGCGAGAt                                                                        <  2:627921/140‑1 (MQ=255)                                                                          ttAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGAc                                                        <  1:587242/140‑1 (MQ=255)                                                                                  gTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGTTTGTTGGGGAGATGACGTTAAATATCGACCAttgttt                                                >  1:258971/1‑140 (MQ=255)                                                                                  gTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCAttgttt                                                >  2:450638/1‑140 (MQ=255)                                                                                      aaaCCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTttgttg                                            <  2:735973/140‑1 (MQ=255)                                                                                      aaaCCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTttgttg                                            <  2:121352/140‑1 (MQ=255)                                                                                           ggCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGtctc                                       <  2:552758/140‑1 (MQ=255)                                                                                            gCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTcc                                      >  1:382581/1‑140 (MQ=255)                                                                                              cAGGCCGACTATGCCGTCGTACGGATTGAGAATACCACCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCgg                                    <  1:257993/140‑1 (MQ=255)                                                                                                   cGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCacta                               <  1:513686/140‑1 (MQ=255)                                                                                                   cGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCacta                               <  1:62747/140‑1 (MQ=255)                                                                                                          gCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCacta                               >  2:657299/1‑133 (MQ=255)                                                                                                          gCCGTCGTACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCacta                               <  1:657299/133‑1 (MQ=255)                                                                                                             gtcgtACCGATTGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATc                     <  2:681698/140‑1 (MQ=255)                                                                                                                gtACCGATTGAAAATGCCAGCTACGGCGCTATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATccac                  <  1:330932/140‑1 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCGTATACGGCGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATa          <  1:329895/140‑1 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATcc                    <  1:89405/130‑1 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATcc                    >  2:89405/1‑130 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATa          <  2:553355/140‑1 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATa          <  2:27842/140‑1 (MQ=255)                                                                                                                        tGAAAATACCAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATa          <  2:444647/140‑1 (MQ=255)                                                                                                                              tACCAGCTCCGGCGCCATAAACGACGTTTACTATCAGCTACAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCt    <  2:206682/140‑1 (MQ=255)                                                                                                                                ccAGCTCCGGCGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGATTGTTGGTCTCCCGCAATACTGATTCATCCGCCATCAATACGGTCTAc  >  2:243966/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  GCGGCGCGCCAGTATGCTGCCCGTCACTTTGAGCAATTCATTGAAAGTGGCTGCGCCAAATTTGCCGATATTTTTAATCAGGTGGAAACCGGCCAGGCCGACTATGCCGTCGTACCGATTGAAAATACCAGCTCCGGTGCCATAAACGACGTTTACGATCTGCTGCAACATACCAGCTTGTCGATTGTTGGCGAGATGACGTTAACTATCGACCATTGTTTGTTGGTCTCCGGCACTACTGATTTATCCACCATCAATACGGTCTAC  >  NC_000913/2738099‑2738365

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