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.6% 158.9 / ‑5.0 42R20C (CGC→TGC) soxRredox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/19);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267367‑4267538
                                                                                       |                                                                                    
tCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGc                                                                                    <  1:590631/90‑1 (MQ=255)
 ccTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCa                                                                                   >  2:747481/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:1374469/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 >  1:138127/1‑90 (MQ=255)
   tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                 <  2:465878/90‑1 (MQ=255)
         tAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGtgtg                                                                           <  1:624435/90‑1 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcg                                                               >  2:483468/1‑90 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATTGgcg                                                               >  1:1009798/1‑90 (MQ=255)
                     aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGgcg                                                               >  1:1510355/1‑90 (MQ=255)
                      aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCTGTGTGGCGGTATCGgcgc                                                              <  1:1391773/90‑1 (MQ=255)
                            ttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                        >  2:1156477/1‑90 (MQ=255)
                                   aaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGa                                                 >  2:70289/1‑90 (MQ=255)
                                    aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                <  2:614021/90‑1 (MQ=255)
                                            aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                        <  1:483468/90‑1 (MQ=255)
                                              cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                      <  2:447867/90‑1 (MQ=255)
                                               ccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTTTGAAAGGAAAGGGtt                                     >  2:396377/1‑90 (MQ=255)
                                                cGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTaa                                           >  1:1150482/1‑83 (MQ=255)
                                                cGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTaa                                           <  2:1150482/83‑1 (MQ=255)
                                                 gCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGa                                   >  2:178930/1‑90 (MQ=255)
                                                  cATTAAAGCGCTGGTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAt                                  <  1:791212/90‑1 (MQ=255)
                                                  cATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAt                                  <  2:1107665/90‑1 (MQ=255)
                                                   aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                 <  2:1124488/90‑1 (MQ=255)
                                                   aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                 <  1:178930/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                              <  2:933473/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                              <  2:829132/90‑1 (MQ=255)
                                                      aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                              <  2:442431/90‑1 (MQ=255)
                                                        aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                            <  2:1234708/90‑1 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                         >  1:92357/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                         >  2:509546/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTAt                         >  2:879382/1‑90 (MQ=255)
                                                           gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCCGCATTTCTATGAAAGTAAAGGGTTGATTACCAGGAt                         >  2:815966/1‑90 (MQ=255)
                                                            ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATTGGCGCTGCATTTTTTTGAAAGGAAAGGGGTGGTTACCAGTATc                        >  1:1606515/1‑90 (MQ=255)
                                                            ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTTTATGAAAGTAAAGGGGTGATTACCAGTATc                        >  1:1354148/1‑90 (MQ=255)
                                                                      ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCg              >  1:993169/1‑90 (MQ=255)
                                                                      ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTAACGGCGCTGCATTTTTATGAAAGAAAAGGGGTGATTACCAGTATCCGTAACAGCg              >  2:1602885/1‑90 (MQ=255)
                                                                       cGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCgg             <  2:1086028/90‑1 (MQ=255)
                                                                        ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGGAAAGGGTTGATTACCAGTATCCGTAACAGCGGc            >  2:659556/1‑90 (MQ=255)
                                                                         gCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                              >  2:1351506/1‑39 (MQ=39)
                                                                         gCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                              <  1:1351506/39‑1 (MQ=39)
                                                                             aGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAGGAAAGGGGTGATTACCAGTATCCGTAACAGCGGCAATCa       >  1:55024/1‑90 (MQ=255)
                                                                               tGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGGAAAAAAAGGGTTGGTTTCCAGGATCCGTAACAGCGGCCATCAgc     >  1:418131/1‑90 (MQ=255)
                                                                                  cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc  <  2:461294/90‑1 (MQ=255)
                                                                                       |                                                                                    
TCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGC  >  NZ_CP009273/4267367‑4267538

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