Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,730,340 C→T 12.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,3400CT12.9% 169.8 / 14.1 62T75T (ACC→ACTbglGtranscriptional antiterminator of the bgl operon
Reads supporting (aligned to +/- strand):  ref base C (21/33);  new base T (6/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 6.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GAGCAGTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGC  >  W3110S.gb/3730271‑3730407
                                                                     |                                                                   
gAGCAGTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCt                                                                    >  1:819869/1‑71 (MQ=255)
        aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCg                                                             >  1:3245951/1‑70 (MQ=255)
        aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             >  1:606752/1‑70 (MQ=255)
        aTGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             >  1:4573120/1‑70 (MQ=255)
          gaactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtat                                                          >  1:912507/1‑71 (MQ=255)
          gaactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtat                                                          >  1:587250/1‑71 (MQ=255)
           aactgaacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatt                                                         <  1:30448/71‑1 (MQ=255)
                aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatta                                                        >  1:126559/1‑67 (MQ=255)
                aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTAtct                                                     >  1:1649558/1‑70 (MQ=255)
                aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTAtct                                                     >  1:1689451/1‑70 (MQ=255)
                aacGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTAtct                                                     >  1:1955049/1‑70 (MQ=255)
                    ggCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTATCTCTTTa                                                <  1:2808628/71‑1 (MQ=255)
                       gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTTATGGCAACCTGTGATCGTATTATCTCttt                                                 <  1:3533774/67‑1 (MQ=255)
                       gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGta                                                           <  1:1718205/57‑1 (MQ=255)
                       gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGta                                                           <  1:1405156/57‑1 (MQ=255)
                       gATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgc                                              <  1:2990398/70‑1 (MQ=255)
                        aTTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGtatta                                                        <  1:2002674/59‑1 (MQ=255)
                         ttAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTATCTCTTTAgcg                                             >  1:1789417/1‑69 (MQ=255)
                         ttAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             <  1:37837/53‑1 (MQ=255)
                         ttAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             <  1:2682076/53‑1 (MQ=255)
                          tAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGta                                                           >  1:3030619/1‑54 (MQ=255)
                            aGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAgg                                         <  1:1160887/70‑1 (MQ=255)
                            aGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGa                                        <  1:4290961/71‑1 (MQ=255)
                              cGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgc                                      <  1:1896972/71‑1 (MQ=255)
                               gAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgcg                                     <  1:4504729/71‑1 (MQ=255)
                                aaCTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgcgc                                    >  1:4635955/1‑71 (MQ=255)
                                 aCTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCt                                   <  1:3089046/71‑1 (MQ=255)
                                   tctTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGa                                                                <  1:1770382/40‑1 (MQ=255)
                                   tctTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGt                                                            <  1:2767226/44‑1 (MQ=255)
                                   tctTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCa                                           <  1:3536429/61‑1 (MQ=255)
                                     ttAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             <  1:4429698/41‑1 (MQ=255)
                                     ttAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCg                                                             <  1:2023458/41‑1 (MQ=255)
                                        aGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAg                                          >  1:324512/1‑57 (MQ=255)
                                        aGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAg                                          >  1:2529611/1‑57 (MQ=255)
                                         gTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgc                                      <  1:4619278/60‑1 (MQ=255)
                                          tCATATTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTAtc                                                      <  1:4486011/43‑1 (MQ=255)
                                          tCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAAt                          <  1:447604/71‑1 (MQ=255)
                                            atatTCCTCTTGAGGTGATGTCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        >  1:1827317/1‑71 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTAtctc                                                    >  1:1489145/1‑43 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTATCTCTTTAGCGCa                                           >  1:307050/1‑52 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACTTGTGATCGTATTATCTCTTTAGCGCa                                           >  1:2477390/1‑52 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        <  1:284717/71‑1 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        >  1:3498366/1‑71 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        <  1:50590/71‑1 (MQ=255)
                                            atatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        <  1:2605880/71‑1 (MQ=255)
                                             tatTCCTCTTGAGGTGATGGCAACCTGTGATCTTATTATCTCTTTAGCGCAGGAGCGCTTGGGaa                             <  1:295484/65‑1 (MQ=255)
                                             tatTCCTCTTGAGGTGATGGCAACCTGTGATCGta                                                           <  1:4307132/35‑1 (MQ=255)
                                             tatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAgcgc                                            <  1:4331156/50‑1 (MQ=255)
                                             tatTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAgcg                                     <  1:620182/57‑1 (MQ=255)
                                               ttCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAg                     <  1:381914/71‑1 (MQ=255)
                                                 ccTCTTGAGGTGATGGCAACCTGTGATCGTATTAtc                                                      >  1:2609298/1‑36 (MQ=255)
                                                    ctTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAg                >  1:2595446/1‑71 (MQ=255)
                                                       gAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCATGAGCGCTTGGGAAAATTACAGGACAGTa              <  1:1006190/70‑1 (MQ=255)
                                                          gTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTa                        <  1:4480762/57‑1 (MQ=255)
                                                           tGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTtat         <  1:1231718/71‑1 (MQ=255)
                                                             aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAgg                                         >  1:1169803/1‑37 (MQ=255)
                                                             aTGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGaa                             <  1:4047101/49‑1 (MQ=255)
                                                              tGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTAtt            <  1:3410516/65‑1 (MQ=255)
                                                                gCAACTTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTg                                 >  1:3980560/1‑42 (MQ=255)
                                                                gCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATAtctc    >  1:4483004/1‑71 (MQ=255)
                                                                  aaCCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGc  >  1:878939/1‑71 (MQ=255)
                                                                  aaCCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGc  >  1:2198386/1‑71 (MQ=255)
                                                                     |                                                                   
GAGCAGTCATGAACTGAACGGGCGATTAAGCGAACTCTTAAGTCATATTCCTCTTGAGGTGATGGCAACCTGTGATCGTATTATCTCTTTAGCGCAGGAGCGCTTGGGAAAATTACAGGACAGTATTTATATCTCGC  >  W3110S.gb/3730271‑3730407

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: