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,4540CT98.4% 237.6 / ‑3.9 62R20C (CGC→TGC) soxRredox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/37);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCG  >  NZ_CP009273/4267370‑4267539
                                                                                    |                                                                                     
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1061041/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  2:661340/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  2:1831952/90‑1 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1362925/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  2:1128150/90‑1 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1510968/1‑90 (MQ=255)
       aaCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTgg                                                                           <  2:1535764/90‑1 (MQ=255)
        aCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGc                                                                          >  1:807032/1‑90 (MQ=255)
         cTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         <  1:1877216/90‑1 (MQ=255)
           tGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGt                                                                       >  2:213664/1‑90 (MQ=255)
            gAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTa                                                                      <  1:1596196/90‑1 (MQ=255)
              ggTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATc                                                                    <  2:966009/90‑1 (MQ=255)
                 aaaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgc                                                                 >  1:768171/1‑90 (MQ=255)
                  aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcg                                                                >  2:579902/1‑90 (MQ=255)
                   aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                               <  2:1061041/90‑1 (MQ=255)
                   aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                               <  1:1479828/90‑1 (MQ=255)
                      gATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGc                                                            >  2:482309/1‑90 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  2:772234/1‑66 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  1:772234/66‑1 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAt                                                          >  2:561134/1‑90 (MQ=255)
                               aaaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATg                                                   <  1:723555/90‑1 (MQ=255)
                                aaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGa                                                  <  1:1192497/90‑1 (MQ=255)
                                 aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 <  1:215850/90‑1 (MQ=255)
                                         aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                         <  1:1390967/90‑1 (MQ=255)
                                         aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                         <  1:212868/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:1891118/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:733992/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:603027/90‑1 (MQ=255)
                                               cATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAt                                   <  1:1773663/90‑1 (MQ=255)
                                                aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                  <  1:482309/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:415898/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:1095572/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:55347/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               >  2:633684/1‑90 (MQ=255)
                                                     aGCGCTGCTACCCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  2:202943/90‑1 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  1:1070996/90‑1 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  2:1629233/90‑1 (MQ=255)
                                                       cgcTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                           <  2:347622/90‑1 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                         >  1:19534/1‑59 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                         <  2:19534/59‑1 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAGGGGTGATTACCAGTAt                          >  2:1065872/1‑90 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTAt                          >  2:1278567/1‑90 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         <  1:1043357/42‑1 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         >  2:1043357/1‑42 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTATc                         >  2:1259976/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCg               >  1:876652/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAGTAAAGGGTTGGTTACCAGTATCCGTAACAGCg               >  2:1832309/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGGAAAGGGGTGATTACCAGTATCCGTAACAGCg               >  1:763475/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCg               >  1:258136/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGGTTACCAGTATCCGTAACAGCg               >  1:19462/1‑90 (MQ=255)
                                                                    cGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCgg              <  1:1763336/90‑1 (MQ=255)
                                                                    cGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCgg              <  1:599227/90‑1 (MQ=255)
                                                                     ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGGCGCCGCATTTCTCTGAAAGGAAAGGGGTGGTTTCCAGTATCCGTAACAGCggg             >  1:1418614/1‑89 (MQ=255)
                                                                             ggCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 <  2:1506337/46‑1 (MQ=255)
                                                                             ggCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 >  1:1506337/1‑46 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:480271/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:1198779/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  1:27322/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:276548/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:864470/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:91278/90‑1 (MQ=255)
                                                                                gAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggcg  <  1:1278567/90‑1 (MQ=255)
                                                                                    |                                                                                     
TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCG  >  NZ_CP009273/4267370‑4267539

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