Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,482,946 A→G 5.0% intergenic (‑83/+185) mdtL ← / ← trmE multidrug efflux system protein/GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,482,9460AG 5.0% 440.3 / 11.0 158intergenic (‑83/+185)mdtL/trmEmultidrug efflux system protein/GTPase
Reads supporting (aligned to +/- strand):  ref base A (71/79);  new base G (3/5);  total (74/84)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 172 total aligned reads displayed.

CTGATAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCTTTTTTTA  >  minE/2482877‑2483013
                                                                     |                                                                   
cTGATAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     <  1:1212865/70‑1 (MQ=255)
  gATAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTc                                                                  >  1:2169848/1‑71 (MQ=255)
   aTAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     <  1:2328222/67‑1 (MQ=255)
   aTAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     <  1:1593620/67‑1 (MQ=255)
   aTAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCt                                                                 >  1:1970944/1‑71 (MQ=255)
     aGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTg                                                                >  1:4327219/1‑70 (MQ=255)
      gTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTg                                                                >  1:1798469/1‑69 (MQ=255)
       tAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGaaa                                                             <  1:495838/71‑1 (MQ=255)
        aaGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     >  1:574248/1‑62 (MQ=255)
        aaGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAAt                                                            >  1:2975641/1‑71 (MQ=255)
        aaGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAAt                                                            >  1:1759198/1‑71 (MQ=255)
         aGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAAtt                                                           >  1:3940606/1‑71 (MQ=255)
          ggAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTg                                                          <  1:2524554/71‑1 (MQ=255)
          ggAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTg                                                          <  1:3499846/71‑1 (MQ=255)
              tGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGa                                                               <  1:1485655/62‑1 (MQ=255)
              tGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGttt                                                       >  1:2179466/1‑70 (MQ=255)
              tGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGttt                                                       >  1:2870892/1‑70 (MQ=255)
                cACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTaaat                                                             >  1:3448693/1‑61 (MQ=255)
                cACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTaaa                                                    >  1:4523437/1‑71 (MQ=255)
                cACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTaa                                                     <  1:741733/70‑1 (MQ=255)
                    tATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     >  1:4624953/1‑50 (MQ=255)
                    tATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAt                                                   >  1:1402404/1‑68 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAtgtg                                                >  1:1459167/1‑70 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAtgtg                                                >  1:412197/1‑70 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAtgtg                                                >  1:4218191/1‑70 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAtgtg                                                >  1:4375421/1‑70 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAtgtg                                                <  1:364835/70‑1 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAt                                                   >  1:2269181/1‑67 (MQ=255)
                     aTCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAt                                                   >  1:713262/1‑67 (MQ=255)
                      tCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTGTTCTGAAATTGTTTAAATGTgaa                                              >  1:1519062/1‑71 (MQ=255)
                      tCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAAt                                                   >  1:2856775/1‑66 (MQ=255)
                      tCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTga                                               <  1:2610937/70‑1 (MQ=255)
                      tCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTga                                               <  1:769786/70‑1 (MQ=255)
                        cTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTg                                                                >  1:4025987/1‑51 (MQ=255)
                        cTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTGAATGTgaat                                             <  1:826778/70‑1 (MQ=255)
                         tCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTaa                                                     <  1:4669770/61‑1 (MQ=255)
                          cGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     >  1:1421764/1‑44 (MQ=255)
                          cGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTa                                                                     >  1:1074701/1‑44 (MQ=255)
                          cGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCt                                                                 >  1:2902471/1‑48 (MQ=255)
                           gCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTgaatcgaa                                         >  1:263555/1‑71 (MQ=255)
                           gCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTgaatcg                                           <  1:3074661/69‑1 (MQ=255)
                            cAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGa                                                               <  1:3369048/48‑1 (MQ=255)
                            cAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTgaatcgaat                                        <  1:678669/71‑1 (MQ=255)
                              aGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAAtt                                                           >  1:3798242/1‑50 (MQ=255)
                              aGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTga                                               <  1:4635564/62‑1 (MQ=255)
                              aGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATca                                      <  1:2242043/71‑1 (MQ=255)
                               gACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATcac                                     <  1:1093199/71‑1 (MQ=255)
                                acacGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTATATGTGAATCGAATcac                                     >  1:2596702/1‑70 (MQ=255)
                                acacGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATcaca                                    <  1:2486703/71‑1 (MQ=255)
                                  acGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGaaa                                                             >  1:1784541/1‑44 (MQ=255)
                                  acGAGGCTACATTATTACACTCAGAGTTTTTGTATATTCTGAAATTGTTTAAAt                                                   <  1:1069490/54‑1 (MQ=255)
                                   cGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATc                                 <  1:2080440/71‑1 (MQ=255)
                                   cGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATc                                 <  1:3824381/71‑1 (MQ=255)
                                      ggCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGtt                              >  1:2655443/1‑71 (MQ=255)
                                      ggCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGtt                              <  1:191391/71‑1 (MQ=255)
                                       gCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATca                                      <  1:3426979/62‑1 (MQ=255)
                                       gCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTc                             <  1:1011145/71‑1 (MQ=255)
                                        cTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTgaatcgaatga                                      <  1:873842/61‑3 (MQ=255)
                                        cTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCg                            >  1:2612393/1‑71 (MQ=255)
                                         tACATTATTACACTCAGAGTTTTTGTCTATTCTGAAAtt                                                           >  1:3723515/1‑39 (MQ=255)
                                         tACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCgg                           >  1:785882/1‑71 (MQ=255)
                                          agattattACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTa                                                      <  1:1005800/41‑1 (MQ=255)
                                           cATTATTACACTCAGAGTTTTTGTCTATTCTGaaa                                                             <  1:3474198/35‑1 (MQ=255)
                                           cATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTaaa                                                    >  1:1524632/1‑44 (MQ=255)
                                           cATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCggg                          >  1:4470211/1‑70 (MQ=255)
                                           cATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGATTAAATGTGAATCGAATCACAATCGTTCggg                          >  1:2973548/1‑70 (MQ=255)
                                            attattACCCTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCTTTCgg                           >  1:644141/1‑68 (MQ=255)
                                            attattACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCggg                          >  1:1630032/1‑69 (MQ=255)
                                             ttattaCACTCAGAGTTTTTGTCTATTCTGAAATTGttt                                                       >  1:1118040/1‑39 (MQ=255)
                                             ttattaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGt                               <  1:2572134/63‑1 (MQ=255)
                                              tattaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCggg                          >  1:2719984/1‑67 (MQ=255)
                                              tattaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAg                       >  1:3972142/1‑70 (MQ=255)
                                               attaCACTCAGAGTTTTTGTCTGTTCTGAAATTGTTTaaa                                                    >  1:2785136/1‑40 (MQ=255)
                                               attaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGGGc                      <  1:2011230/70‑1 (MQ=255)
                                               attaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGc                      >  1:2340480/1‑70 (MQ=255)
                                                ttaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTa                                                      <  1:1257205/37‑1 (MQ=255)
                                                ttaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTa                                                      <  1:3772862/37‑1 (MQ=255)
                                                ttaCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCaa                    <  1:609735/71‑1 (MQ=255)
                                                 taCACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCaa                    >  1:1182048/1‑70 (MQ=255)
                                                    acTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATAtt                <  1:3249286/71‑1 (MQ=255)
                                                    acTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATAtt                <  1:2359358/71‑1 (MQ=255)
                                                     cTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTa               <  1:3246143/71‑1 (MQ=255)
                                                        agagTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCaa                    >  1:734843/1‑63 (MQ=255)
                                                        agagTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAAt            >  1:3608370/1‑71 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTgaatcg                                           <  1:447933/36‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATca                                      <  1:2116635/41‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGtt                              <  1:1631959/49‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGt                               <  1:1200473/48‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTa               <  1:1840088/64‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTa               <  1:2187551/64‑1 (MQ=255)
                                                            tttttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCtt        <  1:4197101/71‑1 (MQ=255)
                                                              tttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATcaca                                    >  1:2923164/1‑41 (MQ=255)
                                                              tttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTa               <  1:155979/62‑1 (MQ=255)
                                                              tttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCtttt      >  1:4076386/1‑71 (MQ=255)
                                                               ttGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACaa                                   <  1:2102781/41‑1 (MQ=255)
                                                                  tCTGTTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCTTTTTTTa  <  1:2693167/71‑1 (MQ=255)
                                                                  tCTGTTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCTTTTTTTa  <  1:1034149/71‑1 (MQ=255)
                                                                  tCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCTTTTTTTa  <  1:1360686/71‑1 (MQ=255)
                                                                   cTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCggga                         >  1:1287073/1‑46 (MQ=255)
                                                                    tATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCg                                <  1:1975782/39‑1 (MQ=255)
                                                                     |                                                                   
CTGATAGTAAGGAATGCACTTATCCTCGCAAGACACGAGGCTACATTATTACACTCAGAGTTTTTGTCTATTCTGAAATTGTTTAAATGTGAATCGAATCACAATCGTTCGGGGAGCAATATTAAATGCTTTTTTTA  >  minE/2482877‑2483013

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: