Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,621,307	A→G	E37E (GAA→GAG)	purM →	phosphoribosylaminoimidazole synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,621,307	0	A	G	100.0%	139.6 / NA	39	E37E (GAA→GAG)	purM	phosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/17);  total (22/17)

ACGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGAC  >  NC_000913/2621178‑2621439                                                                                                                                  |                                                                                                                                     aCGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTAGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGc                                                                                                                             <  2:149934/139‑1 (MQ=255)                        ccGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCt                                                                                                      <  2:641888/139‑1 (MQ=255)                           aTAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                        >  2:46242/1‑134 (MQ=255)                           aTAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                        <  1:46242/134‑1 (MQ=255)                               aaCCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGC‑GGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTc                                                                                                              <  2:404145/123‑1 (MQ=255)                               aaCCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGC‑GGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTc                                                                                                              >  1:404145/1‑123 (MQ=255)                                 ccTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGACGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCAtt                                                                                             >  1:54932/1‑139 (MQ=255)                                   tctctTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGc                                                                                           >  1:453798/1‑139 (MQ=255)                                               aaaGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATc                                                                               >  1:737424/1‑139 (MQ=255)                                                    tGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGaa                                                                          >  2:401303/1‑139 (MQ=255)                                                     gCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAAc                                                                         >  1:73641/1‑139 (MQ=255)                                                                     aCGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTgg                                                         >  1:442119/1‑139 (MQ=255)                                                                      cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCg                                                                                                             <  2:325886/86‑1 (MQ=255)                                                                      cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCg                                                                                                             >  1:325886/1‑86 (MQ=255)                                                                      cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGc                                                        >  2:388930/1‑139 (MQ=255)                                                                      cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGc                                                        >  1:291707/1‑139 (MQ=255)                                                                       gcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCa                                                       <  2:497230/139‑1 (MQ=255)                                                                         gggTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACt                                                     <  2:54932/139‑1 (MQ=255)                                                                          ggTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTg                                                    >  2:621451/1‑139 (MQ=255)                                                                           gTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGa                                                   <  2:73641/139‑1 (MQ=255)                                                                              aTGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACgg                                                >  2:687031/1‑139 (MQ=255)                                                                                     ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGGTGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGc                                               <  1:233758/133‑1 (MQ=255)                                                                                     ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGc                                               >  2:233758/1‑133 (MQ=255)                                                                                      gTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTa                                        >  2:233640/1‑139 (MQ=255)                                                                                      gTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTa                                        <  2:310063/139‑1 (MQ=255)                                                                                      gTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTa                                        >  2:638331/1‑139 (MQ=255)                                                                                         ggAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGGTGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCa                                     >  2:238211/1‑139 (MQ=255)                                                                                           aagaaTCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAg                                   <  2:303515/139‑1 (MQ=255)                                                                                           aagaaTCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAg                                   >  1:593553/1‑139 (MQ=255)                                                                                           aagaaTCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAg                                   <  1:388930/139‑1 (MQ=255)                                                                                           aagaaTCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAg                                   >  1:303515/1‑139 (MQ=255)                                                                                                    aGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTg                          >  1:134823/1‑139 (MQ=255)                                                                                                         tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAAt                     <  2:442119/139‑1 (MQ=255)                                                                                                                gaaAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTa              >  2:10230/1‑139 (MQ=255)                                                                                                                 aaaaCGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTaa             <  1:401303/139‑1 (MQ=255)                                                                                                                          cgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      <  2:162726/121‑1 (MQ=255)                                                                                                                          cgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCaa                      >  1:162726/1‑121 (MQ=255)                                                                                                                           gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGAGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGa   <  2:379429/139‑1 (MQ=255)                                                                                                                            tcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGac  <  1:638331/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                     ACGCGTGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGAC  >  NC_000913/2621178‑2621439

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