Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,558,459 C→G 11.1% L167F (TTG→TTC emrK ← EmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,558,4590CG11.1% 175.5 / 11.9 63L167F (TTG→TTCemrKEmrKY‑TolC multidrug resistance efflux pump, membrane fusion protein component
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base G (2/5);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCC  >  Exported/1558390‑1558524
                                                                     |                                                                 
tgGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTc                                                                   <  1:1576657/70‑1 (MQ=255)
  gCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCaa                                                                 <  1:2557678/70‑1 (MQ=255)
    tGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGc                                                              >  1:1696945/1‑71 (MQ=255)
     gACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGcc                                                             >  1:2298832/1‑71 (MQ=255)
      aCGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCg                                                            >  1:2184358/1‑71 (MQ=255)
        ggTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCt                                                          <  1:3484873/71‑1 (MQ=255)
        ggTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCt                                                          <  1:1693657/71‑1 (MQ=255)
        ggTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCaa                                                                 <  1:410078/64‑1 (MQ=255)
        ggTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCt                                                          <  1:3067635/71‑1 (MQ=255)
         gTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCtt                                                         <  1:2272922/71‑1 (MQ=255)
         gTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCtt                                                         <  1:2125725/71‑1 (MQ=255)
          tttAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTGAATGCCGCt                                                          <  1:2372768/69‑1 (MQ=255)
          tttAATGGCGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCttt                                                        <  1:3105348/71‑1 (MQ=255)
             aaTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTc                                                                   <  1:2587212/57‑1 (MQ=255)
             aaTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTg                                                       >  1:1806221/1‑69 (MQ=255)
             aaTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGc                                                      >  1:3426908/1‑70 (MQ=255)
             aaTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGc                                                      <  1:855210/70‑1 (MQ=255)
              aTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCt                                                     <  1:2555295/70‑1 (MQ=255)
              aTGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTa                                                    >  1:3957/1‑71 (MQ=255)
               tGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTAc                                                   <  1:746524/71‑1 (MQ=255)
                     gtTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttatt                                              >  1:3193629/1‑70 (MQ=255)
                     gtTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttatt                                              >  1:2128047/1‑70 (MQ=255)
                      ttCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttatta                                             <  1:981952/70‑1 (MQ=255)
                      ttCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttatta                                             <  1:868970/70‑1 (MQ=255)
                         aTTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCg                                                            <  1:183500/52‑1 (MQ=255)
                         aTTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCg                                                            <  1:704629/52‑1 (MQ=255)
                          ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCa                                                                  >  1:852964/1‑45 (MQ=255)
                          ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCa                                                                  >  1:2298234/1‑45 (MQ=255)
                          ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTa                                        <  1:847411/71‑1 (MQ=255)
                          ttACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTa                                        <  1:3221911/71‑1 (MQ=255)
                            aCTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATc                                      >  1:911705/1‑71 (MQ=255)
                              tAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCtt                                    >  1:2542462/1‑71 (MQ=255)
                                  gCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCaa                                                                 >  1:821412/1‑38 (MQ=255)
                                  gCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTgg                                 >  1:1953539/1‑70 (MQ=255)
                                    tttATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCtt                                    <  1:3539623/65‑1 (MQ=255)
                                     ttATTCGCTTTATAAGCCTGGATAGCGGCATTCaa                                                                 <  1:2053904/35‑1 (MQ=255)
                                     ttATTCGCTTTATAAGCCTGGATAGCGGCATTCaa                                                                 <  1:1075979/35‑1 (MQ=255)
                                      tATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGc                            >  1:2479374/1‑71 (MQ=255)
                                       aTTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                               >  1:333308/1‑67 (MQ=255)
                                        ttCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACttatt                                              <  1:190179/51‑1 (MQ=255)
                                           gCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCa                        <  1:460049/70‑1 (MQ=255)
                                           gCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAg                       >  1:789015/1‑71 (MQ=255)
                                              ttATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTAt                                       >  1:2126084/1‑52 (MQ=255)
                                               tataAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGtt                   >  1:1125227/1‑71 (MQ=255)
                                                  aaGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTcc                         >  1:1621259/1‑62 (MQ=255)
                                                     ccTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                               >  1:1883668/1‑53 (MQ=255)
                                                      cTGGATAGCGGCATTCAATGCCGCTTTGCTACttatt                                              >  1:1164856/1‑37 (MQ=255)
                                                      cTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCtttt             >  1:2502317/1‑70 (MQ=255)
                                                      cTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCtttt             <  1:2382057/70‑1 (MQ=255)
                                                        ggATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGa           >  1:519582/1‑70 (MQ=255)
                                                         gATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaaa         >  1:2373961/1‑71 (MQ=255)
                                                         gATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaa          <  1:1227254/70‑1 (MQ=255)
                                                         gATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaa          <  1:404497/70‑1 (MQ=255)
                                                         gATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaa          <  1:848674/70‑1 (MQ=255)
                                                          aTAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaaa         >  1:1520340/1‑70 (MQ=255)
                                                          aTAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaaa         >  1:2254248/1‑70 (MQ=255)
                                                          aTAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGaa          >  1:3396480/1‑69 (MQ=255)
                                                           tAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATa       <  1:3604188/71‑1 (MQ=255)
                                                           aaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATc                                      <  1:819442/39‑1 (MQ=37)
                                                           aaGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATc                                      <  1:894284/39‑1 (MQ=37)
                                                            aGCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATaa      >  1:503942/1‑71 (MQ=255)
                                                             gCGGCATTGAATGCCGCTTTGCTACTTATTAACGTATCtt                                    >  1:3423886/1‑40 (MQ=255)
                                                                gCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAAcccc  >  1:2523100/1‑71 (MQ=255)
                                                                gCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAAcc    >  1:723309/1‑69 (MQ=255)
                                                                gCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAAcc    >  1:85075/1‑69 (MQ=255)
                                                                     cAATGCCGCTTTGCTACTTATTAACGTATCTTTgg                                 >  1:2795683/1‑35 (MQ=255)
                                                                     cAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                               >  1:767130/1‑37 (MQ=255)
                                                                     cAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                               >  1:472153/1‑37 (MQ=255)
                                                                     cAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTa                               >  1:1300975/1‑37 (MQ=255)
                                                                     cAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTc                          <  1:1821272/42‑1 (MQ=255)
                                                                     |                                                                 
TGGCTGACGGTTTAATGGTGTGTTCATTACTAAAGCTTTATTCGCTTTATAAGCCTGGATAGCGGCATTCAATGCCGCTTTGCTACTTATTAACGTATCTTTGGTATGCTCCAGCGTTTCTTTTGAAATAACCCC  >  Exported/1558390‑1558524

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: