Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,923 C→T 7.0% V127I (GTA→ATA) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9230CT 7.0% 276.0 / 6.8 85V127I (GTA→ATA) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (39/40);  new base T (4/2);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

GCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGG  >  minE/1439853‑1439993
                                                                      |                                                                      
gcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTAc                                                                        <  1:2274427/71‑1 (MQ=255)
  gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt                                                                      >  1:1132895/1‑71 (MQ=255)
  gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt                                                                      <  1:2873946/71‑1 (MQ=255)
  gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt                                                                      <  1:94569/71‑1 (MQ=255)
    aGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt                                                                      >  1:2782432/1‑69 (MQ=255)
       cGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTc                                                                 <  1:2611029/71‑1 (MQ=255)
        gTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTcc                                                                <  1:1345636/71‑1 (MQ=255)
        gTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTc                                                                 >  1:2690092/1‑70 (MQ=255)
        gTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTc                                                                 >  1:3316357/1‑70 (MQ=255)
          tACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCg                                                                    >  1:3235536/1‑65 (MQ=255)
          tACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAg                                                              >  1:2783518/1‑71 (MQ=255)
          tACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAg                                                              >  1:2858440/1‑71 (MQ=255)
           aCCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAg                                                              <  1:133951/70‑1 (MQ=255)
           aCCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAg                                                              >  1:1542061/1‑70 (MQ=255)
           aCCACGCTGCTGCAGGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAg                                                              <  1:3584891/70‑1 (MQ=255)
            ccACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGa                                                            <  1:2140960/71‑1 (MQ=255)
             cACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAgg                                                             >  1:3779828/1‑69 (MQ=255)
             cACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAgg                                                             >  1:46729/1‑69 (MQ=255)
             cACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGCTCCAgg                                                             >  1:2193691/1‑69 (MQ=255)
               cGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGt                                                         <  1:2317743/71‑1 (MQ=255)
               cGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGt                                                         >  1:1492648/1‑71 (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:2826385/1‑71 (MQ=255)
                          gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                              >  1:3104139/1‑71 (MQ=255)
                          gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                              >  1:188509/1‑71 (MQ=255)
                          gAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGAACTACGTCGGTcc                                                                >  1:2863230/1‑53 (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:1647549/70‑1 (MQ=255)
                                   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                      >  1:2843222/1‑70 (MQ=255)
                                   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                      <  1:2137852/70‑1 (MQ=255)
                                   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                      <  1:2153638/70‑1 (MQ=255)
                                   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCg                                      >  1:315302/1‑70 (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: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)
                                                  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:2066189/62‑1 (MQ=255)
                                                      ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGAcc                           <  1:1667860/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)
                                                       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)
                                                                      |                                                                      
GCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGG  >  minE/1439853‑1439993

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: