Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,260,284	A→G	D372G (GAC→GGC)	hemA →	glutamyl‑tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,260,284	0	A	G	100.0%	56.3 / NA	18	D372G (GAC→GGC)	hemA	glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/12);  total (6/12)

AAACCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGC  >  NZ_CP009273/1260152‑1260418                                                                                                                                     |                                                                                                                                       aaaCCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCa                                                                                                                                  <  1:7474/139‑1 (MQ=255)   aCCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTAGCGCAGCCAGGCAGGGCAAGTTGGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAg                                                                                                                                <  1:267640/139‑1 (MQ=255)                gggCGTGGATGCGAGCAAAAAGCGCCAGCGAAACCTTGCGCGATTATCGCAGCCAGGCAAGGAAAGTTCGCATGGCGTTAACCGCCACAGCGTTAGCGGCCCTTGAGCAGGGGGGGGGCGCGCAAGCCATTATGCAgg                                                                                                                    <  2:125520/137‑1 (MQ=255)                                        ccAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCg                                                                                           <  2:194887/139‑1 (MQ=255)                                                              tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAATCTGACTAACCGCTTGATCCATGCGCCAACGaaa                                                                     >  2:419062/1‑139 (MQ=255)                                                              tCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGaa                                                                      <  2:358154/138‑1 (MQ=255)                                                                gCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAAt                                                                    <  1:359439/138‑1 (MQ=255)                                                                     cAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACtt                                                               >  2:185774/1‑138 (MQ=255)                                                                       ggCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCaa                                                            <  1:45653/139‑1 (MQ=255)                                                                            agCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAAcctcc                                                       >  2:161780/1‑135 (MQ=255)                                                                                 gTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCa                                                                                                                                  <  1:79518/59‑1 (MQ=255)                                                                                 gTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCa                                                                                                                                  >  2:79518/1‑59 (MQ=255)                                                                                     gcgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGa                                              <  1:185774/139‑1 (MQ=255)                                                                                                aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacg                                   <  1:11303/139‑1 (MQ=255)                                                                                                     cAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAg                                                                                                                                <  2:167454/41‑1 (MQ=255)                                                                                                     cAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAg                                                                                                                                >  1:167454/1‑41 (MQ=255)                                                                                                       aaGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACt                                                                                                >  2:21874/1‑71 (MQ=255)                                                                                                       aaGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACt                                                                                                <  1:21874/71‑1 (MQ=255)                                                                                                                                 ggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAATCTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGc  <  1:419062/139‑1 (MQ=255)                                                                                                                                 ggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGc  <  1:310714/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       AAACCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGC  >  NZ_CP009273/1260152‑1260418

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