Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 71,175 A→G 43.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,1750AG43.6% 11.6 / 48.6 39R263R (CGA→CGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (11/11);  new base G (12/5);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71120‑71241
                                                       |                                                                  
cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGc                                                     <  1:1377600/71‑1 (MQ=255)
cTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGc                                                     <  1:485740/71‑1 (MQ=255)
 tATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                          <  1:2308309/65‑1 (MQ=255)
  aTCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAc                                                   >  1:1606547/1‑71 (MQ=255)
   tCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCAcc                                                  >  1:551478/1‑71 (MQ=255)
     gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTaaaaaa                                                        >  1:2478229/1‑63 (MQ=255)
     gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                >  1:1106114/1‑71 (MQ=255)
     gCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCgg                                                >  1:1706246/1‑71 (MQ=255)
                 cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                 >  1:1062619/1‑42 (MQ=255)
                 cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                 >  1:397928/1‑42 (MQ=255)
                 cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                 >  1:1858321/1‑42 (MQ=255)
                 cgCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTa                                                                 >  1:1584059/1‑42 (MQ=255)
                      tGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                               <  1:1502223/71‑1 (MQ=255)
                      tGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGt                               <  1:751863/71‑1 (MQ=255)
                        ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:2510283/42‑1 (MQ=255)
                        ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCa                                           >  1:1479666/1‑57 (MQ=255)
                        ttGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                             >  1:2467546/1‑71 (MQ=255)
                        ttGGTTTTGACGATAAACTCTATTTCTCGAGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGttt                             >  1:1660433/1‑71 (MQ=255)
                         tGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTaaaaaa                                                        >  1:195452/1‑43 (MQ=255)
                          ggTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCg                           >  1:1277143/1‑71 (MQ=255)
                            ttttGACGATCAACTCTATTTCTCGCGGGTATTTaaaa                                                          <  1:2395074/38‑1 (MQ=255)
                             tttGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCa                                           >  1:2553499/1‑52 (MQ=255)
                              ttGACGATCAACTCTATTTCTCTCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                 >  1:49579/1‑61 (MQ=255)
                              ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCcgagcga                                 >  1:1369304/1‑61 (MQ=255)
                              ttGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCAcgagcga                                 >  1:781460/1‑61 (MQ=255)
                              ttGACGATCAACTCTATTTCTCGCGAGTATTTaaaa                                                          <  1:2198747/36‑1 (MQ=255)
                                      cAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTtgtg                >  1:214346/1‑70 (MQ=255)
                                              tttCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGaaaaa        <  1:491442/70‑1 (MQ=255)
                                              tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcga                                     <  1:938420/41‑1 (MQ=255)
                                              tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa         <  1:2254790/69‑1 (MQ=255)
                                              tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa         <  1:1703738/69‑1 (MQ=255)
                                              tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaa         <  1:1298185/69‑1 (MQ=255)
                                              tttCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaa          <  1:103623/68‑1 (MQ=255)
                                                tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGc                        >  1:575059/1‑52 (MQ=255)
                                                tctcGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGt      <  1:51593/70‑1 (MQ=255)
                                                 ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     >  1:478242/1‑70 (MQ=255)
                                                 ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     >  1:2154272/1‑70 (MQ=255)
                                                 ctcGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTg     >  1:2011509/1‑70 (MQ=255)
                                                    gcgAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat  <  1:38205/70‑1 (MQ=255)
                                                       |                                                                  
CTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAAT  >  W3110S.gb/71120‑71241

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: