Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,180,586	(A)7→6	intergenic (+333/+77)	yjaZ → / ← thiH	heat‑shock protein C/2‑iminoacetate synthase ThiH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,180,580	0	A	.	100.0%	72.0 / NA	18	intergenic (+327/+83)	yjaZ/thiH	heat‑shock protein C/2‑iminoacetate synthase ThiH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (6/12);  total (6/12)

TTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTC  >  NZ_CP009273/4180491‑4180664                                                                                          |                                                                                     tttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATa                                                                                      >  1:314680/1‑90 (MQ=255) tttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATa                                                                                      >  2:219330/1‑90 (MQ=255)  ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                     >  1:173751/1‑90 (MQ=255)  ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                     >  2:61462/1‑90 (MQ=255)  ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                     >  2:522265/1‑90 (MQ=255)  ttGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATaa                                                                                     >  1:435259/1‑90 (MQ=255)           ggAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACaa                                                                           <  2:85559/90‑1 (MQ=255)              aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc                                                                        <  1:141027/90‑1 (MQ=255)              aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc                                                                        <  2:435259/90‑1 (MQ=255)              aTTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGc                                                                        <  2:314680/90‑1 (MQ=255)                      gAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTcaaca                                                                <  2:34216/90‑1 (MQ=255)                                        aGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCt                                              >  1:388052/1‑90 (MQ=255)                                                           tGGACAACATTCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATa                           >  1:316057/1‑90 (MQ=255)                                                                    ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                  <  1:61462/90‑1 (MQ=255)                                                                    ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                  <  1:457353/90‑1 (MQ=255)                                                                    ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                  <  2:436136/90‑1 (MQ=255)                                                                    ttCCTGCTGACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGt                  >  2:526285/1‑90 (MQ=255)                                                                            gACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt          >  2:160374/1‑90 (MQ=255)                                                                            gACGGCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACAt          >  2:487777/1‑90 (MQ=255)                                                                                gCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg      <  2:343771/90‑1 (MQ=255)                                                                                gCACTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCg      <  1:263511/90‑1 (MQ=255)                                                                                   cTACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGtct   >  1:104326/1‑90 (MQ=255)                                                                                    tACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGtctc  <  2:150245/90‑1 (MQ=255)                                                                                    tACCAT‑AAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGtctc  <  2:158653/90‑1 (MQ=255)                                                                                          |                                                                                     TTTGCATTACGGAATTTTTAAGAAAGGCAGGGCGAAACGAGGAAGAAGCTTTTCGCCCTGGACAACATTCCTGCTGACGGCACTACCATAAAAAAACACAACGCTTTCAACAACCTTGCGTTTTTCATCTTCCGTAGCGGCTCGGAATAAATTTACGTTCTTACATTCCGTCTC  >  NZ_CP009273/4180491‑4180664

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