Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 309,788 C→T 11.5% V1029I (GTT→ATT) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE309,7880CT11.5% 137.2 / 7.3 52V1029I (GTT→ATT) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (16/30);  new base T (3/3);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACC  >  minE/309718‑309857
                                                                      |                                                                     
gtttATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAac                                                                       >  1:1055826/4‑71 (MQ=255)
       aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                <  1:372150/71‑1 (MQ=255)
       aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                <  1:1622046/71‑1 (MQ=255)
       aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                <  1:3270433/71‑1 (MQ=255)
       aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                <  1:3263763/71‑1 (MQ=255)
       aaTGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                <  1:938070/71‑1 (MQ=255)
         tgatgaTCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACcaca                                                                >  1:5110062/1‑69 (MQ=255)
           atgatCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAg                                                             >  1:2303031/1‑70 (MQ=255)
               tCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATAc                                                        >  1:2835169/1‑71 (MQ=255)
                          gCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaa                                               >  1:4735029/1‑69 (MQ=255)
                              tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                              <  1:3756945/50‑1 (MQ=255)
                              tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                              <  1:417756/50‑1 (MQ=255)
                              tgCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACaaa                                                              <  1:2417350/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:1446701/2‑71 (MQ=16)
                                                      ttAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt                 >  1:221571/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:701302/70‑1 (MQ=255)
                                                       tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt                 <  1:2154671/70‑1 (MQ=255)
                                                       tAAAGCAACGGCGGATCACCACAAAGAACACCGGTACGAAGAAGATTGCCAGCAACGTTGCAGAGACCAt                 <  1:1254904/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:3234260/69‑1 (MQ=255)
                                                         aaGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATc                <  1:2545637/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:3891948/53‑1 (MQ=255)
                                                               cggcgAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc                          <  1:2851986/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)
                                                                     accaccACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACAcc  >  1:4442604/1‑71 (MQ=255)
                                                                      ccaccaCAAAGAATACCGGAACGAAGAAGATTGCCAg                                   <  1:2936979/37‑1 (MQ=255)
                                                                      ccaccaCAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCa         >  1:5409611/1‑63 (MQ=255)
                                                                      |                                                                     
TTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCCCGCCCATTACACC  >  minE/309718‑309857

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: