Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 311,562 C→T 20.7% Q437Q (CAG→CAA acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE311,5620CT20.7% 179.2 / 38.0 87Q437Q (CAG→CAAacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base T (7/11);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTT  >  minE/311493‑311630
                                                                     |                                                                    
gCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCt                                                                     <  1:2642158/71‑1 (MQ=255)
  cgaagaatGCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGc                                                                   <  1:593814/69‑1 (MQ=255)
   aaaGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGcc                                                                  >  1:1816424/1‑71 (MQ=255)
       aagGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTg                                                                    <  1:3090746/65‑1 (MQ=255)
         tgCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCg                                                            <  1:4761962/70‑1 (MQ=255)
         tgCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCg                                                            <  1:457399/70‑1 (MQ=255)
         tgCCATCGGAATAAATACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCg                                                            <  1:2156899/70‑1 (MQ=255)
         ggCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCg                                                            <  1:5518337/71‑1 (MQ=255)
         ggCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCg                                                            <  1:5272128/71‑1 (MQ=255)
            cATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCTTCGACt                                                         >  1:1866184/1‑71 (MQ=255)
             aTCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACt                                                         >  1:1423079/1‑70 (MQ=255)
              tCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                                                       >  1:3063381/1‑71 (MQ=255)
               cGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc                                                                      >  1:5411668/1‑55 (MQ=255)
               cGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATc                                                                      >  1:4732004/1‑55 (MQ=255)
                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                                                     >  1:4652197/1‑71 (MQ=255)
                ggTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCAACTTACg                                                     >  1:5107319/1‑71 (MQ=255)
                  tACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGcccc                                                                <  1:3768275/58‑1 (MQ=255)
                  tACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACgg                                                    <  1:3657130/70‑1 (MQ=255)
                   aCGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTAc                                                      <  1:775211/67‑1 (MQ=255)
                    cGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACggg                                                   >  1:370842/1‑69 (MQ=255)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGATTTTGcga                                                                 >  1:3464653/1‑51 (MQ=19)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGa                                                           >  1:4848865/1‑59 (MQ=255)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGa                                                           >  1:3019711/1‑59 (MQ=255)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGa                                                           >  1:3216126/1‑59 (MQ=255)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGa                                                           >  1:5636199/1‑59 (MQ=255)
                      aaTACCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGa                                                           >  1:1861955/1‑59 (MQ=255)
                       aTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACt                                                         <  1:3574320/60‑1 (MQ=255)
                       aTACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACt                                                         <  1:4281949/60‑1 (MQ=255)
                        tACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCt                                                                     <  1:2928484/47‑1 (MQ=255)
                        tACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCt                                                                     <  1:4381134/47‑1 (MQ=255)
                         aCCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGAttt                                                        <  1:3507797/59‑1 (MQ=255)
                         aCCGCTGACAGCACCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGAttt                                                        <  1:1629952/59‑1 (MQ=255)
                         aCCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttc                                            >  1:1685943/1‑71 (MQ=255)
                         aCCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttc                                            >  1:433074/1‑71 (MQ=255)
                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGc                                                                   >  1:2158641/1‑47 (MQ=255)
                          ccgccgACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttct                                           >  1:1320568/1‑71 (MQ=255)
                               gACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCt                                                                     <  1:4342965/40‑1 (MQ=255)
                                 cAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGcttc                                            >  1:3349190/1‑63 (MQ=255)
                                  aGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTggcgg                                   >  1:246316/1‑71 (MQ=255)
                                   gcaCCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGATTTTTCCGTCGCTTCTTTTggc                                     <  1:4178348/66‑1 (MQ=255)
                                     aCCATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGATTTTTCCGTCGCTTCTTTTggcgggag                                >  1:2750095/1‑68 (MQ=255)
                                      ccATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGAc                                                          >  1:5577918/1‑44 (MQ=255)
                                      ccATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGATTTTTCCGTCGCTTCttt                                         <  1:714095/61‑1 (MQ=255)
                                       cATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATc                                                             <  1:3814667/40‑1 (MQ=16)
                                       cATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGatt                                                         <  1:2529745/44‑3 (MQ=16)
                                       cATCGCGATACCCACCAGTGCGCCCTGAATTTGCGACATCGatt                                                         <  1:4122318/44‑3 (MQ=16)
                                        aTCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCaa                                <  1:3859758/68‑1 (MQ=255)
                                        aTCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAAcc                             <  1:1911963/71‑1 (MQ=255)
                                         tCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACt                                                         <  1:2434777/42‑1 (MQ=255)
                                           gcgATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTa                                                       >  1:1941244/1‑42 (MQ=255)
                                           gcgATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACcttc                          >  1:914193/1‑71 (MQ=255)
                                            cgATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                                                     <  1:2847855/43‑1 (MQ=255)
                                                aCCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                                                     <  1:799649/39‑1 (MQ=255)
                                                aCCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACg                                                     <  1:1161061/39‑1 (MQ=255)
                                                 ccGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCtttt                                        >  1:5651878/1‑51 (MQ=255)
                                                 ccGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCtttt                                        >  1:4673531/1‑51 (MQ=255)
                                                   gACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCa                  >  1:954246/1‑71 (MQ=255)
                                                    aCCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCttt                                         <  1:4650574/47‑1 (MQ=255)
                                                    aCCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCg                     <  1:2506615/67‑1 (MQ=255)
                                                    aCCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCAt                 <  1:3774947/71‑1 (MQ=255)
                                                    aCCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCAt                 <  1:72699/71‑1 (MQ=255)
                                                     ccAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCAt                 <  1:2099848/70‑1 (MQ=255)
                                                     ccAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCAt                 >  1:2593359/1‑70 (MQ=255)
                                                        gagCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTggcggc                                  <  1:3583454/50‑1 (MQ=255)
                                                        gagCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACcttct                         <  1:953740/59‑1 (MQ=255)
                                                          gcgcCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACg           >  1:2943218/1‑71 (MQ=255)
                                                           cgcCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGc          <  1:3676133/71‑1 (MQ=255)
                                                              ccTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAaca             <  1:111315/65‑1 (MQ=255)
                                                              ccTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAaca             <  1:1490332/65‑1 (MQ=255)
                                                               cTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCAt                 <  1:4912134/60‑1 (MQ=255)
                                                                tGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCaa                                <  1:508183/44‑1 (MQ=255)
                                                                tGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATaa               >  1:4982303/1‑61 (MQ=255)
                                                                tGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAAc     >  1:1377659/1‑71 (MQ=255)
                                                                tGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAAc     >  1:4760678/1‑71 (MQ=255)
                                                                tGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAAc     >  1:3517069/1‑71 (MQ=255)
                                                                 gAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACctt                           >  1:1940471/1‑48 (MQ=255)
                                                                 gAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAAc     >  1:1619275/1‑70 (MQ=255)
                                                                 gAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACg    >  1:3677104/1‑71 (MQ=255)
                                                                 gAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACg    >  1:4923046/1‑71 (MQ=255)
                                                                 gAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACg    >  1:340200/1‑71 (MQ=255)
                                                                 gAATCTGCACCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACg    >  1:1657118/1‑71 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACcttc                          >  1:324872/1‑48 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACcttc                          >  1:4306377/1‑48 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAAc     <  1:5082804/69‑1 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGt   >  1:1580229/1‑71 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGt   >  1:722958/1‑71 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGt   >  1:5678957/1‑71 (MQ=255)
                                                                  aaTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGt   <  1:2954064/71‑1 (MQ=255)
                                                                   aTCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGtt  <  1:5282395/71‑1 (MQ=255)
                                                                     |                                                                    
GCCAAAGAAGGCCATCGGTACGAATACCGCCGACAGTACCATCGCGATACCGACCAGAGCGCCCTGAATCTGCCCCATCGACTTACGGGTAGCTTCTTTTGGCGGCAAACCTTCTTCCGCCATAACACGCTCAACGTT  >  minE/311493‑311630

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: