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%	37.7 / NA	10	coding (930/3744 nt)	narG	nitrate reductase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (4/6);  total (4/6)

GTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGC  >  NZ_CP009273/1276160‑1276313                                                                                          |                                                                 gTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATg                                                                  <  1:194717/90‑1 (MQ=255)   cacaCCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATggc                                                                <  2:346483/90‑2 (MQ=255)              aCGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgct                                                    <  1:570079/90‑1 (MQ=255)              aCGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAAtgct                                                    >  2:340082/1‑90 (MQ=255)                  tGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGt                                                <  1:76992/90‑1 (MQ=255)                                 gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                 >  1:195276/1‑90 (MQ=255)                                 gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                 >  1:97451/1‑90 (MQ=255)                                 gtgCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGAc                                 <  2:584526/90‑1 (MQ=255)                                                           aGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGt                                                >  1:173541/1‑49 (MQ=39)                                                           aGGGCACCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGt                                                <  2:173541/49‑1 (MQ=39)                                                                cccGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCCCCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc  <  1:426070/89‑1 (MQ=255)                                                                aCCGATGCGGCAATGGCGCTGGCGAT‑GGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGc  <  1:340082/90‑1 (MQ=255)                                                                                          |                                                                 GTCACACCAGACTACGCTGAAATCGCCAAACTGTGCGATCTGTGGCTGGCACCGAAACAGGGCACCGATGCGGCAATGGCGCTGGCGATGGGCCACGTAATGCTGCGTGAATTCCACCTCGACAACCCAAGCCAGTATTTCACCGACTATGTGC  >  NZ_CP009273/1276160‑1276313

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