Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 48.9% 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,1750AG48.9% ‑3.1 / 65.4 45R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (8/15);  new base G (11/11);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71108‑71242
                                                                   |                                                                   
cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                  >  1:1348924/1‑71 (MQ=255)
  aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                <  1:2568233/71‑1 (MQ=255)
  aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:715793/1‑71 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:2568440/1‑70 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                >  1:1342329/1‑70 (MQ=255)
    ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                              <  1:1384311/71‑1 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaa                                                             <  1:836530/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                          <  1:2449282/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           >  1:1847379/1‑70 (MQ=255)
           ccTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAtt                                                                >  1:64026/1‑62 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                      <  1:607031/71‑1 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                      <  1:1875249/71‑1 (MQ=255)
              aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:2461463/1‑66 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                     <  1:1272736/69‑1 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                   >  1:133495/1‑71 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGcc                                             <  1:899810/71‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGcc                                             <  1:1335800/71‑1 (MQ=255)
                      cGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                            <  1:15483/55‑1 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAg                                           >  1:1032502/1‑71 (MQ=255)
                       gtcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGAACCGGGGCCAg                                           >  1:2337561/1‑71 (MQ=255)
                          ggtcgCAAAGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCa                                                     >  1:792026/1‑58 (MQ=255)
                             cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                     >  1:84341/1‑71 (MQ=255)
                             cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgag                                     >  1:2639616/1‑71 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                         >  1:605388/1‑50 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    <  1:1031310/71‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    <  1:403397/70‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                    <  1:1327787/70‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                  <  1:1422240/37‑1 (MQ=255)
                                  tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                  >  1:1250789/1‑53 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:1251881/42‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                   <  1:1099216/66‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTATCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                        <  1:581808/61‑1 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                            <  1:2170787/39‑1 (MQ=255)
                                       gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCa                                            >  1:1580476/1‑54 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                           <  1:144830/38‑1 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                         <  1:110243/70‑1 (MQ=255)
                                           tGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                       <  1:16514/39‑1 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGt                     <  1:2572035/71‑1 (MQ=255)
                                                aTCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcga                                      >  1:145331/1‑51 (MQ=255)
                                                  cAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGtt                    >  1:2286597/1‑67 (MQ=255)
                                                        tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa            <  1:884505/69‑1 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCTTGCCGGTTGTGAAGaaa           <  1:2050914/68‑1 (MQ=255)
                                                            tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt       >  1:1323446/1‑70 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGTCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg      <  1:1328153/70‑1 (MQ=255)
                                                                gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatg  <  1:813416/71‑1 (MQ=255)
                                                                   |                                                                   
CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATG  >  W3110S.gb/71108‑71242

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: