Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 93.9% 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,1750AG93.9% 1271.3 / 45.3 407R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base G (194/187);  total (206/201)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 408 total aligned reads displayed.

CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71108‑71242
                                                                   |                                                                   
ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  >  1:1642174/3‑71 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                     >  1:1749868/1‑65 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               <  1:354067/71‑1 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:2470111/1‑70 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:1175113/1‑70 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:3340159/1‑70 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:2165757/1‑70 (MQ=255)
      ggATGCCTATCGCCACCGTCGGTCGTAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:858870/71‑1 (MQ=255)
      ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:3095193/71‑1 (MQ=255)
      ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTCACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            >  1:18433/1‑70 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2924364/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:3561062/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2353674/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:734910/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:3462633/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:1070987/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaag                                                          <  1:2261499/71‑2 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:2008652/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:1862051/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:1196797/1‑70 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1832495/66‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       <  1:2958749/70‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                      <  1:923395/71‑1 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:1193234/1‑59 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                        >  1:1988843/1‑67 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2564347/1‑71 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:654956/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:480359/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:582954/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2442559/1‑70 (MQ=255)
                cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:2743254/61‑1 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:1874005/1‑61 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:643377/1‑61 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             >  1:842943/1‑59 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:1793248/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  <  1:2816219/70‑1 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGTTATTTAAAAAATGCACCGGGGc                                              <  1:449040/71‑1 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              >  1:2052919/1‑71 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              <  1:784997/71‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:536391/1‑57 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       >  1:3179374/1‑59 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCa                                            >  1:3187573/1‑70 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTAACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             >  1:1339056/1‑53 (MQ=39)
                        tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                          >  1:1251397/1‑71 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                        <  1:3133702/71‑1 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGcc                                         >  1:3667240/1‑70 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                        >  1:3484772/1‑71 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  <  1:11492/60‑1 (MQ=255)
                             cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2686663/49‑1 (MQ=255)
                             cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                     >  1:1394939/1‑71 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:2867299/1‑71 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCctagc                                    >  1:2901891/1‑70 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                   <  1:3390733/71‑1 (MQ=255)
                                aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                        >  1:2936512/1‑49 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:1365203/52‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                                >  1:427025/1‑55 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                <  1:360383/71‑1 (MQ=255)
                                   gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                               >  1:1483713/1‑71 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTTGCGGGTATTTAAAAAATGc                                                      <  1:3317465/47‑1 (MQ=39)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  >  1:2308880/1‑35 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:1273360/1‑42 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2322217/42‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:626966/42‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                     >  1:1459464/1‑48 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:678377/1‑52 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcga                                      >  1:1308722/1‑63 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                  >  1:2070529/1‑65 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                            >  1:1538189/1‑71 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                            >  1:1682149/1‑71 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:2185400/1‑36 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATTCACCGGGGCCAGCCcgagcg                                   >  1:1061764/1‑60 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                      <  1:1802465/41‑1 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2788326/1‑43 (MQ=255)
                                            gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCg                       <  1:3243077/70‑1 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                      <  1:2503959/70‑1 (MQ=255)
                                              cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                    >  1:3360494/1‑71 (MQ=255)
                                                 tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCAGAGCGAGTTCCGTGCCGGTtgtg                 <  1:2412550/71‑1 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:3226805/1‑35 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2642306/1‑35 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGcc                                         >  1:340761/1‑46 (MQ=255)
                                                  cAACTCTATTTCTAGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 >  1:100058/1‑70 (MQ=255)
                                                     ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaag              >  1:89943/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaag              >  1:1917494/1‑68 (MQ=255)
                                                        tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                        <  1:2698212/41‑1 (MQ=255)
                                                        tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           >  1:2981779/1‑70 (MQ=255)
                                                        tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa           >  1:1620957/1‑70 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                              <  1:2798621/49‑2 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           >  1:1931122/1‑68 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtg                   <  1:3673667/60‑1 (MQ=255)
                                                           ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:3066081/1‑71 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:3018719/1‑70 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg      <  1:2487365/70‑1 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg      <  1:1572970/70‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCCGCCCGAGCGAGTTCCGTGCCGGTGGTGAAGAAAAAGTGAat   >  1:237288/1‑70 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa         >  1:2209747/1‑64 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa         >  1:3116233/1‑64 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa     <  1:3586373/68‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:2388727/71‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:1424809/71‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat   >  1:1560502/1‑70 (MQ=255)
                                                                   |                                                                   
CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71108‑71242

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