Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,954 G→A 6.8% G116G (GGC→GGT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9540GA 6.8% 285.7 / 6.8 88G116G (GGC→GGTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (47/35);  new base A (4/2);  total (51/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.99e-01

TCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCC  >  minE/1439886‑1440024
                                                                    |                                                                      
tCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       <  1:96884/70‑1 (MQ=255)
tCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       <  1:3312659/70‑1 (MQ=255)
  aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                                     >  1:2843222/1‑70 (MQ=255)
  aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                                     <  1:2153638/70‑1 (MQ=255)
  aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                                     <  1:2137852/70‑1 (MQ=255)
  aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                                     >  1:315302/1‑70 (MQ=255)
  aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                                                     <  1:1647549/70‑1 (MQ=255)
   gAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTa                                                                   <  1:1576057/71‑1 (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)
              gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                                                <  1:3880346/63‑1 (MQ=255)
              gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                                                <  1:1654544/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)
                  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:600831/1‑70 (MQ=255)
                     ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                  >  1:3657604/1‑70 (MQ=255)
                       gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                                                 >  1:3515643/1‑53 (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)
                                                   gTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCa                          <  1:3419489/64‑1 (MQ=255)
                                                     tACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCt                                   >  1:1605084/1‑53 (MQ=255)
                                                      acttCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGa                 <  1:2208939/67‑1 (MQ=255)
                                                      aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg                  >  1:2903333/1‑69 (MQ=255)
                                                      aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGaa                >  1:3371628/1‑71 (MQ=255)
                                                      aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGaa                >  1:1440381/1‑71 (MQ=255)
                                                       cGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAAt               <  1:1304788/71‑1 (MQ=255)
                                                         tCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCt                     <  1:1617650/63‑1 (MQ=255)
                                                         tCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGc             <  1:1454130/71‑1 (MQ=255)
                                                          cATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:3697915/36‑1 (MQ=255)
                                                          cATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:3279062/36‑1 (MQ=255)
                                                          cATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGc             >  1:3269406/1‑70 (MQ=255)
                                                          cATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTCGAAtt                             >  1:2486308/1‑54 (MQ=255)
                                                           aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca           >  1:3716298/1‑71 (MQ=255)
                                                           aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca           >  1:828758/1‑71 (MQ=255)
                                                           aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca           >  1:1768496/1‑71 (MQ=255)
                                                            tAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca           >  1:3879485/1‑70 (MQ=255)
                                                             aaaCAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATc                                     >  1:1055595/1‑43 (MQ=255)
                                                             aaaCAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATc                                     >  1:919574/1‑43 (MQ=255)
                                                              aacaacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCcaca         >  1:1512988/1‑70 (MQ=255)
                                                                caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg      >  1:2098678/1‑71 (MQ=255)
                                                                caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGAGAATGCCACACgg      >  1:1062754/1‑71 (MQ=255)
                                                                 aacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACggg     >  1:600698/1‑71 (MQ=255)
                                                                 aacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACggg     >  1:2098816/1‑71 (MQ=255)
                                                                  acGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCt                                   >  1:1917987/1‑40 (MQ=255)
                                                                  acGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGt    >  1:583513/1‑71 (MQ=255)
                                                                  acGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGt    >  1:1725414/1‑71 (MQ=255)
                                                                   cGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGaa                >  1:2689896/1‑58 (MQ=255)
                                                                   cGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTc   >  1:2328904/1‑71 (MQ=255)
                                                                    gCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTcc  <  1:1377654/71‑1 (MQ=255)
                                                                    |                                                                      
TCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCC  >  minE/1439886‑1440024

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: