Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 894,163 C→T 49.0% A88T (GCG→ACG)  sapB ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE894,1630CT49.0% 7.5 / 144.4 102A88T (GCG→ACG) sapBpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base T (22/28);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 105 total aligned reads displayed.

AAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCG  >  minE/894096‑894231
                                                                   |                                                                    
aaGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTaa                                                                    >  1:4110294/1‑70 (MQ=255)
aaGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTaa                                                                    >  1:414031/1‑70 (MQ=255)
aaGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaa                                                                    <  1:2470509/70‑1 (MQ=255)
  gTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAAc                                                                 >  1:4247872/1‑71 (MQ=255)
   ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAAcc                                                                <  1:4000247/71‑1 (MQ=255)
   ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAAcc                                                                <  1:1922759/71‑1 (MQ=255)
   ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAAcc                                                                <  1:2489794/71‑1 (MQ=255)
   ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAAcc                                                                <  1:191499/71‑1 (MQ=255)
   ttATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCCCatcg                                                                  <  1:2785539/69‑4 (MQ=255)
    tATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAAcc                                                                >  1:3398239/1‑70 (MQ=255)
    tATCCTGCCACTTATGGCGCGTAAACCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCCCatcg                                                                  <  1:3758185/68‑4 (MQ=255)
     aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                               <  1:3525114/70‑1 (MQ=255)
     aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                               <  1:3530522/70‑1 (MQ=255)
     aTCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                               <  1:2336689/70‑1 (MQ=255)
       ccTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAg                                                            >  1:422845/1‑71 (MQ=255)
           ccACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                               <  1:720430/64‑1 (MQ=255)
           ccACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGc                                                          <  1:907327/69‑1 (MQ=255)
            cACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAg                                                       <  1:4320638/71‑1 (MQ=255)
             aCTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAg                                                       <  1:1040865/70‑1 (MQ=255)
              cTTATGGCGCGTAATCCCTGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGa                                                     <  1:2479451/71‑1 (MQ=255)
              cTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAg                                                       <  1:3874820/69‑1 (MQ=255)
              cTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGa                                                     <  1:1816766/71‑1 (MQ=255)
               ttATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGTGAGGAt                                                    >  1:4171861/1‑71 (MQ=255)
                tATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATg                                                   >  1:4260959/1‑71 (MQ=255)
                 aTGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATg                                                   >  1:3442116/1‑70 (MQ=255)
                    gcgcGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGcaca                                               <  1:1372980/71‑1 (MQ=255)
                    gcgcGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGcaca                                               <  1:11934/71‑1 (MQ=255)
                     cgcgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAg                                              <  1:820275/71‑1 (MQ=255)
                       cgTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGtt                                            <  1:1830411/71‑1 (MQ=255)
                         tAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTcc                                          <  1:3647064/71‑1 (MQ=255)
                         tAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTcc                                          <  1:948388/71‑1 (MQ=255)
                           aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAt                                        >  1:4205774/1‑71 (MQ=255)
                           aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAt                                        >  1:3860693/1‑71 (MQ=255)
                           aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAt                                        >  1:172556/1‑71 (MQ=255)
                           aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCg                                                               <  1:4634210/48‑1 (MQ=255)
                           aTCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAg                                                       >  1:1199477/1‑56 (MQ=255)
                            tCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAtc                                       <  1:4284635/71‑1 (MQ=255)
                              ccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAtcgt                                     >  1:2559257/1‑71 (MQ=255)
                              ccGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcg                                      <  1:2687230/70‑1 (MQ=255)
                               cGGCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGaa                                                             <  1:3229486/46‑1 (MQ=255)
                                ggCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcgtcg                                   >  1:1754645/1‑71 (MQ=255)
                                 gCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAAcc                                                                <  1:3865571/41‑1 (MQ=255)
                                 gCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAg                                                       <  1:3854539/50‑1 (MQ=255)
                                 gCAATCATCCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAg                                                       <  1:114437/50‑1 (MQ=255)
                                   aaTCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCg                                <  1:4175406/71‑1 (MQ=255)
                                   aaTCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCg                                >  1:3329252/1‑71 (MQ=255)
                                    atcatcCCCACGGGGATGCCAACAATCAGCGTAAAACCg                                                               <  1:1091519/39‑1 (MQ=255)
                                    atcatcCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAg                                                       <  1:1900479/47‑1 (MQ=255)
                                    atcatcCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGca                                                 <  1:3891798/53‑1 (MQ=255)
                                      catcCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCggg                              >  1:3735932/1‑70 (MQ=255)
                                      catcCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCg                                >  1:4416975/1‑68 (MQ=255)
                                       atcCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCggg                              >  1:1346054/1‑69 (MQ=255)
                                        tcCCCACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGcc                                 <  1:4236049/65‑1 (MQ=255)
                                         ccccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                          >  1:3268836/1‑71 (MQ=255)
                                         ccccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                          >  1:1057593/1‑71 (MQ=255)
                                          cccACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                            >  1:3887829/1‑68 (MQ=255)
                                          cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATg                                                   >  1:2873538/1‑45 (MQ=255)
                                          cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtc                                       >  1:594414/1‑57 (MQ=255)
                                          cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                          >  1:1045520/1‑70 (MQ=255)
                                          cccACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGAc                         >  1:1230709/1‑71 (MQ=255)
                                            cACGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                            <  1:2455753/66‑1 (MQ=255)
                                             aCGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAg                                              >  1:1539669/1‑47 (MQ=255)
                                             aCGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctc                      >  1:4437066/1‑71 (MQ=255)
                                             aCGGGGATGCCAACAATCAGAGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctc                      >  1:2858348/1‑71 (MQ=255)
                                              cGGGGATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctct                     >  1:4554262/1‑71 (MQ=255)
                                              cGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGc                                  >  1:1293607/1‑58 (MQ=255)
                                               ggggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAg                           <  1:2137241/64‑1 (MQ=255)
                                               ggggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtt                    >  1:2237324/1‑71 (MQ=255)
                                                gggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAtcgt                                     >  1:4600159/1‑53 (MQ=255)
                                                gggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAtcgt                                     >  1:2679478/1‑53 (MQ=255)
                                                gggATGCCAACAATCAGCGCAAAACCGAAGGCGAg                                                       <  1:2934539/35‑1 (MQ=255)
                                                gggATGCCAACAATCAGCGCAAAACCGAAGGCGAg                                                       <  1:4028823/35‑1 (MQ=255)
                                                gggATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtt                    >  1:2992384/1‑70 (MQ=255)
                                                 ggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                  <  1:2710929/71‑1 (MQ=255)
                                                 ggATGCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                  <  1:2308541/71‑1 (MQ=255)
                                                 ggATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                  <  1:2082085/71‑1 (MQ=255)
                                                   aTGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTaa                <  1:2550779/71‑1 (MQ=255)
                                                    tGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAAc               <  1:2456857/71‑1 (MQ=255)
                                                     gCCAACAATCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt              >  1:2472086/1‑71 (MQ=255)
                                                     gCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt              >  1:3606880/1‑71 (MQ=255)
                                                        aacaaTCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCggg                              <  1:4142661/52‑1 (MQ=255)
                                                        aacaaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAAc               >  1:3373568/1‑67 (MQ=255)
                                                         acaaTCAGCGTAAAACCGAAGGCGAGGATGCACAGTTcc                                          <  1:1791623/39‑1 (MQ=255)
                                                         acaaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACt              >  1:1155931/1‑67 (MQ=255)
                                                          caaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCg         >  1:1888281/1‑71 (MQ=255)
                                                          caaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCg         >  1:308267/1‑71 (MQ=255)
                                                           aaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                            >  1:4193730/1‑51 (MQ=255)
                                                           aaTCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCgg        >  1:445381/1‑71 (MQ=255)
                                                           aaTCAGCGCAAAACCGAAGGCGAAGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTaa                >  1:3636717/1‑63 (MQ=255)
                                                            aTCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGc       <  1:1641259/71‑1 (MQ=255)
                                                             tCAGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCtttt                  <  1:2935650/59‑1 (MQ=255)
                                                             tCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAg                           <  1:3050396/50‑1 (MQ=255)
                                                              cAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                            >  1:2479246/1‑48 (MQ=255)
                                                              cAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGaa                            >  1:2762161/1‑48 (MQ=255)
                                                               aGCGTAAAACCGAAGGCGAGGATGCACAGTTCCAt                                        >  1:3153213/1‑35 (MQ=255)
                                                               aGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCggg                              <  1:1543485/45‑1 (MQ=255)
                                                               aGCGTAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAAt    >  1:4130247/1‑71 (MQ=255)
                                                                gcgcAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATc   <  1:4300900/71‑1 (MQ=255)
                                                                gcgaAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACctct                     <  1:3955644/53‑1 (MQ=255)
                                                                 cgcAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCg  <  1:3048898/71‑1 (MQ=255)
                                                                   |                                                                    
AAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCG  >  minE/894096‑894231

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: