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_CP0092734,267,4540CT84.2% 47.0 / 3.2 19R20C (CGC→TGC) soxRredox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (11/5);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267370‑4267538
                                                                                    |                                                                                    
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:208141/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:43543/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:122659/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  2:306375/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGc                                                                                 >  1:449188/1‑90 (MQ=255)
                  aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGgcg                                                               >  1:334368/1‑90 (MQ=255)
                         ttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                        >  2:34931/1‑90 (MQ=255)
                                 aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                <  2:323374/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                      <  1:495950/90‑1 (MQ=255)
                                             cGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGGGGCATTTCTATGAAAGTAAAAGGTTg                                    >  2:594392/1‑90 (MQ=255)
                                              gCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTTTGAAAGGAAAGGGTTGa                                   >  1:258946/1‑90 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                            <  1:95361/90‑1 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                            <  2:43543/90‑1 (MQ=255)
                                                       cgcTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                          <  2:115000/90‑1 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAAGGGTGATTACCAGTAt                         >  2:95361/1‑90 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGAGCGGCGCGGCCTTTCTATGGAAGGAAAGGGGGGGTTACCAGGAt                         >  2:450663/1‑90 (MQ=255)
                                                                cccccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                          <  1:258359/81‑1 (MQ=255)
                                                                cccccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTa                          >  2:258359/1‑81 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc  <  2:92030/90‑1 (MQ=255)
                                                                                    |                                                                                    
TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267370‑4267538

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