Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,892 G→A 6.3% Y327Y (TAC→TAT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8920GA 6.3% 205.0 / 2.5 63Y327Y (TAC→TATacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/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

CGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATT  >  minE/311822‑311958
                                                                      |                                                                  
cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                    <  1:3546563/71‑1 (MQ=255)
cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                    <  1:3830764/71‑1 (MQ=255)
cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                    <  1:1569925/71‑1 (MQ=255)
 gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                    <  1:2843283/70‑1 (MQ=255)
  aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                  <  1:3447892/71‑1 (MQ=255)
      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)
                        ttttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                                                                    <  1:2675163/47‑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:463632/1‑71 (MQ=255)
                          ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                          >  1:2213810/1‑71 (MQ=255)
                          ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                          >  1:2175611/1‑71 (MQ=255)
                          ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa                                          >  1:187607/1‑71 (MQ=255)
                                aCTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                  >  1:2766196/1‑41 (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)
                                   tCGTGAATAGAGATTTTCACGATCGGCGTGGTGTCGTa                                                                  >  1:61539/1‑38 (MQ=255)
                                   tCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACaa                                                       >  1:651076/1‑49 (MQ=255)
                                   tCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTGGAGAACCTTCATTCCCTGCGGgaagaa                                 >  1:3759953/1‑71 (MQ=255)
                                    cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                  >  1:3795534/1‑37 (MQ=255)
                                    cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                  >  1:1043407/1‑37 (MQ=255)
                                    cGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                                                                  >  1:1077709/1‑37 (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:2620911/43‑1 (MQ=255)
                                             agaTTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAAtttt                                                   <  1:1985641/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:3865418/1‑61 (MQ=255)
                                                                      gTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCa        <  1:3163284/61‑1 (MQ=255)
                                                                      |                                                                  
CGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGACGGGAAGAACGGTTCCATCTTCGCCAGTTCAGCACGGATT  >  minE/311822‑311958

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: