Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,850 A→T 7.0% V341V (GTT→GTA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,8500AT 7.0% 290.6 / 11.2 100V341V (GTT→GTAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (40/53);  new base T (2/5);  total (42/58)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 103 total aligned reads displayed.

CGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGA  >  minE/311781‑311918
                                                                     |                                                                    
cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                                     >  1:706391/1‑71 (MQ=255)
cgcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                                     >  1:3432227/1‑71 (MQ=255)
 gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTaa                                                                     >  1:3644790/1‑70 (MQ=255)
 gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                                    <  1:2472535/71‑1 (MQ=255)
 gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                                    >  1:2901924/1‑71 (MQ=255)
 gcgcgGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                                    >  1:3917369/1‑71 (MQ=255)
    cgGAAGTTCTGCAGGAACAGATACATAACCCGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                                 <  1:3786173/71‑1 (MQ=255)
     ggAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACt                                                                >  1:2842000/1‑71 (MQ=255)
       atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                                                              <  1:4132534/68‑1 (MQ=16)
       atattCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTc                                                              <  1:1257570/68‑1 (MQ=16)
         attCTGCAAGAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGt                                                            >  1:3806065/2‑71 (MQ=18)
          ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                           <  1:3151805/71‑1 (MQ=255)
          ttCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                           <  1:3051962/71‑1 (MQ=255)
            cTGCAGGAACAGAAACATAACCAGGAACACGAGGAAGATCGCTTCGACCAGCGTTTTAACCACt                                                                >  1:2082536/1‑64 (MQ=255)
               cAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                         >  1:3926964/1‑68 (MQ=255)
                aagAACAGATACATCACCAGGAACACCAGCATAATGGCTTCGAACAGCGTTTTTACCACTTCGTGaa                                                         <  1:3641581/65‑1 (MQ=21)
                 ggAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                                    <  1:3727808/71‑1 (MQ=255)
                     cAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTAACCACTTCGTGAATAGAGAtttt                                               <  1:282269/71‑1 (MQ=255)
                      aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                               <  1:1906737/71‑1 (MQ=255)
                      aGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAttt                                                >  1:3224490/1‑70 (MQ=255)
                       gATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTc                                              <  1:1770351/71‑1 (MQ=255)
                        aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                                       >  1:3316468/1‑61 (MQ=255)
                        aTACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                               <  1:360248/69‑1 (MQ=255)
                         tACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                            <  1:3823302/71‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                              <  1:1019695/52‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                              <  1:2870579/52‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTc                                                              <  1:2960730/52‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                                            <  1:3108121/54‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGt                                                            <  1:683765/54‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                           <  1:2335233/55‑1 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                               >  1:2972066/1‑67 (MQ=255)
                          aCATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                            >  1:4121586/1‑70 (MQ=255)
                            aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTg                                                           <  1:2657994/53‑1 (MQ=255)
                            aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                         >  1:66945/1‑71 (MQ=255)
                            aTAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                         >  1:442460/1‑71 (MQ=255)
                             tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                                        >  1:2195688/1‑71 (MQ=255)
                             tAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAAc                                        >  1:2112705/1‑71 (MQ=255)
                              aaCCAGGAACACGTGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                                       >  1:1061898/1‑71 (MQ=255)
                               aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                                    >  1:217126/1‑41 (MQ=255)
                               aCCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGAtttt                                               >  1:3868383/1‑62 (MQ=255)
                                 cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAAc                                                                    <  1:618998/39‑1 (MQ=255)
                                 cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATa                                                       >  1:2920808/1‑52 (MQ=255)
                                 cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                           <  1:1361033/64‑1 (MQ=255)
                                 cAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                                    <  1:2200701/71‑1 (MQ=255)
                                    gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATagag                                                    <  1:86175/52‑1 (MQ=255)
                                    gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCAc                                            >  1:2062484/1‑60 (MQ=255)
                                    gAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtgg                                 >  1:2425559/1‑71 (MQ=255)
                                     aaCACGAGGATGATCGCTTCGACCAGCGTTTTAACTACTTCGTGAATAGAGATTTTCACGAACGTCgtgg                                 <  1:1653084/70‑1 (MQ=255)
                                      acacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                               <  1:2719486/71‑1 (MQ=255)
                                       cacGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                                 >  1:2161213/1‑36 (MQ=255)
                                       cacGAGGATGATCGCTTCGACCAGCGTTTTAACCAc                                                                 >  1:2148393/1‑36 (MQ=255)
                                       cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                         >  1:3436375/1‑44 (MQ=255)
                                       cacGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggtg                               <  1:770770/70‑1 (MQ=255)
                                         cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGaa                                                         <  1:2813876/42‑1 (MQ=255)
                                         cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                            <  1:1569925/71‑1 (MQ=255)
                                         cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                            <  1:3546563/71‑1 (MQ=255)
                                         cGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                            <  1:3830764/71‑1 (MQ=255)
                                          gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCg                                    <  1:3222918/62‑1 (MQ=255)
                                          gAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCg                            <  1:2843283/70‑1 (MQ=255)
                                           aGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTa                          <  1:3447892/71‑1 (MQ=255)
                                             gatgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATag                                                      <  1:1701625/41‑1 (MQ=255)
                                               tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGaa                                         >  1:1592630/1‑52 (MQ=255)
                                               tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACg                                       <  1:3144583/54‑1 (MQ=255)
                                               tgatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATggg                      >  1:2882058/1‑71 (MQ=255)
                                                gatCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                  <  1:2270954/58‑1 (MQ=255)
                                                 atggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGt                     <  1:2409392/67‑1 (MQ=255)
                                                 atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                  <  1:2593067/57‑1 (MQ=255)
                                                 atCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGt                     >  1:72631/1‑70 (MQ=255)
                                                   ggCTTCGAACAGCGTTTTTACCACTTCGTGAATAGAAAGCTGGACGAATGGCGTGGTGTCATAAGGGTa                    >  1:1326359/2‑69 (MQ=255)
                                                   cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAgaga                                                   >  1:3200662/1‑38 (MQ=255)
                                                   cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTGACGAACGGCgtggt                                >  1:3754235/1‑57 (MQ=255)
                                                   cGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                >  1:2418844/1‑57 (MQ=255)
                                                    gCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                   <  1:2217799/53‑1 (MQ=255)
                                                    gCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                   <  1:3901708/53‑1 (MQ=255)
                                                     cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACg                                           <  1:474611/44‑1 (MQ=255)
                                                     cTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                  <  1:1477670/53‑1 (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)
                                                             aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtg                                  <  1:1484704/45‑1 (MQ=255)
                                                             aGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTaaa                  >  1:3857942/1‑61 (MQ=255)
                                                              gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                   <  1:327933/43‑1 (MQ=255)
                                                              gCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCgt                                   <  1:1639191/43‑1 (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)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                >  1:3785071/1‑41 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCgtggt                                >  1:4012222/1‑41 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTc          >  1:2098151/1‑63 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa  >  1:463632/1‑71 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa  >  1:187607/1‑71 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa  >  1:2213810/1‑71 (MQ=255)
                                                                   ttAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGa  >  1:2175611/1‑71 (MQ=255)
                                                                     |                                                                    
CGCGCGGAAGTTCTGCAGGAACAGATACATAACCAGGAACACGAGGATGATCGCTTCGACCAGCGTTTTAACCACTTCGTGAATAGAGATTTTCACGAACGGCGTGGTGTCGTATGGGTAAACAATTTTCAGACCCGA  >  minE/311781‑311918

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: