Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,586,050 G→C 13.7% A409A (GCC→GCG mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,586,0500GC13.7% 285.1 / 27.3 102A409A (GCC→GCGmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base G (34/54);  new base C (8/6);  total (42/60)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 105 total aligned reads displayed.

ATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCG  >  minE/2585981‑2586119
                                                                     |                                                                     
aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                                                       <  1:3945938/70‑1 (MQ=255)
aTGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATg                                                                       <  1:1426324/70‑1 (MQ=255)
 tGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATgg                                                                      <  1:274622/70‑1 (MQ=255)
        cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                              <  1:2360006/71‑1 (MQ=255)
        cttccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATcc                                                              <  1:576957/71‑1 (MQ=255)
          tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                            >  1:625026/1‑71 (MQ=255)
          tccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                            >  1:1285240/1‑71 (MQ=255)
           ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTTATccac                                                            <  1:511122/70‑1 (MQ=255)
           ccttcGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccac                                                            <  1:2412808/70‑1 (MQ=255)
              ttggcggGAGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAg                                                         <  1:2026811/68‑1 (MQ=12)
               cggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCGCCAGt                                                        >  1:1686169/1‑70 (MQ=255)
                ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                      >  1:3302608/1‑71 (MQ=255)
                ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                       <  1:1024032/70‑1 (MQ=255)
                ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                       <  1:3740850/70‑1 (MQ=255)
                ggcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTa                                                       <  1:1783259/70‑1 (MQ=255)
                 gcggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAc                                                     >  1:3736467/1‑71 (MQ=255)
                   ggTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATccacca                                                          <  1:1971060/64‑1 (MQ=255)
                    ggaGCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGccc                                                   >  1:1364286/3‑70 (MQ=255)
                     tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                 <  1:3595756/71‑1 (MQ=255)
                     tAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                 <  1:290303/71‑1 (MQ=255)
                      aGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                >  1:3489502/1‑71 (MQ=255)
                       gCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                               <  1:1975963/71‑1 (MQ=255)
                       gCTTATCCTCCATCATCACGCGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGAtg                                               >  1:3295347/1‑70 (MQ=255)
                         ttATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                             >  1:2878324/1‑71 (MQ=255)
                          tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                             >  1:3234193/1‑70 (MQ=255)
                          tATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGc                                             >  1:1274685/1‑70 (MQ=255)
                           aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                   <  1:1540280/63‑1 (MQ=255)
                           aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCg                                              <  1:2119080/68‑1 (MQ=255)
                           aTCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGcc                                            <  1:3554745/70‑1 (MQ=255)
                            tCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCa                                           <  1:2884450/70‑1 (MQ=255)
                              ttCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATACACCAGTAACCCTATCGCCAGCa                                        >  1:3488866/1‑71 (MQ=255)
                               tCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                       >  1:521560/1‑71 (MQ=255)
                                ccGCAATGACACTCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                        <  1:1871226/69‑1 (MQ=255)
                                ccGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCa                                        <  1:3080823/69‑1 (MQ=255)
                                 cGCAATTACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                      <  1:2019901/70‑1 (MQ=255)
                                 cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAcc                                      <  1:59734/70‑1 (MQ=255)
                                 cGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCa                                     <  1:3776617/71‑1 (MQ=255)
                                  gCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTAt                                                >  1:179953/1‑59 (MQ=255)
                                    aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTaa                                                      >  1:2727636/1‑51 (MQ=255)
                                    aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTa                                                 <  1:1279990/56‑1 (MQ=255)
                                    aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                  >  1:3364156/1‑71 (MQ=255)
                                    aaTGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATcc                                  >  1:253318/1‑71 (MQ=255)
                                      tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                               >  1:1467828/1‑56 (MQ=255)
                                      tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATc                                               >  1:3748009/1‑56 (MQ=255)
                                      tGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATccc                                 >  1:416224/1‑70 (MQ=255)
                                       gACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCAc                                       <  1:786780/63‑1 (MQ=255)
                                        acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATGCCGaa                              >  1:2101319/4‑71 (MQ=2)
                                        acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                              >  1:1284789/4‑71 (MQ=2)
                                        acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                              >  1:1833409/4‑71 (MQ=2)
                                        acgcGCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGaa                              >  1:3790593/4‑71 (MQ=2)
                                            gCTCGACGTTCTCCACCACCGCGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCg                                >  1:3353477/1‑65 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                        <  1:3711481/57‑1 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                        <  1:1200652/57‑1 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCa                                        <  1:148099/57‑1 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGt                                                        <  1:1123023/41‑1 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                            >  1:1872219/1‑69 (MQ=255)
                                            gCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACa                            >  1:2023796/1‑69 (MQ=255)
                                             cTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAAccc                                                   <  1:2353192/45‑1 (MQ=255)
                                              tCGACGTTCTCCACCACCACGATGGCGTCATccacc                                                           <  1:2801586/36‑1 (MQ=40)
                                              tCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                        >  1:187113/1‑71 (MQ=255)
                                               cGACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGt                        >  1:3416428/1‑70 (MQ=2)
                                               cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                             <  1:928251/65‑1 (MQ=255)
                                               cGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                             <  1:413724/65‑1 (MQ=255)
                                                gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                      <  1:85629/71‑3 (MQ=2)
                                                gACGTTCTCCACCACCACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGtta                      <  1:2084177/71‑3 (MQ=2)
                                                gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGt                        <  1:2651753/69‑1 (MQ=255)
                                                gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                       <  1:2405175/70‑1 (MQ=255)
                                                gACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCa                      <  1:2225889/71‑1 (MQ=255)
                                                 aCGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAg                         >  1:2130414/1‑67 (MQ=255)
                                                  cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAAc                             >  1:2417964/1‑62 (MQ=255)
                                                  cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                    >  1:417058/1‑71 (MQ=255)
                                                  cGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                    <  1:3045170/71‑1 (MQ=255)
                                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                   <  1:2033255/71‑1 (MQ=255)
                                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                   <  1:2210775/71‑1 (MQ=255)
                                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                   >  1:1140588/1‑71 (MQ=255)
                                                   gTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACg                   <  1:2441673/71‑1 (MQ=255)
                                                    ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAAc                    <  1:1396246/69‑1 (MQ=255)
                                                    ttCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                  >  1:1738287/1‑71 (MQ=255)
                                                     tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAg                                          <  1:3877786/46‑1 (MQ=255)
                                                     tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtg                 >  1:142341/1‑71 (MQ=255)
                                                     tctcCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                  <  1:526951/70‑1 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                           <  1:836523/56‑1 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACAt                           <  1:1956328/56‑1 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt            <  1:916688/71‑1 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt            >  1:3757272/1‑71 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt            >  1:3481322/1‑71 (MQ=255)
                                                          accaccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGAt            >  1:2685255/1‑71 (MQ=255)
                                                            caccacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgt                  <  1:1421996/63‑1 (MQ=255)
                                                               cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa       <  1:3618446/71‑1 (MQ=255)
                                                               cacGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGa       <  1:2792971/71‑1 (MQ=255)
                                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa      <  1:40117/71‑1 (MQ=255)
                                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa      <  1:3900925/71‑1 (MQ=255)
                                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa      <  1:1450063/71‑1 (MQ=255)
                                                                acGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaa      <  1:2680586/71‑1 (MQ=255)
                                                                 cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACgtgt                >  1:1537788/1‑60 (MQ=255)
                                                                 cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa     <  1:484202/71‑1 (MQ=255)
                                                                 cGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGaaa     <  1:2915695/71‑1 (MQ=255)
                                                                  gATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTc                       >  1:124831/1‑52 (MQ=255)
                                                                   aTGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAacc   <  1:3423260/71‑1 (MQ=255)
                                                                    tGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAaccg  <  1:3864657/71‑1 (MQ=255)
                                                                     |                                                                     
ATGCGTCGCTTCCTTCGGCGGTAGCTTATCTTCCGCAATGACACGCTCGACGTTCTCCACCACCACGATGGCGTCATCCACCAGTAACCCTATCGCCAGCACCATCCCGAACATAGTCAACGTGTTGATGGTGAAACCG  >  minE/2585981‑2586119

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: