Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,540 C→G 13.5% L111L (CTG→CTC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5400CG13.5% 267.2 / 26.8 96L111L (CTG→CTCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (47/36);  new base G (6/7);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 101 total aligned reads displayed.

TGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTG  >  minE/312477‑312608
                                                               |                                                                    
tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                                <  1:2469988/70‑1 (MQ=255)
tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                                <  1:46484/70‑1 (MQ=255)
tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                               <  1:2443771/71‑1 (MQ=255)
tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                               <  1:263968/71‑1 (MQ=255)
  atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                                              <  1:3458726/70‑1 (MQ=255)
     tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                                                           <  1:3633194/70‑1 (MQ=255)
      ctcAACGCTCACCCCTTGCTGCTGGACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                                                         >  1:3768052/1‑71 (MQ=255)
      ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                                          <  1:1285477/70‑1 (MQ=255)
      ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                                          <  1:2309462/70‑1 (MQ=255)
      ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCa                                                         >  1:3244722/1‑71 (MQ=255)
         aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCt                                                      >  1:737035/1‑71 (MQ=255)
          aCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTg                                                     <  1:2484024/71‑1 (MQ=255)
           cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGa                                                    <  1:3679459/71‑1 (MQ=255)
               gatccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCAcc                                                       >  1:2370368/4‑64 (MQ=255)
               cACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                                              >  1:2061855/1‑57 (MQ=255)
                 ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                              <  1:3909276/71‑1 (MQ=255)
                 ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                              >  1:3584377/1‑71 (MQ=255)
                 ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                              >  1:2161056/1‑71 (MQ=255)
                    cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc                                            >  1:144054/1‑70 (MQ=255)
                     ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt                                                                 <  1:2385447/48‑1 (MQ=255)
                      tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTg                                                                <  1:2691196/48‑1 (MQ=19)
                       gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                                                            >  1:2949376/1‑51 (MQ=255)
                       gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                         <  1:1439362/70‑1 (MQ=255)
                       gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACATGCGCGata                                        >  1:4121072/1‑71 (MQ=255)
                       gctgctGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                                                   >  1:4055741/1‑44 (MQ=16)
                       gctgctGAACATCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGt                                                                   >  1:1610271/1‑44 (MQ=16)
                        ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGat                                         <  1:3652408/69‑1 (MQ=255)
                          gctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATcc                                     >  1:3016428/1‑71 (MQ=255)
                           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                  >  1:3629191/1‑57 (MQ=255)
                           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                         <  1:882196/66‑1 (MQ=255)
                           ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg                                    >  1:2171676/1‑71 (MQ=255)
                             gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                                               <  1:510482/58‑1 (MQ=255)
                              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc                                   <  1:1552607/69‑1 (MQ=255)
                              aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                 <  1:3180559/71‑1 (MQ=255)
                               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAc                                                        >  1:3997405/1‑47 (MQ=255)
                               aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt                                 >  1:4026884/1‑70 (MQ=255)
                                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                         >  1:2182012/1‑61 (MQ=255)
                                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                         >  1:1338591/1‑61 (MQ=255)
                                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                         >  1:1272280/1‑61 (MQ=255)
                                 ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                                                  <  1:3660103/35‑1 (MQ=255)
                                  tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                              >  1:35390/1‑70 (MQ=255)
                                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                  >  1:1240221/1‑49 (MQ=255)
                                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                                  >  1:1235457/1‑49 (MQ=255)
                                   cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTa                            >  1:92574/1‑71 (MQ=255)
                                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                           >  1:2360078/1‑69 (MQ=255)
                                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                           >  1:1210264/1‑69 (MQ=255)
                                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                           >  1:3028572/1‑69 (MQ=255)
                                       cGGCATCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTcccg                         <  1:737645/70‑2 (MQ=12)
                                       cGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCa                         <  1:1260947/70‑1 (MQ=255)
                                         gcagcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGAc                      <  1:304415/71‑1 (MQ=255)
                                            gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                    >  1:3535914/1‑69 (MQ=12)
                                            gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                    >  1:3531102/1‑69 (MQ=12)
                                            gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                    >  1:2110872/1‑69 (MQ=12)
                                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata                                        <  1:1416572/50‑1 (MQ=255)
                                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                   <  1:2788957/71‑1 (MQ=255)
                                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa                   <  1:3570354/71‑1 (MQ=255)
                                             caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa                  <  1:3336376/71‑1 (MQ=255)
                                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                >  1:2376829/1‑71 (MQ=255)
                                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg                >  1:2351372/1‑71 (MQ=255)
                                                cGGCATCGCCAGCTGCAGTTTGTTCTGCACCTTAACCTGCGCGATATCCGCAt                                 <  1:791388/53‑1 (MQ=255)
                                                 gggaTCGACAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt              <  1:3347239/68‑1 (MQ=255)
                                                 ggCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGa                       <  1:2990753/62‑1 (MQ=255)
                                                   cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt              >  1:3728487/1‑69 (MQ=255)
                                                   cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt              >  1:2012322/1‑69 (MQ=255)
                                                    aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg           <  1:1578154/71‑1 (MQ=255)
                                                    aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg           <  1:2178774/71‑1 (MQ=255)
                                                    aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg           <  1:841145/71‑1 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg         >  1:1257401/1‑70 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg          <  1:3389683/69‑1 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:3957681/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:978908/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:3519374/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:1535932/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:3393173/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:3445840/1‑71 (MQ=255)
                                                       gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt        >  1:2122566/1‑71 (MQ=255)
                                                        ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg       <  1:1073045/71‑1 (MQ=255)
                                                        ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg       <  1:3038001/71‑1 (MQ=255)
                                                         cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                       >  1:1442397/1‑38 (MQ=255)
                                                         cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa      >  1:3511589/1‑71 (MQ=255)
                                                         cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa      >  1:2646272/1‑71 (MQ=255)
                                                         aaGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa      >  1:3665785/2‑71 (MQ=255)
                                                          aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt     >  1:2523077/1‑71 (MQ=255)
                                                          aGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAGAGGTCAGGGTGAt     >  1:3466175/1‑71 (MQ=255)
                                                          aGCTGCAGTTTGTTATGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt     >  1:3491305/1‑71 (MQ=255)
                                                            cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                       <  1:3258854/51‑1 (MQ=12)
                                                            cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                       <  1:1655155/51‑1 (MQ=12)
                                                            cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                       <  1:1887770/51‑1 (MQ=12)
                                                            cTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGa                       <  1:734976/51‑1 (MQ=12)
                                                            cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg       <  1:1203974/67‑1 (MQ=255)
                                                            cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGAt     >  1:3189015/1‑69 (MQ=255)
                                                            cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGAATCAGTACCAGACTCAAAGGTCAGGGTGATCt   >  1:2820966/1‑71 (MQ=255)
                                                            cTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCAGCATCAGTACCAGACTCAAAGGTCAGGGTGAt     >  1:318973/1‑69 (MQ=255)
                                                             tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAtg                                <  1:3866301/41‑2 (MQ=255)
                                                             tGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg  >  1:3588475/1‑71 (MQ=255)
                                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg           <  1:3910271/60‑1 (MQ=255)
                                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg  <  1:2848358/69‑1 (MQ=255)
                                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg  <  1:2009352/69‑1 (MQ=255)
                                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg  <  1:3822765/69‑1 (MQ=255)
                                                               cAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTg  <  1:3802011/69‑1 (MQ=255)
                                                               |                                                                    
TGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGATCTG  >  minE/312477‑312608

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: