Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,664,529	G→A	I231I (ATC→ATT)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,529	0	G	A	100.0%	335.0 / NA	97	I231I (ATC→ATT)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/55);  total (42/55)

GTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGC  >  minE/1664478‑1664591                                                    |                                                               gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCt                                              >  1:234302/1‑70 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:551464/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:1726653/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:3191312/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:1604389/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:473962/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:1270921/1‑71 (MQ=255) gTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTa                                             >  1:2833107/1‑71 (MQ=255)   gcagcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGCGGCCGAGAATGACTTTATCGTCGCCTAcc                                           >  1:2216304/1‑71 (MQ=255)   gcagcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTAcc                                           >  1:598254/1‑71 (MQ=255)   gcagcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTAcc                                           >  1:1752323/1‑71 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCGCACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2731611/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:616182/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2226414/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2216824/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2181901/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2093584/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2090261/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:3003445/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:3054639/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2723127/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2699772/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2719612/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1132874/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:475143/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1263913/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:12590/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1374327/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1258840/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:382404/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:930142/71‑1 (MQ=255)     agcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:36379/71‑1 (MQ=255)      gcagcaTTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1324176/70‑1 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCt                                              >  1:2909994/1‑54 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTcg                                     >  1:1804118/1‑63 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:1194484/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:123675/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:2479864/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:2437575/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:2217076/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:3127281/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:2958754/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:1452737/1‑70 (MQ=255)                 gTTACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGc                              >  1:1937333/1‑70 (MQ=255)                  gtACGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:2077335/57‑1 (MQ=255)                     cagaggaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2834220/69‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:404873/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2912370/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:403759/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2838736/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:3047550/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2769625/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:549927/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:3246394/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:569833/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2668763/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:934440/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1385347/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1527361/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:165369/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1747375/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1747588/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:968637/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1171988/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:216247/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2179921/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:715865/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:234940/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2395176/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:706908/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2484642/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:630665/71‑1 (MQ=255)                     cGGAGGAATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:622417/71‑1 (MQ=255)                      ggaggaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:1336708/70‑1 (MQ=255)                          gaATCCACACCACCAGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCt                         <  1:2838102/66‑1 (MQ=255)                                      ccaGAGAGGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:277079/38‑1 (MQ=255)                                         gagagGCCGAGAATGACTTTATCGTCGCCTACCTg                                         <  1:1484798/35‑1 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:971961/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1023154/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:926220/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:898662/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:806381/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:750877/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:3073969/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2869212/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2626722/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2525268/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2475953/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2209898/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:2012857/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1984842/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1806865/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1794949/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1485080/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1335298/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1268129/1‑71 (MQ=255)                                            agGCCGAGAATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGc  >  1:1033816/1‑71 (MQ=255)                                                    |                                                               GTGCAGCAGCATTGCGGTTACGGAGGAATCCACACCACCAGAGAGGCCGAGGATGACTTTATCGTCGCCTACCTGCTCGCGGATGCGAGCTACAGCATCGTCGATAATTTTCGC  >  minE/1664478‑1664591

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑