Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,531 C→G 11.5% A114A (GCG→GCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5310CG11.5% 219.0 / 14.7 78A114A (GCG→GCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (6/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAG  >  minE/312461‑312599
                                                                      |                                                                    
tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                                      <  1:4006811/71‑1 (MQ=255)
tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                                      <  1:339736/71‑1 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:1051993/1‑69 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:3920429/1‑69 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2656583/1‑69 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2661605/1‑69 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:3842665/1‑69 (MQ=255)
    gAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                                                    >  1:2771831/1‑69 (MQ=255)
      aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGccct                                                                 >  1:2937306/1‑68 (MQ=255)
        ctgctgTATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                              <  1:2824893/71‑1 (MQ=255)
        ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                                              <  1:1574427/71‑1 (MQ=255)
                tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                       <  1:2469988/70‑1 (MQ=255)
                tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTg                                                       <  1:46484/70‑1 (MQ=255)
                tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                      <  1:263968/71‑1 (MQ=255)
                tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                                      <  1:2443771/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:2309462/70‑1 (MQ=255)
                      ctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                                                 <  1:1285477/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:2161056/1‑71 (MQ=255)
                                 ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                     >  1:3584377/1‑71 (MQ=255)
                                 ccccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                                     <  1:3909276/71‑1 (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:1338591/1‑61 (MQ=255)
                                                cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat                                >  1:2182012/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:1210264/1‑69 (MQ=255)
                                                      gCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTAc                  >  1:2360078/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:3570354/71‑1 (MQ=255)
                                                            gcaACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCa          <  1:2788957/71‑1 (MQ=255)
                                                             caACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCaa         <  1:3336376/71‑1 (MQ=255)
                                                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg       >  1:2351372/1‑71 (MQ=255)
                                                               aCGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAg       >  1:2376829/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:2012322/1‑69 (MQ=255)
                                                                   cATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGt     >  1:3728487/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)
                                                                      |                                                                    
TCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATCAGTACCAGACTCAAAGGTCAG  >  minE/312461‑312599

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: