Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,267,830	G→A	G145D (GGC→GAC)	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,830	0	G	A	100.0%	103.9 / NA	29	G145D (GGC→GAC)	soxR	redox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)

GGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCC  >  NZ_CP009273/4267757‑4267916                                                                          |                                                                                       ggTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGGACCGACGCACGCTTGCTGGaa                                                                        >  2:511837/1‑90 (MQ=255)   ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGa                                                                      >  2:471733/1‑90 (MQ=255)   ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGa                                                                      >  2:329360/1‑90 (MQ=255)      tgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGa                                                                   <  1:297683/90‑1 (MQ=255)        gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAAc                                                                 <  1:591278/90‑1 (MQ=255)        gCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAAc                                                                 <  2:171850/90‑1 (MQ=255)                 cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTaaa                                                        <  1:373243/90‑1 (MQ=255)                   cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAgc                                                      <  2:664264/90‑1 (MQ=255)                   cgcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAgc                                                      <  1:92141/90‑1 (MQ=255)                              cccGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                           <  1:256901/90‑1 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCt                                                                                <  2:237148/50‑1 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCt                                                                                >  1:237148/1‑50 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGa                                                                      >  1:470709/1‑60 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGa                                                                      <  2:470709/60‑1 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAg                                              >  1:400864/1‑84 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAg                                              <  2:400864/84‑1 (MQ=255)                                 gTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCt                                        >  2:177739/1‑90 (MQ=255)                                            cGGGCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgttt                             >  2:272084/1‑90 (MQ=255)                                              ggCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTc                           >  1:292014/1‑90 (MQ=255)                                               gCGACCGCTTAGGAGAAGAAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTcc                          <  2:691591/90‑1 (MQ=255)                                                            agaagaAGGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtcttt             <  2:527957/90‑1 (MQ=255)                                                               agaagGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc          <  1:302862/90‑1 (MQ=255)                                                               agaagGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtctc          >  2:73930/1‑90 (MQ=255)                                                                gaagGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTa                                     >  1:408680/1‑62 (MQ=255)                                                                gaagGTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTa                                     <  2:408680/62‑1 (MQ=255)                                                                    gTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTAt     >  2:472256/1‑90 (MQ=255)                                                                    gTACCGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTAt     >  2:644558/1‑90 (MQ=255)                                                                       ccGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATccc  <  1:73930/90‑1 (MQ=255)                                                                       ccGACGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATccc  <  2:50788/90‑1 (MQ=255)                                                                          |                                                                                       GGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCC  >  NZ_CP009273/4267757‑4267916

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: ‑