Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,267,454	C→T	R20C (CGC→TGC)	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,454	0	C	T	91.7%	76.5 / ‑0.7	24	R20C (CGC→TGC)	soxR	redox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (11/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267370‑4267538                                                                                     |                                                                                     tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:318599/1‑90 (MQ=255) tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGc                                                                                 >  1:491439/1‑90 (MQ=255)         aCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGc                                                                         <  2:431375/90‑1 (MQ=255)             gAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTa                                                                     >  2:105528/1‑90 (MQ=255)               ggTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATc                                                                   <  1:194584/90‑1 (MQ=255)                    aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGgcgc                                                              <  1:611376/90‑1 (MQ=255)                        aTTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGGcgccgcc                                                          >  2:737674/1‑89 (MQ=255)                          ttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                        >  1:576127/1‑90 (MQ=255)                                 aaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGa                                                 >  1:221189/1‑90 (MQ=255)                                  aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                <  1:648895/90‑1 (MQ=255)                                    gAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAg                                              >  1:711298/1‑90 (MQ=255)                                    gAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAg                                              <  2:532841/90‑1 (MQ=255)                                            cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                      <  2:318599/90‑1 (MQ=255)                                               gCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAGGAAAAGGGTGa                                   >  1:465846/1‑90 (MQ=255)                                                 aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                 <  1:170898/90‑1 (MQ=255)                                                 aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                 <  1:412166/90‑1 (MQ=255)                                                      aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                            <  1:42017/90‑1 (MQ=255)                                                                aCCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAc                  <  2:681290/90‑1 (MQ=255)                                                                    ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCg              >  1:463244/1‑90 (MQ=255)                                                                    ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGAAAAGGGGTGATTACCAGTATCCGTAACAGCg              >  2:398768/1‑90 (MQ=255)                                                                      ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGc            >  2:635142/1‑90 (MQ=255)                                                                      ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGGAAAGGGTTGATTACCAGTATCCGTAACAGCGGc            >  2:276620/1‑90 (MQ=255)                                                                                cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc  <  2:268861/90‑1 (MQ=255)                                                                                cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc  <  1:722581/90‑1 (MQ=255)                                                                                     |                                                                                     TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267370‑4267538

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