Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,462 G→A 6.3% F280F (TTC→TTT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,4620GA 6.3% 262.6 / 5.2 79F280F (TTC→TTTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (42/32);  new base A (2/3);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTA  >  minE/1439393‑1439527
                                                                     |                                                                 
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:1658874/1‑71 (MQ=255)
                 aaGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTc                                                 >  1:3070839/1‑71 (MQ=255)
                    ccAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGtt                                                  <  1:4109220/67‑1 (MQ=255)
                     cAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGt                                                               <  1:3210596/53‑1 (MQ=255)
                     cAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGt                                                               <  1:1909043/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: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: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:3648082/70‑1 (MQ=255)
                                                   gTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTg                <  1:681393/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:686870/1‑41 (MQ=255)
                                                           cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAgg                                     >  1:597374/1‑41 (MQ=255)
                                                           cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCt         >  1:3734296/1‑69 (MQ=255)
                                                           cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg       <  1:3340915/71‑1 (MQ=255)
                                                           cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg       <  1:2002852/71‑1 (MQ=255)
                                                           cAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTgg       <  1:1491548/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:4096178/1‑59 (MQ=255)
                                                                gCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTa              >  1:1884437/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)
                                                                     |                                                                 
ACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCGAACTGGTACTGAGCAACAGCTTCGAAGTTCTGTGCTTTGTTCGCCCAACCCAGGGAACCTACGCGAGTTGCGTTGTAGGTCTGGGTGTA  >  minE/1439393‑1439527

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: