Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,621,163 A→G 12.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,1630AG12.9% 71.5 / 5.1 31Q110Q (CAA→CAGmarRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (2/2);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGG  >  W3110S.gb/1621100‑1621230
                                                               |                                                                   
cccGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCAt                                                               >  1:2311622/1‑70 (MQ=255)
   gAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAAt                                                           <  1:2267933/71‑1 (MQ=255)
     aTGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTa                                                         <  1:302045/71‑1 (MQ=255)
        aCAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGtt                                                      >  1:2774926/1‑71 (MQ=255)
        aCAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGtt                                                      >  1:2576427/1‑71 (MQ=255)
         cAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTg                                                     >  1:514244/1‑71 (MQ=255)
                   gTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCt                                           >  1:261744/1‑71 (MQ=255)
                    tACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTa                                                         <  1:84449/56‑1 (MQ=255)
                    tACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc                                                  <  1:626374/63‑1 (MQ=255)
                      cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCa                                        >  1:1735967/1‑71 (MQ=255)
                        ggTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCAc                                       >  1:1053997/1‑70 (MQ=255)
                         gTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAGTGCCATCAATTAGTTGGc                                                   >  1:232614/1‑57 (MQ=255)
                         gTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCAcc                                      >  1:117774/1‑70 (MQ=255)
                           aaaaCTTACCACCGGCGGCGCGGCAATATGTGAACaa                                                                     >  1:2331299/1‑37 (MQ=255)
                            aaaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAg                                   >  1:2262725/1‑70 (MQ=255)
                             aaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc                                                  <  1:781891/54‑1 (MQ=255)
                             aaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGcc                                                  <  1:2719136/54‑1 (MQ=255)
                               cTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCa                                                                >  1:2637712/1‑38 (MQ=255)
                               cTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAAtt                               <  1:2387855/71‑1 (MQ=255)
                                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTaa                             >  1:671374/1‑71 (MQ=255)
                                      ccGGCTGCGCGGCAATATGTGAACAGTGCCATCAATTAGTTGGCCAgg                                               <  1:1155141/48‑1 (MQ=255)
                                      ccGGCGGCGCGGCAATATGTGAACAGTGCCATCAATTAGTTGGCCAgg                                               <  1:1673354/48‑1 (MQ=255)
                                      ccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaa                         >  1:2553015/1‑70 (MQ=255)
                                        ggcggcGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaaa                      <  1:2382362/71‑1 (MQ=255)
                                                   aaTATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGAcggcgg           <  1:404158/71‑1 (MQ=255)
                                                     tatGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGAcggcg            <  1:677474/68‑1 (MQ=255)
                                                       tgtgAACAGTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACg        >  1:257971/1‑70 (MQ=255)
                                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACg        >  1:1878272/1‑70 (MQ=255)
                                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGa       <  1:926929/71‑1 (MQ=255)
                                                        gtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGa                 >  1:630141/1‑60 (MQ=255)
                                                         tgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAg     >  1:1701422/1‑71 (MQ=255)
                                                             cAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTgg  >  1:1796883/1‑70 (MQ=255)
                                                               |                                                                   
CCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGG  >  W3110S.gb/1621100‑1621230

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: