Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,457 G→A 6.4% A282V (GCT→GTT)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4570GA 6.4% 260.2 / 5.3 78A282V (GCT→GTT) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base A (2/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTA  >  minE/1439387‑1439527
                                                                      |                                                                      
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)
       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)
                     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: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:782062/1‑46 (MQ=255)
                                acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                                 >  1:123236/1‑46 (MQ=255)
                                acggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGaa                                                                 >  1:3875356/1‑46 (MQ=255)
                                  ggacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCg                                                                   <  1:1778484/42‑1 (MQ=255)
                                   gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                     >  1:3998221/1‑71 (MQ=255)
                                   gacgCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                     <  1:171568/71‑1 (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: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: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:2111609/58‑1 (MQ=255)
                                                        ggTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTAc                             <  1:379509/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:594248/1‑71 (MQ=255)
                                                               aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt         >  1:1684308/1‑71 (MQ=255)
                                                               aGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt         >  1:493184/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:2002852/71‑1 (MQ=255)
                                                                 cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg       <  1:3340915/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:3738859/71‑1 (MQ=255)
                                                                    cAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGgtg    <  1:1594091/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:4096178/1‑59 (MQ=255)
                                                                      gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTa              >  1:1884437/1‑59 (MQ=255)
                                                                      gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTa  <  1:2900689/71‑1 (MQ=255)
                                                                      |                                                                      
GTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTA  >  minE/1439387‑1439527

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: