Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,276,251	Δ1 bp	coding (932/3744 nt)	narG →	nitrate reductase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,276,249	0	G	.	100.0%	77.1 / NA	19	coding (930/3744 nt)	narG	nitrate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/7);  total (12/7)

CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCT  >  NZ_CP009273/1276162‑1276319                                                                                        |                                                                       cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                      <  1:472115/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                      <  1:66516/90‑2 (MQ=255) cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                      <  2:385075/90‑2 (MQ=255)        gACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTaa                                                              >  2:396976/1‑90 (MQ=255)            aCGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgct                                                          >  1:239862/1‑90 (MQ=255)             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgctg                                                         >  2:429364/1‑90 (MQ=255)             cGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgctg                                                         >  2:114132/1‑90 (MQ=255)                      cGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTc                                                >  1:26461/1‑90 (MQ=255)                       gCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTcc                                               >  2:76322/1‑90 (MQ=255)                               gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                       >  1:397785/1‑90 (MQ=255)                                         tggcttgcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                             <  2:461243/90‑1 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                             <  2:79545/90‑1 (MQ=255)                                         tggctggcACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                             <  2:19842/90‑1 (MQ=255)                                                  ccGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCa                    >  1:54075/1‑90 (MQ=255)                                                    gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                             <  2:375215/79‑1 (MQ=255)                                                    gAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGcc                             >  1:375215/1‑79 (MQ=255)                                                         aGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa             >  1:325843/1‑90 (MQ=255)                                                              aCCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc        <  2:485858/90‑1 (MQ=255)                                                               ccGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa             >  1:75561/1‑84 (MQ=255)                                                               ccGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTa             <  2:75561/84‑1 (MQ=255)                                                                  aTGCGGCAATGGCGCTGGCAAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCg    <  1:114132/90‑1 (MQ=255)                                                                    gCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCt  >  1:211670/1‑90 (MQ=255)                                                                                        |                                                                       CACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGCGTCGCT  >  NZ_CP009273/1276162‑1276319

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