Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,835 T→C 8.0% D156G (GAC→GGC) ‡ ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8350TC 8.0% 338.6 / 16.4 124D156G (GAC→GGC) ‡ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base T (65/49);  new base C (2/8);  total (67/57)
Fisher's exact test for biased strand distribution p-value = 4.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 125 total aligned reads displayed.

CCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACC  >  minE/1439766‑1439905
                                                                     |                                                                      
ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                                                       <  1:2850079/71‑1 (MQ=255)
     tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                                  <  1:1767127/71‑1 (MQ=255)
        gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                                                               >  1:2273931/1‑71 (MQ=255)
        gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                                                               >  1:3023141/1‑71 (MQ=255)
          tGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                             >  1:1771873/1‑71 (MQ=255)
           gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                             >  1:2088991/1‑70 (MQ=255)
            tcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga                                                             >  1:3865664/2‑68 (MQ=2)
             cGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                            <  1:3047435/69‑1 (MQ=255)
             cGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAg                                                          >  1:629630/1‑71 (MQ=255)
             cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                                                            >  1:2134175/1‑67 (MQ=255)
             cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACAATTTACCAGACCGAAGAAATCAGTATTACgat                                                            >  1:3357153/1‑67 (MQ=255)
              gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                                                         >  1:2133154/1‑71 (MQ=255)
              gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                                                         >  1:2364191/1‑71 (MQ=255)
               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                            <  1:2949539/67‑1 (MQ=255)
               tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                        >  1:1662705/1‑71 (MQ=255)
                ttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                                  >  1:2405916/1‑60 (MQ=255)
                 gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                                                       <  1:2924173/69‑1 (MQ=255)
                   tACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                                                     <  1:2593206/70‑1 (MQ=255)
                    aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg                                                    <  1:3891873/70‑1 (MQ=255)
                    aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg                                                    <  1:1926335/70‑1 (MQ=255)
                    aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                            <  1:1806140/62‑1 (MQ=255)
                     cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                    <  1:2765595/69‑1 (MQ=255)
                     cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                    >  1:1247060/1‑69 (MQ=255)
                     cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                    >  1:202365/1‑69 (MQ=255)
                     cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGa                                                   <  1:1979738/70‑1 (MQ=255)
                     cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa                                                   <  1:212156/70‑3 (MQ=255)
                     cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa                                                   <  1:3432106/70‑3 (MQ=255)
                     cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgca                                                    <  1:3448962/69‑2 (MQ=255)
                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                            >  1:1132752/1‑60 (MQ=255)
                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                            >  1:320875/1‑60 (MQ=255)
                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                 <  1:2413742/71‑1 (MQ=255)
                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                 <  1:3639299/71‑1 (MQ=255)
                      ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                 <  1:3569742/71‑1 (MQ=255)
                      ccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt                                                        <  1:2829234/64‑1 (MQ=255)
                        tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                                               >  1:1207240/1‑71 (MQ=255)
                        tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                                               >  1:3453501/1‑71 (MQ=255)
                         ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg                                              >  1:2106683/1‑71 (MQ=255)
                         gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca                                              >  1:2334103/3‑70 (MQ=255)
                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                 >  1:2475107/1‑67 (MQ=255)
                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                             >  1:3323222/1‑71 (MQ=255)
                          gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                             <  1:2767405/71‑1 (MQ=255)
                          atACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca                                              <  1:2362503/69‑2 (MQ=255)
                           tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                            <  1:165070/71‑1 (MQ=255)
                           tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                            <  1:294004/71‑1 (MQ=255)
                            aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                           >  1:1549251/1‑70 (MQ=255)
                            aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                           >  1:1124359/1‑70 (MQ=255)
                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                                >  1:3378310/1‑49 (MQ=255)
                             cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                          <  1:1444877/71‑1 (MQ=255)
                                aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                                       >  1:3771385/1‑71 (MQ=255)
                                 aCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACg                                      <  1:2900722/71‑1 (MQ=255)
                                  cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg                                              >  1:3010504/1‑62 (MQ=255)
                                   aGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTTGCAACACCATTAGCACgac                                    >  1:1383453/1‑69 (MQ=255)
                                    gCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                                                <  1:2088577/42‑1 (MQ=255)
                                     cgaaGTTCAGACCATTTACCAGACCGAAGAAATCAGTAtt                                                                 <  1:876569/38‑1 (MQ=12)
                                      aaaGTTCAGGCCGTCATCCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                        <  1:2793821/48‑1 (MQ=255)
                                      aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                        <  1:3204583/48‑1 (MQ=255)
                                      aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                 >  1:288432/1‑71 (MQ=255)
                                      aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                 >  1:459440/1‑71 (MQ=255)
                                       aaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                             >  1:1433708/1‑42 (MQ=255)
                                        aGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGctgctg                                >  1:3007767/1‑70 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                                           >  1:1529600/1‑58 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                          >  1:2840793/1‑59 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                          >  1:650303/1‑59 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCtgctgc                               >  1:3128980/1‑70 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                              <  1:1545053/71‑1 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                              >  1:3775932/1‑71 (MQ=255)
                                         gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                              <  1:3667624/71‑1 (MQ=255)
                                         gTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                                                     >  1:129690/1‑48 (MQ=255)
                                          ttCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                             <  1:84361/55‑1 (MQ=255)
                                           tCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCAt                             >  1:3433500/1‑70 (MQ=255)
                                            cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                             >  1:3801564/1‑53 (MQ=255)
                                            cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACg                                      >  1:2870202/1‑60 (MQ=255)
                                            cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATg                            <  1:1066064/70‑1 (MQ=255)
                                            cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAACACGCTGCTGCATg                            <  1:2090055/70‑1 (MQ=255)
                                            agggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGCTCGCGAAGCCGt                                             >  1:847743/3‑53 (MQ=255)
                                             aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                                      >  1:419455/1‑43 (MQ=255)
                                              ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                         >  1:649939/1‑71 (MQ=255)
                                              ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                         <  1:2603702/71‑1 (MQ=255)
                                               gCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                    <  1:1387196/43‑1 (MQ=255)
                                               gCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGt                        >  1:3373374/1‑71 (MQ=255)
                                                ccGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTTCATGAAGtt                       <  1:1082344/71‑1 (MQ=255)
                                                 cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                      >  1:2465204/1‑71 (MQ=255)
                                                 cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                      >  1:1218480/1‑71 (MQ=255)
                                                 cGTCAACCAGACCGAAGAAGCCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                      >  1:401130/1‑71 (MQ=255)
                                                    cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                 >  1:3532626/1‑57 (MQ=255)
                                                    cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACg                                      >  1:290745/1‑52 (MQ=255)
                                                    cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt                     >  1:3202437/1‑69 (MQ=255)
                                                     aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                                           >  1:1595158/1‑46 (MQ=255)
                                                     aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGtt                       <  1:2642369/66‑1 (MQ=255)
                                                     aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAg                  <  1:2543525/71‑1 (MQ=255)
                                                      aCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGa                 <  1:2030951/71‑1 (MQ=255)
                                                       ccAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGaa                <  1:260037/71‑1 (MQ=255)
                                                           aCCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                             <  1:2410964/38‑1 (MQ=255)
                                                             cGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                         >  1:761707/1‑56 (MQ=255)
                                                             cGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAg          <  1:2579079/71‑1 (MQ=255)
                                                               aagaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgt        >  1:1657937/1‑71 (MQ=255)
                                                                 gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt      <  1:2064845/71‑1 (MQ=255)
                                                                 gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCCGCTGCATGAAGTTGTCAGAACCGTAGgtgt      >  1:1653934/1‑71 (MQ=255)
                                                                    gTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac   <  1:2843482/71‑1 (MQ=255)
                                                                     tCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcacc  <  1:582000/71‑1 (MQ=255)
                                                                     |                                                                      
CCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACC  >  minE/1439766‑1439905

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: