Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,823 G→A 9.5% L684L (CTG→TTG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8230GA 9.5% 244.9 / 10.8 84L684L (CTG→TTG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (4/4);  total (39/45)
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

GCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATC  >  minE/310753‑310890
                                                                      |                                                                   
gctgctTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                                                                     <  1:646128/71‑1 (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:1314353/70‑1 (MQ=255)
         ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                                                            >  1:2573173/2‑71 (MQ=255)
         ccaagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtc                                                            >  1:2114190/2‑71 (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:3703668/1‑70 (MQ=255)
          caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc                                                            >  1:2867036/1‑70 (MQ=255)
           aagcaaCTGGTTACGGGCCTGGGTTAGGGCATCGTGACCCAGCCCAGCCTGATCAATTAACtcaaagtcga                                                          <  1:396480/71‑3 (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)
               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)
                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                                   <  1:1244384/69‑1 (MQ=255)
                    gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                                   <  1:3880562/69‑1 (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)
                           gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                                 <  1:2684378/48‑1 (MQ=255)
                           gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                                 <  1:1557665/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:3116239/69‑1 (MQ=255)
                                aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAg                                       <  1:534599/69‑1 (MQ=255)
                                aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt                                     >  1:3010988/1‑71 (MQ=255)
                                aGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGtt                                     >  1:2056128/1‑71 (MQ=255)
                                  tCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTc                                    >  1:2475754/1‑70 (MQ=255)
                                   cAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                                 <  1:216579/56‑1 (MQ=255)
                                    aGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa                                                         <  1:2454866/47‑1 (MQ=255)
                                    aGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca                                                           >  1:615354/1‑45 (MQ=255)
                                      tttttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                                   <  1:1105130/51‑1 (MQ=255)
                                        tttCGTGACCAAGGCCCGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGGTGCAGTACCCAGTTCCa                                  >  1:425849/1‑66 (MQ=255)
                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                             >  1:3010677/1‑71 (MQ=255)
                                        tttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTACACGAtt                             >  1:2486565/1‑71 (MQ=255)
                                         ttCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                             <  1:468957/70‑1 (MQ=255)
                                            gTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAg                                       <  1:1794714/57‑1 (MQ=255)
                                              gACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCgg                                                    >  1:910805/1‑42 (MQ=255)
                                              gACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGt                                      <  1:3443286/56‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACg                                <  1:872887/61‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACg                                <  1:1040173/61‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGAtt                             <  1:2579081/64‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGc                       <  1:1876877/70‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGc                       <  1:82898/70‑1 (MQ=255)
                                               aCCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCa                      <  1:3509312/71‑1 (MQ=255)
                                                ccAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAg                     <  1:2932140/71‑1 (MQ=255)
                                                 cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTcc                                   >  1:3937777/1‑56 (MQ=255)
                                                 cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTcc                                   >  1:362563/1‑56 (MQ=255)
                                                 cAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAgg                    >  1:3119178/1‑71 (MQ=255)
                                                   aGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGtt                  <  1:2155527/71‑1 (MQ=255)
                                                   aGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGt                   >  1:2484394/1‑70 (MQ=255)
                                                       cAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTaaa               <  1:3499212/70‑1 (MQ=255)
                                                              gTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAg                                       >  1:1623697/1‑39 (MQ=255)
                                                               tCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCa                      <  1:670322/55‑1 (MQ=255)
                                                               tCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAgg                    >  1:3269946/1‑57 (MQ=255)
                                                               tCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAAc      <  1:3302760/71‑1 (MQ=255)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTaa                <  1:3855605/59‑1 (MQ=255)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTGAAGGCGAAAACCa    >  1:2266626/1‑71 (MQ=255)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCa    >  1:82497/1‑71 (MQ=255)
                                                                 aaTCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCa    <  1:1426591/71‑1 (MQ=255)
                                                                   tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCggg                        <  1:459417/49‑1 (MQ=255)
                                                                   tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTaaa               <  1:2942743/58‑1 (MQ=255)
                                                                   tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATc  >  1:95307/1‑71 (MQ=255)
                                                                   tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATc  >  1:2998070/1‑71 (MQ=255)
                                                                     aGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGaaaa       >  1:1044796/1‑64 (MQ=255)
                                                                      gCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTaaa               <  1:2234498/55‑1 (MQ=255)
                                                                      |                                                                   
GCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATC  >  minE/310753‑310890

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: