Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,471 T→G 6.1% A277A (GCA→GCC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4710TG 6.1% 225.2 / 9.0 82A277A (GCA→GCCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base T (44/33);  new base G (2/3);  total (46/36)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGG  >  minE/1439403‑1439540
                                                                    |                                                                     
tGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGc                                                                     >  1:1274405/1‑71 (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:1419433/62‑1 (MQ=255)
             gggACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                 <  1:1160238/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:192512/1‑47 (MQ=255)
                              aaGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                               >  1:4106496/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:3846794/49‑1 (MQ=16)
                               aGTCGAACTGATACTGAGCAACAACTTCAAAGTTTTGGGTTTTGTTAGc                                                            <  1:859923/49‑1 (MQ=16)
                                 tCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTg                                                                 <  1:2165320/42‑1 (MQ=255)
                                   gAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTg                                  >  1:669401/1‑71 (MQ=255)
                                   gAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTg                                  >  1:2052527/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:1684308/1‑71 (MQ=255)
                                               aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt                      >  1:493184/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:3340915/71‑1 (MQ=255)
                                                 cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                    <  1:1491548/71‑1 (MQ=255)
                                                 cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg                    <  1:2002852/71‑1 (MQ=255)
                                                  aaCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAg                                                   >  1:2134187/1‑39 (MQ=255)
                                                   aCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgt                  >  1:445605/1‑71 (MQ=255)
                                                   aCAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgt                  >  1:1253464/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:1880926/1‑69 (MQ=255)
                                                             aGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGag          >  1:3187340/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)
                                                                    tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGt                            <  1:1048692/44‑1 (MQ=255)
                                                                    tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAgg  >  1:1566981/1‑70 (MQ=255)
                                                                    tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAgg  >  1:4073929/1‑70 (MQ=255)
                                                                    tgCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAgg  <  1:1850245/70‑1 (MQ=255)
                                                                    |                                                                     
TGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTACTGAGCAGCCAGG  >  minE/1439403‑1439540

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: