Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,969 G→A 7.4% Y111Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9690GA 7.4% 261.0 / 7.3 81Y111Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (46/29);  new base A (4/2);  total (50/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCG  >  minE/1439900‑1440038
                                                                     |                                                                     
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)
      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)
       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                                                                 >  1:2167440/1‑69 (MQ=255)
       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                >  1:3657604/1‑70 (MQ=255)
       ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                >  1:600831/1‑70 (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)
                ggCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCca                                                       >  1:79342/1‑70 (MQ=255)
                 gCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCcaca                                                     <  1:375199/71‑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:199946/1‑70 (MQ=255)
                                  gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc                                     >  1:3669581/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:1440381/1‑71 (MQ=255)
                                        aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGaa                              >  1:3371628/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:1768496/1‑71 (MQ=255)
                                             aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca                         >  1:3716298/1‑71 (MQ=255)
                                             aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca                         >  1:828758/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:1725414/1‑71 (MQ=255)
                                                    acGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGt                  >  1:583513/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)
                                                        cGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCa               >  1:2223510/1‑70 (MQ=255)
                                                         gTAGTTACGACCGTAGTCGAAAGAACCCACATCCTgg                                               <  1:3572317/37‑1 (MQ=255)
                                                          tAGTTACGACCGTAGTCGAAAGAACCCACATCATGGAATTTCAGACCTGCGAATGCCACACGGGTCCAg              >  1:1290762/1‑69 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAg           <  1:485373/70‑1 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAg           <  1:1200416/70‑1 (MQ=255)
                                                              tACGACCGTAGTCTAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgtt         >  1:5562/1‑70 (MQ=255)
                                                              tACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttg        >  1:712264/1‑71 (MQ=255)
                                                              tACGACCGTAGTCGAAAGAACCAACATCCTGGAATTTc                                         >  1:1831646/1‑38 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttg        <  1:3569287/70‑1 (MQ=255)
                                                                cGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttgtt      <  1:1126039/71‑1 (MQ=255)
                                                                cGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGCAgttgt       <  1:510866/70‑1 (MQ=255)
                                                                 gACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttt     <  1:2377203/71‑1 (MQ=255)
                                                                  aCCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttt     >  1:2421305/1‑70 (MQ=255)
                                                                   ccGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttc   >  1:1866773/1‑71 (MQ=255)
                                                                    cGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttcg  >  1:1010611/1‑71 (MQ=255)
                                                                     |                                                                     
GTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCG  >  minE/1439900‑1440038

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: