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%	45.0 / NA	12	coding (1733/2055 nt)	helD	DNA helicase IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/2);  total (10/2)

TCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCG  >  NZ_CP009273/1021586‑1021732                                                                          |                                                                          tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                          >  2:316189/1‑90 (MQ=255) tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt                                                          >  2:423631/1‑90 (MQ=255)        aTGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGaa                                                   >  1:272710/1‑90 (MQ=255)          gCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGt                                                 >  1:310355/1‑90 (MQ=255)                   gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg                                        <  1:337486/90‑1 (MQ=255)                        gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                   >  1:20647/1‑90 (MQ=255)                        gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt                                   >  1:433514/1‑90 (MQ=255)                             ttCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGAc                              >  1:233310/1‑90 (MQ=255)                             ttCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGAc                              <  2:22568/90‑1 (MQ=255)                                       ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACATGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                    >  1:356792/1‑90 (MQ=255)                                       ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGcc                    >  2:425721/1‑90 (MQ=255)                                                         gAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCg  >  1:156718/1‑90 (MQ=255)                                                                          |                                                                          TCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCG  >  NZ_CP009273/1021586‑1021732

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