Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,708 T→C 17.2% K55K (AAA→AAG) ‡ acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,7080TC17.2% 136.7 / 20.9 64K55K (AAA→AAG) ‡acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base C (2/9);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

ACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCC  >  minE/312638‑312777
                                                                      |                                                                     
aCATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGtt                                                                       <  1:3322075/71‑1 (MQ=255)
 cATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                      <  1:2185003/71‑1 (MQ=255)
 cATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                      <  1:3734500/71‑1 (MQ=255)
  aTGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGtttt                                                                     >  1:1843486/1‑71 (MQ=255)
     tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCg                                                                   <  1:3021689/70‑1 (MQ=255)
     tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCg                                                                   <  1:4256768/70‑1 (MQ=255)
     tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCg                                                                   <  1:4250009/70‑1 (MQ=255)
     tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGGCCTGCACTGtttt                                                                     <  1:2766564/68‑1 (MQ=255)
       cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                                <  1:1173121/71‑1 (MQ=255)
       cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                                <  1:2673768/71‑1 (MQ=255)
       cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                                >  1:35407/1‑71 (MQ=255)
       cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                                <  1:1243728/71‑1 (MQ=255)
         tCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATc                                                               <  1:4368823/70‑1 (MQ=37)
           aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAg                                                             <  1:3614084/70‑1 (MQ=255)
           aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAg                                                             <  1:3307584/70‑1 (MQ=255)
           aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAg                                                             <  1:5294707/70‑1 (MQ=255)
             gTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAgcgc                                                          <  1:3808628/71‑1 (MQ=255)
               tATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCg                                                        >  1:2702726/1‑71 (MQ=255)
               tATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCg                                                        >  1:5689768/1‑71 (MQ=255)
               tATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCg                                                        >  1:979089/1‑71 (MQ=255)
                  cGATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCa                                                                 >  1:248889/1‑59 (MQ=25)
                  agATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     >  1:3245035/2‑71 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     <  1:1161302/70‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     <  1:125038/70‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     <  1:5188206/70‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     <  1:269216/70‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     <  1:4475126/70‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGt                                                    <  1:388676/71‑1 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGTTCTGCGCATCAGCGCCcgga                                                     <  1:5384541/70‑2 (MQ=255)
                   gATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCcggt                                                     <  1:1290593/70‑2 (MQ=255)
                    aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                      >  1:45095/1‑52 (MQ=255)
                    aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                      >  1:4070017/1‑52 (MQ=255)
                    aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                      >  1:2680137/1‑52 (MQ=255)
                    aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                     >  1:2838786/1‑69 (MQ=255)
                      aCCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTag                                                  >  1:1469681/1‑70 (MQ=255)
                       ccGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTc                                                                    >  1:1595539/1‑51 (MQ=255)
                        cGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagttt                                               <  1:3017300/71‑4 (MQ=255)
                         gTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTaggagg                                              >  1:3180705/1‑71 (MQ=255)
                          ttCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGcc                                                         <  1:3628965/59‑1 (MQ=255)
                             atatTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCgg                                           <  1:2095546/70‑1 (MQ=255)
                              tatTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagt                                                 <  1:3632931/63‑2 (MQ=255)
                                ttCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCa                                                              >  1:4439182/1‑48 (MQ=255)
                                ttCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTagga                                                >  1:3892778/1‑62 (MQ=255)
                                ttCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATc                                                               <  1:3153770/47‑1 (MQ=255)
                                 tCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagtt                                                >  1:1346194/1‑59 (MQ=255)
                                  cTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATc                                      <  1:180377/70‑1 (MQ=255)
                                  cTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCg                                     >  1:2623153/1‑71 (MQ=255)
                                       cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGgag                                         <  1:2043722/62‑1 (MQ=255)
                                       cGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAg                                                             <  1:1227074/42‑1 (MQ=255)
                                        gATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTa                                  >  1:356855/1‑68 (MQ=255)
                                         aTAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATag                                                  <  1:1605184/51‑1 (MQ=255)
                                         aTAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATag                                                  <  1:2063184/51‑1 (MQ=255)
                                             ccTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGt                                    <  1:3274199/61‑1 (MQ=255)
                                                gtgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcggcgg                        >  1:2464280/1‑70 (MQ=255)
                                                 tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcggc                          <  1:758756/67‑1 (MQ=255)
                                                 tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGt                       <  1:3920667/70‑1 (MQ=255)
                                                 tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATAGTTACTGCCGGCGGt                       <  1:4952196/70‑1 (MQ=255)
                                                     aCCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAAt                  >  1:1350858/1‑71 (MQ=255)
                                                     aCCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAAt                  >  1:4572112/1‑71 (MQ=255)
                                                        gtgtCCTGCACTGTTTTCGCATCAGCGCCGGTGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGt               <  1:4952369/71‑1 (MQ=255)
                                                        gtgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCg                >  1:1962058/1‑70 (MQ=255)
                                                         tgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTa                                  >  1:2581657/1‑51 (MQ=255)
                                                             cTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATc                                      >  1:3605734/1‑43 (MQ=255)
                                                             cTGCACTGTTTTCGCATCAGCGCCGGGGTAGCAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGa           <  1:2107639/70‑1 (MQ=255)
                                                               gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGt                                    <  1:4156464/43‑1 (MQ=255)
                                                               gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGatat        >  1:2573190/1‑71 (MQ=255)
                                                                cACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATAtt       <  1:2573076/71‑1 (MQ=255)
                                                                cACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATAtt       >  1:1729493/1‑71 (MQ=255)
                                                                  cTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTgc     >  1:969068/1‑71 (MQ=255)
                                                                    gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTgcgc   >  1:4582172/1‑71 (MQ=255)
                                                                    gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTgcgc   >  1:4916248/1‑71 (MQ=255)
                                                                    gTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTgcgc   >  1:1742671/1‑71 (MQ=255)
                                                                     ttttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAg             <  1:198815/60‑1 (MQ=255)
                                                                     ttttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAg             <  1:1977454/60‑1 (MQ=255)
                                                                     ttttCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGcc  <  1:3193246/71‑1 (MQ=255)
                                                                      |                                                                     
ACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCGCC  >  minE/312638‑312777

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