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,4540CT97.7% 316.3 / ‑1.6 85R20C (CGC→TGC) soxRredox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (34/49);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATAT  >  NZ_CP009273/4267367‑4267543
                                                                                       |                                                                                         
tCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGc                                                                                         >  1:1612384/1‑90 (MQ=255)
tCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGc                                                                                         >  1:981255/1‑90 (MQ=255)
 ccTCAAGTTAACTTGAGGTAAAGCGATTTTTGGAAAAGAAATTACCCCGCATTAAAGCGTGGTTAACCCCCGGCGAGGGGGCGAAATGCa                                                                                        >  1:22339/1‑90 (MQ=255)
 ccTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCa                                                                                        >  1:624337/1‑90 (MQ=255)
 ccTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCa                                                                                        >  2:1108504/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:1081311/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:422673/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:1545619/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:603136/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      <  1:606360/90‑1 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  2:1286509/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      <  2:1362375/90‑1 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      <  2:1403352/90‑1 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  2:399135/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:113311/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:11178/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  2:82272/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      <  1:1094143/90‑1 (MQ=255)
         tAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGtgtg                                                                                >  2:1263341/1‑90 (MQ=255)
         tAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGtgtg                                                                                >  2:1034156/1‑90 (MQ=255)
          aaCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTgg                                                                               <  2:1610216/90‑1 (MQ=255)
            cTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                             <  2:853882/90‑1 (MQ=255)
            cTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                             <  2:180822/90‑1 (MQ=255)
            cTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                             <  1:482394/90‑1 (MQ=255)
                 ggTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATc                                                                        <  2:213383/90‑1 (MQ=255)
                    aaaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgc                                                                     <  2:397742/90‑1 (MQ=255)
                    aaaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgc                                                                     <  2:222825/90‑1 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      <  2:1509997/72‑1 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                      >  1:1509997/1‑72 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGgcg                                                                    >  2:301139/1‑90 (MQ=255)
                      aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                                   <  2:1369307/90‑1 (MQ=255)
                      aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                                   <  1:996094/90‑1 (MQ=255)
                      aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                                   <  1:793424/90‑1 (MQ=255)
                       gCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCt                                                                  <  2:1542008/90‑1 (MQ=255)
                          gattATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGGCGCTGCa                                                               >  2:1135290/3‑90 (MQ=255)
                            ttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                             <  1:1351190/90‑1 (MQ=255)
                                 gAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTAt                                                        <  2:732907/90‑1 (MQ=255)
                                   aaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGa                                                      <  1:1108504/90‑1 (MQ=255)
                                    aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                     <  1:859306/90‑1 (MQ=255)
                                           tACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAg                                              >  1:702295/1‑90 (MQ=255)
                                            aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                             <  1:795978/90‑1 (MQ=255)
                                            aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                             <  1:1466072/90‑1 (MQ=255)
                                            aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                             <  2:1562919/90‑1 (MQ=255)
                                              cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                           <  2:1612384/90‑1 (MQ=255)
                                              cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGAGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                           <  1:1183484/90‑1 (MQ=255)
                                                cGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTg                                         <  1:1023579/90‑1 (MQ=255)
                                                cGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTTTGAAAGGAAAGGGGTg                                         >  1:582745/1‑90 (MQ=255)
                                                 gCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAGATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGa                                        >  2:290028/1‑90 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGTGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                                   <  1:374653/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAGATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                                   <  1:290028/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                                   <  2:1623341/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                                   <  2:724733/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                                   <  2:981255/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                                 <  1:518560/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                                 <  1:1319466/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                                 <  2:867758/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                                 <  2:216510/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                                 <  1:969208/90‑1 (MQ=255)
                                                          cgcTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                               <  2:1620482/90‑1 (MQ=255)
                                                          cgcTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTAGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                               <  1:82272/90‑1 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAGGGGTGGTTACCAGTAt                              >  2:792214/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAGGGGTGATTACCAGTAt                              >  1:388899/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTAt                              >  2:1379715/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                              >  1:164541/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                              >  2:1240978/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGAAAAGGGGTGATTACCAGTAt                              >  1:6783/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCCGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                              >  1:1539933/1‑90 (MQ=255)
                                                            ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                             <  2:1159053/74‑1 (MQ=255)
                                                            ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                             >  1:1159053/1‑74 (MQ=255)
                                                            ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATc                             >  2:1643314/1‑90 (MQ=255)
                                                                  aCCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAc                       <  2:22339/90‑1 (MQ=255)
                                                                      ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCg                   >  1:140467/1‑90 (MQ=255)
                                                                      ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCg                   >  2:795197/1‑90 (MQ=255)
                                                                      ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCg                   >  1:631749/1‑90 (MQ=255)
                                                                        ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                           >  2:1525782/1‑64 (MQ=255)
                                                                        ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                           <  1:1525782/64‑1 (MQ=255)
                                                                        ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGc                 >  1:272083/1‑90 (MQ=255)
                                                                          cGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCaa               <  2:1530256/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  2:631749/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  1:446606/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  2:272083/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  1:795197/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  2:11178/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  2:1545619/90‑1 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc       <  1:1301529/90‑1 (MQ=255)
                                                                                       cGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGatat  <  2:582745/90‑1 (MQ=255)
                                                                                       |                                                                                         
TCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATAT  >  NZ_CP009273/4267367‑4267543

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