Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,780 G→A 9.6% R1031R (CGC→CGT acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,7800GA 9.6% 144.0 / 5.1 52R1031R (CGC→CGTacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base A (2/3);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCC  >  minE/309710‑309848
                                                                      |                                                                    
aTTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                                                                      >  1:213003/1‑71 (MQ=255)
  tACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                                                                    >  1:4857399/1‑71 (MQ=255)
  tACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                                                                    >  1:1765905/1‑71 (MQ=255)
   acacGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg                                                                   >  1:5509266/1‑71 (MQ=255)
   acacGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcgg                                                                   >  1:5186827/1‑71 (MQ=255)
   acacGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcg                                                                    >  1:4655138/1‑70 (MQ=255)
       gTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                               >  1:3241305/1‑71 (MQ=255)
        gtttATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAac                                                              >  1:1055826/4‑71 (MQ=255)
             tCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                                                                     <  1:1343530/59‑1 (MQ=255)
             tCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                                                                      >  1:918399/1‑58 (MQ=255)
             tCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                                                                      >  1:722164/1‑58 (MQ=255)
             tCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcgg                                                                      >  1:4774689/1‑58 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       <  1:938070/71‑1 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       <  1:372150/71‑1 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       <  1:3270433/71‑1 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       <  1:3263763/71‑1 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       <  1:1622046/71‑1 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCCGCGGCGa                                                                >  1:3959701/1‑62 (MQ=255)
               aaTGATGATCGACAGTATGGCTGTGATCGATATCTTCATTCTTGCGGCTAAAgcggcgg                                                                   >  1:1550917/1‑59 (MQ=255)
                 tgatgaTCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggcggcg                                                                 <  1:3217467/59‑1 (MQ=255)
                 tgatgaTCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                       >  1:5110062/1‑69 (MQ=255)
                   atgatCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                    >  1:2303031/1‑70 (MQ=255)
                       tCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                               >  1:2835169/1‑71 (MQ=255)
                        cGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAgcggc                                                                     <  1:4057310/48‑1 (MQ=255)
                        cGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                               <  1:4874872/54‑1 (MQ=255)
                            aGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGa                                                                <  1:3452607/49‑1 (MQ=255)
                            aGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGa                                                                <  1:3586751/49‑1 (MQ=255)
                                  gCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                      >  1:4735029/1‑69 (MQ=255)
                                      tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                     <  1:417756/50‑1 (MQ=255)
                                      tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                     <  1:2417350/50‑1 (MQ=255)
                                      tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                     <  1:3756945/50‑1 (MQ=255)
                                      tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                    <  1:3083854/51‑1 (MQ=255)
                                       gCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                    <  1:3002248/50‑1 (MQ=255)
                                       gCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAtt                               >  1:3632136/1‑71 (MQ=255)
                                       gCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACAACAAAGAATACCGGAACgaag                                     >  1:5068710/1‑65 (MQ=255)
                                        cTCGATATCTTCATTCTTGCGGCTATAGCGGCGGCGAACCACCACAAAGAATAGCGGAACgtagaa                                   <  1:2436516/66‑1 (MQ=38)
                                         tCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca                                                            >  1:5370036/1‑40 (MQ=255)
                                         tCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca                                                            >  1:4044873/1‑40 (MQ=255)
                                         tCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca                                                            >  1:5692974/1‑40 (MQ=255)
                                            atatCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                    <  1:3299348/45‑1 (MQ=255)
                                               tCTTCATTCTTGCGGCTAAAGCGGCGGCGAACTACCACAAAGAATACCGGAAc                                         <  1:3705274/53‑1 (MQ=255)
                                                 ttCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                     >  1:3118460/1‑71 (MQ=255)
                                                  tCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGt                    >  1:3062001/1‑71 (MQ=255)
                                                   cATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                   >  1:726645/1‑71 (MQ=255)
                                                    aTTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                   <  1:5632875/70‑1 (MQ=255)
                                                    aTTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGtt                   <  1:5639346/70‑1 (MQ=255)
                                                       cTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGa                                                   >  1:1629771/1‑35 (MQ=255)
                                                       cTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAt                                >  1:2279661/1‑54 (MQ=255)
                                                          gCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGt              <  1:1037312/69‑1 (MQ=255)
                                                          gCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGt              <  1:4802089/69‑1 (MQ=255)
                                                              ttAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt        >  1:221571/2‑71 (MQ=16)
                                                              ttAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt        >  1:1446701/2‑71 (MQ=16)
                                                              cTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                     <  1:5569396/58‑1 (MQ=255)
                                                               tAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg                     <  1:3829365/57‑1 (MQ=255)
                                                               tAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg                  <  1:1550459/60‑1 (MQ=255)
                                                               tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt        <  1:2154671/70‑1 (MQ=255)
                                                               tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt        <  1:1254904/70‑1 (MQ=255)
                                                               tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt        <  1:701302/70‑1 (MQ=255)
                                                                aaaGCGGCGGCGAACCACCACAAAGAATACCGGAACgaagaa                                   <  1:611669/42‑1 (MQ=255)
                                                                aaaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATcc      >  1:3255651/1‑71 (MQ=255)
                                                                 aaGCGGCGGCGAACCACCACAAAGAATACCGGAACg                                        >  1:3643413/1‑36 (MQ=255)
                                                                 aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATc       <  1:2545637/69‑1 (MQ=255)
                                                                 aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATc       <  1:3879147/69‑1 (MQ=255)
                                                                 aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATc       <  1:3234260/69‑1 (MQ=255)
                                                                  aGCGGCGGCGAACCAACACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCg    >  1:3422099/1‑71 (MQ=255)
                                                                    cggcggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTCCAGTACCGTTGCGGTCACCATcc      >  1:4777737/1‑66 (MQ=255)
                                                                     ggcggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGcc  <  1:3265306/70‑1 (MQ=255)
                                                                      |                                                                    
ATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCC  >  minE/309710‑309848

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: