Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,584,343 C→G 9.5% L978L (CTG→CTC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,584,3430CG 9.5% 352.9 / 20.1 116L978L (CTG→CTCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base C (42/63);  new base G (3/8);  total (45/71)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 120 total aligned reads displayed.

CCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCC  >  minE/2584276‑2584413
                                                                   |                                                                      
ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     <  1:3681596/71‑1 (MQ=255)
ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     <  1:3885301/71‑1 (MQ=255)
ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     <  1:1428751/71‑1 (MQ=255)
ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                                       >  1:2490451/1‑69 (MQ=255)
ccGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCa                                                                       <  1:3667066/69‑1 (MQ=255)
  gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagag                                                                    <  1:1834526/70‑1 (MQ=255)
  gcgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     >  1:2647922/1‑69 (MQ=255)
   cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                                                  >  1:1297359/1‑71 (MQ=255)
   cgTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                                                  >  1:4116409/1‑71 (MQ=255)
    gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                                                 >  1:492515/1‑71 (MQ=255)
    gTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAgg                                                                 >  1:3922321/1‑71 (MQ=255)
       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     >  1:2634427/1‑64 (MQ=255)
       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                                              >  1:3018587/1‑71 (MQ=255)
       ttGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtc                                                               <  1:2590390/70‑1 (MQ=255)
       ctGCGCGCCACTGCCGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGtc                                                               >  1:3042924/2‑70 (MQ=255)
        tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                                              <  1:3297886/70‑1 (MQ=255)
        tGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                                              <  1:3691148/70‑1 (MQ=255)
         gcgcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtcat                                                             <  1:3285327/70‑1 (MQ=255)
           gcgcCAGAACCGGCACCATGACTTATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                                          <  1:73211/71‑1 (MQ=255)
           gcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                                          <  1:487228/71‑1 (MQ=255)
           gcgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAg                                                          <  1:3677551/71‑1 (MQ=255)
            cgcCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                                         <  1:2766735/71‑1 (MQ=255)
               cAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                                       <  1:902317/70‑1 (MQ=255)
                 gAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTgg                                                    <  1:828765/71‑1 (MQ=255)
                   aCCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGc                                                  >  1:1064371/1‑71 (MQ=255)
                    ccGGCACCGTTACTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATCAGGATAGGACg                                                 <  1:1490653/71‑1 (MQ=9)
                    ccGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCg                                                 >  1:119297/1‑71 (MQ=255)
                     cGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAgg                                                         <  1:3487752/62‑1 (MQ=255)
                      ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                                                >  1:3151704/1‑70 (MQ=255)
                      ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                                               <  1:1305174/71‑1 (MQ=255)
                      ggCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                                               <  1:3973594/71‑1 (MQ=255)
                       gCACCATGACTGATAACCAGCGTCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaa                                              <  1:1738442/71‑1 (MQ=255)
                         aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                             <  1:3193860/70‑1 (MQ=255)
                         aCCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAAc                                            <  1:2938334/71‑1 (MQ=255)
                          ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                                                <  1:1491381/66‑1 (MQ=255)
                          ccATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                                               <  1:325839/67‑1 (MQ=255)
                           cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGt                                                                <  1:2112082/49‑1 (MQ=255)
                           cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGt                                                                <  1:4049633/49‑1 (MQ=255)
                           cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                          <  1:1497753/71‑1 (MQ=255)
                           cATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                          <  1:460122/71‑1 (MQ=255)
                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                                       >  1:1367422/1‑57 (MQ=255)
                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                                       >  1:2989982/1‑57 (MQ=255)
                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTggg                                                   >  1:4004351/1‑61 (MQ=255)
                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                         >  1:2226244/1‑71 (MQ=255)
                            aTGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                         >  1:3108791/1‑71 (MQ=255)
                             tGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGtca                                                              <  1:1936658/49‑1 (MQ=255)
                             tGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                          <  1:3213435/69‑1 (MQ=255)
                               aCTGATAGCTAGCGGTAATACGCCGAGAATAAAGGCGAGAGAGGTCATAGGATAGGACGCAGACGCATAcg                                     >  1:1001754/1‑71 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCagagag                                                                  >  1:3195361/1‑43 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGa                                                        >  1:4127276/1‑53 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                                       >  1:155850/1‑54 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                             >  1:2387376/1‑64 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                             >  1:27258/1‑64 (MQ=255)
                               aCTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATcc                                      >  1:102082/1‑71 (MQ=255)
                                 tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAt                                                       <  1:2689535/52‑1 (MQ=255)
                                 tGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCcg                                     <  1:3504817/70‑1 (MQ=255)
                                  gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCaga                                                                     <  1:1826979/37‑1 (MQ=255)
                                  gATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGCGTAAACGCATCcgc                                    <  1:2455645/69‑1 (MQ=255)
                                   aTAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGAtt                                                      >  1:1133546/1‑51 (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:3313112/71‑1 (MQ=255)
                                         aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc                            <  1:656111/71‑1 (MQ=255)
                                         aGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTc                            <  1:2644862/71‑1 (MQ=255)
                                          gaggCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCg                           <  1:487447/69‑1 (MQ=255)
                                          gCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                             >  1:2294722/1‑53 (MQ=255)
                                              cagcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCgatga                       <  1:1594923/71‑1 (MQ=255)
                                               agcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCAt                                        >  1:3310816/1‑53 (MQ=255)
                                                gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTa                                               <  1:2500439/45‑1 (MQ=255)
                                                gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCa                                         <  1:111352/51‑1 (MQ=255)
                                                gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt                     <  1:2357021/71‑1 (MQ=255)
                                                gcaCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGAtt                     >  1:3735974/1‑71 (MQ=255)
                                                 caCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGt                                                <  1:2186735/43‑1 (MQ=255)
                                                  aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                   >  1:15852/1‑71 (MQ=255)
                                                  aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                   <  1:794632/71‑1 (MQ=255)
                                                  aCGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGc                   >  1:899600/1‑71 (MQ=255)
                                                    gCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCct                 >  1:1602426/1‑71 (MQ=255)
                                                      cGAGAATAAAGGCCAGAGAGGTCATCAGGATTggg                                                   <  1:2754577/35‑1 (MQ=255)
                                                      cGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCctct               <  1:2207389/71‑1 (MQ=255)
                                                      cGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCctct               <  1:821019/71‑1 (MQ=255)
                                                       gagaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTaaa                                             >  1:601585/1‑40 (MQ=255)
                                                        agaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                          <  1:2082660/42‑1 (MQ=255)
                                                        agaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGc                                          <  1:3571546/42‑1 (MQ=255)
                                                         gaATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATc            >  1:1749552/1‑71 (MQ=255)
                                                          aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                    <  1:1654566/46‑2 (MQ=9)
                                                          aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                    <  1:2101599/46‑2 (MQ=9)
                                                          aaTAAAGGCGAGAGAGGTCATCAGGATAGGACGCAGACGCATAcgt                                    <  1:2553341/46‑2 (MQ=9)
                                                          aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                               >  1:1976319/1‑51 (MQ=255)
                                                          aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                               >  1:3843247/1‑51 (MQ=255)
                                                          aaTAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGcgccgc                               >  1:3178103/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)
                                                                   cAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTcc  >  1:1346506/1‑71 (MQ=255)
                                                                   |                                                                      
CCGCGTTTTGCGCGCCAGAACCGGCACCATGACTGATAACCAGCGGCAGCACGCCGAGAATAAAGGCCAGAGAGGTCATCAGGATTGGGCGTAAACGCATCCGCGCCGCTTCGATGATTGCCTCTATCGGCGTTTTCC  >  minE/2584276‑2584413

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: