Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,895 T→A 6.4% P326P (CCA→CCT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8950TA 6.4% 169.6 / 2.1 60P326P (CCA→CCTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base A (2/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311828‑311962
                                                                   |                                                                   
tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg                                                                  >  1:2882058/1‑71 (MQ=255)
  atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGt                                                                 <  1:2409392/67‑1 (MQ=255)
  atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGt                                                                 >  1:72631/1‑70 (MQ=255)
    ggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                                                                >  1:1326359/2‑69 (MQ=255)
       ttCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTaaa                                                              >  1:437715/1‑68 (MQ=255)
         cGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAAGGCGTGGTGTCGTATGGGTAAACa                                                            >  1:2040389/1‑68 (MQ=255)
          gACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtt                                                         <  1:3340499/70‑1 (MQ=255)
           aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:1053984/71‑1 (MQ=255)
           aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:3171169/71‑1 (MQ=255)
           aCCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtt                                                         <  1:1015302/69‑1 (MQ=255)
             cAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCa                                                     <  1:2790105/71‑1 (MQ=255)
              aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTaaa                                                              >  1:3857942/1‑61 (MQ=255)
               gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:75179/71‑1 (MQ=255)
                cGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATgg                                                                   >  1:3838680/1‑54 (MQ=255)
                  ttttAACCACTTCGTGAATATAGATTTTCACTAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGa                                                   <  1:3202357/68‑1 (MQ=255)
                  ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAg                                                    <  1:1429467/67‑1 (MQ=255)
                   tttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:1285712/63‑1 (MQ=255)
                    ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc                                                      >  1:2098151/1‑63 (MQ=255)
                    ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                              >  1:2213810/1‑71 (MQ=255)
                    ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                              >  1:463632/1‑71 (MQ=255)
                    ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                              >  1:2175611/1‑71 (MQ=255)
                    ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                              >  1:187607/1‑71 (MQ=255)
                          aCTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGAc                                                  <  1:420580/61‑1 (MQ=255)
                           cTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                                       <  1:646312/71‑1 (MQ=255)
                           cTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTAGAGAACCTTCATTCCCTGCGGgaag                                       <  1:2744248/71‑1 (MQ=255)
                            ttCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:888904/54‑1 (MQ=255)
                             tCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                                           >  1:651076/1‑49 (MQ=255)
                             tCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTGGAGAACCTTCATTCCCTGCGGgaagaa                                     >  1:3759953/1‑71 (MQ=255)
                              cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc                                                      >  1:1099397/1‑53 (MQ=255)
                              cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaag                                       <  1:3793707/68‑1 (MQ=255)
                              cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAAc                                    >  1:110222/1‑71 (MQ=255)
                                taaaTAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAAc                                    <  1:1177759/67‑1 (MQ=255)
                                 tattAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACgg                                  <  1:4104252/67‑1 (MQ=255)
                                  aaTAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAt                                                          <  1:3500427/45‑1 (MQ=255)
                                  aaTAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGgaagaa                                     >  1:2007828/1‑66 (MQ=255)
                                    tAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAttt                                                        >  1:2761428/1‑45 (MQ=255)
                                    tAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTcc                              <  1:467296/71‑1 (MQ=255)
                                      gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                                           >  1:2458169/1‑40 (MQ=255)
                                      gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                                           >  1:3675489/1‑40 (MQ=255)
                                      gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCAt                            <  1:985052/71‑1 (MQ=255)
                                      gagaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCAt                            <  1:1438278/71‑1 (MQ=255)
                                       agaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:1985641/43‑1 (MQ=255)
                                       agaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                       <  1:2620911/43‑1 (MQ=255)
                                       agaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATc                           <  1:4053385/71‑1 (MQ=255)
                                          ttttCACGAACGGCGTGGTGTCGTATGGGTAAAAAtttt                                                       >  1:401097/1‑39 (MQ=25)
                                              cACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGcc                     <  1:222910/70‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAg                   <  1:3123314/71‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAg                   <  1:2897831/71‑1 (MQ=255)
                                               aCGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAg                   >  1:2083655/1‑71 (MQ=255)
                                                cGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGtt                                >  1:3323121/1‑57 (MQ=255)
                                                        gtggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACg          <  1:2385732/71‑1 (MQ=255)
                                                        gtggtgTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATATTCGCCAGTTCAGCACg          <  1:2857905/71‑1 (MQ=255)
                                                         tggtgTCGTATGTGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACg          <  1:1825211/70‑1 (MQ=255)
                                                            tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCAt                            >  1:1002855/1‑49 (MQ=255)
                                                            tgtCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGAtt      >  1:1493013/1‑71 (MQ=255)
                                                              tCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGc                      <  1:4060135/53‑1 (MQ=255)
                                                               cGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTc                               <  1:2132783/43‑1 (MQ=255)
                                                                gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCa            <  1:3163284/61‑1 (MQ=255)
                                                                gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCa            >  1:3865418/1‑61 (MQ=255)
                                                                  aTGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCg  >  1:2562262/1‑69 (MQ=255)
                                                                   |                                                                   
TGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATTGCCG  >  minE/311828‑311962

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: