Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 61.1% 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,1750AG61.1% 105.6 / 148.7 126R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (13/36);  new base G (46/30);  total (59/67)
Fisher's exact test for biased strand distribution p-value = 2.35e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 127 total aligned reads displayed.

ACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  W3110S.gb/71107‑71245
                                                                    |                                                                      
aCTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                       <  1:4117087/70‑1 (MQ=255)
   aCCCGGATGCCTATCGCCACCGTGGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                     >  1:5046623/1‑69 (MQ=255)
       ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                               <  1:2699782/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2121468/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:194548/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:185860/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:1719826/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:1445837/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:5380408/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:3214709/71‑1 (MQ=255)
        gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:3550743/71‑1 (MQ=255)
         aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              >  1:3376547/1‑70 (MQ=255)
         aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              >  1:705352/1‑70 (MQ=255)
          tGCCTATCTCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1492432/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1520471/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:1392189/1‑71 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:4049776/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:410651/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:872935/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1239591/71‑1 (MQ=255)
          tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:205569/71‑1 (MQ=255)
            ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:4566705/1‑71 (MQ=255)
            ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:843431/1‑71 (MQ=255)
             cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:2542487/71‑1 (MQ=255)
             cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:4585508/71‑1 (MQ=255)
             cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:1504248/71‑1 (MQ=255)
             cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                          >  1:295877/1‑70 (MQ=255)
             cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                         <  1:3252881/71‑1 (MQ=255)
              tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:2957741/70‑1 (MQ=255)
               aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaagaaa                                                            >  1:4074347/1‑66 (MQ=255)
               aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaagaaa                                                            >  1:1742492/1‑66 (MQ=255)
               aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaagaaa                                                            >  1:5336626/1‑66 (MQ=255)
                tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      <  1:2461944/71‑1 (MQ=255)
                tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      >  1:192361/1‑71 (MQ=255)
                tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       <  1:715873/70‑1 (MQ=255)
                tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       <  1:2564827/70‑1 (MQ=255)
                tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       <  1:4213662/70‑1 (MQ=255)
                  gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:4152228/1‑63 (MQ=255)
                     aCCGTCTGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 <  1:3825119/71‑1 (MQ=255)
                     aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                 <  1:2969850/71‑1 (MQ=255)
                      ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:4430054/56‑1 (MQ=255)
                      ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:2868362/56‑1 (MQ=255)
                      ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:3547696/56‑1 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGGACCGGGGCCAg                                              >  1:2342232/1‑71 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       >  1:3168399/1‑62 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       >  1:3175086/1‑62 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                              >  1:3949247/1‑71 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                              >  1:5277390/1‑71 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                              >  1:2137000/1‑71 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                              >  1:487487/1‑71 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                    >  1:1069024/1‑65 (MQ=255)
                        gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                    >  1:4925114/1‑65 (MQ=255)
                           ggtcgCAGTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                           >  1:2403601/1‑71 (MQ=255)
                            gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                          <  1:2367982/71‑1 (MQ=255)
                            gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcg                                          <  1:5583411/71‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTacaa                                                              <  1:4666801/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:3538112/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:1164475/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:4452889/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:3656930/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:2028161/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:1860643/50‑1 (MQ=255)
                             tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:4414149/50‑1 (MQ=255)
                              cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                        >  1:670015/1‑71 (MQ=255)
                              cgCAATGTTCGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                        >  1:2377746/1‑71 (MQ=255)
                                  aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:3200257/1‑71 (MQ=255)
                                  aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:1103428/1‑71 (MQ=255)
                                  aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                              <  1:3895878/45‑1 (MQ=255)
                                   tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                   >  1:4103221/1‑71 (MQ=255)
                                     ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                                 >  1:2092908/1‑70 (MQ=255)
                                     ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                                 >  1:4449536/1‑70 (MQ=255)
                                     ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                                 >  1:2731175/1‑70 (MQ=255)
                                     ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                  <  1:1100549/70‑1 (MQ=255)
                                       ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                                <  1:2214867/70‑1 (MQ=255)
                                         ttttGACGATCAACTCTATTTCTCGCGTGTATTTaa                                                                <  1:273944/36‑1 (MQ=255)
                                           ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                           >  1:1735555/1‑71 (MQ=255)
                                           ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                           >  1:5568315/1‑71 (MQ=255)
                                           ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                           >  1:910648/1‑71 (MQ=255)
                                                   cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                    >  1:3284552/1‑70 (MQ=255)
                                                   cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                    >  1:2801030/1‑70 (MQ=255)
                                                   cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                    >  1:4787991/1‑70 (MQ=255)
                                                     aCTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       <  1:4495271/49‑1 (MQ=255)
                                                       tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaaga                >  1:2724489/1‑70 (MQ=255)
                                                        ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa               >  1:225872/1‑70 (MQ=255)
                                                        ctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa               >  1:2999250/1‑70 (MQ=255)
                                                         tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa              >  1:278649/1‑70 (MQ=255)
                                                           tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaa             <  1:2357573/69‑1 (MQ=255)
                                                             tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                     >  1:138819/1‑43 (MQ=255)
                                                             tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt          >  1:4189318/1‑70 (MQ=255)
                                                              ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa               <  1:990963/64‑1 (MQ=255)
                                                              ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGGGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg         <  1:522545/70‑1 (MQ=255)
                                                                 gcgGGTATTTAAAAAATGCACCGGGGCCAGACCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGaa       >  1:5536501/1‑69 (MQ=255)
                                                                 gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg     <  1:4550850/71‑1 (MQ=255)
                                                                 gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat      <  1:4084637/70‑1 (MQ=255)
                                                                  cgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatga    <  1:3733854/71‑1 (MQ=255)
                                                                  cgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatga    <  1:5332735/71‑1 (MQ=255)
                                                                  cgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatga    <  1:4299688/71‑1 (MQ=255)
                                                                   gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat   <  1:637450/71‑1 (MQ=255)
                                                                   gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat   <  1:1242951/71‑1 (MQ=255)
                                                                    aGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg  <  1:5116838/71‑1 (MQ=255)
                                                                    |                                                                      
ACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  W3110S.gb/71107‑71245

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: