Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,786 A→G 9.2% V1029V (GTT→GTC) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,7860AG 9.2% 144.5 / 5.2 54V1029V (GTT→GTC) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (19/30);  new base G (2/3);  total (21/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 = 5.28e-01

GTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTA  >  minE/309717‑309853
                                                                     |                                                                   
gTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                                    >  1:3241305/1‑71 (MQ=255)
 gtttATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAac                                                                   >  1:1055826/4‑71 (MQ=255)
        aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                            <  1:938070/71‑1 (MQ=255)
        aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                            <  1:1622046/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)
        aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCCGCGGCGa                                                                     >  1:3959701/1‑62 (MQ=255)
          tgatgaTCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                            >  1:5110062/1‑69 (MQ=255)
            atgatCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                         >  1:2303031/1‑70 (MQ=255)
                tCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                    >  1:2835169/1‑71 (MQ=255)
                 cGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGaa                                                                    <  1:4874872/54‑1 (MQ=255)
                     aGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGa                                                                     <  1:3586751/49‑1 (MQ=255)
                     aGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGa                                                                     <  1:3452607/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:4044873/1‑40 (MQ=255)
                                  tCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca                                                                 >  1:5692974/1‑40 (MQ=255)
                                  tCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAacca                                                                 >  1:5370036/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:1254904/70‑1 (MQ=255)
                                                        tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt             <  1:701302/70‑1 (MQ=255)
                                                        tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt             <  1:2154671/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:3879147/69‑1 (MQ=255)
                                                          aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATc            <  1:2545637/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)
                                                              aacggcggATCACCACAAAGAACACCGGTACGAAGAAGATTg                                   >  1:956319/3‑42 (MQ=18)
                                                                cggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                      <  1:4354259/53‑1 (MQ=255)
                                                                cggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                      <  1:2851986/53‑1 (MQ=255)
                                                                cggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                      <  1:3891948/53‑1 (MQ=255)
                                                                  gcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTAGCg                          <  1:1792654/47‑1 (MQ=255)
                                                                   cgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTa  <  1:4164197/70‑1 (MQ=255)
                                                                     aaCCACCACAAAGAATACGGAACGAAGAAGATTGCCAGTAc                            >  1:3598405/1‑41 (MQ=38)
                                                                     |                                                                   
GTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTA  >  minE/309717‑309853

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: