Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,443,410 G→A 31.6% L179L (CTG→TTG)  ytfJ ← predicted transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,443,4100GA31.6% 81.2 / 45.0 57L179L (CTG→TTG) ytfJpredicted transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base A (14/4);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGATAACGATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCA  >  W3110S.gb/4443341‑4443477
                                                                     |                                                                   
tAGATAACGATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTgcag                                                                     >  1:4582448/1‑70 (MQ=255)
tAGATAACGATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTgcag                                                                     >  1:2602940/1‑70 (MQ=255)
  gATAACGATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGcaacag                                                                  >  1:3863894/1‑71 (MQ=255)
  gATAACGATAAGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGcagcag                                                                  >  1:4253407/1‑71 (MQ=255)
        gATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCAt                                                            <  1:1189093/71‑1 (MQ=255)
        gATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCAt                                                            <  1:3418634/71‑1 (MQ=255)
        gATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCAt                                                            >  1:3678353/1‑71 (MQ=255)
         aTAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATc                                                           <  1:3440673/71‑1 (MQ=255)
         aTAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCAt                                                            >  1:4208238/1‑70 (MQ=255)
         aTAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATc                                                           <  1:5025293/71‑1 (MQ=255)
               gAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTtgc                                                     <  1:1344924/71‑1 (MQ=255)
                 aaaTTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCACTtgctg                                                   <  1:3029105/71‑1 (MQ=255)
                  aaTTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTtgctg                                                   >  1:3366899/1‑70 (MQ=255)
                    ttattTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACt                                                        <  1:125499/63‑1 (MQ=255)
                    ttattTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCAc                                                         >  1:4003263/1‑62 (MQ=255)
                    ttattTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCACt                                                        >  1:4657667/1‑63 (MQ=255)
                    ttattTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCACTTGCTGCAc                                                <  1:5125062/71‑1 (MQ=255)
                    ttattTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGAAGCAGGTCCATCACt                                                        <  1:1745639/63‑1 (MQ=255)
                      attTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACt                                                        <  1:2031007/61‑1 (MQ=255)
                      attTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACt                                                        <  1:4388598/61‑1 (MQ=255)
                             cAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTgcag                                                                     <  1:2122106/41‑1 (MQ=255)
                              aGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgt                                      <  1:1700697/71‑1 (MQ=255)
                                aaTTTATGTTATTGCCATTATTTGAGCAACTTTTgcagc                                                                    <  1:2892934/39‑1 (MQ=255)
                                aaTTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgtgt                                    <  1:2845017/71‑1 (MQ=255)
                                aaTTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgtg                                     >  1:5118150/1‑70 (MQ=255)
                                   ttATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAg                                 >  1:1363280/1‑71 (MQ=255)
                                   ttATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAg                                 >  1:4119916/1‑71 (MQ=255)
                                    tATGTTATTGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCAATtgc                                                     >  1:4416050/1‑50 (MQ=255)
                                        ttATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAg                                 >  1:2418824/1‑66 (MQ=255)
                                            tGCCATTATTTGAGCAACTTTTGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAgcg                               >  1:126621/1‑64 (MQ=255)
                                              ccATTATTTGAGCAACTTTTGCAACAGGTCCATCa                                                          >  1:4175008/1‑35 (MQ=255)
                                              ccATTATTTGAGCAACTTTTGCAACAGGTCCATCAc                                                         >  1:1403497/1‑36 (MQ=255)
                                              ccATTATTTGAGCAACTTTTGCAACAGGTCCATCAc                                                         >  1:1371383/1‑36 (MQ=255)
                                              ccATTATTTGAGCAACTTTTGCAACAGGTCCATCAc                                                         >  1:3231219/1‑36 (MQ=255)
                                               cATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCttt                     <  1:3293584/71‑1 (MQ=255)
                                                attattTGAGCAACTTTTGCAGCAGGTCCATCACTtgctg                                                   <  1:3573717/40‑1 (MQ=255)
                                                attattTGAGCAACTTTTGCAGCAGGTCCATCACTtgctg                                                   <  1:620321/40‑1 (MQ=255)
                                                 ttattTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACctct                                            <  1:59518/46‑1 (MQ=255)
                                                 ttattTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACctct                                            <  1:408346/46‑1 (MQ=255)
                                                  tattTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCa                                                 >  1:2768341/1‑40 (MQ=255)
                                                  tattTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCa                                                 >  1:2605365/1‑40 (MQ=255)
                                                   attTGAGCAACTTTTGCAGCAGGTCCATCACTtgct                                                    >  1:796637/1‑36 (MQ=255)
                                                   attTGAGCAACTTTTGCAGCAGGTCCATCACTtgct                                                    >  1:2842641/1‑36 (MQ=255)
                                                   attTGAGCAAATTTTGCAGCAGGTCCATCACTtgctgc                                                  <  1:4447274/38‑1 (MQ=255)
                                                    tttGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGccc                <  1:3063958/71‑1 (MQ=255)
                                                    tttGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGccc                <  1:3895497/71‑1 (MQ=255)
                                                    tttGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGccc                <  1:2809598/71‑1 (MQ=255)
                                                    tttGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGcc                 <  1:3269259/70‑1 (MQ=255)
                                                     ttGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCAc                                                >  1:4668892/1‑38 (MQ=255)
                                                          cAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCa          <  1:687215/71‑1 (MQ=255)
                                                           aaCTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCAc         <  1:2637776/71‑1 (MQ=255)
                                                            aCTTTTGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCAcg        >  1:4915825/1‑71 (MQ=255)
                                                             cTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgt                                      <  1:723034/40‑1 (MQ=255)
                                                               tttGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgtgt                                    >  1:1952984/1‑40 (MQ=255)
                                                               tttGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGgtgt                                    >  1:649190/1‑40 (MQ=255)
                                                               tttGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGa     >  1:4821292/1‑71 (MQ=255)
                                                               tttGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGa     >  1:5119964/1‑71 (MQ=255)
                                                               tttGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCtt                      <  1:5089025/54‑1 (MQ=255)
                                                               tttGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGa     >  1:3063095/1‑71 (MQ=255)
                                                               tttGCAACAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCACGTCTTTGGCCCATTGCa          >  1:1412369/1‑66 (MQ=255)
                                                                 tGCAGCGGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGAc    >  1:658750/1‑70 (MQ=255)
                                                                 tGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGAc    >  1:1643923/1‑70 (MQ=255)
                                                                 tGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGAc    >  1:1220122/1‑70 (MQ=255)
                                                                   cagcagGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCa  <  1:779065/70‑1 (MQ=255)
                                                                     |                                                                   
TAGATAACGATAGGGGAAAATTATTTTTGCAGAATTTATGTTATTGCCATTATTTGAGCAACTTTTGCAGCAGGTCCATCACTTGCTGCACCTCTTCCGGTGTGAGCGCCCCGTCTTTGGCCCATTGCACGCGACCA  >  W3110S.gb/4443341‑4443477

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: