Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 85.7% 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,1750AG85.7% 534.0 / 65.6 196R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base G (82/86);  total (98/98)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 196 total aligned reads displayed.

CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGA  >  W3110S.gb/71108‑71243
                                                                   |                                                                    
ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                   <  1:2015853/69‑1 (MQ=255)
ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                   <  1:861460/69‑1 (MQ=255)
cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                   >  1:3054661/1‑71 (MQ=255)
    ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                               <  1:2961507/71‑1 (MQ=255)
    ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                               <  1:3145270/71‑1 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                              <  1:3566624/71‑1 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                               >  1:3468396/1‑70 (MQ=255)
      ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                             <  1:3071279/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:1597403/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:1223447/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:340112/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                           <  1:1626472/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                             >  1:589783/1‑69 (MQ=255)
         tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                          >  1:3429273/1‑71 (MQ=255)
          gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                         >  1:455169/1‑71 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:2933662/71‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:212606/71‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:3508406/71‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:3783425/71‑1 (MQ=255)
             tATCGCCACCTTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                      >  1:3541086/1‑71 (MQ=255)
             tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                           >  1:2535564/1‑66 (MQ=255)
              aTCGCCACCGTCGGTCGCATTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                           <  1:1337354/65‑1 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                        >  1:3993593/1‑68 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                        >  1:1074115/1‑68 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                     >  1:755174/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                          >  1:397453/1‑65 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                    >  1:295297/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                    <  1:967619/71‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                    >  1:1514066/1‑71 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                      <  1:1734233/69‑1 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCTGGTATTTAAAAAATGCACCgg                                                  >  1:3932766/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                        >  1:2455500/1‑65 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                     >  1:1906670/1‑68 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                  >  1:3319920/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                  >  1:1351929/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                   >  1:1189893/1‑70 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                   <  1:3200050/70‑1 (MQ=255)
                      cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:2112540/56‑1 (MQ=255)
                      cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:1305900/56‑1 (MQ=255)
                      cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:2061724/56‑1 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                            >  1:1447937/1‑71 (MQ=255)
                         cggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                                  >  1:3295111/1‑47 (MQ=255)
                         cggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                           >  1:2707174/1‑70 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            >  1:263041/1‑52 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:84271/51‑1 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcg                                        <  1:2722681/71‑1 (MQ=255)
                            tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                            <  1:3464366/50‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:1384325/52‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                    <  1:2143950/71‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                    <  1:2501603/71‑1 (MQ=255)
                               cAATGTTGGGTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:1039894/52‑1 (MQ=255)
                                 aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                  >  1:3839881/1‑71 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                        <  1:1226953/48‑1 (MQ=255)
                                    ttGGTTTTGCCGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                      >  1:534361/1‑64 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCTGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               >  1:1744669/1‑70 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:3099763/47‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                     >  1:3236169/1‑49 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                              <  1:2640102/56‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               >  1:449095/1‑70 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               >  1:894713/1‑70 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               >  1:1424335/1‑70 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               >  1:3698265/1‑70 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                          >  1:1806821/1‑44 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                               >  1:634727/1‑71 (MQ=255)
                                    ttGGATTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                      >  1:724033/1‑48 (MQ=39)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                          >  1:1294483/1‑42 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                              >  1:124431/1‑70 (MQ=255)
                                       gTTTTGACTATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCTAGCGAGTTCCGt                            <  1:362315/71‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                             <  1:1768653/37‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                       <  1:2681512/43‑1 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                      >  1:4062899/1‑42 (MQ=255)
                                           tGACGATCAACTCTATTTCTCGCTGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                         <  1:1167077/70‑1 (MQ=255)
                                           tGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                  <  1:2785881/61‑1 (MQ=255)
                                            gACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCg                        <  1:362033/70‑1 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGt                      <  1:3022512/71‑1 (MQ=255)
                                              cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCCAGCGAGTTCCGTGCCGGt                      <  1:2605422/70‑1 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                              >  1:3866258/1‑42 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcga                                       >  1:1472888/1‑49 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCTTGCCGGTtgtg                  >  1:2372601/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                  >  1:1573308/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                  >  1:2748298/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                  >  1:3909423/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                  >  1:3946965/1‑70 (MQ=255)
                                                      tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                  >  1:2232170/1‑50 (MQ=255)
                                                      tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaaga              >  1:3800176/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                               <  1:3341984/52‑2 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa             >  1:808606/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa             >  1:3757841/1‑70 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCGCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaa           <  1:1811100/69‑1 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg         >  1:388333/1‑71 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGt                            <  1:1144883/52‑1 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt        >  1:3354269/1‑70 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt        >  1:1703322/1‑70 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt        >  1:1831365/1‑70 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg       <  1:2847083/70‑1 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg       <  1:3264059/70‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg   <  1:3632720/71‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat    <  1:2838239/70‑1 (MQ=255)
                                                                gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat    <  1:1781637/70‑1 (MQ=255)
                                                                 cgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatga  <  1:4039173/71‑1 (MQ=255)
                                                                   |                                                                    
CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGA  >  W3110S.gb/71108‑71243

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