Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 312,705 T→G 12.2% T56T (ACA→ACC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE312,7050TG12.2% 179.0 / 23.6 82T56T (ACA→ACCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (22/50);  new base G (7/3);  total (29/53)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCG  >  minE/312636‑312775
                                                                     |                                                                      
ggACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTg                                                                       <  1:1190008/71‑1 (MQ=255)
   cATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                    <  1:1026371/71‑1 (MQ=255)
   cATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                    <  1:248404/71‑1 (MQ=255)
   cATGTACATCAGCTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGttt                                                                    <  1:1211857/71‑1 (MQ=255)
     tGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTc                                                                  <  1:1535085/71‑1 (MQ=255)
       tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTc                                                                  >  1:1247704/1‑69 (MQ=255)
       tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTc                                                                  >  1:612396/1‑69 (MQ=255)
       tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCg                                                                 >  1:352436/1‑70 (MQ=255)
       tACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGc                                                                <  1:1873206/71‑1 (MQ=255)
        aCATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCa                                                               >  1:2709462/1‑71 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCTCAt                                                              <  1:1813843/71‑1 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              <  1:3837007/71‑1 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              <  1:514134/71‑1 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              <  1:3274708/71‑1 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              <  1:265665/71‑1 (MQ=255)
         cATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              <  1:1594504/71‑1 (MQ=255)
          aTCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCAt                                                              >  1:964414/1‑70 (MQ=255)
          aTCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             <  1:3355219/71‑1 (MQ=255)
          aTCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             <  1:272665/71‑1 (MQ=255)
          aTCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             <  1:1739761/71‑1 (MQ=255)
          aTCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCAt                                                              >  1:1037871/1‑70 (MQ=37)
           tCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCa                                                            <  1:53633/71‑1 (MQ=255)
             aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGc                                                                >  1:2659004/1‑65 (MQ=255)
             aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGTATCAg                                                           <  1:480624/70‑1 (MQ=255)
             aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAg                                                           <  1:405278/70‑1 (MQ=255)
             aGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAgc                                                          >  1:3362124/1‑71 (MQ=37)
                ttATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             >  1:2431728/1‑65 (MQ=255)
                 tATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             >  1:1697515/1‑64 (MQ=255)
                 tATCGATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGccc                                                      >  1:714867/1‑70 (MQ=255)
                  aTCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTc                                                                  <  1:1077446/58‑1 (MQ=255)
                    cGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                   >  1:1807409/1‑71 (MQ=255)
                    cGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                   >  1:2923637/1‑71 (MQ=255)
                     gATACCGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCcgga                                                   <  1:3615062/70‑2 (MQ=255)
                      aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCg                                                      <  1:2667872/66‑1 (MQ=255)
                      aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCg                                                      <  1:24974/66‑1 (MQ=255)
                      aTACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTa                                                 <  1:858137/71‑1 (MQ=255)
                         ccGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTagga                                              <  1:1470941/71‑1 (MQ=255)
                         ccGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCa                                                            >  1:893231/1‑57 (MQ=255)
                         ccGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCa                                                            >  1:683653/1‑57 (MQ=255)
                         ccGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagtt                                              >  1:3706071/1‑69 (MQ=25)
                          cGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagttt                                             <  1:1412321/71‑4 (MQ=255)
                          cGTTCATATTCTGTTCGATAACCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCCGGATagttt                                             <  1:173351/71‑4 (MQ=255)
                           gTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             <  1:3919259/54‑1 (MQ=255)
                              cATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             >  1:1537787/1‑51 (MQ=255)
                               atatTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGcc                                                       <  1:1086841/56‑1 (MQ=255)
                                 atTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATc                                                             <  1:43547/48‑1 (MQ=255)
                                  ttCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGAt                                     <  1:3928209/71‑1 (MQ=255)
                                    cTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGtttt                                                                   >  1:575913/1‑39 (MQ=255)
                                    cTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGtttt                                                                   >  1:346902/1‑39 (MQ=255)
                                    cTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCTGGGTAGGAGGCGGAGATCg                                   >  1:2437973/1‑71 (MQ=255)
                                       ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGcc                                                       <  1:1009101/48‑1 (MQ=255)
                                       ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGcc                                                       <  1:1310219/48‑1 (MQ=255)
                                       ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCg                                   <  1:2771027/68‑1 (MQ=255)
                                       ttCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGt                                  <  1:315028/69‑1 (MQ=255)
                                        tCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTAc                               >  1:3030298/1‑71 (MQ=255)
                                         cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGTGCCGGGGTAGGAGGCGGAGATCGTTAc                               >  1:1280558/1‑70 (MQ=255)
                                         cGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACt                              >  1:348681/1‑71 (MQ=255)
                                           aTAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTagga                                              <  1:3662500/53‑1 (MQ=255)
                                           aTAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGc                            <  1:3503577/71‑1 (MQ=255)
                                             aaCCTGCGTCACCGTATCCTGCACGGTCTGCGCATCAGCGCCcggat                                                  >  1:1763248/1‑45 (MQ=16)
                                                cTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcggcg                       >  1:3112881/1‑71 (MQ=255)
                                                   tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcggcgg                      <  1:56683/69‑1 (MQ=255)
                                                   tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGt                     <  1:163203/70‑1 (MQ=255)
                                                   tgtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGt                     <  1:575948/70‑1 (MQ=255)
                                                    gtCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGt                                  <  1:1452694/56‑1 (MQ=255)
                                                       aCCGTGTCCTGCACTGTTTTCGCATCAGCGCCgggg                                                   >  1:1440289/1‑36 (MQ=255)
                                                         cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCg              <  1:3818560/71‑1 (MQ=255)
                                                         cGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCg              <  1:386248/71‑1 (MQ=255)
                                                          gtgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGt             <  1:3424427/71‑1 (MQ=255)
                                                          gtgtCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGt             <  1:3062016/71‑1 (MQ=255)
                                                               cTGCACTGTTTTCGCATCAGCGCCGGGGTAGTAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGat        <  1:2128348/71‑1 (MQ=255)
                                                                 gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGgaga                                      >  1:872266/1‑39 (MQ=255)
                                                                 gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGtt                                 >  1:2321990/1‑44 (MQ=255)
                                                                 gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcgg                         <  1:1920837/52‑1 (MQ=255)
                                                                 gCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcgg                         <  1:1913564/52‑1 (MQ=255)
                                                                  cACTGTTTTCGCATCAGCGCCGGGGTAGGATGCGGAGATCGTTACTGCCGGCGGTGCATTCGTAGGATAtt     <  1:1753862/71‑1 (MQ=255)
                                                                  cACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCGGGCGGTGCAATCGTAGGatat      <  1:129791/70‑1 (MQ=255)
                                                                  cACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGatat      <  1:720345/70‑1 (MQ=255)
                                                                  cACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATAtt     <  1:221956/71‑1 (MQ=255)
                                                                   aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCcggcg                       >  1:1228041/1‑52 (MQ=255)
                                                                   aCTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTg    <  1:2905290/71‑1 (MQ=255)
                                                                   aCTGTTTTCGAATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTg    <  1:793099/71‑1 (MQ=255)
                                                                     tGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGAt                                     >  1:2879259/1‑36 (MQ=255)
                                                                     tGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGt                                  <  1:1727680/39‑1 (MQ=255)
                                                                     tGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTgcg  >  1:959724/1‑71 (MQ=255)
                                                                     tGTTTTCGCATCAGCGCCGGGGTAGGAGGAGGAGATCGtt                                 >  1:541529/1‑40 (MQ=255)
                                                                     |                                                                      
GGACATGTACATCAGGTTATCGATACCGTTCATATTCTGTTCGATAACCTGTGTCACCGTGTCCTGCACTGTTTTCGCATCAGCGCCGGGGTAGGAGGCGGAGATCGTTACTGCCGGCGGTGCAATCGTAGGATATTGCG  >  minE/312636‑312775

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: