Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 92.8% 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,1750AG92.8% 973.2 / 40.7 317R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base G (145/147);  total (156/161)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 318 total aligned reads displayed.

CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  W3110S.gb/71108‑71245
                                                                   |                                                                      
ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                      <  1:2219885/68‑1 (MQ=255)
  aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                   <  1:4034736/71‑1 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTCACGATCAACTCTATTTCTCGCGGGTATTTaa                                                                <  1:2836734/71‑1 (MQ=255)
      ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:3996237/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:3893847/71‑1 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               >  1:292595/1‑70 (MQ=255)
        aTGCCTATCGCCAACGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                             <  1:4340413/71‑1 (MQ=255)
         tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:375489/71‑1 (MQ=255)
           ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:2049204/1‑71 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:818475/1‑70 (MQ=255)
             tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:3635116/65‑1 (MQ=255)
             tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:1113864/64‑1 (MQ=255)
             tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        >  1:66413/1‑71 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        >  1:124864/1‑70 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:3430010/1‑65 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:1436389/63‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                         <  1:1505779/68‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        <  1:4309890/69‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      <  1:3567614/71‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      <  1:2159919/71‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                       >  1:4413570/1‑70 (MQ=255)
                atccACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:2700689/3‑64 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                    >  1:373033/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                     <  1:438823/70‑1 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCTCGAGTATTTAAAAAATGCACCGGGGc                                                 <  1:922200/71‑1 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                                 <  1:3419835/71‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:3826680/57‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                                <  1:1322685/71‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                         <  1:485802/62‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGAGGGTATTTaaaa                                                              <  1:200268/57‑1 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                           >  1:2515816/1‑58 (MQ=255)
                        tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1814011/56‑1 (MQ=255)
                        tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1743875/56‑1 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                           >  1:2812244/1‑71 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2535180/51‑1 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                     <  1:2888028/60‑1 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                          <  1:1946960/71‑1 (MQ=255)
                            tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2795158/50‑1 (MQ=255)
                             cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                        >  1:1778498/1‑71 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2576210/48‑1 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                       <  1:761341/71‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                               <  1:3521578/46‑1 (MQ=255)
                                aaTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          >  1:2484767/1‑50 (MQ=255)
                                 aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    <  1:239220/71‑1 (MQ=255)
                                 aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:4059650/1‑71 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1675090/46‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          <  1:1161887/48‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                          <  1:715222/48‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                   >  1:267548/1‑71 (MQ=255)
                                   gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcga                                         >  1:2124559/1‑64 (MQ=255)
                                   gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                  >  1:412864/1‑71 (MQ=255)
                                   gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                                  >  1:2094131/1‑71 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            >  1:476824/1‑44 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                      >  1:3398967/1‑50 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                                 <  1:3508334/71‑2 (MQ=255)
                                    ttGGTATTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                        >  1:1004352/1‑64 (MQ=255)
                                      tgTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                                <  1:3293431/69‑3 (MQ=255)
                                       gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                              <  1:2064615/39‑1 (MQ=255)
                                       gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                              <  1:648355/71‑1 (MQ=255)
                                       gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                              <  1:3359586/71‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        <  1:54213/44‑1 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                        >  1:1539709/1‑42 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                    >  1:2947810/1‑62 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                           >  1:3921085/1‑71 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGAACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                           >  1:30543/1‑71 (MQ=255)
                                            gACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                         >  1:697882/1‑71 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                            <  1:1851055/35‑1 (MQ=255)
                                              cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                       <  1:4187736/71‑1 (MQ=255)
                                              cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                       <  1:1078415/71‑1 (MQ=255)
                                                aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                          <  1:386745/50‑1 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                                >  1:232368/1‑56 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                    >  1:3232922/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                    >  1:436762/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGAGAGTTCCGTGCCGGTtgtg                    >  1:1587056/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                    >  1:1200579/1‑70 (MQ=255)
                                                   aaCTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                     >  1:3582054/1‑52 (MQ=255)
                                                   aaCTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                    >  1:392642/1‑69 (MQ=255)
                                                     ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaag                 >  1:1128548/1‑70 (MQ=255)
                                                     ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaag                 >  1:245565/1‑70 (MQ=255)
                                                     ctctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaag                 >  1:3955407/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa               >  1:872985/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa               >  1:1657035/1‑70 (MQ=255)
                                                       ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa               >  1:2900262/1‑70 (MQ=255)
                                                        tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa              >  1:10581/1‑70 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa            <  1:3743704/70‑1 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa              <  1:3098469/68‑1 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa              <  1:3165538/68‑1 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTATGAAGaaaa             <  1:2995377/69‑1 (MQ=255)
                                                          cttCTCGCGGGTATTTAAAAAATGCACCGGGGCCACCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa              <  1:1453064/67‑1 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGGCGGTTGTGAAGAAAAAGt          >  1:3134267/1‑70 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt          <  1:83004/70‑1 (MQ=255)
                                                            tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt          <  1:3203881/70‑1 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg         >  1:3312766/1‑70 (MQ=255)
                                                              tcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg         <  1:2737738/69‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaa                  >  1:1967654/1‑58 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAACGTGa        <  1:2250224/68‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGaa       <  1:4267918/69‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg     <  1:2025521/71‑1 (MQ=255)
                                                                gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAACGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAat      <  1:1611151/70‑1 (MQ=255)
                                                                   aGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgatg  <  1:4095871/71‑1 (MQ=255)
                                                                   |                                                                      
CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATG  >  W3110S.gb/71108‑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: