Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,496 A→G 22.6% F459F (TTT→TTC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,4960AG22.6% 253.8 / 60.2 115F459F (TTT→TTCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base G (12/14);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 125 total aligned reads displayed.

GTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCC  >  minE/311426‑311566
                                                                      |                                                                      
gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                                        <  1:93049/71‑1 (MQ=255)
gTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCa                                                                        <  1:1905006/71‑1 (MQ=255)
   ccGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      <  1:3573257/70‑1 (MQ=255)
    cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                                    <  1:2455357/71‑1 (MQ=255)
    cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                                    <  1:1989027/71‑1 (MQ=255)
    cGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                                    <  1:669992/71‑1 (MQ=255)
     gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCTATAAATTGCCCCAGTAGAACCGCCgaagaa                                                                   <  1:3210604/68‑1 (MQ=12)
     gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                                                                   <  1:3393074/68‑1 (MQ=12)
     gaaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaa                                                                   <  1:2556783/68‑1 (MQ=12)
     gACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaag                                                                     >  1:1456554/1‑69 (MQ=255)
      aaaGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCgaagaat                                                                  >  1:4834798/3‑70 (MQ=12)
       cAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAgg                                                                 >  1:96034/1‑71 (MQ=255)
        aGCGCCATTGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGc                                                                >  1:454850/1‑71 (MQ=12)
          cgcCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                                                   >  1:4182131/1‑66 (MQ=255)
              aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          <  1:4161754/71‑1 (MQ=255)
              aTTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          <  1:5331884/71‑1 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCaaa                                                                      <  1:3313267/58‑1 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          >  1:1600655/1‑70 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          >  1:486730/1‑70 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                                         <  1:64888/71‑1 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                                         <  1:3403910/71‑1 (MQ=255)
               ttGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGt                                                         <  1:5332686/71‑1 (MQ=255)
                tGCCGAAACGATGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCgga                                                         >  1:1026478/1‑69 (MQ=9)
                 gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:1691408/1‑71 (MQ=255)
                 gCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:5294934/1‑71 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       <  1:4707146/70‑1 (MQ=255)
                  cTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:4232480/71‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaagaa                                                                   <  1:2727620/55‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:4185394/68‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:771359/68‑1 (MQ=255)
                     aaaCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACg                                                      <  1:5077240/68‑1 (MQ=255)
                       aCAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          >  1:2289126/1‑62 (MQ=255)
                        cAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAcc                                                <  1:45062/71‑1 (MQ=255)
                          aTGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATAccgc                                              >  1:4106045/1‑71 (MQ=255)
                           tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                                             <  1:4220291/71‑2 (MQ=9)
                           tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                                             >  1:4012000/1‑70 (MQ=9)
                           tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                                             <  1:115177/71‑2 (MQ=9)
                           tGGTGATAGAGAACTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                                             <  1:4585766/71‑2 (MQ=9)
                            ggTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAaaga                                                                    <  1:2447143/47‑1 (MQ=255)
                             gTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGa                                           <  1:2636719/71‑1 (MQ=255)
                              tAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCgg                                                          <  1:400823/55‑1 (MQ=255)
                                aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGaa                                                    >  1:3368717/1‑59 (MQ=255)
                                aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        >  1:4973468/1‑71 (MQ=255)
                                aTAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAg                                        >  1:5367023/1‑71 (MQ=255)
                                 tAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACa                                         <  1:679377/69‑1 (MQ=255)
                                  agagAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTa                                      >  1:1350901/1‑71 (MQ=255)
                                   gagaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTAc                                                       >  1:5282410/1‑53 (MQ=255)
                                    agaACTGACGATAGATAGCACCACTAGAACCCCCAaag                                                                     <  1:2297252/38‑1 (MQ=38)
                                     gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa                                   <  1:1460261/71‑1 (MQ=255)
                                     gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa                                   <  1:4286376/71‑1 (MQ=255)
                                     gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa                                   <  1:3284158/71‑1 (MQ=255)
                                     gaACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa                                   <  1:2545083/71‑1 (MQ=255)
                                      aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                                  <  1:3447088/71‑1 (MQ=9)
                                      aaCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                                  <  1:950705/71‑1 (MQ=9)
                                       aCTGGCGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCAt                                  <  1:2281518/70‑1 (MQ=9)
                                       aCTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                                        <  1:1229471/48‑1 (MQ=255)
                                       aCTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATc                                 <  1:567990/71‑1 (MQ=255)
                                        cTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcg                                >  1:3512703/1‑71 (MQ=255)
                                          gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCAt                                                             >  1:753224/1‑40 (MQ=255)
                                          gACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCAt                                                             >  1:4791495/1‑40 (MQ=255)
                                           gcGATAAATTGCCCCAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                             <  1:2261693/70‑1 (MQ=12)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                     <  1:4097488/63‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                     <  1:2462337/63‑1 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                             >  1:3048474/1‑71 (MQ=255)
                                           aCGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGa                             >  1:2999704/1‑71 (MQ=255)
                                             gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                                        >  1:3021514/1‑42 (MQ=255)
                                             gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                                        >  1:2907225/1‑42 (MQ=255)
                                             gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTa                                                        >  1:2097921/1‑42 (MQ=255)
                                             gatagataGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGAt                            <  1:5466771/70‑1 (MQ=255)
                                                    agCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGAc                      >  1:2910826/1‑69 (MQ=255)
                                                     gccccAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATa                           >  1:3179630/4‑63 (MQ=12)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa                    <  1:4451474/70‑1 (MQ=255)
                                                     gCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCa                    <  1:1459092/70‑1 (MQ=255)
                                                       cccAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCAcc                                    <  1:5367771/51‑1 (MQ=12)
                                                       cccAGTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCAcc                                    <  1:1448029/51‑1 (MQ=12)
                                                         cAGTAGAACCGCCGAAGAATGCCATCGGAATAAATAccgct                                             <  1:5424972/41‑2 (MQ=9)
                                                         cAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCagag                 >  1:3760587/1‑69 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCa                                   >  1:4201841/1‑50 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:254290/1‑71 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAgcgc              >  1:5664523/1‑71 (MQ=255)
                                                          aGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGAAAGTACCATCGCGATACCGACCAGAgcgc              >  1:3427005/1‑71 (MQ=255)
                                                           gTAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGcc             >  1:5383800/1‑71 (MQ=12)
                                                           gTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGcc             >  1:786918/1‑71 (MQ=255)
                                                            tAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATc                                 <  1:1082062/50‑1 (MQ=9)
                                                            tAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATc                                 <  1:3242049/50‑1 (MQ=9)
                                                            tAGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATc                                 <  1:5592764/50‑1 (MQ=9)
                                                            tAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGccc            >  1:2052512/1‑71 (MQ=255)
                                                             aGAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCt           >  1:1240559/1‑71 (MQ=255)
                                                             aGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:780975/1‑71 (MQ=255)
                                                              gAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                             >  1:5134930/1‑52 (MQ=9)
                                                              gAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                             >  1:1058364/1‑52 (MQ=9)
                                                              gAACCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGa                             >  1:817718/1‑52 (MQ=9)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTAc                                     <  1:2898020/44‑1 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATcgcg                              >  1:1799234/1‑51 (MQ=255)
                                                              gAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCt           >  1:1090433/1‑70 (MQ=255)
                                                               aaCCGCCGAAGAATGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTg          >  1:5073295/1‑70 (MQ=9)
                                                                aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:718203/70‑1 (MQ=255)
                                                                aCCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGa         <  1:260009/70‑1 (MQ=255)
                                                                   gCCAAAGAAGGCCATCGGTACGAATACCGCGGACAGTACCATCGCGATACCGACCaga                  >  1:3391834/1‑58 (MQ=255)
                                                                      aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGcc  >  1:1816424/1‑71 (MQ=255)
                                                                      |                                                                      
GTACCGACAGCGCCATTGCTGAAACAATGGTAATAGAGAACTGACGATAGATAGCACCAGTAGAACCGCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCC  >  minE/311426‑311566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: