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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5340CG11.2% 229.8 / 15.5 80L113L (CTG→CTCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (6/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGA  >  minE/312467‑312604
                                                                   |                                                                      
aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGccct                                                                      >  1:2937306/1‑68 (MQ=255)
  ctgctgTATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                                   <  1:2824893/71‑1 (MQ=255)
  ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                                   <  1:1574427/71‑1 (MQ=255)
          tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                            <  1:46484/70‑1 (MQ=255)
          tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                            <  1:2469988/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)
                     cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                                   >  1:1001038/1‑52 (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)
                                tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                            <  1:983296/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:1272280/1‑61 (MQ=255)
                                          cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                     >  1:2182012/1‑61 (MQ=255)
                                          cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                     >  1:1338591/1‑61 (MQ=255)
                                           ttcttGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                                              <  1:3660103/35‑1 (MQ=255)
                                            tcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAg                          >  1:35390/1‑70 (MQ=255)
                                             cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                              >  1:1235457/1‑49 (MQ=255)
                                             cttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                                              >  1:1240221/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:2110872/1‑69 (MQ=12)
                                                      gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                >  1:3535914/1‑69 (MQ=12)
                                                      gcaACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATCAGGATCGGTCCCGGACtg                >  1:3531102/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:2178774/71‑1 (MQ=255)
                                                              aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg       <  1:841145/71‑1 (MQ=255)
                                                              aTCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAg       <  1:1578154/71‑1 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAggg     >  1:1257401/1‑70 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAgg      <  1:3389683/69‑1 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:3393173/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:3519374/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:3957681/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:3445840/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:1535932/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:978908/1‑71 (MQ=255)
                                                                 gCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCACGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGt    >  1:2122566/1‑71 (MQ=255)
                                                                  ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg   <  1:3038001/71‑1 (MQ=255)
                                                                  ccAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTg   <  1:1073045/71‑1 (MQ=255)
                                                                   cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat                                   >  1:1442397/1‑38 (MQ=255)
                                                                   cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa  >  1:2646272/1‑71 (MQ=255)
                                                                   cAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGa  >  1:3511589/1‑71 (MQ=255)
                                                                   |                                                                      
AGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAGGGTGA  >  minE/312467‑312604

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: