Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,432 G→A 8.6% F290F (TTC→TTT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4320GA 8.6% 286.7 / 11.7 93F290F (TTC→TTTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (54/31);  new base A (5/3);  total (59/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGC  >  minE/1439365‑1439502
                                                                   |                                                                      
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:3584847/1‑70 (MQ=255)
    gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                                                  >  1:3460184/1‑70 (MQ=255)
    gtcgtAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg                                                                  >  1:2758161/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)
             aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                                        <  1:224211/71‑1 (MQ=255)
             aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                                        >  1:1587104/1‑71 (MQ=255)
             aCGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACt                                                        <  1:2695219/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: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:3243455‑M1/1‑70 (MQ=255)
                            ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                          >  1:860289‑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:908807/1‑71 (MQ=255)
                                ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt                                     >  1:1497516/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:434528/1‑71 (MQ=255)
                                        cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                             >  1:1470464/1‑71 (MQ=255)
                                        cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                             >  1:3622126/1‑71 (MQ=255)
                                        cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                             >  1:465393/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: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: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:1629976/1‑70 (MQ=255)
                                                          acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg            >  1:3512753/1‑70 (MQ=255)
                                                            gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGcc                     >  1:2047331/1‑59 (MQ=255)
                                                            gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGaa         >  1:3477011/1‑71 (MQ=255)
                                                             cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc        <  1:1950751/71‑1 (MQ=255)
                                                             cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc        <  1:2047919/71‑1 (MQ=255)
                                                              aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                        <  1:3145725/54‑1 (MQ=255)
                                                              aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAcc       >  1:2904924/1‑71 (MQ=255)
                                                               gACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCt      <  1:1629197/71‑1 (MQ=255)
                                                                aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa     >  1:2061859/1‑71 (MQ=255)
                                                                aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa     >  1:3220294/1‑71 (MQ=255)
                                                                aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTATTCGCCCAACCCAGGGAACCTa     >  1:1071219/1‑71 (MQ=255)
                                                                 ccGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCg                       >  1:2651472/1‑52 (MQ=255)
                                                                 ccGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGaa         >  1:1765266/1‑66 (MQ=255)
                                                                  cGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                        <  1:589999/50‑1 (MQ=255)
                                                                  cGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa     <  1:3089718/69‑1 (MQ=255)
                                                                   gAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAcgc  <  1:1510194/71‑1 (MQ=255)
                                                                   |                                                                      
CGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGC  >  minE/1439365‑1439502

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: