Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,621,163 A→G 22.9% Q110Q (CAA→CAG marR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,621,1630AG22.9% 62.6 / 17.6 35Q110Q (CAA→CAGmarRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base G (3/5);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTG  >  W3110S.gb/1621095‑1621229
                                                                    |                                                                  
ccGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAt                                                                   <  1:3720261/70‑1 (MQ=255)
  gAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGc                                                                 >  1:2399779/1‑70 (MQ=255)
      ccGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAGTGCCATCa                                                            <  1:3302962/71‑1 (MQ=255)
          aTGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAAt                                                          >  1:3385683/1‑69 (MQ=255)
                 gcgcGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc                                                 <  1:3829417/71‑1 (MQ=255)
                           cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAGTGCCATCAAt                                                          <  1:1687423/52‑1 (MQ=255)
                           cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAGTGCCATCAAt                                                          <  1:2583303/52‑1 (MQ=255)
                           cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCa                                       >  1:3970894/1‑71 (MQ=255)
                             ggTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAg                                                       >  1:773399/1‑53 (MQ=255)
                             ggTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCAc                                      >  1:1328126/1‑70 (MQ=255)
                                 aaaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGACCAGGACCTGCACCAAg                                  >  1:1853340/1‑70 (MQ=255)
                                      tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTaa                            >  1:774184/1‑71 (MQ=255)
                                      tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTa                             >  1:102001/1‑70 (MQ=255)
                                            cggcggCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaa                      <  1:1022657/71‑1 (MQ=255)
                                                        aaTATGTGAACAATGCCATCAATTAGTTGGCCAGGa                                             <  1:84579/36‑1 (MQ=255)
                                                        aaTATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCt                  <  1:1136092/63‑1 (MQ=255)
                                                          tatGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAg                                  <  1:2431969/45‑1 (MQ=255)
                                                          tatGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaa                       <  1:97486/56‑1 (MQ=255)
                                                          tatGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaa                       <  1:1421643/56‑1 (MQ=255)
                                                           atGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGAc        >  1:371662/1‑70 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaaa                     >  1:2979526/1‑56 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaaa                     >  1:2793403/1‑56 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaaa                     >  1:2732815/1‑56 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACg       >  1:1200789/1‑70 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACg       >  1:1261261/1‑70 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGa      <  1:1895657/71‑1 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGa      <  1:1202468/71‑1 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGa      <  1:3798351/71‑1 (MQ=255)
                                                            tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGCCGGACGa      <  1:756310/71‑1 (MQ=255)
                                                             gtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACg       >  1:310628/1‑69 (MQ=255)
                                                              tgAACAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAg    >  1:404271/1‑71 (MQ=255)
                                                              tgAACAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAg    >  1:3668750/1‑71 (MQ=255)
                                                              tgAACAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAg    >  1:3177659/1‑71 (MQ=255)
                                                                aaCAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTg  <  1:1296476/71‑1 (MQ=255)
                                                                 aCAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTg  <  1:650595/70‑1 (MQ=255)
                                                                    |                                                                  
CCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTG  >  W3110S.gb/1621095‑1621229

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: