Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,267,806	A→G	D137G (GAC→GGC)	soxR →	redox‑sensitive transcriptional activator SoxR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,267,806	0	A	G	100.0%	57.6 / NA	19	D137G (GAC→GGC)	soxR	redox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/8);  total (11/8)

TAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCT  >  NZ_CP009273/4267722‑4267878                                                                                     |                                                                         tAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCt                                                                     >  2:260952/1‑90 (MQ=255) tAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCt                                                                     >  2:173167/1‑90 (MQ=255)     ggCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAgg                                                                 >  2:439326/1‑90 (MQ=255)      gcgcTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGga                                                                >  2:339582/1‑90 (MQ=255)          tGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaag                                                            <  2:58508/90‑1 (MQ=255)          tGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaag                                                            <  1:339582/90‑1 (MQ=255)          tGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGagaag                                                            >  2:302334/1‑90 (MQ=255)            cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAgaagaa                                                          <  1:344606/90‑1 (MQ=255)            cGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAgaagaa                                                          <  1:439326/90‑1 (MQ=255)                  gAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTAcc                                                    <  1:73434/90‑1 (MQ=255)                                 aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGGGAAGAAGGTACCGGCGCACGCTTGCTg                                     >  2:272012/1‑90 (MQ=255)                                      ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCtt                                                                    >  1:171620/1‑54 (MQ=255)                                      ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCtt                                                                    <  2:171620/54‑1 (MQ=255)                                      ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGa                                <  1:323367/90‑1 (MQ=255)                                         tgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGa                             >  1:165499/1‑90 (MQ=255)                                                    cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTaaa                  <  1:277177/90‑1 (MQ=255)                                                        cAGTGATTGCCCGTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAgcgc              >  1:392453/1‑90 (MQ=255)                                                                    gTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt  >  2:124525/1‑90 (MQ=255)                                                                    gTTGCGTAACCCGGGCGGCCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt  >  2:289272/1‑90 (MQ=255)                                                                                     |                                                                         TAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCT  >  NZ_CP009273/4267722‑4267878

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