Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,444 G→A 8.7% Y286Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4440GA 8.7% 181.1 / 5.0 57Y286Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTG  >  minE/1439375‑1439513
                                                                     |                                                                     
gCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTgg                                                                       <  1:1418986/70‑1 (MQ=255)
gCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTgg                                                                       <  1:2790721/70‑1 (MQ=255)
       cccAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGc                                                               <  1:186307/71‑1 (MQ=255)
         cAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCaa                                                             <  1:2900696/71‑1 (MQ=255)
          aGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCaa                                                             >  1:2467765/1‑70 (MQ=255)
           ggTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCaa                                                             >  1:2653260/1‑69 (MQ=255)
            gTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAg                                                          >  1:1684937/1‑71 (MQ=255)
              ttttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCt                                                        >  1:3009260/1‑71 (MQ=255)
               tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGcg                                                        <  1:1836276/70‑2 (MQ=255)
               tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCtt                                                       >  1:357376/1‑71 (MQ=255)
               tttACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCt                                                        <  1:855093/70‑1 (MQ=255)
                  aCCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                                     <  1:4000594/70‑1 (MQ=255)
                   ccTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCaa                                                             >  1:3401434/1‑61 (MQ=255)
                   ccTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGa                                                    <  1:1426935/70‑1 (MQ=255)
                       tAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAg                                                  >  1:1024019/1‑68 (MQ=255)
                        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)
                                ggTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGt                                                 <  1:261651/60‑1 (MQ=255)
                                 gTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCaa                                                             <  1:276747/47‑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)
                                   aaGCCAGGGACGGACGCAGACCGAAGTCGAACTgg                                                                       >  1:833724/1‑35 (MQ=255)
                                   aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAGGTTCTGTGCTTTGTTc                                   >  1:1056659/1‑71 (MQ=255)
                                   aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                   >  1:3070839/1‑71 (MQ=255)
                                   aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                   >  1:1658874/1‑71 (MQ=255)
                                      ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                    <  1:4109220/67‑1 (MQ=255)
                                       cAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGt                                                 <  1:1909043/53‑1 (MQ=255)
                                       cAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGt                                                 <  1:3210596/53‑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)
                                            acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                   >  1:123236/1‑46 (MQ=255)
                                            acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                   >  1:782062/1‑46 (MQ=255)
                                            acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                   >  1:3875356/1‑46 (MQ=255)
                                              ggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                                     <  1:1778484/42‑1 (MQ=255)
                                               gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                       <  1:171568/71‑1 (MQ=255)
                                               gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                       >  1:3998221/1‑71 (MQ=255)
                                                 cgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGa                     <  1:1212527/71‑1 (MQ=255)
                                                  gCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                   <  1:2218482/40‑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: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)
                                                                     |                                                                     
GCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTG  >  minE/1439375‑1439513

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: