Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 470,764 C→A 6.0% T302N (ACC→AAC)  mdlB → fused predicted multidrug transporter subunits and ATP binding components of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb470,7640CA 6.0% 335.0 / 3.5 99T302N (ACC→AAC) mdlBfused predicted multidrug transporter subunits and ATP binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base C (51/42);  new base A (3/3);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

TTGAAGTGGGCGTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACC  >  W3110S.gb/470695‑470833
                                                                     |                                                                     
ttGAAGTGGGCGTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGAc                                                                       >  1:4479112/1‑70 (MQ=255)
    aGTGGGCGTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACg                                                                   >  1:490249/1‑70 (MQ=255)
     gTGGGCGTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGc                                                                  >  1:2926448/1‑70 (MQ=255)
           gTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGc                                                           >  1:1462264/1‑71 (MQ=255)
           gTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGc                                                           <  1:854514/71‑1 (MQ=255)
            tGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCaa                                                                <  1:4430217/65‑1 (MQ=255)
            tGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGc                                                           >  1:2608202/1‑70 (MQ=255)
             gCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGc                                                                  <  1:5049740/62‑1 (MQ=255)
              cTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGa                                                         <  1:1351211/70‑1 (MQ=255)
              cTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGa                                                         <  1:4976851/70‑1 (MQ=255)
              cTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAt                                                        <  1:411781/71‑1 (MQ=255)
               tGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGa                                                         <  1:4933455/69‑1 (MQ=255)
                gTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAtgc                                                      >  1:2160454/1‑71 (MQ=255)
                 tATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAtgc                                                      >  1:4114573/1‑70 (MQ=255)
                   tGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAtgctgc                                                   <  1:4581115/71‑1 (MQ=255)
                   tGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAtgctg                                                    >  1:1157132/1‑70 (MQ=255)
                     cGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCaa                                                 <  1:5099265/71‑1 (MQ=255)
                     cGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCaa                                                 >  1:1123190/1‑71 (MQ=255)
                     cGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCaa                                                 >  1:2010770/1‑71 (MQ=255)
                     cGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCa                                                  <  1:4454499/70‑1 (MQ=255)
                     cGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGAACACGCAACAGGCGATGCTGCaa                                                 <  1:1054835/71‑1 (MQ=255)
                       tttATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACa                                               <  1:5283673/71‑1 (MQ=255)
                       tttATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGAACACGCAACAg                                                             <  1:2043066/57‑1 (MQ=255)
                       tttATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGAACACGCAACAg                                                             <  1:1103705/57‑1 (MQ=255)
                         tATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAgg                                             >  1:4079815/1‑71 (MQ=255)
                          aTCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAgg                                             <  1:4699013/70‑1 (MQ=255)
                          aTCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGc                                            >  1:1809252/1‑71 (MQ=255)
                           tCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGc                                            >  1:159084/1‑70 (MQ=255)
                            cAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCtg                                          >  1:4081054/1‑71 (MQ=255)
                             aGCTATCTTGGGCTACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATTCTGCAACAGGCtg                                          <  1:2970163/70‑1 (MQ=255)
                             aGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGc                                                                  <  1:2432247/46‑1 (MQ=255)
                               cTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCaa                                                                <  1:4983637/46‑1 (MQ=255)
                               cTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCtgttg                                       >  1:444967/1‑71 (MQ=255)
                                  tCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGc                                                           <  1:1728251/48‑1 (MQ=255)
                                  tCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCa                                                  <  1:3670446/57‑1 (MQ=255)
                                    ttGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAAc                                                               >  1:2870443/1‑42 (MQ=255)
                                    ttGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGa                                                         >  1:3172818/1‑48 (MQ=255)
                                    ttGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCt                                  >  1:3477621/1‑71 (MQ=255)
                                    ttGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCt                                  >  1:2222367/1‑71 (MQ=255)
                                      gggCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGttgttg                                    <  1:449348/67‑1 (MQ=255)
                                      gggCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCt                                  <  1:5223941/69‑1 (MQ=255)
                                      gggCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTgg                                >  1:3512496/1‑71 (MQ=255)
                                       ggCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGt                               <  1:2698200/71‑1 (MQ=255)
                                         cGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGt                               >  1:831828/1‑69 (MQ=255)
                                          gACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCt                                           <  1:966068/56‑1 (MQ=255)
                                          gACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGt                               >  1:784388/1‑68 (MQ=255)
                                          gACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGt                               >  1:778063/1‑68 (MQ=255)
                                           aCTTAACGAACCATTAATCGAACTGAACACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAgc                           >  1:3040113/1‑71 (MQ=255)
                                            cTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGAt                                                        <  1:2059609/41‑1 (MQ=255)
                                             ttAACGAACCATTAATCGAACTGACCACGCAACAGGc                                                           <  1:3236133/37‑1 (MQ=255)
                                             ttAACGAACCATTAATCGAACTGACCACGCAACAGGCGa                                                         <  1:4376521/39‑1 (MQ=255)
                                             ttAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTg                                 >  1:5117903/1‑63 (MQ=255)
                                             ttAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAgcgc                         >  1:370801/1‑71 (MQ=255)
                                              tAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGc                                   >  1:3178794/1‑60 (MQ=255)
                                               aaCGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgt                       >  1:1176945/1‑71 (MQ=255)
                                               aaCGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgt                       >  1:4591758/1‑71 (MQ=255)
                                               aaCGAACCATTAATCGAACTGAACACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgt                       >  1:1540814/1‑71 (MQ=255)
                                               aaCGAACCATTAATCGAACTGAACACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgt                       >  1:1534663/1‑71 (MQ=255)
                                                aCGAACCATTAATCGAACTGACAACGCAACAGGCGATGCTGCAACAgg                                             >  1:3707391/1‑48 (MQ=255)
                                                 cGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCtgttg                                       >  1:4710677/1‑53 (MQ=255)
                                                 cGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgtgt                     <  1:894814/71‑1 (MQ=255)
                                                 cGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgtgt                     >  1:1246829/1‑71 (MQ=255)
                                                  gAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGc                                            <  1:3119734/47‑1 (MQ=255)
                                                  gAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGc                                            <  1:2065642/47‑1 (MQ=255)
                                                  gAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCgt                       <  1:4011838/68‑1 (MQ=255)
                                                   aaCCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGttt                   >  1:379027/1‑71 (MQ=255)
                                                    aCCATTACTCGAACTGACCACGCAACAGGCGAtgc                                                      <  1:928656/35‑1 (MQ=255)
                                                    aCCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAAc                                                >  1:2242102/1‑41 (MQ=255)
                                                    aCCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGttt                   >  1:2296976/1‑70 (MQ=255)
                                                    aCCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTg                  <  1:2756278/71‑1 (MQ=255)
                                                    aCCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTg                  <  1:3943878/71‑1 (MQ=255)
                                                     ccATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGa                 <  1:5558/71‑1 (MQ=255)
                                                      cATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGaa                <  1:1538818/71‑1 (MQ=255)
                                                       aTTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAAc               >  1:2603304/1‑71 (MQ=255)
                                                        ttAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCtgttgt                                      >  1:4926870/1‑47 (MQ=255)
                                                        ttAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACt              <  1:956882/71‑1 (MQ=255)
                                                         tAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGc                                   <  1:242556/49‑1 (MQ=255)
                                                         tAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACt              >  1:1428828/1‑70 (MQ=255)
                                                         tAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACt              >  1:5016913/1‑70 (MQ=255)
                                                          aaTCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTg             <  1:4910033/70‑1 (MQ=255)
                                                          aaTCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGa            >  1:4566753/1‑71 (MQ=255)
                                                          aaTCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGa            >  1:1405608/1‑71 (MQ=255)
                                                          aaTCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGa            >  1:3666645/1‑71 (MQ=255)
                                                           aTCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAAcag             <  1:4283510/69‑3 (MQ=255)
                                                            tCGAACTGACCACGCAACAGGCGATGCTGCAACAgg                                             >  1:4858422/1‑36 (MQ=255)
                                                             cGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGcgcg                        >  1:403507/1‑56 (MQ=255)
                                                             cGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGa                 >  1:3795824/1‑63 (MQ=255)
                                                             cGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGa            >  1:4223125/1‑68 (MQ=255)
                                                              gAACTGACCACGCAACAGGCGATGCTGCAACAGTCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATgga        >  1:4057215/1‑71 (MQ=255)
                                                              gAACTGACCACGCAACAGGCGATGCTGCAACAGGCtgttg                                       >  1:2399396/1‑40 (MQ=255)
                                                              gAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATgga        >  1:1862179/1‑71 (MQ=255)
                                                              gAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATgga        <  1:2135012/71‑1 (MQ=255)
                                                                aCTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATggacg      >  1:2132368/1‑71 (MQ=255)
                                                                 cTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATggacg      >  1:1364953/1‑70 (MQ=255)
                                                                 cTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATggacg      >  1:5015088/1‑70 (MQ=255)
                                                                    aCCACGCATCAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGAcc  <  1:1624240/71‑1 (MQ=255)
                                                                    aCCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATggacggac   >  1:1264471/1‑70 (MQ=255)
                                                                    aCCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGAcc  <  1:370011/71‑1 (MQ=255)
                                                                    aCCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGAcc  <  1:2145382/71‑1 (MQ=255)
                                                                    aCCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGAcc  <  1:4926521/71‑1 (MQ=255)
                                                                     |                                                                     
TTGAAGTGGGCGTGCTGTATGCGTTTATCAGCTATCTTGGGCGACTTAACGAACCATTAATCGAACTGACCACGCAACAGGCGATGCTGCAACAGGCTGTTGTTGCTGGTGAGCGCGTGTTTGAACTGATGGACGGACC  >  W3110S.gb/470695‑470833

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: