Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 38.6% 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,1750AG38.6% 52.9 / 80.9 70R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (28/15);  new base G (13/14);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTG  >  W3110S.gb/71109‑71238
                                                                  |                                                               
caCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                             >  1:2150668/2‑71 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                            >  1:980952/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                            >  1:3345712/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                >  1:4132440/1‑66 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                >  1:1855544/1‑66 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                          >  1:3346138/1‑70 (MQ=255)
       aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGCTCACCTCTATTTCTCGCGGGTATTTaaaa                                                       >  1:1205066/1‑70 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       <  1:4226235/69‑1 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                 >  1:3702672/1‑58 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                 >  1:2494988/1‑58 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                     >  1:59106/1‑70 (MQ=255)
           cTATCGGCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                     <  1:3224269/68‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       >  1:1121588/1‑66 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                       <  1:1374349/65‑1 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                             <  1:883257/58‑1 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                     >  1:1678217/1‑66 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAg                                                                >  1:3470380/1‑54 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                    >  1:1933549/1‑66 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                               >  1:2990095/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                               >  1:3312486/1‑71 (MQ=255)
               cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAttt                                                           >  1:2056664/1‑58 (MQ=255)
                gCCTCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                              <  1:2853105/70‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                              <  1:799921/70‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       <  1:4057528/61‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                             >  1:1419535/1‑71 (MQ=255)
                     cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                 >  1:578197/1‑46 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                       <  1:3795373/71‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                       <  1:250615/54‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       <  1:802674/54‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       <  1:1653351/54‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                   >  1:3339722/1‑58 (MQ=255)
                        cggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                 >  1:2813030/1‑43 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                    <  1:556268/71‑1 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                    <  1:3878045/71‑1 (MQ=255)
                            cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGGCAGCCcgag                                 >  1:2708356/1‑71 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                               <  1:2218992/71‑1 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                               <  1:1579799/71‑1 (MQ=255)
                                aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                 >  1:2599589/1‑51 (MQ=255)
                                aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                             >  1:2392995/1‑71 (MQ=255)
                                aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                       >  1:502042/1‑45 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                              >  1:1865853/1‑36 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                              >  1:1140077/1‑36 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAt                                                             >  1:3873959/1‑37 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                 >  1:4510186/1‑49 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                 >  1:3856555/1‑49 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                          <  1:588841/71‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                          <  1:3411186/71‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                          <  1:2581233/71‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                           <  1:4067587/70‑1 (MQ=255)
                                    tGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                             >  1:2756746/1‑67 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                        >  1:258159/1‑71 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                           >  1:2408461/1‑68 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                         <  1:3586305/70‑1 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGAGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtttg                         >  1:1062828/1‑69 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGt                       <  1:1303982/71‑1 (MQ=255)
                                       ttttGACGATCAACTCTATATCTCGCGGGTATTTAAAAAATGCACCg                                              <  1:2312371/47‑1 (MQ=38)
                                        tttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATTCACCGGGGc                                          >  1:634302/1‑50 (MQ=255)
                                        tttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                     >  1:203092/1‑71 (MQ=255)
                                            aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                              >  1:763498/1‑58 (MQ=255)
                                            aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGt                 >  1:4163312/1‑71 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGc                                      <  1:3243773/47‑1 (MQ=255)
                                                tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgggg                                           >  1:2425586/1‑41 (MQ=255)
                                                tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgt              >  1:1750223/1‑70 (MQ=255)
                                                 cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGtt                >  1:911298/1‑67 (MQ=255)
                                                   aCTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaag          <  1:760634/71‑1 (MQ=255)
                                                     tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt              <  1:2807304/65‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa       <  1:2696112/68‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa       <  1:1183091/68‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg    >  1:2353664/1‑71 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa      <  1:14887/69‑1 (MQ=255)
                                                          ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaaa     >  1:1238171/1‑69 (MQ=255)
                                                          ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaaa     >  1:3435122/1‑69 (MQ=255)
                                                          ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt   >  1:1918284/1‑71 (MQ=255)
                                                            ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg  <  1:1275977/70‑1 (MQ=255)
                                                                  |                                                               
TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTG  >  W3110S.gb/71109‑71238

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: