Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,649 G→A 8.0% D408D (GAC→GAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,6490GA 8.0% 274.9 / 9.3 87D408D (GAC→GATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base A (4/3);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

GGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAG  >  minE/311580‑311718
                                                                     |                                                                     
ggTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGt                                                                      <  1:3891864/71‑1 (MQ=255)
  tAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCa                                                                    >  1:810034/1‑71 (MQ=255)
  tAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCa                                                                    <  1:3209438/71‑1 (MQ=255)
    gCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATc                                                                  <  1:333877/71‑1 (MQ=255)
     cttcttTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                 <  1:2068252/71‑1 (MQ=255)
     cttcttTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATcc                                                                 <  1:794742/71‑1 (MQ=255)
       tcttTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccac                                                               >  1:711409/1‑71 (MQ=255)
         ttttGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacc                                                              >  1:2165208/1‑70 (MQ=255)
          tttGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            >  1:1891190/1‑71 (MQ=255)
            tGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                             >  1:655654/1‑68 (MQ=255)
            tGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            >  1:120840/1‑69 (MQ=255)
            tGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            >  1:516705/1‑69 (MQ=255)
            tGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            >  1:702067/1‑69 (MQ=255)
               cggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAgg                                                        >  1:96511/1‑70 (MQ=255)
               cggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                       >  1:2992948/1‑71 (MQ=255)
                ggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGt                                                                      <  1:3644213/55‑1 (MQ=255)
                ggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACa                                                          <  1:3951450/67‑1 (MQ=255)
                ggcAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACa                                                          <  1:941101/67‑1 (MQ=255)
                   aaaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                             <  1:371832/61‑1 (MQ=255)
                   aaaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                             <  1:2918154/61‑1 (MQ=255)
                   aaaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATccacca                                                             <  1:1578762/61‑1 (MQ=255)
                   aaaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            >  1:1315617/1‑62 (MQ=255)
                   aaaCCTTCTTCCGCCATAACACGCTAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                  >  1:3383171/1‑71 (MQ=255)
                    aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                  <  1:290721/71‑1 (MQ=255)
                    aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                  >  1:1955805/1‑71 (MQ=255)
                    aaCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATg                                                  <  1:2664213/71‑1 (MQ=255)
                     aCCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATgg                                                 >  1:1359380/1‑71 (MQ=255)
                      ccTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAAc                                                           >  1:1686300/1‑60 (MQ=255)
                       cttcttCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCg                                               <  1:123282/71‑1 (MQ=255)
                        ttcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGc                                                >  1:2761354/1‑69 (MQ=255)
                        ttcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                              >  1:3515405/1‑71 (MQ=255)
                         tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                       >  1:1427650/1‑61 (MQ=255)
                         tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                             >  1:488728/1‑71 (MQ=255)
                         tcttcCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAg                                             >  1:3675685/1‑71 (MQ=255)
                             ccGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCa                                           <  1:1651762/69‑1 (MQ=255)
                              cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                         <  1:903897/70‑1 (MQ=255)
                              cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                         <  1:3713951/70‑1 (MQ=255)
                              cGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCAcc                                         <  1:274408/70‑1 (MQ=255)
                                ccAAAACACGCTCAACGTTTTCTACCACAACGATGCCGTCATc                                                                  >  1:2657470/1‑43 (MQ=38)
                                 cATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATcc                                     >  1:333152/1‑71 (MQ=255)
                                     acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                 >  1:2133125/1‑71 (MQ=255)
                                     acacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGa                                  >  1:3355033/1‑70 (MQ=255)
                                      cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            <  1:1747085/43‑1 (MQ=255)
                                      cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCaa                                                            <  1:3498776/43‑1 (MQ=255)
                                      cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATgg                                                 <  1:1873226/54‑1 (MQ=255)
                                      cacGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGaa                                 >  1:3239778/1‑70 (MQ=255)
                                       acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                      <  1:3904518/64‑1 (MQ=255)
                                       acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATc                                      <  1:3175583/64‑1 (MQ=255)
                                       acGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                               <  1:1348260/71‑1 (MQ=255)
                                        cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACa                                                          >  1:3380469/1‑43 (MQ=255)
                                        cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                              >  1:3003479/1‑71 (MQ=255)
                                        cGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAt                              <  1:1507946/71‑1 (MQ=255)
                                          cTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACAtt                             <  1:1341490/70‑1 (MQ=255)
                                           tCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                       >  1:2432467/1‑43 (MQ=255)
                                             aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGa                                              >  1:1735466/1‑50 (MQ=255)
                                             aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGtt                          <  1:2023695/70‑1 (MQ=255)
                                             aaCGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTa                         <  1:692119/71‑1 (MQ=255)
                                               cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                       <  1:848997/71‑1 (MQ=255)
                                               cGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGc                       >  1:1053744/1‑71 (MQ=255)
                                                gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGc                                                       >  1:3020226/1‑38 (MQ=255)
                                                gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                      >  1:992196/1‑71 (MQ=255)
                                                gTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCg                      <  1:2724515/71‑1 (MQ=255)
                                                 ttttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAg                        <  1:1369796/68‑1 (MQ=255)
                                                   ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                   >  1:282127/1‑71 (MQ=255)
                                                   ttCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCgtgt                   >  1:2238548/1‑71 (MQ=255)
                                                    tCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGtt                  >  1:3860039/1‑71 (MQ=255)
                                                      tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                >  1:1699142/1‑71 (MQ=255)
                                                      tACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                >  1:3897047/1‑71 (MQ=255)
                                                      caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                >  1:3350701/2‑69 (MQ=9)
                                                      caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                >  1:482721/2‑69 (MQ=9)
                                                      caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttga                >  1:93297/2‑69 (MQ=9)
                                                      caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                 <  1:1540903/69‑2 (MQ=9)
                                                      caccaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGttg                 <  1:193640/69‑2 (MQ=9)
                                                       accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt               <  1:1366210/71‑1 (MQ=9)
                                                       accaccACTATCGCATCATCGACGAGCAGCCCGATGGCAAGCACCATCCCGAACATCGTTAGTGTGTTGAt               >  1:2544236/1‑71 (MQ=9)
                                                       aCCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt               >  1:3452488/1‑71 (MQ=255)
                                                        ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTAt               <  1:1304492/70‑1 (MQ=255)
                                                        ccACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATc              <  1:1727597/71‑1 (MQ=255)
                                                         cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACa                               <  1:3109005/53‑1 (MQ=255)
                                                         cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTa                <  1:1515133/68‑1 (MQ=255)
                                                         cacaACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg             <  1:564412/71‑1 (MQ=255)
                                                         cacaACGATCGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCg             <  1:1458078/71‑1 (MQ=255)
                                                           caACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgag           <  1:2228932/71‑1 (MQ=255)
                                                            aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATccc                                    >  1:1344730/1‑45 (MQ=255)
                                                            aaCGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCgaga          <  1:2189356/71‑1 (MQ=255)
                                                                 tGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGcc      <  1:2971964/70‑1 (MQ=255)
                                                                    cGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAg  <  1:3797414/71‑1 (MQ=255)
                                                                     |                                                                     
GGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTTTTCTACCACAACGATGGCGTCATCCACCAACAGGCCGATGGCGAGCACCATCCCGAACATTGTTAGCGTGTTTATCGAGAAGCCAAAG  >  minE/311580‑311718

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: