Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,976,020 C→T 59.4% L256F (CTC→TTC)  galR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,976,0200CT59.4% 14.1 / 36.3 32L256F (CTC→TTC) galRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (7/12);  total (13/19)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GACGAGGAAGAAATTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGATGATGTGCTGGTGTCACG  >  W3110S.gb/2975952‑2976090
                                                                    |                                                                      
gACGAGGAAGAAATTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTtctt                                                                      <  1:1427890/71‑2 (MQ=255)
 aCGAGGAAGAAATTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCa                                                                     <  1:592092/71‑1 (MQ=255)
      gaagaaATTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTtctc                                                                      <  1:1956317/65‑1 (MQ=255)
            aTTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGTGTGCGATGGGCGTTTTCAATGATAATgg                                                           >  1:2759429/1‑70 (MQ=255)
            aTTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATgg                                                           >  1:2630362/1‑70 (MQ=255)
               tCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTAtt                                                       >  1:758067/1‑71 (MQ=255)
               tCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTAtt                                                       >  1:754817/1‑71 (MQ=255)
                cACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTAtt                                                       <  1:2564899/70‑1 (MQ=255)
                         aGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCgg                                             >  1:596035/1‑71 (MQ=255)
                               ttATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTga                                         <  1:1381148/69‑1 (MQ=255)
                                 ataACGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTa                                                 >  1:2264783/1‑59 (MQ=255)
                                   aaCGATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAAt                                                                   <  1:16990/39‑1 (MQ=255)
                                     cGATTCAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGt                                 <  1:1504643/71‑1 (MQ=255)
                                      gATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATaa                                                              >  1:1026863/1‑41 (MQ=255)
                                      gATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGtt                                <  1:2373586/71‑1 (MQ=255)
                                       aTTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATgg                                                           <  1:2150466/43‑1 (MQ=255)
                                         tCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAAt                             >  1:1697190/1‑71 (MQ=255)
                                         tCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAAt                             >  1:2382638/1‑71 (MQ=255)
                                          cAATGGCGGCGGGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAAt                             >  1:980899/1‑70 (MQ=255)
                                                cggcggGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCttt                      <  1:2575609/71‑1 (MQ=255)
                                                   cggGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:1014062/71‑1 (MQ=255)
                                                   cggGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:959805/71‑1 (MQ=255)
                                                   cggGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:918640/71‑1 (MQ=255)
                                                   cggGTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:1122561/71‑1 (MQ=255)
                                                    gggTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:2331589/70‑1 (MQ=255)
                                                    gggTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:1098064/70‑1 (MQ=255)
                                                    gggTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGAGAGATTTCGTTAATTGGCTTtgatg                  <  1:707897/71‑1 (MQ=255)
                                                    gggTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgat                   <  1:1484050/70‑1 (MQ=255)
                                                      gTGCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTga                                         <  1:503633/46‑1 (MQ=255)
                                                      gTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGtt                                <  1:1279458/55‑1 (MQ=255)
                                                       tGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTtgatga                 <  1:536633/69‑1 (MQ=255)
                                                        gCGATGGGCGTTTTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAAt                             >  1:172032/1‑56 (MQ=255)
                                                        gCGATGGGAGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGAttt                                    >  1:1980371/1‑49 (MQ=255)
                                                         cGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGatgatg               >  1:330601/1‑69 (MQ=255)
                                                                    ctcAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGATGATGTGCTGGTGTCACg  >  1:143164/1‑71 (MQ=255)
                                                                    ctcAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGATGATGTGCTGGTGTCACg  >  1:2545795/1‑71 (MQ=255)
                                                                    ctcAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGATGATGTGCTGGTGTCACg  >  1:2494455/1‑71 (MQ=255)
                                                                    |                                                                      
GACGAGGAAGAAATTTCACTGCGGTAGCCTGTTATAACGATTCAATGGCGGCGGGTGCGATGGGCGTTCTCAATGATAATGGTATTGATGTACCGGGTGAGATTTCGTTAATTGGCTTTGATGATGTGCTGGTGTCACG  >  W3110S.gb/2975952‑2976090

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: