Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,779 C→G 10.1% A698A (GCG→GCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,7790CG10.1% 221.6 / 11.6 79A698A (GCG→GCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (4/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTA  >  minE/310709‑310848
                                                                      |                                                                     
ccAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                       >  1:13703/1‑71 (MQ=255)
   gACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTAGCTGCTTCTGCAAGCAACTGGTTACGCGcc                                                                    >  1:1639666/1‑71 (MQ=255)
   gACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTGCGCTGCTTCTGCAAGCAACTGGTTACGCGcc                                                                    <  1:14495/71‑1 (MQ=255)
    aCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGcc                                                                    <  1:2369559/70‑1 (MQ=255)
    aCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                                   >  1:3424576/1‑71 (MQ=255)
     ccGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                                   >  1:377911/1‑70 (MQ=255)
      cGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                                 <  1:1890409/71‑1 (MQ=255)
       gTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAg                                                                >  1:3103774/1‑71 (MQ=255)
        tttGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                               <  1:2958128/71‑1 (MQ=255)
        tttGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                               <  1:162940/71‑1 (MQ=255)
          tGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                               >  1:2213759/1‑69 (MQ=255)
          tGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                               >  1:439135/1‑69 (MQ=255)
           ggACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                                     <  1:634323/62‑1 (MQ=255)
           ggACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                            >  1:3130390/1‑71 (MQ=255)
             aCGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                          <  1:2614093/71‑1 (MQ=255)
              cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                         >  1:2941970/1‑71 (MQ=255)
                tACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAg                                                                <  1:2650051/62‑1 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                       >  1:1071104/1‑53 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                      <  1:2918122/70‑1 (MQ=255)
                  cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                     >  1:3228633/1‑71 (MQ=255)
                          aaCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                             >  1:3483735/1‑71 (MQ=255)
                           aCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGGCTGAGt                                                               <  1:2706700/52‑1 (MQ=255)
                             atatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                           >  1:2914019/1‑70 (MQ=255)
                                 cAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                              <  1:996557/63‑1 (MQ=255)
                                    ggTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                     <  1:2500475/53‑1 (MQ=255)
                                    ggTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                     <  1:827407/53‑1 (MQ=255)
                                         ttCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                              <  1:1884911/55‑1 (MQ=255)
                                         ttCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                <  1:1746145/69‑1 (MQ=255)
                                         ttCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                              >  1:2655032/1‑71 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                             >  1:3455283/1‑38 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                             >  1:2907207/1‑38 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                                  >  1:711949/1‑49 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                              <  1:1078909/69‑1 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                              >  1:3141402/1‑69 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                              >  1:2062369/1‑69 (MQ=255)
                                           cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                             <  1:2051395/70‑1 (MQ=255)
                                            gctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                     <  1:48507/61‑1 (MQ=255)
                                            gctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                           <  1:646128/71‑1 (MQ=255)
                                             ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                     <  1:575452/44‑1 (MQ=255)
                                             ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                     <  1:1094045/44‑1 (MQ=255)
                                               gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                          >  1:813155/1‑37 (MQ=255)
                                               gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                          >  1:3344638/1‑37 (MQ=255)
                                               gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                         >  1:3763047/1‑70 (MQ=255)
                                                ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                         <  1:901874/69‑1 (MQ=255)
                                                ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                         <  1:471647/69‑1 (MQ=255)
                                                  tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                        >  1:1954965/1‑68 (MQ=255)
                                                    tGCAAGCAACTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt                   >  1:316966/1‑71 (MQ=255)
                                                    accaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagt                   >  1:3579898/3‑71 (MQ=255)
                                                    accaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagt                   >  1:351289/3‑71 (MQ=255)
                                                     ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                  >  1:2114190/2‑71 (MQ=255)
                                                     ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                  >  1:2573173/2‑71 (MQ=255)
                                                     ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                  <  1:1314353/70‑1 (MQ=255)
                                                     ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagt                   <  1:2486406/69‑1 (MQ=255)
                                                      caagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcg                 <  1:2758862/71‑2 (MQ=255)
                                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                 <  1:2711687/71‑1 (MQ=255)
                                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                 <  1:3896191/71‑1 (MQ=255)
                                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc                  >  1:2867036/1‑70 (MQ=255)
                                                      caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc                  >  1:3703668/1‑70 (MQ=255)
                                                       aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                <  1:396480/71‑3 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTg                                                   <  1:3422555/36‑1 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                    <  1:3277726/51‑1 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                           >  1:3301772/1‑60 (MQ=255)
                                                       aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                           >  1:92915/1‑60 (MQ=255)
                                                        agcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaac               >  1:2726355/1‑70 (MQ=255)
                                                        agcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                 >  1:905232/1‑69 (MQ=255)
                                                         gcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaag              >  1:829777/1‑71 (MQ=255)
                                                           aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                              <  1:1002894/37‑1 (MQ=255)
                                                           aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc            <  1:2924668/71‑1 (MQ=255)
                                                           aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc            <  1:1383708/71‑1 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACTAAGGCCAGCCTGGTCAATCAGc                          >  1:2284306/1‑55 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGc                          >  1:480434/1‑55 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg          <  1:3873835/71‑1 (MQ=255)
                                                             cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg          <  1:3215487/71‑1 (MQ=255)
                                                              tGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                <  1:2366504/48‑1 (MQ=255)
                                                               ggTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                         >  1:399158/1‑38 (MQ=255)
                                                               ggTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                         >  1:983871/1‑38 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                         >  1:2292176/1‑37 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt         <  1:3880562/69‑1 (MQ=255)
                                                                gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt         <  1:1244384/69‑1 (MQ=255)
                                                                  tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                      >  1:1734052/1‑38 (MQ=255)
                                                                  tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg       >  1:2920819/1‑69 (MQ=255)
                                                                  tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg       >  1:637619/1‑69 (MQ=255)
                                                                     gcgcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt                   <  1:913837/54‑1 (MQ=255)
                                                                     gcgcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTa  <  1:365271/71‑1 (MQ=255)
                                                                      cgcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                             >  1:3863507/1‑43 (MQ=255)
                                                                      |                                                                     
CCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTA  >  minE/310709‑310848

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: