Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,936 A→T 7.5% T122T (ACT→ACA) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9360AT 7.5% 226.6 / 8.7 80T122T (ACT→ACA) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base T (4/2);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACC  >  minE/1439868‑1440004
                                                                    |                                                                    
cGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGt                                                                    >  1:1492648/1‑71 (MQ=255)
cGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGt                                                                    <  1:2317743/71‑1 (MQ=255)
 gctgctGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGtt                                                                   <  1:1950039/71‑1 (MQ=255)
     ctgcATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACg                                                                <  1:853179/70‑1 (MQ=255)
      tgcATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGt                                                               >  1:3697373/1‑70 (MQ=255)
        cATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGaa                                                                      >  1:3909906/1‑61 (MQ=255)
        cATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCAt                                                            <  1:3415665/71‑1 (MQ=255)
         aTGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATa                                                           <  1:2642394/71‑1 (MQ=255)
          tGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGt                                                               >  1:1276129/1‑66 (MQ=255)
          tGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                          <  1:2359947/71‑1 (MQ=255)
           gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                                         >  1:3104139/1‑71 (MQ=255)
           gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                                         >  1:188509/1‑71 (MQ=255)
           gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                                         >  1:2826385/1‑71 (MQ=255)
            aaGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGaa                                                                      >  1:1486456/1‑57 (MQ=255)
            aaGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                                         <  1:939866/70‑1 (MQ=255)
            aaGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                        <  1:2376953/71‑1 (MQ=255)
             aGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaaca                                                       <  1:848054/71‑1 (MQ=255)
             aGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaaca                                                       <  1:3177711/71‑1 (MQ=255)
              gTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaacaa                                                      <  1:2699585/71‑1 (MQ=255)
                  tCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                   <  1:96884/70‑1 (MQ=255)
                  tCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                   <  1:3312659/70‑1 (MQ=255)
                    aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                 >  1:315302/1‑70 (MQ=255)
                    aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                 <  1:2137852/70‑1 (MQ=255)
                    aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                 <  1:2153638/70‑1 (MQ=255)
                    aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                 >  1:2843222/1‑70 (MQ=255)
                    aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                 <  1:1647549/70‑1 (MQ=255)
                     gAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTa                                               <  1:1576057/71‑1 (MQ=255)
                      aaCCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaacaac                                                     >  1:1993026/1‑64 (MQ=255)
                       aCCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                             >  1:2833090/1‑71 (MQ=255)
                          gTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGctg                                                 <  1:1797124/64‑3 (MQ=255)
                          gTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                          <  1:798849/71‑1 (MQ=255)
                           tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                           >  1:2214039/1‑69 (MQ=255)
                                gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCAt                                                            <  1:498840/47‑1 (MQ=255)
                                gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                            <  1:1654544/63‑1 (MQ=255)
                                gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                            <  1:3880346/63‑1 (MQ=255)
                                gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCg                                     >  1:3335823/1‑70 (MQ=255)
                                 tCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGt                                    <  1:3634502/70‑1 (MQ=255)
                                   gccaccaAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGt                                 >  1:1083262/2‑71 (MQ=2)
                                   accaccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCaa                                                            >  1:289273/1‑43 (MQ=255)
                                    ccaccaAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGt                                 <  1:1694418/70‑1 (MQ=2)
                                     caccGAATTCTGGCAGTACGTCGGTCAAGGAAGTTACGTCATAAACAACGCCGTa                                               >  1:616193/1‑55 (MQ=255)
                                      accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAg                                  <  1:1279486/67‑1 (MQ=255)
                                      accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                              <  1:2271558/71‑1 (MQ=255)
                                      accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                              <  1:762587/71‑1 (MQ=255)
                                       ccaaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                              >  1:1949291/4‑70 (MQ=2)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                            >  1:310308/1‑56 (MQ=255)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGAcc                                      <  1:1667860/62‑1 (MQ=255)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGAcc                                      <  1:2066189/62‑1 (MQ=255)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                               >  1:2167440/1‑69 (MQ=255)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                              >  1:3657604/1‑70 (MQ=255)
                                       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                              >  1:600831/1‑70 (MQ=255)
                                        cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGt                                                               <  1:2898704/36‑1 (MQ=255)
                                         gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                             >  1:3515643/1‑53 (MQ=255)
                                          aaTTCTGGCAGTACGTCGGTCCAGGAAGTTACGTc                                                              >  1:35051/1‑35 (MQ=255)
                                          aaTTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                           >  1:1611697/1‑70 (MQ=255)
                                          aaTTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGa                          <  1:2656700/71‑1 (MQ=255)
                                           aTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGaaa                            >  1:2394127/1‑68 (MQ=9)
                                           aTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGAAAgta                         >  1:155944/1‑69 (MQ=9)
                                            ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGAAAGTAc                        <  1:1918023/71‑1 (MQ=9)
                                              cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAcc                       >  1:291460/1‑70 (MQ=255)
                                               tGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCca                     >  1:1408694/1‑71 (MQ=255)
                                                ggCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACg                                         <  1:3208009/50‑1 (MQ=255)
                                                ggCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCca                     >  1:79342/1‑70 (MQ=255)
                                                 gCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCcaca                   <  1:375199/71‑1 (MQ=255)
                                                      aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAg                                              >  1:3292891/1‑39 (MQ=255)
                                                      aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                             <  1:2661232/40‑1 (MQ=255)
                                                       cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGTACCCACATCCTgg             <  1:2976083/71‑1 (MQ=255)
                                                       cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATcc                >  1:2994251/1‑68 (MQ=255)
                                                        gTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                           <  1:2125259/56‑1 (MQ=255)
                                                        gTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAccc                      <  1:995554/61‑1 (MQ=255)
                                                         tCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGt                                 <  1:753038/49‑1 (MQ=255)
                                                          cGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAt          >  1:782519/1‑71 (MQ=255)
                                                           ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAtt         <  1:2887498/71‑1 (MQ=255)
                                                           ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAt          <  1:2630664/70‑1 (MQ=255)
                                                               cAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCa      <  1:874016/70‑1 (MQ=255)
                                                                 ggAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc   <  1:2276580/71‑1 (MQ=255)
                                                                 ggAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc   >  1:1523235/1‑71 (MQ=255)
                                                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaaa                            >  1:585564/1‑45 (MQ=255)
                                                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTc       >  1:2799608/1‑66 (MQ=255)
                                                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAcc  >  1:1255727/1‑71 (MQ=255)
                                                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc   >  1:3669581/1‑70 (MQ=255)
                                                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc   >  1:199946/1‑70 (MQ=255)
                                                                    aGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGa    >  1:2132747/1‑67 (MQ=255)
                                                                    |                                                                    
CGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACC  >  minE/1439868‑1440004

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