Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,429 A→G 8.8% G291G (GGT→GGC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4290AG 8.8% 239.7 / 12.0 91G291G (GGT→GGCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (51/32);  new base G (5/3);  total (56/35)
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.37e-01

TATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTA  >  minE/1439359‑1439499
                                                                      |                                                                      
tatCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGa                                                                        >  1:1032517/1‑71 (MQ=255)
 atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                                       >  1:2298242/1‑71 (MQ=255)
 atCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGAc                                                                       >  1:1199024/1‑71 (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: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:2695219/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)
                          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:955654/1‑71 (MQ=255)
                                      ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt                                  >  1:1497516/1‑71 (MQ=255)
                                      ttAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGtt                                  >  1:908807/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:3622126/1‑71 (MQ=255)
                                              cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                          >  1:1470464/1‑71 (MQ=255)
                                              cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCtt                          >  1:465393/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:3472081‑M2/5‑71 (MQ=255)
                                                 gtaaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                       >  1:1647389‑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)
                                                        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:1629976/1‑70 (MQ=255)
                                                                acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg         >  1:3512753/1‑70 (MQ=255)
                                                                acgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg         >  1:3947181/1‑70 (MQ=255)
                                                                  gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGcc                  >  1:2047331/1‑59 (MQ=255)
                                                                  gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGaa      >  1:3477011/1‑71 (MQ=255)
                                                                   cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc     <  1:2047919/71‑1 (MQ=255)
                                                                   cAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc     <  1:1950751/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:3220294/1‑71 (MQ=255)
                                                                      aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTa  >  1:2061859/1‑71 (MQ=255)
                                                                      aCCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTATTCGCCCAACCCAGGGAACCTa  >  1:1071219/1‑71 (MQ=255)
                                                                      |                                                                      
TATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTA  >  minE/1439359‑1439499

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: