Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 893,899 T→C 50.4% I176V (ATC→GTC)  sapB ← predicted antimicrobial peptide transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE893,8990TC50.4% ‑2.6 / 119.3 125I176V (ATC→GTC) sapBpredicted antimicrobial peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (44/18);  new base C (44/19);  total (88/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 125 total aligned reads displayed.

TCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTC  >  minE/893835‑893966
                                                                |                                                                   
cctgttgttGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:1590106/66‑1 (MQ=255)
   gttgttGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCAttt                                                                <  1:86173/67‑1 (MQ=255)
   gttgttGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCAttt                                                                <  1:2865179/67‑1 (MQ=255)
      tttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2385534/62‑1 (MQ=255)
      gttgGCGCAACCTACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2061645/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATCTCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2202555/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:1035360/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2108576/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:1914043/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2499725/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:172091/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2746171/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:2749409/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:1473505/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:3259341/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:133529/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCAtt                                                                 <  1:742303/63‑1 (MQ=255)
      gttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTAATGACCAtt                                                                 <  1:1024714/63‑1 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2677619/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:903438/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2153844/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:805457/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2292300/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2352340/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2381332/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3397760/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3319837/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3109508/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:2737836/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3310527/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3273006/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3261346/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3018005/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3108336/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:3147769/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1038755/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1674944/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1702951/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1539645/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1801049/1‑68 (MQ=255)
        tgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATccc                                                          >  1:1804948/1‑68 (MQ=255)
           cgcAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGcc                                                       <  1:3344023/68‑1 (MQ=255)
           cgcAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGcc                                                       <  1:1448055/68‑1 (MQ=255)
           cgcAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGcc                                                       <  1:292789/68‑1 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCa                                                      >  1:1882774/1‑62 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:236484/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:1929654/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:2536554/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:3335169/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:1416945/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:2671738/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:1799268/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:2739618/1‑66 (MQ=255)
                  gACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:704596/1‑66 (MQ=255)
                  gACAGGCTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGt                                                  >  1:439922/1‑66 (MQ=255)
                             cACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGcc                                     >  1:2062204/1‑68 (MQ=255)
                             cACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCa                                    >  1:2142137/1‑69 (MQ=255)
                             cACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCa                                    >  1:2581015/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGt                                                  <  1:133284/44‑1 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGt                                                  <  1:3362054/44‑1 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGt                                                  <  1:713077/44‑1 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:1200808/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:2634654/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:1163594/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:303768/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:3019453/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:193666/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:1043354/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:2774593/1‑69 (MQ=255)
                                        tCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCatcaatca                         >  1:894504/1‑69 (MQ=255)
                                          atatGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGTGAGTTCGAAAGCCAGTCAtcaatcaa                        >  1:701557/1‑68 (MQ=255)
                                          atatGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGTGAGTTCGAAAGCCAGGCAtcaatcaa                        >  1:567733/1‑68 (MQ=255)
                                          atatGGCGGATTGCGCTCATGACCATTTCATCCCGCCACGGTGAGTTCGAAAGCCAGGCAtcaatcaa                        >  1:1316430/1‑68 (MQ=255)
                                            atGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACg                      >  1:1142080/1‑68 (MQ=255)
                                               gCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaaa                  >  1:728915/1‑69 (MQ=255)
                                               gCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaaa                  >  1:2771453/1‑69 (MQ=255)
                                               gCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCaaa                  >  1:2691/1‑69 (MQ=255)
                                                    ttGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCg              >  1:2576839/1‑68 (MQ=255)
                                                     tGCGCTCATGATCATTTCATCCCTCCATGGTGAGTCCGAAAGCCAGGCAtcaatcaa                        >  1:511678/1‑57 (MQ=255)
                                                     tGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGt            >  1:2389373/1‑69 (MQ=255)
                                                     tGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGt            <  1:36504/69‑1 (MQ=255)
                                                     tGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGt            <  1:45394/69‑1 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:2760830/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:1729880/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:801663/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:2917642/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:1935004/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:1512577/1‑69 (MQ=255)
                                                           cATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:1392215/1‑69 (MQ=255)
                                                           cATGATCATTACATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCgg      >  1:639368/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  <  1:586701/69‑1 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:67333/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  <  1:3062136/69‑1 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:1523240/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:3106104/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:1157948/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:859061/1‑69 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  <  1:2806729/69‑1 (MQ=255)
                                                               aTCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:1266246/1‑69 (MQ=255)
                                                               aTCATTTCATCCAGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTc  >  1:2773909/1‑69 (MQ=255)
                                                                |                                                                   
TCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTC  >  minE/893835‑893966

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: