Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 313,780 C→T 6.2% A103A (GCG→GCA acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE313,7800CT 6.2% 253.2 / 4.7 81A103A (GCG→GCAacrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (3/2);  total (42/39)
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.85e-01

TCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAA  >  minE/313710‑313849
                                                                      |                                                                     
tCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTc                                                                       >  1:185794/1‑71 (MQ=255)
 cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      >  1:1357740/1‑71 (MQ=255)
 cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      >  1:3302138/1‑71 (MQ=255)
 cACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      >  1:3297553/1‑71 (MQ=255)
  aCCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGc                                                                     <  1:355407/71‑1 (MQ=255)
   ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      <  1:1414649/69‑1 (MQ=255)
   ccGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      <  1:1466063/69‑1 (MQ=255)
    cGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCt                                                                   >  1:2733143/1‑71 (MQ=255)
     gTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGAACACCTTTCGCACTGTCGTATGTCGCCTg                                                                  <  1:3397587/71‑1 (MQ=255)
      tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                                 <  1:2672726/71‑1 (MQ=255)
      tCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                                 <  1:842687/71‑1 (MQ=255)
        aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGTCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                                <  1:786936/70‑1 (MQ=255)
        aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGa                                                                 >  1:3890704/1‑69 (MQ=255)
        aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                                <  1:2552917/70‑1 (MQ=255)
        aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                                <  1:2664489/70‑1 (MQ=255)
        aaTTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGAt                                                                <  1:588666/70‑1 (MQ=255)
          ttGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAgg                                                             >  1:383311/1‑71 (MQ=255)
           tGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAg                                                              >  1:2139751/1‑69 (MQ=255)
                  tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAgg                                                      <  1:2881273/70‑1 (MQ=255)
                  tatTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGa                                                     >  1:3806933/1‑71 (MQ=255)
                    ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCg                                                                      <  1:3325655/52‑1 (MQ=255)
                    ttGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                   <  1:2134782/71‑1 (MQ=255)
                      ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTc                                                                       >  1:699980/1‑49 (MQ=255)
                      ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGt                                                            <  1:2684919/60‑1 (MQ=255)
                      ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                   >  1:21546/1‑69 (MQ=255)
                      ggCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                  <  1:1700155/70‑1 (MQ=255)
                      ggCTGAAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                   >  1:556766/1‑69 (MQ=255)
                        cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                <  1:2781175/70‑1 (MQ=255)
                        cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                <  1:3621205/70‑1 (MQ=255)
                        cTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                               >  1:3065537/1‑71 (MQ=255)
                           cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                               <  1:3575031/68‑1 (MQ=255)
                           cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCt                                              <  1:3543853/69‑1 (MQ=255)
                           cAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                            <  1:685667/71‑1 (MQ=255)
                            aGCCTGGGCTTTCGACAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                            >  1:3426661/1‑70 (MQ=255)
                             gCCTGGGTTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                          <  1:657635/71‑1 (MQ=255)
                             gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                            >  1:2569564/1‑69 (MQ=255)
                             gCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                            >  1:620324/1‑69 (MQ=255)
                                tGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATa                                          >  1:3057524/1‑68 (MQ=255)
                                 gggCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCa                                                  <  1:458518/59‑1 (MQ=255)
                                    ttttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGt                                                            <  1:2940451/45‑1 (MQ=255)
                                    cTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATag                                         <  1:1896816/65‑1 (MQ=255)
                                     tttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAg                                                       <  1:3934072/50‑1 (MQ=255)
                                      ttCGCCAGTTCGCCTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTGGTacagggac                                   <  1:1106380/69‑2 (MQ=255)
                                      ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAt                                                    >  1:2659880/1‑52 (MQ=255)
                                      ttCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagagag                                   >  1:2727577/1‑69 (MQ=255)
                                       tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATaga                                        >  1:1044064/1‑63 (MQ=255)
                                       tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                  <  1:2579346/69‑1 (MQ=255)
                                       tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                 <  1:1487958/70‑1 (MQ=255)
                                       tCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGac                                 <  1:3551790/70‑1 (MQ=255)
                                        cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                                           >  1:137180/1‑43 (MQ=255)
                                        cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTc                                                           >  1:2077633/1‑43 (MQ=255)
                                        cGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGaca                                <  1:3874872/70‑1 (MQ=255)
                                         gCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATc                                               >  1:1605930/1‑54 (MQ=255)
                                          ccAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGAt                                           <  1:198998/57‑1 (MQ=255)
                                           cAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTg                            <  1:3451331/71‑1 (MQ=255)
                                             gATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTg                                             >  1:1723289/1‑52 (MQ=255)
                                               tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                       >  1:2331124/1‑56 (MQ=255)
                                               tCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                        <  1:1729312/71‑1 (MQ=255)
                                                cACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                        >  1:241718/1‑70 (MQ=255)
                                                 aCCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCg                       <  1:848699/70‑1 (MQ=255)
                                                 aCCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGa                      <  1:3254427/71‑1 (MQ=255)
                                                  ccTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCtt                         >  1:992920/1‑67 (MQ=255)
                                                   cTTTCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGATCTg                                             >  1:911433/1‑46 (MQ=255)
                                                    tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                   >  1:3695322/1‑71 (MQ=255)
                                                    tttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                   >  1:3272140/1‑71 (MQ=255)
                                                     ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                >  1:3783257/1‑41 (MQ=255)
                                                     ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                >  1:920251/1‑41 (MQ=255)
                                                     ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                >  1:2717618/1‑41 (MQ=255)
                                                     ttCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAAt                                                >  1:2754584/1‑41 (MQ=255)
                                                     ttCGCAATGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTc                        >  1:3799638/1‑65 (MQ=255)
                                                      tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                >  1:2542452/1‑40 (MQ=255)
                                                      tCGCGCTGTCATAATTTGCCTGATAGGTCGCGGGATCGAt                                                >  1:3475039/1‑40 (MQ=255)
                                                            tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGa                                            >  1:552886/1‑38 (MQ=255)
                                                            tGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagag                                       <  1:957234/43‑1 (MQ=255)
                                                                gTATGTCGCCTGATAGGTCGCAGGATCAATCTGATagagaga                                    >  1:1043863/1‑42 (MQ=255)
                                                                 tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt      >  1:2322824/1‑71 (MQ=255)
                                                                 tATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt      <  1:1184875/71‑1 (MQ=255)
                                                                  aTGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCttt      >  1:2728984/1‑70 (MQ=255)
                                                                   tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCt               <  1:2553512/60‑1 (MQ=255)
                                                                   tGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTg     >  1:867376/1‑70 (MQ=255)
                                                                    gTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGc                           >  1:1868462/1‑47 (MQ=255)
                                                                     tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa  <  1:2958293/71‑1 (MQ=255)
                                                                     tCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGaaa  >  1:3076080/1‑71 (MQ=255)
                                                                      |                                                                     
TCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAA  >  minE/313710‑313849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: