Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,369 A→G 10.7% L970L (TTA→CTA) ‡ mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,3690AG10.7% 229.8 / 17.4 84L970L (TTA→CTA) ‡mdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base G (2/7);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGCATCATCTCAACGGCA  >  minE/2584299‑2584439
                                                                      |                                                                      
gCACCATGACTGATAACCAGCGTCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaa                                                                        <  1:1738442/71‑1 (MQ=255)
  aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       <  1:3193860/70‑1 (MQ=255)
  aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAAc                                                                      <  1:2938334/71‑1 (MQ=255)
   ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACg                                                                     <  1:479912/71‑1 (MQ=255)
    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:460122/71‑1 (MQ=255)
    cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:1497753/71‑1 (MQ=255)
     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:3738684/70‑1 (MQ=255)
     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                                                   <  1:1152603/71‑1 (MQ=255)
     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                                                   >  1:3108791/1‑71 (MQ=255)
     aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                                                   >  1:2226244/1‑71 (MQ=255)
      tGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:3213435/69‑1 (MQ=255)
       taCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCAta                                                                 >  1:4021534/2‑70 (MQ=2)
        aCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATAGGATAGGACGCAGACGCATAcg                                                               >  1:1001754/1‑71 (MQ=255)
        aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:27258/1‑64 (MQ=255)
        aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:2387376/1‑64 (MQ=255)
        aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                                                                >  1:102082/1‑71 (MQ=255)
          tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcg                                                               <  1:3504817/70‑1 (MQ=255)
           gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGCGTAAACGCATCcgc                                                              <  1:2455645/69‑1 (MQ=255)
            aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACg                                                                     >  1:1974473/1‑62 (MQ=255)
            aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcgc                                                              <  1:616185/69‑1 (MQ=255)
             tAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgcc                                                           >  1:2377341/1‑71 (MQ=255)
              aaCCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                                                   >  1:2717271/1‑62 (MQ=255)
               aCCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                                                         <  1:342431/71‑1 (MQ=255)
                  aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCtt                                                       >  1:3384175/1‑70 (MQ=255)
                  aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc                                                      <  1:2644862/71‑1 (MQ=255)
                  aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc                                                      <  1:3313112/71‑1 (MQ=255)
                  aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc                                                      <  1:656111/71‑1 (MQ=255)
                   gaggCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCg                                                     <  1:487447/69‑1 (MQ=255)
                   gCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:2294722/1‑53 (MQ=255)
                    cGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    >  1:3607859/1‑55 (MQ=255)
                       cagcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga                                                 <  1:1594923/71‑1 (MQ=255)
                        agcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCAt                                                                  >  1:3310816/1‑53 (MQ=255)
                         gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                                                   <  1:111352/51‑1 (MQ=255)
                         gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt                                               >  1:3735974/1‑71 (MQ=255)
                         gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt                                               <  1:2357021/71‑1 (MQ=255)
                           aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                                             >  1:15852/1‑71 (MQ=255)
                           aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                                             >  1:899600/1‑71 (MQ=255)
                           aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                                             <  1:794632/71‑1 (MQ=255)
                             gCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCct                                           >  1:1602426/1‑71 (MQ=255)
                               cGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCctct                                         <  1:821019/71‑1 (MQ=255)
                               cGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCctct                                         <  1:2207389/71‑1 (MQ=255)
                                gagaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:601585/1‑40 (MQ=255)
                                gagaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:818340/1‑40 (MQ=255)
                                gagaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                                                       >  1:2787207/1‑40 (MQ=255)
                                 agaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:3571546/42‑1 (MQ=255)
                                 agaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                                                    <  1:2082660/42‑1 (MQ=255)
                                  gaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATc                                      >  1:1749552/1‑71 (MQ=255)
                                   aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                                              <  1:2553341/46‑2 (MQ=9)
                                   aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                                              <  1:252030/46‑2 (MQ=9)
                                   aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                                              <  1:2101599/46‑2 (MQ=9)
                                   aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                                              <  1:1654566/46‑2 (MQ=9)
                                   aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                                                         >  1:3178103/1‑51 (MQ=255)
                                   aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                                                         >  1:1976319/1‑51 (MQ=255)
                                   aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                                                         >  1:3843247/1‑51 (MQ=255)
                                     tAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACg                                                                     <  1:788433/37‑1 (MQ=9)
                                     tAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACg                                                                     <  1:693222/37‑1 (MQ=9)
                                      aaaGGCGAGAGAGGTCATCAGGATAGGACGCAGACg                                                                     <  1:118483/36‑1 (MQ=9)
                                           ccAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga                                                 >  1:4035219/1‑51 (MQ=255)
                                           ccAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgat                                                   <  1:915931/49‑1 (MQ=255)
                                           ccAGAGAGGTCATCAGGATTGGGCGTAAACGCATCAGCGCCGCTTCgatga                                                 >  1:2185675/1‑51 (MQ=255)
                                            cAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTcc                            >  1:1346506/1‑71 (MQ=255)
                                                 agGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga                                                 <  1:3615588/45‑1 (MQ=255)
                                                    tcatcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga                                                 <  1:2204305/42‑1 (MQ=255)
                                                     cgtcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                                             <  1:1547348/43‑1 (MQ=255)
                                                     catcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTcc                            >  1:1755320/1‑62 (MQ=255)
                                                     catcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTcc                            >  1:965729/1‑62 (MQ=255)
                                                     catcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCttcttt                     >  1:1995259/1‑69 (MQ=255)
                                                      atcaGGATTGGGCGTAAACGCATCCGCGCCGCTTc                                                      <  1:1341294/35‑1 (MQ=255)
                                                      atcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt                                               <  1:811881/42‑1 (MQ=255)
                                                      atcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTcc                            <  1:2948180/61‑1 (MQ=255)
                                                       tcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCga                                                    <  1:4039603/36‑1 (MQ=255)
                                                       tcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCg                                                     <  1:923050/35‑1 (MQ=255)
                                                       tcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                                             <  1:295196/43‑1 (MQ=255)
                                                       tcaGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATc                                      <  1:1580106/50‑1 (MQ=255)
                                                        caGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatgat                                                >  1:3615218/1‑39 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCg                                     >  1:2782770/1‑49 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCTGCGTTTTCCcttc                        >  1:2304155/1‑62 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTccc                           >  1:2646466/1‑59 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTccc                           >  1:2180471/1‑59 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCttctttct                   >  1:733781/1‑67 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcat               <  1:1719309/71‑1 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcat               <  1:1967454/71‑1 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGc                 >  1:3687074/1‑69 (MQ=255)
                                                         aGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGc                 >  1:491991/1‑69 (MQ=255)
                                                          ggATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcatc              <  1:3633268/71‑1 (MQ=255)
                                                          ggATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGca                >  1:2501980/1‑69 (MQ=255)
                                                           gATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCcttct                       >  1:1474700/1‑61 (MQ=255)
                                                           gATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCcttct                       >  1:3139256/1‑61 (MQ=255)
                                                                ggCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCttctttct                   >  1:3719333/1‑60 (MQ=255)
                                                                ggCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGCATCATCTCa        <  1:383683/71‑1 (MQ=255)
                                                                    tAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCct                          <  1:59303/49‑1 (MQ=255)
                                                                    tAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcatc              <  1:3743386/61‑1 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCttctttc                    >  1:2577034/1‑53 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTg                  >  1:2772580/1‑55 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcatca             >  1:2219211/1‑60 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGcatca             >  1:1136793/1‑60 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGCATCATCTCAACGGc   >  1:2415120/1‑70 (MQ=255)
                                                                      aaCGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGCATCATCTCAACGGCa  <  1:2893217/71‑1 (MQ=255)
                                                                      |                                                                      
GCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCCCTTCTTTCTGCATCATCTCAACGGCA  >  minE/2584299‑2584439

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: