Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,482,204 A→G 52.1% T220T (ACT→ACC mdtL ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,482,2040AG52.1% 7.0 / 85.7 98T220T (ACT→ACCmdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base G (7/44);  total (28/70)
Fisher's exact test for biased strand distribution p-value = 8.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

Only 100 of 101 total aligned reads displayed.

GTGGCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGAAAAAACGGTTAAGCAGCGA  >  minE/2482139‑2482272
                                                                 |                                                                    
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:3190999/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:600268/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:1960309/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:2462341/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:1915383/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:1863139/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:748463/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:1690509/69‑1 (MQ=255)
gtggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:3228511/69‑1 (MQ=255)
 tggCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTg                                                                   <  1:862877/68‑1 (MQ=255)
   gCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:1313567/69‑1 (MQ=255)
   gCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:1338316/69‑1 (MQ=255)
   gCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:2437200/69‑1 (MQ=255)
   gCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:473788/69‑1 (MQ=255)
   gCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:3154711/69‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:3228614/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:2143932/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:2114962/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:949479/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:3086857/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAgg                                                                <  1:3388549/62‑1 (MQ=255)
          cACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCAcc                                                          >  1:587092/1‑68 (MQ=255)
            ccGCGCTCATAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1721741/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:263152/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:265459/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:3094456/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:2740871/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:2548953/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:2806124/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:486938/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:2440114/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:291700/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:3382354/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:2378098/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:643789/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:55474/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1380101/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1128461/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:857700/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1360103/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1378847/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:598624/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1700754/68‑1 (MQ=255)
            ccGCGCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:687953/68‑1 (MQ=255)
             cgcgCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1461750/67‑1 (MQ=255)
             cgcgCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:743279/67‑1 (MQ=255)
             cgcgCTCAAAACCCATGATTTCCATCAGCAATACTGGCGACGTGGTGACGAAGGTGAGGATCACCGa                                                        <  1:2031473/67‑1 (MQ=255)
                             gATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGAt                                      >  1:800277/1‑69 (MQ=255)
                             gATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGAt                                      >  1:1230940/1‑69 (MQ=255)
                              aTTTCCATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1198781/50‑1 (MQ=255)
                                ttCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATaac                                   <  1:511183/69‑1 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:2673121/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:978401/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:3288348/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:1350549/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:3421305/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 >  1:277464/1‑69 (MQ=255)
                                  ccATCAGCAATACTGGCGAAGTGTTGACGAAGGTGAGGATCACCGaa                                                       >  1:3296321/1‑47 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACTCTGAGGGTGGTGATaacaa                                 <  1:808052/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:3338045/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:3280749/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:81172/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:694649/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:765972/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:1156413/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:322371/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:204735/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:1923694/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:2336841/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:1919110/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:1719173/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:1595561/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:3033054/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:117549/68‑1 (MQ=255)
                                   cATCAGCAATACTGGCGACGTGTTGACGAAGGTGAGCATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:2932546/68‑1 (MQ=255)
                                   cATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:2849904/68‑1 (MQ=255)
                                   cATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATaacaa                                 <  1:3083725/68‑1 (MQ=255)
                                          aaTACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTg                         >  1:3006054/1‑69 (MQ=255)
                                          aaTACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTg                         >  1:1270520/1‑69 (MQ=255)
                                          aaTACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTg                         >  1:3394240/1‑69 (MQ=255)
                                          aaTACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTg                         >  1:2246377/1‑69 (MQ=255)
                                          aaTACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTg                         >  1:217861/1‑69 (MQ=255)
                                           aTACTGGCGACGTGTTGACGAAGGTGAGGATCACCGa                                                        <  1:1817353/37‑1 (MQ=25)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAacaaca                               >  1:3262022/1‑59 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:1289407/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:1175759/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:3168240/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:2159130/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:545469/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:2903043/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:2930369/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:2931141/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:3095023/1‑68 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAg                       >  1:2868784/1‑67 (MQ=255)
                                              ccGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAg                       >  1:1431849/1‑67 (MQ=255)
                                              ccGATGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAgg                      >  1:3414252/1‑68 (MQ=255)
                                                 gTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGaaaa                  <  1:751565/69‑1 (MQ=255)
                                                                 aGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGAAAAAACGGTTAAGCAGCGa  >  1:2332324/1‑69 (MQ=255)
                                                                 aGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGAAAAAACGGTTAAGCAGCGa  >  1:3063741/1‑69 (MQ=255)
                                                                 aGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGAAAAAACGGTTAAGCAGCGa  >  1:1175383/1‑69 (MQ=255)
                                                                 |                                                                    
GTGGCGTATTCACCGCGCTCAAACCCCATGATTTCCATCAGCAATACCGGTGACGTGTTGACGAAAGTGAGGATCACCGAAACGCTGAGGGTGGTGATAACAACACGGCTGAGGAAAAAACGGTTAAGCAGCGA  >  minE/2482139‑2482272

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: