Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,028,476	A→G	E170G (GAA→GGA)	hemG →	menaquinone‑dependent protoporphyrinogen IX dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,028,476	0	A	G	100.0%	96.8 / NA	27	E170G (GAA→GGA)	hemG	menaquinone‑dependent protoporphyrinogen IX dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/10);  total (17/10)

GAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTT  >  NZ_CP009273/4028396‑4028558                                                                                 |                                                                                   gAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCa                                                                           <  1:356362/90‑1 (MQ=255)  aaGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCAt                                                                          >  2:315117/1‑90 (MQ=255)                tGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAAcc                                                            >  2:421734/1‑90 (MQ=255)                      ggATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCt                                                      <  2:203426/90‑1 (MQ=255)                         tACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGaa                                                   >  1:149206/1‑90 (MQ=255)                                  agaagTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATa                                          <  1:315117/90‑1 (MQ=255)                                    aagTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATaaa                                        <  2:519888/90‑1 (MQ=255)                                    aagTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATaaa                                        <  2:449016/90‑1 (MQ=255)                                    aagTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATaaa                                        <  2:414951/90‑1 (MQ=255)                                         gTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATa                                   >  1:20664/1‑90 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATa                                   <  1:154513/81‑1 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATa                                   >  2:154513/1‑81 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg                          >  2:55056/1‑90 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg                          >  2:468345/1‑90 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg                          >  2:438838/1‑90 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg                          >  2:146142/1‑90 (MQ=255)                                                  gATTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcg                          >  1:395845/1‑90 (MQ=255)                                                   aTTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATa                                                >  2:448459/1‑67 (MQ=255)                                                   aTTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATa                                                <  1:448459/67‑1 (MQ=255)                                                   aTTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcgc                         >  2:423333/1‑90 (MQ=255)                                                   aTTGGGAGCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATgcgc                         <  2:129680/90‑1 (MQ=255)                                                          gCAGGTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGa                  >  2:435271/1‑90 (MQ=255)                                                             ggTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGaaaa               >  2:162128/1‑90 (MQ=255)                                                             ggTGGCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGaaaa               >  1:301547/1‑90 (MQ=255)                                                                ggCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAat              <  2:422243/88‑1 (MQ=255)                                                                ggCGAATTTCGCCCGAGGAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAat              >  1:422243/1‑88 (MQ=255)                                                                          gCCCGAGGAATCGCCCATTCAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAAttt  >  2:160941/1‑90 (MQ=255)                                                                                 |                                                                                   GAAGATGTCAGGCGGTGAAACGGATACGCGCAAAGAAGTTGTCTATACCGATTGGGAGCAGGTGGCGAATTTCGCCCGAGAAATCGCCCATTTAACCGACAAACCGACGCTGAAATAAGCATAAAGAATAAAAAATGCGCGGTCAGAAAATTATTTTAAATTT  >  NZ_CP009273/4028396‑4028558

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