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%	152.6 / NA	44	E37E (GAA→GAG)	purM	phosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/20);  total (24/20)

TGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGACAC  >  NC_000913/2621183‑2621441                                                                                                                             |                                                                                                                                       tGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCt                                                                                                                          >  2:313863/1‑139 (MQ=255)                  aCCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcgc                                                                                                         >  2:852557/1‑139 (MQ=255)                            ccTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCAtt                                                                                               <  2:179965/139‑1 (MQ=255)                                 ctTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTg                                                                                                            >  2:701791/1‑121 (MQ=255)                                 ctTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTg                                                                                                            <  1:701791/121‑1 (MQ=255)                                       tACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                          <  1:1000807/117‑1 (MQ=255)                                       tACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                          >  2:1000807/1‑117 (MQ=255)                                            aGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGt                                                                               >  2:831215/1‑139 (MQ=255)                                              aTGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGCAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGa                                                                             >  1:1127072/1‑139 (MQ=255)                                                gCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAAc                                                                           >  1:604226/1‑139 (MQ=255)                                                            atTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGttt                                                               >  2:585336/1‑139 (MQ=255)                                                               gACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                             <  2:423369/90‑1 (MQ=255)                                                               gACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                             >  1:423369/1‑90 (MQ=255)                                                               gACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTg                                                            >  2:31731/1‑139 (MQ=255)                                                                 cgcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGc                                                          >  2:561349/1‑139 (MQ=255)                                                                  gcgGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCa                                                         <  2:313771/139‑1 (MQ=255)                                                                     ggTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTg                                                      >  2:127415/1‑139 (MQ=255)                                                                      gTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGa                                                     <  2:604226/139‑1 (MQ=255)                                                                         aTGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACgg                                                  <  1:585336/139‑1 (MQ=255)                                                                          tGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGc                                                 >  2:1217397/1‑139 (MQ=255)                                                                          tGCTCTGGTTGGAAGAATC‑AAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTg                                                                   <  1:716956/120‑1 (MQ=255)                                                                          tGCTCTGGTTGGAAGAATC‑AAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTg                                                                   >  2:716956/1‑120 (MQ=255)                                                                                ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                             >  1:1094353/1‑73 (MQ=255)                                                                                ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGt                                                                                                             <  2:1094353/73‑1 (MQ=255)                                                                                ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGt                                           >  2:900774/1‑139 (MQ=255)                                                                                ggTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGt                                           >  2:820117/1‑139 (MQ=255)                                                                                 gTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTa                                          >  2:153774/1‑139 (MQ=255)                                                                                    ggAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCaaaaa                                                                                     >  1:380642/1‑93 (MQ=255)                                                                                    ggAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCaaaaa                                                                                     <  2:380642/93‑1 (MQ=255)                                                                                                ggCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTg                                                      <  2:285622/112‑1 (MQ=255)                                                                                                ggCGTAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTg                                                      >  1:285622/1‑112 (MQ=255)                                                                                                    tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                          >  1:955610/1‑56 (MQ=255)                                                                                                    tAGTGAAGAAAACGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTgcg                                                                                                          <  2:955610/56‑1 (MQ=255)                                                                                                              aaCGCGTCGTCCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTaaaa             <  1:852557/139‑1 (MQ=255)                                                                                                                   gtcgtcCGTAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCg                                                                                >  2:790535/1‑67 (MQ=255)                                                                                                                   gtcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCg                                                                                <  1:790535/67‑1 (MQ=255)                                                                                                                   gtcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCa        >  1:477841/1‑139 (MQ=255)                                                                                                                    tcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCAc       <  1:31731/139‑1 (MQ=255)                                                                                                                    tcgtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCAc       <  1:1217397/139‑1 (MQ=255)                                                                                                                      gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGc                                                                                             >  1:168679/1‑51 (MQ=255)                                                                                                                      gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGc                                                                                             <  2:168679/51‑1 (MQ=255)                                                                                                                      gtcCGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGa     <  1:313863/139‑1 (MQ=255)                                                                                                                         cGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGacac  <  1:900774/139‑1 (MQ=255)                                                                                                                         cGGAGGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGacac  <  1:153774/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                       TGGGGACCCAAGCAGTGACCGATAAAACCTCTCTTAGCTACAAAGATGCCGGTGTTGATATTGACGCGGGTAATGCTCTGGTTGGAAGAATCAAAGGCGTAGTGAAGAAAACGCGTCGTCCGGAAGTGATGGGCGGTCTGGGCGGCTTCGGTGCGCTGTGTGCATTGCCGCAAAAATATCGTGAACCCGTGCTGGTTTCTGGCACTGACGGCGTAGGTACCAAGCTGCGTCTGGCAATGGACTTAAAACGTCACGACAC  >  NC_000913/2621183‑2621441

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