Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 51.0% 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,1750AG51.0% 7.4 / 151.0 100R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (28/23);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGA  >  W3110S.gb/71109‑71239
                                                                  |                                                                
tACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                 <  1:4696067/68‑1 (MQ=255)
 aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                             >  1:5079582/1‑71 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                             >  1:4932606/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAttt                                                            <  1:4729385/71‑1 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                             >  1:347565/1‑70 (MQ=255)
   ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTCGACGATCAACTCTATTTCTCGCGAGTATTTa                                                           <  1:2074955/71‑1 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGGAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                          <  1:864814/71‑1 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                           >  1:4358680/1‑70 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                           >  1:4452925/1‑70 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                          <  1:1231211/71‑1 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                          <  1:3402968/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                         <  1:3290660/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                          >  1:26944/1‑70 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                      <  1:4580726/71‑1 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                          >  1:2979554/1‑67 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                     >  1:2496396/1‑71 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                     >  1:46387/1‑71 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                     >  1:1509206/1‑71 (MQ=255)
          ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                     <  1:848042/70‑1 (MQ=255)
          ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                     <  1:2854518/70‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                       >  1:1375812/1‑67 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                   <  1:749223/71‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                        <  1:2073162/66‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                    >  1:2381933/1‑70 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                   <  1:1610402/71‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                   <  1:2938396/71‑1 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                        <  1:5042605/65‑1 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                 >  1:2987378/1‑71 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                 >  1:1565511/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                            <  1:774211/59‑1 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                >  1:3711458/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                 <  1:2572815/70‑1 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                               >  1:1041618/1‑56 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                      >  1:4437815/1‑65 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                >  1:3685345/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                >  1:2872669/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                >  1:2878560/1‑71 (MQ=255)
               cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                               >  1:4212238/1‑71 (MQ=255)
               cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                        >  1:2666086/1‑62 (MQ=255)
               cGCCACCGTCGGGCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                        <  1:1197116/62‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                       >  1:2304616/1‑62 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                  >  1:776072/1‑67 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                              >  1:2710250/1‑71 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                               <  1:1207484/70‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                               <  1:3375048/70‑1 (MQ=255)
                   aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                           <  1:4242944/71‑1 (MQ=255)
                    ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                           >  1:2831086/1‑70 (MQ=255)
                    ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                           >  1:2050604/1‑70 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                          >  1:3568744/1‑53 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                    <  1:3126727/58‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                             >  1:2635148/1‑65 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTACTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGc                                       >  1:3658814/1‑71 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                     <  1:55688/71‑1 (MQ=255)
                            cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                  >  1:4817024/1‑71 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                 <  1:333993/71‑1 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                      >  1:1636590/1‑50 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                         <  1:1137093/47‑1 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                 <  1:4352816/71‑1 (MQ=255)
                              cAATTTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                <  1:4712182/71‑1 (MQ=255)
                                 tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                           <  1:746002/41‑1 (MQ=255)
                                 tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                           <  1:3335422/41‑1 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                    <  1:4642783/63‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                         >  1:1815307/1‑41 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcga                                   >  1:3411731/1‑63 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcg                                    <  1:1745770/62‑1 (MQ=255)
                                    tGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                          >  1:3894022/1‑69 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                        >  1:3111832/1‑40 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGTCCAGCCCGAGCGAGTTCCGt                        <  1:3408110/71‑1 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                        <  1:4404721/71‑1 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                         <  1:1756262/70‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                        <  1:653005/70‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                    <  1:1857762/42‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                 <  1:2350642/45‑1 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                      >  1:3739742/1‑38 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                        <  1:4004445/36‑1 (MQ=255)
                                             cGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGtt                 <  1:2076762/71‑1 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt               <  1:1750837/71‑1 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt               <  1:191678/71‑1 (MQ=255)
                                                tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg              <  1:2906562/71‑1 (MQ=255)
                                                tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt               >  1:2548475/1‑70 (MQ=255)
                                                tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt               >  1:372287/1‑70 (MQ=255)
                                                tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                <  1:2359736/37‑1 (MQ=255)
                                                 cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg              >  1:3411358/1‑70 (MQ=255)
                                                  aaCTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                              >  1:1527815/1‑37 (MQ=255)
                                                     tctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaaga          >  1:4533805/1‑70 (MQ=255)
                                                       tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa         <  1:1437815/69‑1 (MQ=255)
                                                       tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa         >  1:1612485/1‑69 (MQ=255)
                                                       tATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa        >  1:4989020/1‑70 (MQ=255)
                                                       tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa        >  1:1500544/1‑70 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                  >  1:2164353/1‑42 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                  >  1:381522/1‑42 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg     >  1:802016/1‑71 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa       <  1:1067804/69‑1 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg     >  1:3928155/1‑71 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg     >  1:3457057/1‑71 (MQ=255)
                                                          ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt    >  1:3157002/1‑71 (MQ=255)
                                                           tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt    >  1:1572377/1‑70 (MQ=255)
                                                             tcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg   <  1:4548727/69‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa  <  1:4469619/68‑1 (MQ=255)
                                                               gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGa  <  1:5098316/68‑1 (MQ=255)
                                                                  |                                                                
TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGA  >  W3110S.gb/71109‑71239

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: