Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,585,696 A→G 46.6% D527D (GAT→GAC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,585,6960AG46.6% 15.2 / 115.4 131D527D (GAT→GACmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (19/51);  new base G (44/17);  total (63/68)
Fisher's exact test for biased strand distribution p-value = 2.86e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01

Only 100 of 142 total aligned reads displayed.

CCGGCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTG  >  minE/2585631‑2585764
                                                                 |                                                                    
ccGGCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATc                                                                    >  1:519432/1‑68 (MQ=255)
  ggCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCgg                                                                  <  1:3338369/68‑1 (MQ=255)
  ggCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCgg                                                                  <  1:2435536/68‑1 (MQ=255)
  ggCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCgg                                                                  <  1:2422421/68‑1 (MQ=255)
  ggCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCgg                                                                  <  1:2326128/68‑1 (MQ=255)
  ggCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCgg                                                                  <  1:2226538/68‑1 (MQ=255)
          tCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAgt                                                            <  1:2303448/66‑1 (MQ=255)
           cagcagGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        >  1:1039542/1‑69 (MQ=255)
           cagcagGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        >  1:1757085/1‑69 (MQ=255)
                   aGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAtt                                                >  1:1429714/1‑69 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:1357741/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:3313990/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:800338/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:241305/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:3378498/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:1222446/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:2252855/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:2048115/69‑1 (MQ=255)
                    gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:1854992/69‑1 (MQ=255)
                      ccaccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:263789/67‑1 (MQ=255)
                       caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                             >  1:314046/1‑68 (MQ=255)
                       caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                             >  1:2589061/1‑68 (MQ=255)
                       caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                             >  1:2228773/1‑68 (MQ=255)
                       caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                             >  1:1019832/1‑68 (MQ=255)
                       caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                             >  1:15275/1‑68 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:3175962/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:149501/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:3233319/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2385676/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2345369/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:3207797/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2254895/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:1257877/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2205019/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2165988/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:565885/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2925731/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:167018/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:2881634/1‑69 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGATGAGTTGATTTTTc                                            >  1:2816175/1‑68 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           <  1:1417815/69‑1 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           <  1:1894270/69‑1 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           <  1:695390/69‑1 (MQ=255)
                        accaTGTAGCGACCGGTACAACGCAACAGAGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                           >  1:3287023/1‑69 (MQ=255)
                           aTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:2619586/62‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:556049/39‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:1204608/39‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:2753631/39‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:1558680/39‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:49541/39‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:660097/69‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:1495794/69‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:378777/69‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:2119240/69‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:1678190/69‑1 (MQ=255)
                                         aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                          <  1:670443/69‑1 (MQ=255)
                                          cAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                        <  1:2435034/38‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGtt                                >  1:1415422/1‑59 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGc                           >  1:2683423/1‑64 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGc                       >  1:3207072/1‑68 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:114965/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:771170/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:3412493/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:375722/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:865728/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:503809/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:249753/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:3048431/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:2190047/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:2046242/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:2236030/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:1361021/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:1307588/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      >  1:3038474/1‑69 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:1478004/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:1572454/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:1274401/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:1772197/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:2242731/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                      <  1:967414/69‑1 (MQ=255)
                                             cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGAGTGTTGAAGCGTGCg                      <  1:1194/69‑1 (MQ=255)
                                                aaCAGCGAGCGGGGGCTGTCGGTATAGTGCTGAGTTGAttt                                               <  1:1458253/41‑1 (MQ=37)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:94436/69‑1 (MQ=255)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:2019577/69‑1 (MQ=255)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:3159669/69‑1 (MQ=255)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:2393244/69‑1 (MQ=255)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:1173501/69‑1 (MQ=255)
                                                  cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                 <  1:2480772/69‑1 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:412173/1‑37 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:1783466/1‑37 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:3311346/1‑37 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:851045/1‑37 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:945777/1‑37 (MQ=255)
                                                        gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                           >  1:2770416/1‑37 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:3425835/69‑1 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:2002765/69‑1 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:2404141/69‑1 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:71769/69‑1 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:3246127/69‑1 (MQ=255)
                                                                 aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg  <  1:2536573/69‑1 (MQ=255)
                                                                 |                                                                    
CCGGCGCAAATCAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTG  >  minE/2585631‑2585764

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: