Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 75.3% 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,1750AG75.3% 150.0 / 48.9 77R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (12/7);  new base G (21/37);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71109‑71242
                                                                  |                                                                   
caCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGt                                                                   <  1:683602/68‑1 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               <  1:1011274/71‑1 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               <  1:2003332/71‑1 (MQ=255)
  cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               <  1:3644601/71‑1 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:2007196/1‑70 (MQ=255)
    cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:2975444/1‑70 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:2143462/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:698088/71‑1 (MQ=255)
     ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                             <  1:1000115/70‑1 (MQ=255)
      gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2915968/71‑1 (MQ=255)
      gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1851331/71‑1 (MQ=255)
      gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1806617/71‑1 (MQ=255)
      gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1846603/71‑1 (MQ=255)
       aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                          <  1:935978/71‑1 (MQ=255)
        tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         <  1:3463076/71‑1 (MQ=255)
           cTATCGCCACCGTCGGTTGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2741365/66‑1 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       >  1:186507/1‑70 (MQ=255)
           cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       >  1:2860176/1‑70 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:179505/64‑1 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:3000558/64‑1 (MQ=255)
            tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              <  1:2599234/62‑1 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:1322673/1‑66 (MQ=255)
             aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:1318406/1‑66 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:1138793/71‑1 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAcc                                                   >  1:276433/1‑71 (MQ=255)
              tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                    >  1:2426338/1‑70 (MQ=255)
               cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                          <  1:318883/63‑1 (MQ=255)
                gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                      >  1:2262510/1‑66 (MQ=255)
                   aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       >  1:1964231/1‑62 (MQ=255)
                    ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1277349/57‑1 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:1931622/1‑57 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          >  1:1349374/1‑56 (MQ=255)
                      gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAg                                           <  1:3435155/71‑1 (MQ=255)
                       tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCgg                                                 >  1:3462380/1‑64 (MQ=255)
                         ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                            >  1:323676/1‑51 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:337621/52‑1 (MQ=255)
                          gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:2238903/52‑1 (MQ=255)
                             gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:3583244/1‑71 (MQ=255)
                              cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                   <  1:1406912/71‑1 (MQ=255)
                                aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:2206397/1‑68 (MQ=255)
                                 tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                      <  1:2849970/49‑1 (MQ=255)
                                  gTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATg                                                       >  1:1305592/1‑47 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:2765588/42‑1 (MQ=255)
                                   ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                               <  1:3074633/70‑1 (MQ=255)
                                    tGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                               <  1:711855/69‑1 (MQ=255)
                                     ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:473870/1‑42 (MQ=255)
                                      gTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCg                            <  1:2835787/70‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:3090550/39‑1 (MQ=255)
                                       ttttGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              <  1:2274297/35‑1 (MQ=255)
                                        tttGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                          >  1:2876739/1‑38 (MQ=255)
                                         ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGcc                        >  1:2735359/1‑71 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                  >  1:3674975/1‑71 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgt                  <  1:998583/71‑1 (MQ=255)
                                               aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                   >  1:2946964/1‑70 (MQ=255)
                                               aTCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:683004/50‑1 (MQ=255)
                                               aTCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:681067/50‑1 (MQ=255)
                                                 cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 >  1:1500891/1‑70 (MQ=255)
                                                    ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTtaaga             <  1:464448/71‑1 (MQ=255)
                                                     tctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcg                                       <  1:837801/44‑1 (MQ=255)
                                                       tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa           >  1:2183814/1‑70 (MQ=255)
                                                       tATTTCTCGCCGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTgaagaa            <  1:2071436/69‑1 (MQ=255)
                                                        aTTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:2304029/69‑1 (MQ=255)
                                                         tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaa           <  1:1556879/68‑1 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                  >  1:2651276/1‑45 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg        >  1:60251/1‑71 (MQ=255)
                                                         tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAg        >  1:2216414/1‑71 (MQ=255)
                                                          ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                               >  1:2671130/1‑47 (MQ=255)
                                                          ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:1911295/1‑71 (MQ=255)
                                                          ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:157549/1‑71 (MQ=255)
                                                          ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:1919365/1‑71 (MQ=255)
                                                           tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAg        >  1:1745237/1‑69 (MQ=255)
                                                            ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg      <  1:2220067/70‑1 (MQ=255)
                                                            ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg      >  1:1833388/1‑70 (MQ=255)
                                                            ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCAGGTTGTGAAGAAAAAGTg      <  1:1561440/70‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa     <  1:3505669/68‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGa     <  1:1041830/68‑1 (MQ=255)
                                                               gcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:1842287/71‑1 (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: