Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,336,529 C→T 32.3% L321L (CTG→TTG)  cysB → DNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,336,5290CT32.3% 41.2 / 26.8 34L321L (CTG→TTG) cysBDNA‑binding transcriptional dual regulator, O‑acetyl‑L‑serine‑binding
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base T (8/3);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTCGTTGATGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCGGAATATAT  >  W3110S.gb/1336461‑1336595
                                                                    |                                                                  
aTGTCGTTGATGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAATTg                                                                  >  1:606800/1‑71 (MQ=255)
 tGTCGTTGATGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTg                                                                  <  1:342663/70‑1 (MQ=255)
     gTTGATGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCgg                                                              >  1:1582981/1‑70 (MQ=255)
          tGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAATTGCCGGAAAaat                                                        <  1:1428693/71‑1 (MQ=255)
              gCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAAtt                                                    <  1:3387854/71‑1 (MQ=255)
               cTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAAttt                                                   >  1:1156688/1‑71 (MQ=255)
                    gCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGc                                              >  1:2657482/1‑71 (MQ=255)
                       ttGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGtt                                           <  1:1803679/71‑1 (MQ=255)
                       ttGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGtt                                           <  1:1953961/71‑1 (MQ=255)
                                 aTGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGc                                  <  1:1503161/70‑1 (MQ=255)
                                          ttGAGGTCATGTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATg                        >  1:1584315/1‑71 (MQ=255)
                                            gAGGTCATGTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTa                          >  1:427022/1‑67 (MQ=255)
                                            gAGGTCATGTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTa                          >  1:4059962/1‑67 (MQ=255)
                                            gAGGTCATGTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTa                          >  1:1026796/1‑67 (MQ=255)
                                            gAGGTCATGTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACATTTTa                          >  1:1655436/1‑67 (MQ=255)
                                              ggTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGc                    >  1:963652/1‑71 (MQ=255)
                                                tCATGTTTAAAGATATAAAACTGCCGGAAAAATAAttt                                                   <  1:2289761/38‑1 (MQ=255)
                                                  aTGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTATCGaa                 <  1:1131376/70‑1 (MQ=255)
                                                  aTGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGaa                 <  1:408097/70‑1 (MQ=255)
                                                  aTGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAg                >  1:498091/1‑71 (MQ=255)
                                                    gTTTAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTa                          >  1:840073/1‑59 (MQ=255)
                                                     tttAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCttt                            <  1:2411388/56‑1 (MQ=255)
                                                     tttAAAGATATAAAATTGCCGGAAAAATAATTTCTTGCGTTATTTACGGCACCttt                            <  1:143539/56‑1 (MQ=255)
                                                     tttAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGt              <  1:3157589/70‑1 (MQ=255)
                                                        aaaGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGaa                 <  1:718790/64‑1 (MQ=255)
                                                         aaGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCgg         <  1:3943949/71‑1 (MQ=255)
                                                            atatAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCttt                            <  1:3066239/49‑1 (MQ=255)
                                                              ataAAATTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAg                >  1:162515/1‑59 (MQ=255)
                                                               taAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCtttt                           >  1:2654601/1‑47 (MQ=255)
                                                               taAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCGGAatata   >  1:1274478/1‑71 (MQ=255)
                                                                aaaaCTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAg                     >  1:29962/1‑52 (MQ=255)
                                                                aaaaCTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCGGAatatat  >  1:3525451/1‑71 (MQ=255)
                                                                aaaaCTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCGGAatata   >  1:4061850/1‑70 (MQ=255)
                                                                 aaaCTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAg                     >  1:3418251/1‑51 (MQ=255)
                                                                    |                                                                  
ATGTCGTTGATGCGGCTGTCGCATTGCGCTCTAATGAAGAAATTGAGGTCATGTTTAAAGATATAAAACTGCCGGAAAAATAATTTCTTGCGTTATTTTCGGCACCTTTTATGTAGCGAAGGTGCCGGAATATAT  >  W3110S.gb/1336461‑1336595

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: