Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 82.4% 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,1750AG82.4% 176.7 / 28.4 74R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (8/5);  new base G (30/31);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

GCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71105‑71241
                                                                      |                                                                  
gCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                    <  1:1815916/71‑1 (MQ=255)
   cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                 >  1:450453/1‑71 (MQ=255)
     aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                               >  1:358461/1‑71 (MQ=255)
      cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAttt                                                              <  1:1897393/71‑1 (MQ=255)
       ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                               >  1:555337/1‑69 (MQ=255)
       ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                             <  1:239753/71‑1 (MQ=255)
         ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                            >  1:376259/1‑70 (MQ=255)
         ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                            >  1:1177035/1‑70 (MQ=255)
         ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                            >  1:69064/1‑70 (MQ=255)
          gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:1902630/71‑1 (MQ=255)
           aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                         <  1:1306988/71‑1 (MQ=255)
           aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                         <  1:540607/71‑1 (MQ=255)
            tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1906604/1‑71 (MQ=255)
              ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                          <  1:174777/67‑1 (MQ=255)
              ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                       <  1:831059/70‑1 (MQ=255)
               cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                      >  1:1714721/1‑70 (MQ=255)
                 aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:1692178/1‑71 (MQ=255)
                 aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:2006375/1‑71 (MQ=255)
                 aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:388667/1‑71 (MQ=255)
                 aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                    >  1:1869750/1‑70 (MQ=255)
                  tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   <  1:1643373/70‑1 (MQ=255)
                  tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                  <  1:1219658/71‑1 (MQ=255)
                    gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTAGCGGGTATTTaaaaaa                                                        >  1:1289806/1‑63 (MQ=255)
                     ccACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                               <  1:1825424/71‑1 (MQ=255)
                      cACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgggg                                              >  1:269532/1‑71 (MQ=255)
                      cACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGCTATTTa                                                             <  1:269792/56‑1 (MQ=255)
                       aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                          >  1:977052/1‑58 (MQ=255)
                       aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                             >  1:313702/1‑71 (MQ=255)
                        ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                               >  1:1350398/1‑52 (MQ=255)
                        ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:228606/1‑64 (MQ=255)
                         cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                             >  1:245249/1‑69 (MQ=255)
                          gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                          >  1:630529/1‑71 (MQ=255)
                           tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                         <  1:1094223/71‑1 (MQ=255)
                           tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGc                                         <  1:1934411/71‑1 (MQ=255)
                               tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                           <  1:1710045/49‑1 (MQ=255)
                               tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                  <  1:2056138/58‑1 (MQ=255)
                               tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                      <  1:483256/70‑1 (MQ=255)
                                cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                 <  1:1276120/58‑1 (MQ=255)
                                 gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                 <  1:489596/57‑1 (MQ=255)
                                 gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                    <  1:687198/70‑1 (MQ=255)
                                  cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                  <  1:1693745/71‑1 (MQ=255)
                                   aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:1743920/1‑48 (MQ=255)
                                    aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                             >  1:311731/1‑42 (MQ=255)
                                     tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                <  1:1573475/70‑1 (MQ=255)
                                     tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                      >  1:810371/1‑64 (MQ=255)
                                     tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                               <  1:785870/71‑1 (MQ=255)
                                     tGTAGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                 >  1:401655/1‑53 (MQ=39)
                                      gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                             <  1:853212/40‑1 (MQ=255)
                                       ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                         <  1:768900/43‑1 (MQ=255)
                                       ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             <  1:1222028/71‑2 (MQ=255)
                                       ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                             >  1:1458797/1‑70 (MQ=255)
                                       ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                             >  1:1935980/1‑39 (MQ=255)
                                         ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgggg                                              >  1:459373/1‑52 (MQ=255)
                                         ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCa                                           >  1:1225591/1‑55 (MQ=255)
                                         ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                           >  1:1302041/1‑71 (MQ=255)
                                             ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                >  1:1902921/1‑46 (MQ=255)
                                                 cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                   <  1:706791/71‑1 (MQ=255)
                                                 cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                   >  1:317255/1‑71 (MQ=255)
                                                 cGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGtt                   >  1:1428177/1‑71 (MQ=255)
                                                 ataTCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGtt                   <  1:1951380/69‑1 (MQ=255)
                                                  gATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                               <  1:1096505/42‑1 (MQ=255)
                                                    tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                        >  1:1757723/1‑63 (MQ=255)
                                                     cAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:1553699/1‑35 (MQ=255)
                                                     cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                               >  1:933781/1‑39 (MQ=255)
                                                      aaCTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTga               <  1:1596120/70‑1 (MQ=255)
                                                        ctctATTTCTCGCGGGTATTTAAAAAATGCACCggg                                               <  1:1862314/36‑1 (MQ=255)
                                                           tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                          >  1:1355350/1‑54 (MQ=39)
                                                           tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa           <  1:244994/69‑1 (MQ=255)
                                                           tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa          >  1:1631793/1‑70 (MQ=255)
                                                             tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg       <  1:730556/71‑1 (MQ=255)
                                                                ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGATCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:1731177/70‑1 (MQ=255)
                                                                ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     <  1:618662/70‑1 (MQ=255)
                                                                  cgcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat  >  1:924580/1‑71 (MQ=255)
                                                                   gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa    <  1:1656277/68‑1 (MQ=255)
                                                                   gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat  <  1:1436858/70‑1 (MQ=255)
                                                                      |                                                                  
GCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71105‑71241

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: