Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,816,759:1 +G 48.1% coding (417/441 nt) soxR → DNA‑binding transcriptional dual regulator, Fe‑S center for redox‑sensing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,816,7591.G48.1% 38.1 / 68.2 27coding (417/441 nt)soxRDNA‑binding transcriptional dual regulator, Fe‑S center for redox‑sensing
Reads supporting (aligned to +/- strand):  ref base . (11/3);  new base G (7/6);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAA  >  Exported/2816691‑2816828
                                                                     |                                                                     
gACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTaga                                                                     <  1:633351/71‑1 (MQ=255)
      tGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGagaaga                                                                <  1:3084124/71‑1 (MQ=255)
         aTTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAgg                                                             >  1:1235420/1‑70 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTaga                                                                     >  1:2469353/1‑61 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAgaagaa                                                               >  1:3360926/1‑68 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGt                                                            <  1:3301379/71‑1 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGt                                                            <  1:330405/71‑1 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGt                                                            <  1:114862/71‑1 (MQ=255)
          ttGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGt                                                            <  1:4163965/71‑1 (MQ=255)
              ttGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCg                                                        >  1:2987377/1‑71 (MQ=255)
                gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCgg                                                       >  1:3159841/1‑70 (MQ=255)
                gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGgc                                                      >  1:3133796/1‑71 (MQ=255)
                  ggCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTa                                                                        <  1:1357467/51‑1 (MQ=255)
                    cTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCAc                                                  >  1:1460765/1‑71 (MQ=255)
                    cTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCAc                                                  >  1:2044613/1‑71 (MQ=255)
                           tCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTagga                                                                    >  1:2078559/1‑43 (MQ=255)
                           tCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGc                                           >  1:5490880/1‑70 (MQ=255)
                            cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTaga                                                                     >  1:5518853/1‑43 (MQ=255)
                            cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTaga                                                                     >  1:1083296/1‑43 (MQ=255)
                            cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAgg                                                             >  1:872420/1‑51 (MQ=255)
                            cgcgCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAgg                                                             >  1:1182128/1‑51 (MQ=255)
                                cAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCAc                                                  <  1:399633/59‑1 (MQ=255)
                                 aGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAAt                                    >  1:1961746/1‑71 (MQ=255)
                                     aTTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCATTGt                                  <  1:2844756/69‑1 (MQ=255)
                                       tGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTAcc                                                         >  1:3681458/1‑44 (MQ=255)
                                               gCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAgtg                        >  1:2026325/1‑69 (MQ=255)
                                                cGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTa                      <  1:3984659/70‑1 (MQ=255)
                                                 gTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAgg                    <  1:4382118/71‑1 (MQ=255)
                                                  tAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGt                   >  1:1017079/1‑71 (MQ=255)
                                                   aaCCCGGGCGACCGCTTAGGAGAAGAAGCTACCGgcgc                                                    >  1:3559808/1‑38 (MQ=25)
                                                   aaCCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGt                   >  1:2020063/1‑70 (MQ=255)
                                                                  ttAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTg    >  1:166431/1‑70 (MQ=255)
                                                                  ttAGAGAAGAAGGTACCGGCGCACTAGGGGAGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTg    >  1:1358416/1‑70 (MQ=255)
                                                                    agagAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGa                             >  1:1012375/1‑43 (MQ=255)
                                                                    agagAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGaa  <  1:146302/70‑1 (MQ=255)
                                                                    agagAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGa   <  1:1139932/69‑1 (MQ=255)
                                                                     |                                                                     
GACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGAGAAGAAGGTACCGGCGCACTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAA  >  Exported/2816691‑2816828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: