Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,585,708 T→C 50.1% Q523Q (CAA→CAG mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,585,7080TC50.1% ‑1.4 / 111.0 121Q523Q (CAA→CAGmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base T (15/45);  new base C (43/17);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 5.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

Only 100 of 159 total aligned reads displayed.

CAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCGG  >  minE/2585642‑2585776
                                                                  |                                                                    
cagcagGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    >  1:1039542/1‑69 (MQ=255)
cagcagGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    >  1:2038962/1‑69 (MQ=255)
cagcagGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    >  1:1757085/1‑69 (MQ=255)
        aGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAtt                                                            >  1:1429714/1‑69 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:241305/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:2048115/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:1494165/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:3313990/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:1321484/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:332500/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:1222446/69‑1 (MQ=255)
         gACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:1202292/69‑1 (MQ=255)
           ccaccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:263789/67‑1 (MQ=255)
            caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                                         >  1:2228773/1‑68 (MQ=255)
            caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                                         >  1:2589061/1‑68 (MQ=255)
            caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                                         >  1:314046/1‑68 (MQ=255)
            caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                                         >  1:1019832/1‑68 (MQ=255)
            caccaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGAttttt                                                         >  1:15275/1‑68 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTc                                                        >  1:170122/1‑68 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:1553256/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2385676/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2881634/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2120265/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2925731/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:3233319/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:3175962/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:565885/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2254895/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:2420591/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:3207797/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:149501/1‑69 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTGTCGGTATAGTGATGAGTTGATTTTTc                                                        >  1:2816175/1‑68 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       <  1:1894270/69‑1 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       <  1:1417815/69‑1 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       <  1:695390/69‑1 (MQ=255)
             accaTGTAGCGACCGGTACAACGCAACAGAGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCa                                                       >  1:3287023/1‑69 (MQ=255)
                   tAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGc                                                 >  1:1947079/1‑69 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATGGGTATAGTgttg                                                                    <  1:1788130/39‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    <  1:49541/39‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    <  1:2753631/39‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTgttg                                                                    <  1:556049/39‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                                      <  1:670443/69‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                                      <  1:1678190/69‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                                      <  1:1642614/69‑1 (MQ=255)
                              aCAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCg                                      <  1:378777/69‑1 (MQ=255)
                                 aCGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTTAttt                                                           <  1:2509028/45‑1 (MQ=37)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGc                                   >  1:3207072/1‑68 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:1361021/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:510139/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:2265525/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:3136230/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:3038474/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:2236030/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:865728/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:2717638/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:3267713/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:771170/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:2046242/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:2190047/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTGTCGGTATAGTGCTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  >  1:1164187/1‑69 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:344223/69‑1 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:1572454/69‑1 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:2242731/69‑1 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:945016/69‑1 (MQ=255)
                                  cGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:1274401/69‑1 (MQ=255)
                                   gCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCg                                  <  1:2619686/68‑1 (MQ=255)
                                     aaCAGCGAGCGGGGGCTGTCGGTATAGTGCTGAGTTGAttt                                                           <  1:1458253/41‑1 (MQ=37)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:1408516/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:3404442/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:1270360/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:1173501/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:595438/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:2480772/69‑1 (MQ=255)
                                       cAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAg                             <  1:2397561/69‑1 (MQ=255)
                                             gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       >  1:412173/1‑37 (MQ=255)
                                             gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       >  1:851045/1‑37 (MQ=255)
                                             gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       >  1:945777/1‑37 (MQ=255)
                                             gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       >  1:2383697/1‑37 (MQ=255)
                                             gcgGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCa                                                       >  1:1783466/1‑37 (MQ=255)
                                                      aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg              <  1:3425835/69‑1 (MQ=255)
                                                      aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg              <  1:2447277/69‑1 (MQ=255)
                                                      aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg              <  1:2193818/69‑1 (MQ=255)
                                                      aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg              <  1:2002765/69‑1 (MQ=255)
                                                      aTCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTtgtg              <  1:71769/69‑1 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:317746/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:2587682/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:78916/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:1130390/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:1136693/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:655049/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:2148760/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:3042364/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAgg                    >  1:466564/1‑54 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTc     >  1:2772789/1‑69 (MQ=255)
                                                               gtgtTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTc     >  1:3159953/1‑69 (MQ=255)
                                                                  ttgagttgaTTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCgg  >  1:2778449/1‑69 (MQ=255)
                                                                  ttgagttgaTTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCgg  >  1:123901/1‑69 (MQ=255)
                                                                  ttgagttgaTTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCgg  >  1:2693847/1‑69 (MQ=255)
                                                                  ttgagttgaTTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCgg  >  1:873114/1‑69 (MQ=255)
                                                                  ttgagttgaTTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCgg  >  1:985415/1‑69 (MQ=255)
                                                                  |                                                                    
CAGCAGGTAGACCACCATGTAGCGACCGGTACAACGCAACAGCGAGCGGGTGCTATCGGTATAGTGTTGAGTTGATTTTTCAAACAGCGTGTTGAAGCGTGCGAACAGGGCGTTAGGTTTGTGACCGCCTTCCGG  >  minE/2585642‑2585776

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: