Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,414 C→G 9.1% L296L (CTG→CTC) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4140CG 9.1% 267.7 / 12.0 88L296L (CTG→CTC) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (46/34);  new base G (5/3);  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 = 1.00e+00

ACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTC  >  minE/1439346‑1439480
                                                                    |                                                                  
acatATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAg                                                                  >  1:1508271/1‑71 (MQ=255)
acatATTTCAGGATATCTTCGTCGGCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAg                                                                  <  1:2246300/71‑1 (MQ=255)
 catATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAgg                                                                 <  1:3630275/71‑1 (MQ=255)
 catATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAgg                                                                 >  1:3801280/1‑71 (MQ=255)
  atatTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAg                                                                  >  1:2789421/1‑69 (MQ=255)
     tttCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggac                                                              <  1:2195790/70‑1 (MQ=255)
     tttCAGGATATCTTCGTCGGCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGcc                                                                    <  1:1034601/64‑1 (MQ=255)
      ttCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacg                                                             <  1:3149474/70‑1 (MQ=255)
       tCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                           >  1:2263455/1‑71 (MQ=255)
       tCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgg                                                            <  1:1572496/70‑1 (MQ=255)
             tatCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGa                                                     >  1:1032517/1‑71 (MQ=255)
              atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                    >  1:1199024/1‑71 (MQ=255)
              atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                    >  1:2298242/1‑71 (MQ=255)
               tCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                           <  1:253573/63‑1 (MQ=255)
               tCTTCGTCGGCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                           <  1:3305122/63‑1 (MQ=255)
               tCTTAGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAgg                                                                 <  1:1418896/57‑1 (MQ=255)
                 ttCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCa                                                       <  1:425678/65‑1 (MQ=255)
                   cgtcgtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCg                                                  <  1:1498991/68‑1 (MQ=255)
                   cgtcgtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGaa                                                <  1:3224966/70‑1 (MQ=255)
                      cgtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTc                                             >  1:2532630/1‑70 (MQ=255)
                       gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTcta                                           <  1:1558213/71‑3 (MQ=255)
                       gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                            >  1:2758161/1‑70 (MQ=255)
                       gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                            >  1:3460184/1‑70 (MQ=255)
                       gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                            >  1:3584847/1‑70 (MQ=255)
                          gtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACt                                        >  1:490949/1‑71 (MQ=255)
                          gtAGCCACGACCCAGGTGTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACt                                        >  1:1622655/1‑71 (MQ=255)
                            aGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGcc                                                                    <  1:2814247/41‑1 (MQ=255)
                            aGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGcc                                                                    <  1:1069169/41‑1 (MQ=255)
                              ccACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGt                                     <  1:3153461/70‑1 (MQ=255)
                              ccACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTa                                    >  1:1158037/1‑71 (MQ=255)
                                aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGgacggacg                                                         <  1:3453680/48‑1 (MQ=255)
                                aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  <  1:2695219/71‑1 (MQ=255)
                                aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  <  1:224211/71‑1 (MQ=255)
                                aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  >  1:1587104/1‑71 (MQ=255)
                                  gACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAg                                               >  1:2036901/1‑56 (MQ=255)
                                  gACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  >  1:1391016/1‑69 (MQ=255)
                                  gACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGa                                <  1:1939108/71‑1 (MQ=255)
                                    cccAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACt                                        >  1:1017801/1‑61 (MQ=255)
                                     ccAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGa                                                 <  1:2296163/51‑1 (MQ=255)
                                      cAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacg                                                             <  1:2955307/38‑1 (MQ=255)
                                       aGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAAc                           <  1:3687990/71‑1 (MQ=255)
                                        ggTTTTTACCTTTAGACTGCAGGTAAGCCAGggac                                                              <  1:383951/35‑1 (MQ=255)
                                         gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTa                                    >  1:1830990/1‑60 (MQ=255)
                                         gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                         >  1:958831/1‑71 (MQ=255)
                                         gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                         >  1:1384207/1‑71 (MQ=255)
                                          tttttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  <  1:3264392/61‑1 (MQ=255)
                                           ttttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  <  1:2586165/60‑1 (MQ=255)
                                            tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGa                                >  1:2519782/1‑61 (MQ=255)
                                            tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCtt                      >  1:3609629/1‑71 (MQ=255)
                                            tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCt                       <  1:3653053/70‑1 (MQ=255)
                                             ttacctgcagtctaaaggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                  <  1:3106701‑M2/37‑1 (MQ=255)
                                             TTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                            >  1:749762‑M1/1‑48 (MQ=255)
                                              TACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAAc                                         >  1:252680‑M1/1‑50 (MQ=255)
                                              TACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                         >  1:2919814‑M1/1‑66 (MQ=255)
                                               acctgcagtctaaaggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    >  1:2292183‑M2/20‑70 (MQ=255)
                                               aCCTTTAGATTGCAGGTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTca                    <  1:1219581/70‑2 (MQ=18)
                                               ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGa                                                 <  1:1031123‑M1/41‑1 (MQ=255)
                                               ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    >  1:860289‑M1/1‑70 (MQ=255)
                                               ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    >  1:3243455‑M1/1‑70 (MQ=255)
                                                cctgcagtctaaaggtaaGCCAGGGACGGACGCAGACCg                                                  <  1:2206125‑M2/21‑1 (MQ=255)
                                                ccTTTAGATTGCAGGTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCaaa                  <  1:3452313/71‑1 (MQ=16)
                                                CCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    >  1:3684646‑M1/1‑69 (MQ=255)
                                                 cTTTAGATTGCAGGTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAg                 >  1:91448/1‑71 (MQ=16)
                                                   ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGc                              >  1:837870/1‑56 (MQ=255)
                                                   ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt               >  1:1497516/1‑71 (MQ=255)
                                                   ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt               >  1:908807/1‑71 (MQ=255)
                                                   ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt               >  1:955654/1‑71 (MQ=255)
                                                    tAGACTGCAGGTAAGCCAGGGACGGACGCAGACCg                                                  <  1:3612231/35‑1 (MQ=255)
                                                    tAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc              >  1:3518907/1‑71 (MQ=255)
                                                    tAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc              >  1:193162/1‑71 (MQ=255)
                                                    tAGACTGCAGCTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc              >  1:1662945/1‑71 (MQ=255)
                                                      gACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCtg            <  1:2175608/71‑1 (MQ=255)
                                                        cTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCtgt           >  1:3358854/1‑70 (MQ=255)
                                                         tGCAGGTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTtggg          >  1:3539664/1‑68 (MQ=12)
                                                           cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt       >  1:465393/1‑71 (MQ=255)
                                                           cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt       >  1:1470464/1‑71 (MQ=255)
                                                           cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt       >  1:434528/1‑71 (MQ=255)
                                                           cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt       >  1:3622126/1‑71 (MQ=255)
                                                             ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACa                          <  1:3770869‑M2/45‑1 (MQ=255)
                                                             ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg     <  1:1808799‑M2/66‑1 (MQ=255)
                                                             ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg     >  1:2432084‑M2/6‑71 (MQ=255)
                                                             ggTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTg     <  1:671175/71‑1 (MQ=12)
                                                              gtaaGCCATGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt    >  1:2553311‑M2/5‑71 (MQ=255)
                                                              gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt    >  1:1647389‑M2/5‑71 (MQ=255)
                                                              gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt    >  1:3472081‑M2/5‑71 (MQ=255)
                                                               taaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt   <  1:3317240‑M2/68‑1 (MQ=255)
                                                                aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt    >  1:3777063/1‑69 (MQ=12)
                                                                aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt    >  1:1591140/1‑69 (MQ=12)
                                                                aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt    >  1:3566872/1‑69 (MQ=12)
                                                                 aGCCAGGGCCGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc  <  1:1460458‑M2/69‑1 (MQ=255)
                                                                 aGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc  <  1:3164980‑M2/69‑1 (MQ=255)
                                                                   ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    <  1:1980543/50‑1 (MQ=255)
                                                                   ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                    <  1:302289/50‑1 (MQ=255)
                                                                    |                                                                  
ACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTC  >  minE/1439346‑1439480

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