Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,468 C→T 6.4% Q278Q (CAG→CAA ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4680CT 6.4% 261.4 / 5.6 78Q278Q (CAG→CAAompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (43/30);  new base T (2/3);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCC  >  minE/1439399‑1439537
                                                                     |                                                                     
aGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc                                                                       <  1:966320/70‑1 (MQ=255)
aGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTc                                                                       <  1:3420331/70‑1 (MQ=255)
aGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCt                                                                      >  1:1810798/1‑71 (MQ=255)
    tGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGc                                                                  >  1:1274405/1‑71 (MQ=255)
      cAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTAGAAGTTc                                                                       <  1:551192/64‑1 (MQ=255)
         gTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                              >  1:1665663/1‑70 (MQ=255)
         gTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                                                             >  1:3461042/1‑71 (MQ=255)
          tAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                            <  1:2011551/71‑1 (MQ=255)
           aaGCGATGGACGGGCGCAGGCCAAAGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGttttct                                                             >  1:2652852/1‑67 (MQ=9)
           aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAGGTTCTGTGCTTTGTTc                                                           >  1:1056659/1‑71 (MQ=255)
           aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                           >  1:1658874/1‑71 (MQ=255)
           aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                           >  1:3070839/1‑71 (MQ=255)
              ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                            <  1:4109220/67‑1 (MQ=255)
               cAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCt                                                                 >  1:2892149/1‑61 (MQ=255)
                aGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCa                                                      <  1:179870/71‑1 (MQ=255)
                 gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                              <  1:1160238/62‑1 (MQ=255)
                 gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                              <  1:1419433/62‑1 (MQ=255)
                  ggacggacGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAAc                                                    >  1:149465/1‑71 (MQ=255)
                       gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                               >  1:3998221/1‑71 (MQ=255)
                       gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                               <  1:171568/71‑1 (MQ=255)
                         cgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGa                                             <  1:1212527/71‑1 (MQ=255)
                            aGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc                                           >  1:289134/1‑70 (MQ=255)
                             gACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAg                                                >  1:31521/1‑64 (MQ=255)
                                caaaGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGc                                                         >  1:3545070/3‑52 (MQ=255)
                                  aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                            >  1:4106496/1‑47 (MQ=255)
                                  aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                            >  1:192512/1‑47 (MQ=255)
                                  aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGt                                                             >  1:3117034/1‑46 (MQ=255)
                                  aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAcgcg                                    <  1:2940728/71‑1 (MQ=255)
                                   aGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAcc                                          <  1:3986003/64‑1 (MQ=255)
                                   aGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGc                                                         <  1:859923/49‑1 (MQ=16)
                                   aGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGc                                                         <  1:3846794/49‑1 (MQ=16)
                                     tCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                              <  1:2165320/42‑1 (MQ=255)
                                       gAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTg                               >  1:2052527/1‑71 (MQ=255)
                                       gAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTg                               >  1:669401/1‑71 (MQ=255)
                                          cTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGtt                                                            <  1:3270655/39‑1 (MQ=12)
                                            ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAc                                       <  1:379509/58‑1 (MQ=255)
                                            ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAc                                       <  1:2111609/58‑1 (MQ=255)
                                             gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGTGAACCTACGCGAGTTGCGtt                           <  1:2735097/69‑1 (MQ=255)
                                             gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGtt                           <  1:2660799/69‑1 (MQ=255)
                                             gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                          <  1:681393/70‑1 (MQ=255)
                                             gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                          <  1:3648082/70‑1 (MQ=255)
                                              tACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                          >  1:288948/1‑69 (MQ=255)
                                                cTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAgg                      >  1:216314/1‑71 (MQ=255)
                                                  gAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTc                    >  1:1691355/1‑71 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCaa                                                     >  1:2115172/1‑37 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCt                                         >  1:395924/1‑49 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGtt                           <  1:1796790/63‑1 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTc                    >  1:4046781/1‑70 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                   >  1:493184/1‑71 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                   >  1:1684308/1‑71 (MQ=255)
                                                   aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                   >  1:594248/1‑71 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                               >  1:597374/1‑41 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                               >  1:686870/1‑41 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                   >  1:3734296/1‑69 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                 <  1:1491548/71‑1 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                 <  1:3340915/71‑1 (MQ=255)
                                                     cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                 <  1:2002852/71‑1 (MQ=255)
                                                      aaCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAg                                                >  1:2134187/1‑39 (MQ=255)
                                                       aCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgt               >  1:1253464/1‑71 (MQ=255)
                                                       aCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgt               >  1:445605/1‑71 (MQ=255)
                                                        cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg              <  1:1594091/71‑1 (MQ=255)
                                                        cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg              <  1:3738859/71‑1 (MQ=255)
                                                         aGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg              >  1:536009/1‑70 (MQ=255)
                                                          gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc                                           >  1:1534145/1‑40 (MQ=255)
                                                          gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                          <  1:1737911/57‑1 (MQ=255)
                                                          gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTa                        >  1:1884437/1‑59 (MQ=255)
                                                          gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTa                        >  1:4096178/1‑59 (MQ=255)
                                                          gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTa            <  1:2900689/71‑1 (MQ=255)
                                                            ttCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAAc                                           >  1:1849278/1‑38 (MQ=255)
                                                            ttCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTa                        >  1:1084236/1‑57 (MQ=255)
                                                               gAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                          <  1:217462/52‑1 (MQ=255)
                                                                 aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGa                                   >  1:1250944/1‑41 (MQ=255)
                                                                 aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGag       >  1:3187340/1‑69 (MQ=255)
                                                                 aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGag       >  1:1880926/1‑69 (MQ=255)
                                                                 aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGag       >  1:1664102/1‑69 (MQ=255)
                                                                   ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAcgcg                                    <  1:1714729/38‑1 (MQ=255)
                                                                   ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                             <  1:3574368/45‑1 (MQ=255)
                                                                   ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGt                            <  1:3518658/46‑1 (MQ=255)
                                                                   ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGt                            <  1:2793439/46‑1 (MQ=255)
                                                                   ttCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGt                            <  1:1607860/46‑1 (MQ=255)
                                                                    tCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGagcagc   >  1:3747784/1‑70 (MQ=255)
                                                                     cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAg                                  >  1:1456094/1‑38 (MQ=255)
                                                                     cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCg                             >  1:2066152/1‑43 (MQ=255)
                                                                     cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg              <  1:313591/58‑1 (MQ=255)
                                                                     cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg              <  1:276540/58‑1 (MQ=255)
                                                                     cTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGcc  <  1:2799660/70‑1 (MQ=255)
                                                                     |                                                                     
AGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCC  >  minE/1439399‑1439537

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