Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,720,043 A→G 11.5% T479A (ACG→GCG)  yfiQ → fused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,720,0430AG11.5% 154.1 / 12.2 61T479A (ACG→GCG) yfiQfused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (22/32);  new base G (4/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

GCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACC  >  W3110S.gb/2719973‑2720108
                                                                      |                                                                 
gCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTa                                                                   >  1:1684385/1‑71 (MQ=255)
 cccAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAc                                                                  <  1:1910750/71‑1 (MQ=255)
 cccAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAc                                                                  <  1:3731281/71‑1 (MQ=255)
      cAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGc                                                             >  1:2043424/1‑71 (MQ=255)
        aTCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTc                                                           <  1:3805407/71‑1 (MQ=255)
          cTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCCCGa                                                         >  1:2274165/1‑71 (MQ=255)
           tGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTc                                                           <  1:744347/68‑1 (MQ=255)
           tGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCACGa                                                         >  1:422337/1‑70 (MQ=255)
            gACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGa                                                         <  1:2301279/69‑1 (MQ=255)
               ttCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTc                                                           <  1:2673944/64‑1 (MQ=255)
                 ccAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCg                                                          >  1:4648525/1‑63 (MQ=255)
                  cAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGa                                                         >  1:1008324/1‑63 (MQ=255)
                  cAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGAt                                                        >  1:3810550/1‑64 (MQ=255)
                  cAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATg                                                 <  1:2973801/71‑1 (MQ=255)
                     tACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTAc                                                                  <  1:3688472/51‑1 (MQ=255)
                     tACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAg                                              >  1:568335/1‑71 (MQ=255)
                      aCCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGt                                             >  1:4240193/1‑71 (MQ=255)
                      aCCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGt                                             >  1:4354997/1‑71 (MQ=255)
                       ccGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGAt                                                        >  1:38924/1‑59 (MQ=255)
                       ccGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGtt                                            >  1:4490216/1‑71 (MQ=255)
                         gCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCa                                          <  1:3907805/71‑1 (MQ=255)
                         gCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCa                                          <  1:1349805/71‑1 (MQ=255)
                          cAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTGCGTCGCTCGATACCCATGAAGTTCAg                                         >  1:536098/1‑71 (MQ=255)
                                    tcttctGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATAc                                                      <  1:4642788/48‑1 (MQ=255)
                                     cttctGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCGGCCCATCCTGc                               >  1:3844814/1‑70 (MQ=255)
                                       tctGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAg                                         <  1:1553098/58‑1 (MQ=255)
                                        ctGTTGCAACAGGCGATTGCCGAAGGGGCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGc                               <  1:1807382/67‑1 (MQ=255)
                                        ctGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGc                           >  1:2552475/1‑71 (MQ=255)
                                          gTTGCAACAGGCGATTGCCGAAGGGGCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCg                          >  1:1284406/1‑70 (MQ=255)
                                           ttGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATg                                                 <  1:2176277/46‑1 (MQ=255)
                                           ttGCAACAGGCGATTCCCGAAGGGGCTACGTCGCTCGATACCCATg                                                 <  1:3295797/46‑1 (MQ=255)
                                             gCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGaa                                               >  1:311571/1‑46 (MQ=255)
                                             gCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATcc                                  >  1:1363767/1‑59 (MQ=255)
                                              cAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATg                      <  1:2268093/70‑1 (MQ=255)
                                              cAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATg                      <  1:3464279/70‑1 (MQ=255)
                                               aaCAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGc                    <  1:3619103/71‑1 (MQ=255)
                                                aCAGGCGATTGCCGAAGGGGCTACATCGCTCGATACCCATGAAGtt                                            <  1:2518618/46‑1 (MQ=255)
                                                 cAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCg                          <  1:4587959/63‑1 (MQ=255)
                                                  aGGCGATTGCCGAAGGGGCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCAt                  >  1:3163640/1‑70 (MQ=255)
                                                  aGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCa                                          <  1:4204678/46‑1 (MQ=255)
                                                  aGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATg                 <  1:3390237/71‑1 (MQ=255)
                                                  aGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATg                 <  1:2129389/71‑1 (MQ=255)
                                                   ggCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATg                                                 <  1:1097476/38‑1 (MQ=255)
                                                   ggCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAGGCGTATGGCATg                 >  1:3903197/1‑70 (MQ=255)
                                                   ggCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGa                <  1:622255/71‑1 (MQ=255)
                                                      gATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTa                        <  1:3198758/60‑1 (MQ=255)
                                                      gATTGCCAAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAaca             <  1:2807802/71‑1 (MQ=255)
                                                       aTTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCATCCCATCCTGCAAGCGTATGGCATGAacac            <  1:1358034/71‑1 (MQ=255)
                                                        ttGCCGAAGGGGCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACg           >  1:1341782/1‑71 (MQ=255)
                                                        ttGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCg                          <  1:2695554/56‑1 (MQ=255)
                                                        ttGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACg           >  1:1048023/1‑71 (MQ=255)
                                                         tGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTa                        >  1:481275/1‑57 (MQ=255)
                                                              aaGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATc                                   >  1:1772812/1‑41 (MQ=255)
                                                              aaGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCt     <  1:1660949/71‑1 (MQ=255)
                                                                 gggCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATg                 <  1:3473173/56‑1 (MQ=255)
                                                                 gggCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTAcc  >  1:3431056/1‑71 (MQ=255)
                                                                 gggCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTAcc  >  1:3425497/1‑71 (MQ=255)
                                                                   gCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATgg                     <  1:1588925/50‑1 (MQ=255)
                                                                   gCTGCGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATgg                     <  1:784820/50‑1 (MQ=255)
                                                                   gCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGt                         <  1:3918522/46‑1 (MQ=255)
                                                                   gCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCa                   >  1:3080124/1‑52 (MQ=255)
                                                                     tACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAac              <  1:3007492/55‑1 (MQ=255)
                                                                      |                                                                 
GCCCAGCAATCTGACTTCCAATACCGCAGAAGCGCATCTTCTGTTGCAACAGGCGATTGCCGAAGGGGCTACGTCGCTCGATACCCATGAAGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACC  >  W3110S.gb/2719973‑2720108

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: