Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,832 C→G 7.4% A347A (GCG→GCC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8320CG 7.4% 242.3 / 5.4 80A347A (GCG→GCCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (29/45);  new base G (4/2);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTA  >  minE/311762‑311900
                                                                      |                                                                    
cAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                                                      >  1:2326278/1‑71 (MQ=255)
 aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATc                                                                      <  1:1798241/70‑1 (MQ=255)
 aaTGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCg                                                                     <  1:1150541/71‑1 (MQ=255)
   tGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCt                                                                   >  1:207882/1‑71 (MQ=255)
    ggTGGGGATCAGCGTTGCTCGCATATTCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCtt                                                                  >  1:3127576/1‑71 (MQ=2)
      tCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCg                                                                <  1:2622455/71‑1 (MQ=255)
              aaCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                         <  1:3276463/70‑1 (MQ=255)
               aCGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtt                                                       <  1:884341/71‑1 (MQ=255)
                cGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGttt                                                      <  1:3800745/71‑1 (MQ=255)
                  tCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                    <  1:2852763/71‑1 (MQ=255)
                      gcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCg                                                         <  1:1585492/62‑1 (MQ=255)
                         gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                    <  1:568395/64‑1 (MQ=255)
                         gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACtt                                             >  1:2785062/1‑71 (MQ=255)
                         gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACtt                                             <  1:1908208/71‑1 (MQ=255)
                         gAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                                              >  1:243692/1‑70 (MQ=255)
                          atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                                            >  1:688874/4‑71 (MQ=16)
                          aaGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                            <  1:2326430/71‑1 (MQ=255)
                           aGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCg                                           <  1:576963/71‑1 (MQ=255)
                            gTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                          <  1:3177601/71‑1 (MQ=255)
                             ttCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTg                                         <  1:2640254/71‑1 (MQ=21)
                               cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGc                                                                    >  1:3766991/1‑42 (MQ=255)
                               cTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                  >  1:2822574/1‑60 (MQ=255)
                                 gCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                    <  1:2295321/56‑1 (MQ=255)
                                  cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                     <  1:2749706/70‑1 (MQ=255)
                                  cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCCTGAATa                                     <  1:912443/70‑1 (MQ=255)
                                   aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGa                                                               >  1:23756/1‑43 (MQ=255)
                                   aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                                    <  1:998935/70‑1 (MQ=255)
                                   aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                                    <  1:41066/70‑1 (MQ=255)
                                   aGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                                    <  1:1467086/70‑1 (MQ=255)
                                    ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                     >  1:3036528/1‑68 (MQ=255)
                                     gAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGtt                                                       <  1:1705182/49‑1 (MQ=16)
                                     gAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                     >  1:3152570/1‑67 (MQ=255)
                                      aaCATATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaag                                >  1:936893/1‑68 (MQ=255)
                                      aaCAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGAATagaaag                                >  1:2180403/1‑68 (MQ=18)
                                       aCAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTa                                                    <  1:665536/50‑1 (MQ=255)
                                        cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt                              <  1:2916420/71‑1 (MQ=255)
                                         aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                             >  1:2720296/1‑71 (MQ=255)
                                         aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                             <  1:2519929/71‑1 (MQ=255)
                                          gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                             <  1:2376201/70‑1 (MQ=255)
                                          gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                            <  1:16216/71‑1 (MQ=255)
                                           aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                            >  1:1030851/1‑54 (MQ=255)
                                           aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                           >  1:2419244/1‑71 (MQ=255)
                                            tACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGtt                                                       >  1:1949651/1‑42 (MQ=255)
                                              cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                  <  1:1544294/45‑1 (MQ=255)
                                              cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                         >  1:3619781/1‑70 (MQ=255)
                                              cATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                         >  1:3199439/1‑70 (MQ=255)
                                               aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                       <  1:458249/71‑1 (MQ=255)
                                                tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                         <  1:1422026/68‑1 (MQ=255)
                                                 aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                     >  1:3267093/1‑71 (MQ=255)
                                                 aaCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                     >  1:3135286/1‑71 (MQ=255)
                                                   ccAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGc                   <  1:1785394/71‑1 (MQ=255)
                                                    cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACggag                  <  1:1311492/71‑3 (MQ=255)
                                                    cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                  >  1:1722786/1‑71 (MQ=255)
                                                     aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                 <  1:1272663/71‑1 (MQ=255)
                                                     aGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                 <  1:553765/71‑1 (MQ=255)
                                                      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAAcc                                                 <  1:2593737/38‑1 (MQ=255)
                                                      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                  >  1:1449251/1‑37 (MQ=255)
                                                      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                <  1:1442507/71‑1 (MQ=255)
                                                      ggAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                <  1:1604772/71‑1 (MQ=255)
                                                        aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt              <  1:2208850/71‑1 (MQ=255)
                                                          cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt              <  1:2805917/69‑1 (MQ=255)
                                                            cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACgg                    >  1:392723/1‑61 (MQ=255)
                                                            cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg          <  1:1542576/71‑1 (MQ=255)
                                                            cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg          >  1:1150993/1‑71 (MQ=255)
                                                            cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg          >  1:1162280/1‑71 (MQ=255)
                                                             gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGt         <  1:2697793/71‑1 (MQ=255)
                                                             gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGt         <  1:1905565/71‑1 (MQ=255)
                                                              aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa        <  1:2057586/71‑1 (MQ=255)
                                                              aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa        >  1:2392995/1‑71 (MQ=255)
                                                              aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa        >  1:1894549/1‑71 (MQ=255)
                                                              aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa        <  1:2483885/71‑1 (MQ=255)
                                                               ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCa                           >  1:2842723/1‑51 (MQ=255)
                                                               ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt       <  1:3598567/71‑1 (MQ=255)
                                                               ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt       >  1:609365/1‑71 (MQ=255)
                                                               ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTAt       >  1:1509933/1‑71 (MQ=255)
                                                               ggATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAAGGCGTGGTGTCGTAt       >  1:3645292/1‑71 (MQ=255)
                                                                gatgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATg      <  1:3214708/71‑1 (MQ=255)
                                                                    atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGt   <  1:1875796/70‑1 (MQ=255)
                                                                    atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa  <  1:2937811/71‑1 (MQ=255)
                                                                     tCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                  >  1:158440/1‑38 (MQ=255)
                                                                     tCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                  >  1:795887/1‑38 (MQ=255)
                                                                     tCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTa  <  1:85513/70‑1 (MQ=255)
                                                                      |                                                                    
CAATGGTCGGAATCAACGTCGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTA  >  minE/311762‑311900

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: