Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 310,821 C→T 9.3% L684L (CTG→CTA) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE310,8210CT 9.3% 249.8 / 11.3 86L684L (CTG→CTA) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base T (4/4);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATC  >  minE/310752‑310890
                                                                     |                                                                     
cGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATc                                                                       <  1:2051395/70‑1 (MQ=255)
 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:2867036/1‑70 (MQ=255)
           caagcaaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc                                                            >  1:3703668/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:1383708/71‑1 (MQ=255)
                aaCTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc                                                      <  1:2924668/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)
                     gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                                   <  1:3880562/69‑1 (MQ=255)
                     gTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGt                                                   <  1:1244384/69‑1 (MQ=255)
                       tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                                 >  1:637619/1‑69 (MQ=255)
                       tACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTg                                                 >  1:2920819/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:2046649/48‑1 (MQ=255)
                            gcCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtca                                                                 <  1:2684378/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:2998070/1‑71 (MQ=255)
                                                                    tCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATc  >  1:95307/1‑71 (MQ=255)
                                                                     |                                                                     
CGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGCAGTACCCAGTTCCACGATTGCGGGCAGGTTAAAGGCGAAAACCATC  >  minE/310752‑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: