Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 53.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,1750AG53.6% 16.4 / 94.7 69R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (16/16);  new base G (21/16);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

CCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71111‑71242
                                                                |                                                                   
cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGCCGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:4595221/1‑70 (MQ=255)
cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAttt                                                               <  1:4924780/71‑1 (MQ=255)
  cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:2049963/1‑70 (MQ=255)
  cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:142055/1‑70 (MQ=255)
  cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                              >  1:2835920/1‑70 (MQ=255)
   ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:4952911/71‑1 (MQ=255)
   ggATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:1688615/71‑1 (MQ=255)
    gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:590316/71‑1 (MQ=255)
    gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaa                                                            <  1:2216210/70‑1 (MQ=255)
    gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                           <  1:770481/71‑1 (MQ=255)
      tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:3679482/1‑71 (MQ=255)
      tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         <  1:4619317/71‑1 (MQ=255)
       gCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:3623724/1‑70 (MQ=255)
         cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                               <  1:5118188/62‑1 (MQ=255)
         cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                      <  1:639620/71‑1 (MQ=255)
            tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:3245917/1‑71 (MQ=255)
            tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:2939909/1‑65 (MQ=255)
            tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                      <  1:312314/68‑1 (MQ=255)
            tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                     <  1:4338651/69‑1 (MQ=255)
             cGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:2557590/70‑1 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:4378513/1‑71 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:2995813/1‑71 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:3979051/1‑71 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                 >  1:2748996/1‑71 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                  <  1:4116022/70‑1 (MQ=255)
              gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                  <  1:4626912/70‑1 (MQ=255)
                 aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              >  1:1986567/1‑71 (MQ=255)
                 aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              >  1:1480940/1‑71 (MQ=255)
                 aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                              >  1:1378304/1‑71 (MQ=255)
                  ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                <  1:2964134/52‑1 (MQ=255)
                    gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                        >  1:2569794/1‑58 (MQ=255)
                    gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                    >  1:1428579/1‑62 (MQ=255)
                       ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGcc                                         <  1:3573467/70‑1 (MQ=255)
                       ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGcc                                         <  1:1022362/70‑1 (MQ=255)
                       ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                        >  1:979548/1‑71 (MQ=255)
                           gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:268511/1‑71 (MQ=255)
                            cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                         >  1:676465/1‑49 (MQ=255)
                            cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    >  1:3176918/1‑70 (MQ=255)
                              aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:3502374/53‑1 (MQ=255)
                              aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   <  1:2366958/53‑1 (MQ=255)
                              aTGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGcc                                             <  1:2851870/59‑1 (MQ=255)
                                gTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         <  1:3839293/45‑1 (MQ=255)
                                 ttGGTTTTGACGATCAACTCTATTTCTCTCGGGTATTTAAAAAATGCAc                                                    >  1:3520796/1‑49 (MQ=39)
                                 ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                              >  1:127547/1‑39 (MQ=255)
                                 ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                    >  1:4076089/1‑49 (MQ=255)
                                  tGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                     >  1:1105150/1‑47 (MQ=255)
                                    gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                             >  1:916765/1‑69 (MQ=255)
                                      tttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                                <  1:1575808/64‑1 (MQ=255)
                                       ttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:145693/1‑38 (MQ=255)
                                       ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGc                         >  1:3295469/1‑70 (MQ=255)
                                       ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                        >  1:124252/1‑71 (MQ=255)
                                        tGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCAAGCCCGAGCGAGTTTCGTGc                         <  1:3806366/69‑1 (MQ=255)
                                            gATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtg                   <  1:4419078/71‑1 (MQ=255)
                                               cAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:3397762/1‑36 (MQ=255)
                                               cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                 <  1:3854123/70‑1 (MQ=255)
                                               cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                 >  1:1963803/1‑70 (MQ=255)
                                                  ctctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                 >  1:2617323/1‑51 (MQ=255)
                                                    ctATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGtt                               >  1:2928221/1‑51 (MQ=255)
                                                     tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa           >  1:4364241/1‑70 (MQ=255)
                                                      aTTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                              <  1:4595894/50‑1 (MQ=255)
                                                       tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa         <  1:298711/70‑1 (MQ=255)
                                                        ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgt                  <  1:1103441/60‑1 (MQ=255)
                                                        ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgt                  <  1:3043087/60‑1 (MQ=255)
                                                         tctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGt       >  1:3028893/1‑70 (MQ=255)
                                                          ctcGCGAGTATTTAAAAAATGCACCGGGGCCCGCCCGAGCGAGTTTc                             >  1:1290516/1‑47 (MQ=255)
                                                          ctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                             >  1:2421730/1‑47 (MQ=255)
                                                             gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGa     <  1:3685429/68‑1 (MQ=255)
                                                             gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:5034420/71‑1 (MQ=255)
                                                             gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat   <  1:2446440/70‑1 (MQ=255)
                                                                aGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaaga             <  1:1170574/57‑1 (MQ=255)
                                                                |                                                                   
CCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71111‑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: