Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,834 G→A 6.0% D156D (GAC→GAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8340GA 6.0% 389.4 / 9.3 116D156D (GAC→GATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (58/51);  new base A (3/4);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 127 total aligned reads displayed.

CACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCAC  >  minE/1439764‑1439904
                                                                      |                                                                      
 aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                                       >  1:2897085/1‑71 (MQ=255)
 aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                                       <  1:2022812/71‑1 (MQ=255)
  ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                                                      <  1:2533883/71‑1 (MQ=255)
    aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                                    >  1:375926/1‑71 (MQ=255)
    aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAAAAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                                    <  1:1163407/71‑1 (MQ=255)
       tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                                 >  1:856898/1‑71 (MQ=255)
       tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                                 >  1:3487694/1‑71 (MQ=255)
         ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                                                                <  1:1277748/70‑1 (MQ=255)
          gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                                    <  1:3481406/65‑1 (MQ=255)
            cttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACg                                                             <  1:56719/67‑1 (MQ=2)
             gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            >  1:829369/1‑70 (MQ=255)
             gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            >  1:3802273/1‑70 (MQ=255)
             gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            <  1:1424771/70‑1 (MQ=255)
                 tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                               <  1:1841386/63‑1 (MQ=255)
                 tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            >  1:510827/1‑66 (MQ=255)
                 tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                       >  1:3220468/1‑71 (MQ=255)
                  ttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                                   >  1:2695949/1‑58 (MQ=255)
                   tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTc                                                      >  1:214731/1‑70 (MQ=255)
                   attACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGttg                                                     >  1:3924542/2‑69 (MQ=255)
                    ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                                                    <  1:2411618/71‑1 (MQ=9)
                    ttACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                                    <  1:3110226/71‑1 (MQ=255)
                    ttACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                                    <  1:1804202/71‑1 (MQ=255)
                    ttACACTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc                                                    >  1:1130816/1‑71 (MQ=255)
                     tACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc                                                    >  1:1937053/1‑70 (MQ=255)
                     tACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                                    >  1:3849331/1‑70 (MQ=255)
                       cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                                 <  1:382139/55‑1 (MQ=255)
                       cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa                                                 <  1:3749021/71‑1 (MQ=255)
                        ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                <  1:3279985/71‑1 (MQ=255)
                         cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                               >  1:2445327/1‑55 (MQ=255)
                         cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                                               <  1:2598011/71‑1 (MQ=255)
                          tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGt                                                       >  1:3415493/4‑62 (MQ=2)
                          tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                                       >  1:3092619/1‑46 (MQ=255)
                           gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGa                                                  <  1:695181/64‑1 (MQ=255)
                            gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            >  1:2716471/1‑55 (MQ=255)
                            gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            >  1:3478811/1‑55 (MQ=255)
                            gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGGCGCGAAGcc                                              <  1:2959033/69‑1 (MQ=255)
                            gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGAGTTACGGTAGGTCGCGAAGCCGt                                            >  1:1755088/1‑71 (MQ=255)
                              aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                                   >  1:3979937/1‑46 (MQ=255)
                              aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                                    >  1:1692942/1‑45 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         >  1:3807269/1‑71 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         >  1:3251623/1‑71 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         <  1:1489981/71‑1 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         >  1:1429344/1‑71 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         <  1:2004825/71‑1 (MQ=255)
                               cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         <  1:3002901/71‑1 (MQ=255)
                                tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                                    <  1:1233672/43‑1 (MQ=255)
                                tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            <  1:2347834/51‑1 (MQ=255)
                                tGAACAGCAAAGTTCAGCCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            <  1:2017753/51‑1 (MQ=255)
                                 gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                            <  1:1128244/50‑1 (MQ=255)
                                 gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAg                                                         >  1:1760466/1‑53 (MQ=255)
                                 gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                                               >  1:539170/1‑63 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCACCCAGACCGAAGAAGTCAgtg                                                                  >  1:1301783/1‑43 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtg                                                                  >  1:33904/1‑43 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                               >  1:896892/1‑46 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                                     >  1:2840096/1‑56 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                                     >  1:3684352/1‑56 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                                          >  1:2715422/1‑67 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                                       >  1:791053/1‑70 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                                       <  1:610035/70‑1 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                                       <  1:3859286/70‑1 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                                      >  1:1326101/1‑71 (MQ=255)
                                  aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                                      >  1:3444189/1‑71 (MQ=255)
                                    cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                                    >  1:3674322/1‑71 (MQ=255)
                                    cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                                    >  1:3726392/1‑71 (MQ=255)
                                      gCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                               >  1:795070/1‑42 (MQ=255)
                                       cAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                                               <  1:1241078/57‑1 (MQ=255)
                                        aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                >  1:2463292/1‑71 (MQ=255)
                                           gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                                              <  1:4078005/54‑1 (MQ=255)
                                           gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgct                                   <  1:3670763/65‑1 (MQ=255)
                                           gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                             >  1:1941883/1‑71 (MQ=255)
                                            ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTa                                                          >  1:455223/1‑41 (MQ=255)
                                               aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                       <  1:81904/41‑1 (MQ=255)
                                               aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGa                          <  1:1720907/70‑1 (MQ=255)
                                               aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGa                          <  1:3871500/70‑1 (MQ=255)
                                                  ccGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGtt                      <  1:3748706/71‑1 (MQ=255)
                                                   cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                     <  1:1330184/71‑1 (MQ=255)
                                                     tCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTc                   <  1:2383714/71‑1 (MQ=255)
                                                      cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc                                               >  1:2060136/1‑42 (MQ=255)
                                                      cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGt                       >  1:1703648/1‑66 (MQ=255)
                                                      cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt                    >  1:386393/1‑69 (MQ=255)
                                                      cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGAATGAAGt                       >  1:2824373/1‑66 (MQ=255)
                                                       aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                         >  1:2956460/1‑47 (MQ=255)
                                                       aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt                    >  1:1665230/1‑68 (MQ=255)
                                                        aCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                     <  1:4121975/66‑1 (MQ=255)
                                                        aCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                     <  1:1256100/66‑1 (MQ=255)
                                                         ccAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGaa               >  1:3659485/1‑71 (MQ=255)
                                                          cAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                                    >  1:89399/1‑49 (MQ=255)
                                                           aGACCAAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAAcc             <  1:3266978/71‑1 (MQ=255)
                                                            gACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCg            <  1:1317516/71‑1 (MQ=255)
                                                              ccGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTa          <  1:1052370/71‑1 (MQ=255)
                                                                gaagaaGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAg         <  1:2306325/70‑1 (MQ=255)
                                                                 aagaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                                          >  1:2622520/1‑36 (MQ=255)
                                                                 aagaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgt       <  1:246267/71‑1 (MQ=255)
                                                                  agaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCa                  >  1:2705833/1‑59 (MQ=255)
                                                                  agaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCg            >  1:3081415/1‑65 (MQ=255)
                                                                  agaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtg      <  1:3156836/71‑1 (MQ=255)
                                                                   gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt     >  1:293570/1‑71 (MQ=255)
                                                                    aagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTc    <  1:2610395/71‑1 (MQ=255)
                                                                      gTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac  <  1:3990502/71‑1 (MQ=255)
                                                                      gTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac  <  1:330694/71‑1 (MQ=255)
                                                                      |                                                                      
CACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCAC  >  minE/1439764‑1439904

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: