Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,510 T→C 13.5% E121E (GAA→GAG acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,5100TC13.5% 186.2 / 17.6 74E121E (GAA→GAGacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base C (4/6);  total (33/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 = 1.00e+00

GTTGATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATC  >  minE/312441‑312578
                                                                     |                                                                    
gtTGATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActt                                                                     >  1:2580555/1‑71 (MQ=255)
 ttGATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtc                                                                    >  1:2119909/1‑71 (MQ=255)
   gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTTCTGCTGAActtctt                                                                  >  1:1784556/1‑71 (MQ=255)
   gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                  >  1:1052482/1‑71 (MQ=255)
   gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                  >  1:198622/1‑71 (MQ=255)
   gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtct                                                                   <  1:1030666/70‑1 (MQ=255)
   gATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtct                                                                   <  1:1640507/70‑1 (MQ=255)
    aTAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAActtctt                                                                  <  1:3878093/70‑1 (MQ=255)
     ttaCGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                >  1:3361757/3‑71 (MQ=255)
     tAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGc                                                                >  1:3007291/1‑71 (MQ=255)
       aCGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCgg                                                              <  1:2672271/71‑1 (MQ=255)
         gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                 >  1:419610/1‑66 (MQ=255)
         gCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                 >  1:855360/1‑66 (MQ=255)
            gACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagca                                                         <  1:3648748/71‑1 (MQ=255)
             aCAACCATCAGGAAGCTGCTGGATGACTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCaa                                                        <  1:614115/71‑1 (MQ=255)
                  cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgca                                                            >  1:212754/1‑62 (MQ=255)
                  cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                   <  1:3329887/71‑1 (MQ=255)
                  cATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCa                                                   >  1:1869731/1‑71 (MQ=255)
                    tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCAt                                                  >  1:1543853/1‑70 (MQ=255)
                    tCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATc                                                 >  1:2243390/1‑71 (MQ=255)
                      aGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGc                                               <  1:3408095/71‑1 (MQ=255)
                          aGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGc                                           <  1:269802/71‑1 (MQ=255)
                            ctgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGTCATCGCCAGCTg                                         <  1:1460033/71‑1 (MQ=255)
                             tgctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTg                                         >  1:3808758/1‑70 (MQ=255)
                               ctgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAg                                      <  1:3829422/71‑1 (MQ=255)
                                tgGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGgcagc                                                          <  1:2432203/50‑1 (MQ=255)
                                  gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTg                                                                 <  1:3521365/41‑1 (MQ=255)
                                  gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGttt                                   >  1:1551572/1‑71 (MQ=255)
                                  gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                     <  1:482382/69‑1 (MQ=255)
                                  gatgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGt                                     <  1:2680066/69‑1 (MQ=255)
                                    tgatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                 >  1:2612182/1‑71 (MQ=255)
                                     gatTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCg                                                               <  1:3843372/40‑1 (MQ=255)
                                      atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                <  1:3815906/70‑1 (MQ=255)
                                      atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                <  1:3364271/70‑1 (MQ=255)
                                      atTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGtt                                <  1:2798951/70‑1 (MQ=255)
                                       tttCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCg                                                >  1:778843/1‑53 (MQ=255)
                                       tttCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                              <  1:1387447/71‑1 (MQ=255)
                                        ttCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGtt                                    <  1:2403070/64‑1 (MQ=255)
                                         tctcAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                             <  1:1227821/70‑1 (MQ=255)
                                             aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGt                                 <  1:1126872/62‑1 (MQ=255)
                                             aaCGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCAc                          <  1:2075133/69‑1 (MQ=255)
                                               cGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGc                                                    <  1:2075627/41‑1 (MQ=255)
                                                    atccccTGCTGCTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTg                 >  1:3585923/3‑71 (MQ=255)
                                                    aCCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                 >  1:240322/1‑71 (MQ=255)
                                                      cccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgc                <  1:9975/70‑1 (MQ=255)
                                                       ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc              >  1:2540559/1‑71 (MQ=255)
                                                       ccTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTgcgc              <  1:3724069/71‑1 (MQ=255)
                                                        cctgctgcTGAACCTCCTGCGGCAGCAACGTCGTGGCGAGCTGGAGTTTGTTCTGCa                           <  1:980077/55‑1 (MQ=255)
                                                        cctgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTTCa                           <  1:1587391/55‑1 (MQ=255)
                                                        cTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTg                             >  1:2035699/1‑55 (MQ=255)
                                                         ttGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGc                            <  1:604006/55‑1 (MQ=255)
                                                          tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa            >  1:1081387/1‑70 (MQ=255)
                                                          tgctgcTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa            >  1:2814008/1‑70 (MQ=255)
                                                          tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa            >  1:114107/1‑70 (MQ=255)
                                                          tgctgcTGAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGa            >  1:69298/1‑70 (MQ=255)
                                                           gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGaa                     >  1:3295860/1‑60 (MQ=255)
                                                           gctgctGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGa            <  1:923469/69‑1 (MQ=255)
                                                            ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat         <  1:1885283/71‑1 (MQ=255)
                                                            ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGatat         <  1:2940871/71‑1 (MQ=255)
                                                            ctgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGata          <  1:3637246/70‑1 (MQ=255)
                                                            ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCt                        <  1:1467817/56‑1 (MQ=255)
                                                            ctgctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCAGCACCTGCACTTGCGCGatat         <  1:2365335/71‑1 (MQ=255)
                                                             tgctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCt                              >  1:3035932/1‑49 (MQ=255)
                                                              gctgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTgcgc              <  1:2523939/64‑1 (MQ=255)
                                                               ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTg                 >  1:3640926/1‑60 (MQ=255)
                                                               ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg      >  1:1844797/1‑71 (MQ=255)
                                                               ctgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCg      >  1:734831/1‑71 (MQ=255)
                                                                tgAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGat           >  1:1130171/1‑65 (MQ=255)
                                                                tgAACCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcag      >  1:3561527/1‑68 (MQ=255)
                                                                 gAACTTCTTGCGGCAGCAACTGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc     >  1:846843/1‑70 (MQ=255)
                                                                 gAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCa    >  1:265643/1‑71 (MQ=255)
                                                                 gAACTTCTTGCGGCAGCAACGGCATCGCCAGCAGCAGtttc                                  >  1:851633/1‑40 (MQ=255)
                                                                  aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                    <  1:2731992/54‑1 (MQ=255)
                                                                  aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAAc                    <  1:9925/54‑1 (MQ=255)
                                                                  aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc     <  1:849499/69‑1 (MQ=255)
                                                                  aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc     <  1:854367/69‑1 (MQ=255)
                                                                  aaCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGc     <  1:859700/69‑1 (MQ=255)
                                                                  aaCCTCCTGCGGCAGCAACGGCGTGGCGAGCTGGAGTTTGTTCTGCACCTGCACTTGCGCGATATcagg     <  1:598713/69‑4 (MQ=19)
                                                                   aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt   >  1:3482184/1‑70 (MQ=255)
                                                                   aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt   >  1:3351000/1‑70 (MQ=255)
                                                                   aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt   >  1:2828559/1‑70 (MQ=255)
                                                                   aCTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATc  <  1:1409466/71‑1 (MQ=255)
                                                                    cttcttGCGGCAGCAACGGTATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCAt   >  1:2097907/1‑69 (MQ=255)
                                                                    cttcttGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCt                  >  1:2098160/1‑54 (MQ=255)
                                                                     |                                                                    
GTTGATAACGCCGACAACCATCAGGAAGCTGCTGGATGATTTCTCAACGCTCACCCCTTGCTGCTGAACTTCTTGCGGCAGCAACGGCATCGCCAGCTGCAGTTTGTTCTGCACCTGAACCTGCGCGATATCCGCATC  >  minE/312441‑312578

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: