Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,021,663	(A)5→4	coding (1737/2055 nt)	helD →	DNA helicase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,021,659	0	A	.	100.0%	23.9 / NA	7	coding (1733/2055 nt)	helD	DNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (6/1);  total (6/1)

TGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACA  >  NZ_CP009273/1021572‑1021728                                                                                        |                                                                      tGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaa                                                                     >  2:527181/1‑90 (MQ=255)               tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                      >  1:724646/1‑90 (MQ=255)                ctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTg                                                     >  1:473963/1‑90 (MQ=255)                                      gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                               >  1:242096/1‑90 (MQ=255)                                        gcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTAt                             >  1:621983/1‑90 (MQ=255)                                         cATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATg                            <  2:65319/90‑1 (MQ=255)                                                          gTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa           >  2:320310/1‑90 (MQ=255)                                                                   acaTGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACa  >  1:129932/1‑90 (MQ=255)                                                                                        |                                                                      TGCTGGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACA  >  NZ_CP009273/1021572‑1021728

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