Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,730,340 C→T 7.9% T75T (ACC→ACT bglG → transcriptional antiterminator of the bgl operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,730,3400CT 7.9% 275.0 / 10.3 89T75T (ACC→ACTbglGtranscriptional antiterminator of the bgl operon
Reads supporting (aligned to +/- strand):  ref base C (29/53);  new base T (4/3);  total (33/56)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGCTAA  >  W3110S.gb/3730276‑3730410
                                                                |                                                                      
gTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGAt                                                                  >  1:2821896/1‑71 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:3629041/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:3183446/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:3066987/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                >  1:2861177/1‑70 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                >  1:1951302/1‑70 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:1328333/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                >  1:2161289/1‑70 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                >  1:2553602/1‑70 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:3847857/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:2258344/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:533951/70‑1 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                >  1:886890/1‑70 (MQ=255)
   aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                                <  1:2242801/70‑1 (MQ=255)
    tgaactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCg                                                                <  1:1088135/69‑1 (MQ=255)
      aactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                            <  1:3190745/71‑1 (MQ=255)
      aactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                            <  1:2231491/71‑1 (MQ=255)
       actgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                            <  1:3313994/70‑1 (MQ=255)
        ctgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                            <  1:2572826/69‑1 (MQ=255)
         tgaacGGGCGATTAAGCGCACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCg                                                                <  1:3869906/64‑1 (MQ=255)
         tgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCg                                                                <  1:969956/64‑1 (MQ=255)
          gaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtattat                                                          >  1:1425562/1‑69 (MQ=255)
          gaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtattat                                                          >  1:2693296/1‑69 (MQ=255)
          gaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTAtct                                                        >  1:1496563/1‑71 (MQ=255)
           aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGta                                                              >  1:3045493/1‑64 (MQ=255)
           aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGta                                                              >  1:3191334/1‑64 (MQ=255)
           aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGta                                                              >  1:4073039/1‑64 (MQ=255)
            acGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGAt                                                                  >  1:1001725/1‑59 (MQ=255)
            acGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGAt                                                                  >  1:3597504/1‑59 (MQ=255)
             cGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:3621662/71‑1 (MQ=255)
             cGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:2868769/71‑1 (MQ=255)
             cGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:3541183/71‑1 (MQ=255)
                  gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcg                                                <  1:3822147/71‑1 (MQ=255)
                  gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcg                                                <  1:2540526/71‑1 (MQ=255)
                  gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcg                                                <  1:2545385/71‑1 (MQ=255)
                   aTTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGt                                                               <  1:2445616/55‑1 (MQ=255)
                   aTTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTAtct                                                        <  1:1316254/62‑1 (MQ=255)
                    ttAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcg                                                >  1:2441244/1‑69 (MQ=255)
                     tAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGta                                                              >  1:476749/1‑54 (MQ=255)
                     tAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAg                                             >  1:1738522/1‑71 (MQ=255)
                      aaGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAg                                             >  1:922392/1‑70 (MQ=255)
                       aGCGAACTCTTAAGTCATATTGCTCTTGAGGTGATGGCAACCTGTGTTCTTATTATCTCTTTAGCGCAgg                                            <  1:3661659/70‑1 (MQ=255)
                       aGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAgg                                            <  1:1333932/70‑1 (MQ=255)
                          gAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgc                                         <  1:3921845/70‑1 (MQ=255)
                            aCTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCt                                      >  1:3873345/1‑71 (MQ=255)
                            aCTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCt                                      <  1:2735564/71‑1 (MQ=255)
                            aCTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCt                                      <  1:612452/71‑1 (MQ=255)
                             ctctTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCtt                                     <  1:2746701/71‑1 (MQ=255)
                                  aaGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATctct                                                      >  1:4052353/1‑49 (MQ=255)
                                  aaGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGa                                 >  1:2227999/1‑70 (MQ=255)
                                   aGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGaaa                               <  1:1002336/71‑1 (MQ=255)
                                   aGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGaaa                               <  1:97479/71‑1 (MQ=255)
                                   aGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGaaa                               <  1:984488/71‑1 (MQ=255)
                                    gTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:1102783/48‑1 (MQ=255)
                                    gTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:1233742/48‑1 (MQ=255)
                                    gTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:1056535/48‑1 (MQ=255)
                                     tCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                            <  1:2994558/40‑1 (MQ=255)
                                     tCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAAt                             <  1:3286335/71‑1 (MQ=255)
                                       atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAAtt                            <  1:2205193/70‑1 (MQ=255)
                                       atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                           <  1:516506/71‑1 (MQ=255)
                                       atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                           <  1:1682994/71‑1 (MQ=255)
                                       atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                           >  1:442109/1‑71 (MQ=255)
                                        tatTCCTCTTGAGGTGATGGCAACCTGTGATCGtat                                                             <  1:3490807/36‑1 (MQ=255)
                                        tatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCtt                                                     <  1:1087112/44‑1 (MQ=255)
                                        tatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcgc                                               <  1:3419192/50‑1 (MQ=255)
                                         atTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGaaaa                              >  1:2507684/1‑66 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCATGAGCGCTTGGGAAAATTACAg                        <  1:410794/71‑1 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        >  1:821555/1‑71 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        <  1:2092630/71‑1 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        <  1:1237651/71‑1 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        <  1:485568/71‑1 (MQ=255)
                                          ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        <  1:181378/71‑1 (MQ=255)
                                            ccTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGa                      >  1:3941283/1‑71 (MQ=255)
                                             ctctTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGAc                     <  1:2316005/71‑1 (MQ=255)
                                             ctctTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGAc                     <  1:787374/71‑1 (MQ=255)
                                                   aGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTAtt               <  1:995193/71‑1 (MQ=255)
                                                   aGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTAtc               <  1:1000048/71‑2 (MQ=255)
                                                     gTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTta             <  1:3567310/71‑1 (MQ=255)
                                                      tGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtat            <  1:3856561/71‑1 (MQ=255)
                                                      tGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtat            <  1:1979156/71‑1 (MQ=255)
                                                        aTGGCAACCTGTTATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtat            >  1:1596019/1‑69 (MQ=255)
                                                        aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTa                 >  1:1004084/1‑64 (MQ=255)
                                                        aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTa                 >  1:1457044/1‑64 (MQ=255)
                                                        aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTa                 >  1:150741/1‑64 (MQ=255)
                                                        aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtata           >  1:629449/1‑70 (MQ=255)
                                                          ggCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        >  1:3798786/1‑55 (MQ=255)
                                                          ggCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                        >  1:2377855/1‑55 (MQ=255)
                                                           gCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATAtctc       >  1:3138382/1‑71 (MQ=255)
                                                           gCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATAtctc       <  1:3076350/71‑1 (MQ=255)
                                                               ccTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTTCAGGACAGTATTTATATCTCGCt    >  1:3368264/1‑70 (MQ=255)
                                                               ccTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtat            >  1:543852/1‑62 (MQ=255)
                                                               ccTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTta             <  1:810922/61‑1 (MQ=255)
                                                               ccTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTta             <  1:393111/61‑1 (MQ=255)
                                                               ccTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGCt    >  1:717356/1‑70 (MQ=255)
                                                                cTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGCTaa  >  1:3840567/1‑71 (MQ=255)
                                                                |                                                                      
GTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGCTAA  >  W3110S.gb/3730276‑3730410

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: