Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	786,717	C→T	intergenic (+32/+126)	aroG → / ← gpmA	3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive/2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,717	0	C	T	100.0%	124.9 / NA	39	intergenic (+32/+126)	aroG/gpmA	3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive/2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/21);  total (18/21)

ACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG  >  NC_000913/786591‑786850                                                                                                                               |                                                                                                                                      aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAAt                                                                                                                          <  1:632171/140‑1 (MQ=255) aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAAt                                                                                                                          <  1:42528/140‑1 (MQ=255)      aaGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATcacca                                                                                                                     <  2:566821/140‑1 (MQ=255)              cACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATcacca                                                                                                                     <  2:180914/132‑1 (MQ=255)              cACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATcacca                                                                                                                     >  1:180914/1‑132 (MQ=255)              cACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGGGGTGTATCCGATGAATCCCCACACGCCTg                                                                                                             >  1:561469/1‑140 (MQ=255)               aCCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGa                                                                                                            >  1:56460/1‑140 (MQ=255)                   aTGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTGATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAg                                                                                                        >  2:496787/1‑140 (MQ=255)                   aTGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAg                                                                                                        <  2:809911/140‑1 (MQ=255)                   aTGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAg                                                                                                        >  2:797154/1‑140 (MQ=255)                   aTGCCTGCATCGGCTGGGAAGACACCGATGCTCTGGTACTTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAg                                                                                                        <  2:274748/140‑1 (MQ=255)                             cGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCg                                                                                              <  2:162519/140‑1 (MQ=255)                                          aCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCa                                                                                  >  2:173127/1‑139 (MQ=255)                                          aCCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCa                                                                                  <  1:173127/139‑1 (MQ=255)                                             gATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAtgt                                                                              >  2:651989/1‑140 (MQ=255)                                              aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAtgtg                                                                             >  2:644991/1‑140 (MQ=255)                                              aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAtgtg                                                                             <  2:796404/140‑1 (MQ=255)                                              aTGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAtgtg                                                                             <  1:162898/140‑1 (MQ=255)                                              aTGCTCTGGTACGTCAACTGGCTAATGCAGTAAGATCGCCTCGCGGGTTAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAAtgtg                                                                             <  2:404273/140‑1 (MQ=255)                                               tGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGc                                                                            >  1:198634/1‑140 (MQ=255)                                                  tctGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTcc                                                                         >  1:462491/1‑140 (MQ=255)                                                  tctGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTcc                                                                         >  1:184968/1‑140 (MQ=255)                                                        aCGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGtt                                                                   <  2:198634/140‑1 (MQ=255)                                                              aCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGcaac                                                             >  1:123689/1‑140 (MQ=255)                                                              aCTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGcaac                                                             >  1:160419/1‑140 (MQ=255)                                                                 ggCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACaaa                                                          <  1:797154/140‑1 (MQ=255)                                                                  gCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACaaaa                                                         >  1:640651/1‑140 (MQ=255)                                                                         cAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGa                                                  <  2:68817/140‑1 (MQ=255)                                                                            tAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTc                                               >  2:504307/1‑140 (MQ=255)                                                                            tAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTc                                               >  2:105802/1‑140 (MQ=255)                                                                              aaaGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAc                                             <  2:56460/140‑1 (MQ=255)                                                                                    cgTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAg                                       >  1:629460/1‑140 (MQ=255)                                                                                    cgTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAg                                       >  1:197279/1‑140 (MQ=255)                                                                                      tCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTc                                     >  1:465040/1‑140 (MQ=255)                                                                                            gTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCttt                               <  1:651989/140‑1 (MQ=255)                                                                                             tAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTtc                              <  2:197279/140‑1 (MQ=255)                                                                                                     aactGTCGGATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTa                      <  2:490283/137‑1 (MQ=255)                                                                                                            ggATGCGCCGTCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCGATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAAt               <  2:561469/140‑1 (MQ=255)                                                                                                                      tCAGAGTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACt     <  1:685118/140‑1 (MQ=255)                                                                                                                         gagTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCTGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCg  <  1:504307/140‑1 (MQ=255)                                                                                                                         gagTGGTGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCg  <  1:105802/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                      ACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCG  >  NC_000913/786591‑786850

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