Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,423 A→G 10.1% R293R (CGT→CGC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4230AG10.1% 208.1 / 12.7 79R293R (CGT→CGCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (46/25);  new base G (5/3);  total (51/28)
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.02e-01

TCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGG  >  minE/1439353‑1439492
                                                                      |                                                                     
tCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                                       >  1:2263455/1‑71 (MQ=255)
      tatCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGa                                                                 >  1:1032517/1‑71 (MQ=255)
       atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                                >  1:2298242/1‑71 (MQ=255)
       atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                                >  1:1199024/1‑71 (MQ=255)
        tCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                                       <  1:253573/63‑1 (MQ=255)
        tCTTCGTCGGCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGggacgga                                                                       <  1:3305122/63‑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:3584847/1‑70 (MQ=255)
                gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                                        >  1:3460184/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:1587104/1‑71 (MQ=255)
                         aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                              <  1:224211/71‑1 (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)
                                aGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAAc                                       <  1:3687990/71‑1 (MQ=255)
                                  gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTa                                                >  1:1830990/1‑60 (MQ=255)
                                  gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                                     >  1:1384207/1‑71 (MQ=255)
                                  gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                                     >  1:958831/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:955654/1‑71 (MQ=255)
                                            ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt                           >  1:908807/1‑71 (MQ=255)
                                            ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt                           >  1:1497516/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:434528/1‑71 (MQ=255)
                                                    cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                   >  1:1470464/1‑71 (MQ=255)
                                                    cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                   >  1:465393/1‑71 (MQ=255)
                                                    cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                   >  1:3622126/1‑71 (MQ=255)
                                                      ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACa                                      <  1:3770869‑M2/45‑1 (MQ=255)
                                                      ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                 >  1:2432084‑M2/6‑71 (MQ=255)
                                                      ggtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                 <  1:1808799‑M2/66‑1 (MQ=255)
                                                      ggTAAGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTg                 <  1:671175/71‑1 (MQ=12)
                                                       gtaaGCCATGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                >  1:2553311‑M2/5‑71 (MQ=255)
                                                       gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                >  1:3472081‑M2/5‑71 (MQ=255)
                                                       gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                >  1:1647389‑M2/5‑71 (MQ=255)
                                                        taaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt               <  1:3317240‑M2/68‑1 (MQ=255)
                                                         aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                >  1:3566872/1‑69 (MQ=12)
                                                         aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                >  1:3777063/1‑69 (MQ=12)
                                                         aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGt                >  1:1591140/1‑69 (MQ=12)
                                                          aGCCAGGGCCGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc              <  1:1460458‑M2/69‑1 (MQ=255)
                                                          aGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc              <  1:3164980‑M2/69‑1 (MQ=255)
                                                            ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                <  1:302289/50‑1 (MQ=255)
                                                            ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                <  1:1980543/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)
                                                                  acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAAccc     <  1:3771531/71‑1 (MQ=255)
                                                                   cggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                <  1:3892933/43‑1 (MQ=255)
                                                                   cggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAAccc     <  1:79359/70‑1 (MQ=255)
                                                                    ggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCtg                        <  1:1479414/50‑1 (MQ=255)
                                                                      acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg  >  1:3947181/1‑70 (MQ=255)
                                                                      acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg  >  1:3512753/1‑70 (MQ=255)
                                                                      acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg  >  1:1629976/1‑70 (MQ=255)
                                                                      |                                                                     
TCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGG  >  minE/1439353‑1439492

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: