Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 42.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,1750AG42.2% 29.6 / 81.6 66R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base G (15/12);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGAT  >  W3110S.gb/71108‑71244
                                                                   |                                                                     
ccaCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTa                                                                    <  1:3416887/69‑1 (MQ=255)
cTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                    >  1:705401/1‑71 (MQ=255)
 caCCCGAATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGTGTAt                                                                   >  1:547224/2‑71 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTAttt                                                                 <  1:951733/71‑1 (MQ=255)
   cccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                  >  1:3542454/1‑70 (MQ=255)
    ccGGATGCCTATCGCCACCGTCTGTCGCAATGTTGGTTTTCACGATCAACTCTATTTCTCGCGAGTATTTa                                                                <  1:298300/71‑1 (MQ=255)
    ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTa                                                                <  1:3949740/71‑1 (MQ=255)
    ccGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTAtt                                                                  >  1:2751546/1‑69 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                               <  1:3751669/71‑1 (MQ=255)
     cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTa                                                                >  1:3125834/1‑70 (MQ=255)
       gATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                             <  1:751961/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                            <  1:2290105/71‑1 (MQ=255)
        aTGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaa                                                            <  1:3152559/71‑1 (MQ=255)
         tGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                           >  1:1218039/1‑71 (MQ=255)
            cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCTGGTATTTAAAAAATg                                                         >  1:2833781/1‑70 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                             >  1:131515/1‑63 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaa                                                               >  1:491374/1‑61 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                             >  1:2015997/1‑63 (MQ=255)
               tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCAc                                                      <  1:2898709/70‑1 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATg                                                         >  1:3723252/1‑65 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                   >  1:2833011/1‑71 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                          >  1:1753461/1‑64 (MQ=255)
                 gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAAt                                                          >  1:3060263/1‑64 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGc                                                >  1:583418/1‑71 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaa                                                              <  1:2811099/57‑1 (MQ=255)
                    aCCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                <  1:4088816/71‑1 (MQ=255)
                     ccGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGcc                                               <  1:1652661/71‑1 (MQ=255)
                        tcggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                           >  1:2523793/1‑56 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                          >  1:3225115/1‑71 (MQ=255)
                          ggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGccc                                          <  1:1786845/71‑1 (MQ=255)
                           gtcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCg                                                    <  1:768738/60‑1 (MQ=255)
                            tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                            <  1:2131187/51‑1 (MQ=255)
                            tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                             >  1:218950/1‑50 (MQ=255)
                            tcgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                             >  1:3729442/1‑50 (MQ=255)
                              gCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                      <  1:3504127/71‑1 (MQ=255)
                               cAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                     <  1:3022986/71‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                             <  1:2790782/42‑1 (MQ=255)
                                    ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                                <  1:1513424/71‑1 (MQ=255)
                                     aggTTTTGACGATCAACTCTATTTCTCGCGAGAATATAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                               >  1:1286315/2‑71 (MQ=255)
                                      ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                              >  1:2361623/1‑71 (MQ=255)
                                       gTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                              <  1:3033364/70‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCg                                                    <  1:3406336/47‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                             <  1:491639/38‑1 (MQ=255)
                                        ttttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagc                                      >  1:60677/1‑61 (MQ=255)
                                         tttGACGATCAACTCTATTTCTCGCGGGTATTTaaaaa                                                            >  1:4197789/1‑38 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCggg                                                  >  1:261025/1‑47 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGcc                          >  1:2503736/1‑71 (MQ=255)
                                          ttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                           >  1:190643/1‑70 (MQ=255)
                                            gACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                           >  1:2575341/1‑36 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                        <  1:2633316/70‑1 (MQ=255)
                                             aCGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCgg                        <  1:456783/70‑1 (MQ=255)
                                               gATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCa                                                       >  1:3986331/1‑37 (MQ=255)
                                                 tCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                   <  1:580786/71‑1 (MQ=255)
                                                 tCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTtgtg                   <  1:2593089/71‑1 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                   >  1:3086634/1‑70 (MQ=255)
                                                  cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCAGTGCCGGTtgtg                   >  1:3248483/1‑70 (MQ=255)
                                                   aaCTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                     >  1:1864788/1‑51 (MQ=255)
                                                     ctctATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCAAGCCCGAGCGAGTTTCg                              <  1:1725224/56‑1 (MQ=255)
                                                        tATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTgaagaa              <  1:4012599/69‑1 (MQ=255)
                                                          tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa           <  1:707243/70‑1 (MQ=255)
                                                          tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa             <  1:4173866/68‑1 (MQ=255)
                                                           ttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCgg                        <  1:2696449/56‑1 (MQ=255)
                                                           ttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt         >  1:1534808/1‑71 (MQ=255)
                                                             ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa           >  1:1214854/1‑67 (MQ=255)
                                                               cgcgGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTg                            >  1:818054/1‑48 (MQ=38)
                                                                  gAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAatgat  >  1:791121/1‑71 (MQ=255)
                                                                   |                                                                     
CTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGAT  >  W3110S.gb/71108‑71244

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: