Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 86.2% 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,1750AG86.2% 800.5 / 93.3 296R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (26/15);  new base G (93/162);  total (119/177)
Fisher's exact test for biased strand distribution p-value = 1.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 298 total aligned reads displayed.

TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71109‑71242
                                                                  |                                                                   
caCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                                 <  1:4412795/70‑1 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:8503/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:1766545/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:798121/1‑70 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:3176606/1‑70 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:1774357/1‑70 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                             <  1:4000321/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:2011780/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             >  1:364457/1‑70 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             >  1:1495680/1‑70 (MQ=255)
      gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2212090/71‑1 (MQ=255)
       aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:3626290/71‑1 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGa                                                                     >  1:3979999/1‑59 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         <  1:3396810/71‑1 (MQ=255)
         gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAAt                                                        >  1:3451157/1‑71 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       >  1:373750/1‑70 (MQ=255)
           cTATAGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1737590/66‑1 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                                 <  1:468499/59‑1 (MQ=255)
            tATCGCCACCGTCGGT‑‑CAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                     >  1:459116/1‑69 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:4525798/1‑71 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2070789/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  >  1:4281884/1‑56 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAt                                                                 >  1:2172146/1‑57 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:715875/71‑1 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:2496467/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:770405/71‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  <  1:2017638/70‑1 (MQ=255)
                   aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                                <  1:3310922/69‑1 (MQ=255)
                   aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              <  1:4482074/71‑1 (MQ=255)
                   aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                      >  1:1191487/1‑63 (MQ=255)
                    ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               >  1:436522/1‑53 (MQ=255)
                    ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              >  1:3441382/1‑70 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:2106364/1‑57 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:2987830/1‑62 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                           <  1:1619655/71‑1 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                           >  1:353647/1‑71 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         <  1:1361308/56‑1 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1262366/52‑1 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                        <  1:4312587/71‑1 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGcc                                         <  1:384736/70‑1 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:1564231/44‑1 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGccc                                        <  1:2550121/70‑1 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:3912613/71‑1 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:2911084/71‑1 (MQ=255)
                           tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgggg                                               <  1:1968539/62‑1 (MQ=255)
                            cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                     >  1:1696269/1‑71 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:1478275/1‑48 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    <  1:3268933/71‑1 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                   <  1:1229082/71‑1 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                   <  1:1840504/71‑1 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:2969849/1‑70 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                     >  1:2346736/1‑69 (MQ=255)
                                aTGTTGGTTTTTACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:1797409/1‑47 (MQ=38)
                                aTGTTGGTTTTGACGATTAACTCTATTTCTCGCGGGTATTTaa                                                             <  1:1150202/43‑1 (MQ=38)
                                aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         <  1:690187/47‑1 (MQ=255)
                                 tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          >  1:1912171/1‑45 (MQ=255)
                                 tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             <  1:2451573/42‑1 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:284195/1‑45 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                  <  1:3603183/35‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:142154/42‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:620653/42‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             <  1:656911/40‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                             <  1:3210422/40‑1 (MQ=255)
                                    tGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                              >  1:3076664/1‑69 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                             <  1:1311434/70‑3 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 <  1:3464148/49‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                      <  1:2737691/43‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTg                          <  1:4576905/71‑1 (MQ=255)
                                        tttGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                            <  1:1057919/36‑1 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  <  1:3061994/45‑1 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                         <  1:1652480/70‑1 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:3345248/1‑38 (MQ=255)
                                          tGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                         <  1:140839/69‑1 (MQ=255)
                                            aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                      <  1:75419/70‑1 (MQ=255)
                                             cGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGtt                    >  1:277176/1‑71 (MQ=255)
                                             cGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                 >  1:2440084/1‑58 (MQ=255)
                                               aTCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                        >  1:3882219/1‑65 (MQ=255)
                                                tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  >  1:4232697/1‑38 (MQ=255)
                                                 cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:2684622/48‑1 (MQ=255)
                                                 cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttcc                             <  1:4302288/58‑3 (MQ=255)
                                                 cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGCCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 >  1:4114243/1‑70 (MQ=255)
                                                   aCTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              <  1:3136965/39‑1 (MQ=255)
                                                    ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaaga             <  1:971361/71‑1 (MQ=255)
                                                    ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaa               >  1:2310111/1‑69 (MQ=255)
                                                        aTTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           >  1:511064/1‑69 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttc                              >  1:3303763/1‑48 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaa               >  1:4555480/1‑64 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa         <  1:1921483/70‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa         <  1:3770704/70‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:1879701/68‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:4041905/68‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:3566865/68‑1 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa           <  1:543315/68‑1 (MQ=255)
                                                           tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:936084/1‑70 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa     <  1:3420447/68‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:2427052/71‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:2917194/71‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:4208536/71‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:1533707/71‑1 (MQ=255)
                                                                 gAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                  >  1:1029154/1‑37 (MQ=255)
                                                                  |                                                                   
TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71109‑71242

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: