Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 43.5% R263R (CGA→CGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb71,1750AG43.5% 12.3 / 60.6 46R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (10/16);  new base G (7/13);  total (17/29)
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.08e-01

CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGAT  >  W3110S.gb/71106‑71244
                                                                     |                                                                     
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                      >  1:2016630/1‑71 (MQ=255)
cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                      <  1:2612042/71‑1 (MQ=255)
 accaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                     <  1:115213/68‑1 (MQ=255)
 accaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                     <  1:1727710/68‑1 (MQ=255)
   tACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                       >  1:1070536/1‑67 (MQ=255)
   tACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                       >  1:496394/1‑67 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                 <  1:1185505/71‑1 (MQ=255)
     cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                  >  1:1840165/1‑70 (MQ=255)
      ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                                <  1:2499735/71‑1 (MQ=255)
          aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                            <  1:2215116/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                         <  1:354895/71‑1 (MQ=255)
             ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                         <  1:2371329/71‑1 (MQ=255)
              cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                         >  1:1858081/1‑70 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                       >  1:3040216/1‑54 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                                   <  1:2944001/58‑1 (MQ=255)
                aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                      >  1:252249/1‑71 (MQ=255)
                 tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                           >  1:305214/1‑65 (MQ=255)
                   gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                   <  1:1078692/71‑1 (MQ=255)
                            ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                          <  1:379116/71‑1 (MQ=255)
                               cgCAATATTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                    >  1:403118/1‑58 (MQ=255)
                                gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCGcgagc                                      <  1:2709460/71‑1 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                     <  1:2673587/71‑1 (MQ=255)
                                 cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                     <  1:162053/71‑1 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                   <  1:388888/71‑1 (MQ=255)
                                   aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                   <  1:952592/71‑1 (MQ=255)
                                    tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                  <  1:2509365/71‑1 (MQ=255)
                                     gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                             <  1:3081501/43‑1 (MQ=255)
                                     gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                 >  1:2227662/1‑71 (MQ=255)
                                      ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                             >  1:71295/1‑42 (MQ=255)
                                        ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                               <  1:1851105/70‑1 (MQ=255)
                                         gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGt                             >  1:264047/1‑71 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                           <  1:2227702/70‑1 (MQ=255)
                                            ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                          >  1:107809/1‑71 (MQ=255)
                                               aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                        <  1:3089266/70‑1 (MQ=255)
                                                 gATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                     <  1:1861913/71‑1 (MQ=255)
                                                   tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                         >  1:2009393/1‑49 (MQ=255)
                                                       ctctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCGAGCCCGAGCGAGTTTCTTGCCGGTTGTgaaga               <  1:1020578/71‑1 (MQ=255)
                                                        tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                           >  1:345021/1‑58 (MQ=255)
                                                         ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCg                         >  1:1445724/1‑59 (MQ=255)
                                                         ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa              >  1:879287/1‑70 (MQ=255)
                                                          tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                              >  1:1086735/1‑53 (MQ=39)
                                                          tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa             >  1:1178561/1‑70 (MQ=255)
                                                            tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa             <  1:819719/68‑1 (MQ=255)
                                                            tttCTCGAGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa            <  1:1865311/69‑1 (MQ=255)
                                                             ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                 <  1:157868/47‑1 (MQ=255)
                                                                 cgcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTga                  <  1:195370/58‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCAGCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa              <  1:1886725/61‑1 (MQ=255)
                                                                  gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGa       <  1:493843/68‑1 (MQ=255)
                                                                    gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat  >  1:2256777/1‑71 (MQ=255)
                                                                     |                                                                     
CACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGAT  >  W3110S.gb/71106‑71244

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: