Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,788 C→T 9.5% L695L (CTG→CTA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,7880CT 9.5% 240.9 / 12.0 84L695L (CTG→CTAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (4/4);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTC  >  minE/310720‑310856
                                                                    |                                                                    
ggACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                                    >  1:3130390/1‑71 (MQ=255)
  aCGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                                  <  1:2614093/71‑1 (MQ=255)
   cGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                                 >  1:2941970/1‑71 (MQ=255)
       cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                              <  1:2918122/70‑1 (MQ=255)
       cGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                             >  1:3228633/1‑71 (MQ=255)
               aaCATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                                     >  1:3483735/1‑71 (MQ=255)
                  atatCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGc                                                   >  1:2914019/1‑70 (MQ=255)
                      cAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCaa                                                      <  1:996557/63‑1 (MQ=255)
                         ggTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                             <  1:827407/53‑1 (MQ=255)
                         ggTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                             <  1:2500475/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:2062369/1‑69 (MQ=255)
                                cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                      <  1:1078909/69‑1 (MQ=255)
                                cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                                      >  1:3141402/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:1094045/44‑1 (MQ=255)
                                  ctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                             <  1:575452/44‑1 (MQ=255)
                                    gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                                  >  1:3344638/1‑37 (MQ=255)
                                    gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                                  >  1:813155/1‑37 (MQ=255)
                                    gctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                 >  1:3763047/1‑70 (MQ=255)
                                     ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                 <  1:471647/69‑1 (MQ=255)
                                     ctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                                 <  1:901874/69‑1 (MQ=255)
                                       tCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtc                                >  1:1954965/1‑68 (MQ=255)
                                         tGCAAGCAACTTGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagt                           >  1:316966/1‑71 (MQ=255)
                                         accaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagt                           >  1:351289/3‑71 (MQ=255)
                                         accaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagt                           >  1:3579898/3‑71 (MQ=255)
                                          ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                          >  1:2573173/2‑71 (MQ=255)
                                          ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                          >  1:2114190/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:3896191/71‑1 (MQ=255)
                                           caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                         <  1:2711687/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:92915/1‑60 (MQ=255)
                                            aagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                                   >  1:3301772/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:3215487/71‑1 (MQ=255)
                                                  cTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg                  <  1:3873835/71‑1 (MQ=255)
                                                   tGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCa                                        <  1:2366504/48‑1 (MQ=255)
                                                    ggTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                                 >  1:983871/1‑38 (MQ=255)
                                                    ggTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                                 >  1:399158/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)
                                                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                               <  1:1557665/48‑1 (MQ=255)
                                                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                               <  1:2684378/48‑1 (MQ=255)
                                                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                               <  1:2046649/48‑1 (MQ=255)
                                                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAcc        <  1:2932648/71‑1 (MQ=255)
                                                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAcc        <  1:3127684/71‑1 (MQ=255)
                                                             ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt                >  1:3950220/1‑62 (MQ=255)
                                                             ccTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt                >  1:1222993/1‑62 (MQ=255)
                                                              cTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTAcc        >  1:867461/1‑69 (MQ=255)
                                                              cTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCa      <  1:844486/71‑1 (MQ=255)
                                                                 aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCa             >  1:3917331/1‑61 (MQ=255)
                                                                 aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAg     <  1:534599/69‑1 (MQ=255)
                                                                 aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAg     <  1:3116239/69‑1 (MQ=255)
                                                                 aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt   >  1:2056128/1‑71 (MQ=255)
                                                                 aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt   >  1:3010988/1‑71 (MQ=255)
                                                                   tCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc  >  1:2475754/1‑70 (MQ=255)
                                                                    cAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg               <  1:216579/56‑1 (MQ=255)
                                                                    |                                                                    
GGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTC  >  minE/310720‑310856

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: