Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,009,285	T→C	E82G (GAG→GGG)	ysgA ←	dienelactone hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,009,285	0	T	C	100.0%	103.9 / NA	29	E82G (GAG→GGG)	ysgA	dienelactone hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/18);  total (11/18)

AGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGG  >  NZ_CP009273/4009168‑4009403                                                                                                                      |                                                                                                                       aGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACaaa                                                                                                   <  1:342674/139‑1 (MQ=255)         cGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTcgc                                                                                             <  2:230401/137‑1 (MQ=255)           aGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGat                                                                                         <  1:181025/139‑1 (MQ=255)                ggCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGtt                                                                                    >  2:378085/1‑139 (MQ=255)                 gCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGtt                                                                                    >  1:292309/1‑138 (MQ=255)                       ttGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCAc                                                                             >  1:329783/1‑139 (MQ=255)                            aCCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCa                                                                         >  2:464917/1‑138 (MQ=255)                              cAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCaaaa                                                                      <  2:292309/139‑1 (MQ=255)                                 aCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAAtt                                                                   >  1:477392/1‑139 (MQ=255)                                   cGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTc                                                                 <  2:329783/139‑1 (MQ=255)                                   cGCCAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTc                                                                 >  1:299367/1‑139 (MQ=255)                                     cTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCt                                                               >  2:302673/1‑139 (MQ=255)                                        aGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCa                                                            <  2:477392/139‑1 (MQ=255)                                             cGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTAc                                                        >  1:473646/1‑138 (MQ=255)                                               tGGGGTTATCGGCAAAATAATCCGGATCGCCTTCGCGGAAGTAAAGTTAAGGGGCGATAGCAGGAT‑CCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATg                                                    <  1:557239/139‑1 (MQ=255)                                               tGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAAt                                                     >  2:112416/1‑139 (MQ=255)                                                 gggATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATg                                                    <  2:108557/138‑1 (MQ=255)                                                     tatCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCAc                                               <  2:608118/139‑1 (MQ=255)                                                          gCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCa                                          <  1:464917/139‑1 (MQ=255)                                                           cAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAg                                         <  2:433737/139‑1 (MQ=255)                                                           cAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAg                                         <  2:51241/139‑1 (MQ=255)                                                           cAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAg                                         <  2:534026/139‑1 (MQ=255)                                                                tCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGc                                    <  1:491088/139‑1 (MQ=255)                                                                                 gcgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCa                                                                                                                     <  2:359755/41‑1 (MQ=255)                                                                                 gcgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCa                                                                                                                     >  1:359755/1‑41 (MQ=255)                                                                                 gcgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGt                   <  1:302673/139‑1 (MQ=255)                                                                                                  gTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAgg  <  2:389233/139‑1 (MQ=255)                                                                                                               aTACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATg                                                    <  2:188363/76‑1 (MQ=255)                                                                                                               aTACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATg                                                    >  1:188363/1‑76 (MQ=255)                                                                                                                      |                                                                                                                       AGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGG  >  NZ_CP009273/4009168‑4009403

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