Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,416 G→T 10.4% L296M (CTG→ATG) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4160GT10.4% 268.5 / 11.3 86L296M (CTG→ATG) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (44/33);  new base T (6/3);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAA  >  minE/1439346‑1439486
                                                                      |                                                                      
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)
      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)
                              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:193162/1‑71 (MQ=255)
                                                    tAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc                    >  1:3518907/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:1470464/1‑71 (MQ=255)
                                                           cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt             >  1:465393/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:1591140/1‑69 (MQ=12)
                                                                aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt          >  1:3777063/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)
                                                                     aGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGccc    >  1:2357432/1‑70 (MQ=255)
                                                                     aGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGc      >  1:1691457/1‑68 (MQ=255)
                                                                      gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCg       >  1:3157843/1‑66 (MQ=255)
                                                                      gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCaa  >  1:2678961/1‑71 (MQ=255)
                                                                      |                                                                      
ACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAA  >  minE/1439346‑1439486

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: