Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,857 A→C 9.5% F149V (TTC→GTC)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8570AC 9.5% 224.0 / 19.7 95F149V (TTC→GTC) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (47/39);  new base C (6/3);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGT  >  minE/1439788‑1439925
                                                                     |                                                                    
ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTATGGTAGGTCGCGAAg                                                                     <  1:541593/71‑1 (MQ=255)
ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                                     <  1:3569742/71‑1 (MQ=255)
ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                                     <  1:3639299/71‑1 (MQ=255)
ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                                     <  1:2413742/71‑1 (MQ=255)
  tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGc                                                                    >  1:2967278/4‑70 (MQ=9)
  tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                                                                   >  1:3453501/1‑71 (MQ=255)
  tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                                                                   >  1:1207240/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:277325/1‑71 (MQ=255)
    gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                                 >  1:2406348/1‑71 (MQ=255)
    gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                                 >  1:3323222/1‑71 (MQ=255)
    gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                                 <  1:2767405/71‑1 (MQ=255)
    atACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca                                                                  <  1:2362503/69‑2 (MQ=255)
    atACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca                                                                  <  1:956235/69‑2 (MQ=255)
     tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                                                <  1:294004/71‑1 (MQ=255)
     tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                                                <  1:165070/71‑1 (MQ=255)
      aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                                               >  1:1124359/1‑70 (MQ=255)
      aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                                               >  1:1549251/1‑70 (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)
                aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                                     >  1:288432/1‑71 (MQ=255)
                aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                                     >  1:459440/1‑71 (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:3667624/71‑1 (MQ=255)
                   gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                                  >  1:3775932/1‑71 (MQ=255)
                   gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                                  >  1:1738353/1‑71 (MQ=255)
                   gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                                  <  1:1545053/71‑1 (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:2113501/1‑70 (MQ=255)
                      cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATg                                                <  1:1066064/70‑1 (MQ=255)
                      cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAACACGCTGCTGCATg                                                <  1:2090055/70‑1 (MQ=255)
                      agggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGCTCGCGAAGCCGt                                                                 >  1:847743/3‑53 (MQ=255)
                        ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                             >  1:649939/1‑71 (MQ=255)
                        ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                             >  1:1162079/1‑71 (MQ=255)
                        ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                             <  1:2603702/71‑1 (MQ=255)
                         gCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGt                                            >  1:3373374/1‑71 (MQ=255)
                          ccGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTTCATGAAGtt                                           <  1:1082344/71‑1 (MQ=255)
                           cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCa                                                            >  1:3421887/1‑53 (MQ=255)
                           cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                                          >  1:1218480/1‑71 (MQ=255)
                           cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                                          >  1:2465204/1‑71 (MQ=255)
                           cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                                          <  1:3459522/71‑1 (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)
                                  ctaaCCAAAGAAGCCATTGTTACGGTAGGTTGCAACACCATTAGCACGACCGTTCATGAAGtt                                           <  1:3884808/60‑1 (MQ=255)
                                    gACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCg                                 >  1:2087700/1‑71 (MQ=255)
                                     aCCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                                 <  1:2410964/38‑1 (MQ=255)
                                       cGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                             >  1:761707/1‑56 (MQ=255)
                                       cGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAg                              <  1:2579079/71‑1 (MQ=255)
                                       cGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAg                              <  1:2905831/71‑1 (MQ=255)
                                         aagaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgt                            >  1:1657937/1‑71 (MQ=255)
                                           gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt                          <  1:2509336/71‑1 (MQ=255)
                                           gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt                          <  1:2064845/71‑1 (MQ=255)
                                           gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt                          >  1:3738513/1‑71 (MQ=255)
                                           gaagTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCCGCTGCATGAAGTTGTCAGAACCGTAGgtgt                          >  1:1653934/1‑71 (MQ=255)
                                              gTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac                       <  1:2843482/71‑1 (MQ=255)
                                               tCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcacc                      <  1:582000/71‑1 (MQ=255)
                                                   tgtTACGGTAGGTTGCAACACCATTAGCACGACCGTTCATGAAGTTGt                                         >  1:326769/1‑48 (MQ=255)
                                                     ttACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAg                              <  1:823178/57‑1 (MQ=255)
                                                     ttACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGa                 <  1:864994/70‑1 (MQ=255)
                                                      tACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt                                         >  1:1779460/1‑45 (MQ=255)
                                                       aCGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcacca                     >  1:490852/1‑64 (MQ=255)
                                                       aCGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAAtt              >  1:258360/1‑71 (MQ=255)
                                                       aCGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAAt               >  1:3208506/1‑70 (MQ=255)
                                                        cGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAAtt              <  1:2900499/70‑1 (MQ=255)
                                                         ggTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAAtt              <  1:1727357/69‑1 (MQ=255)
                                                           tAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTgg          >  1:3724823/1‑71 (MQ=255)
                                                              gTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAg       <  1:219735/71‑1 (MQ=255)
                                                              gTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAAACGTAGGTGTCACCACCGAATTCTGGCAg       <  1:3893317/71‑1 (MQ=255)
                                                                cgcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcacc                      >  1:632700/1‑54 (MQ=255)
                                                                cgcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTa     >  1:1800746/1‑71 (MQ=255)
                                                                 gcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtg                           <  1:2852761/48‑1 (MQ=255)
                                                                 gcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTg           <  1:585853/64‑1 (MQ=255)
                                                                 gcgAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTAc    <  1:2274427/71‑1 (MQ=255)
                                                                   gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt  <  1:2873946/71‑1 (MQ=255)
                                                                   gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt  <  1:94569/71‑1 (MQ=255)
                                                                   gAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt  >  1:1132895/1‑71 (MQ=255)
                                                                     aGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGgtgt                          >  1:2183563/1‑45 (MQ=255)
                                                                     aGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac                       >  1:1729746/1‑48 (MQ=255)
                                                                     aGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTcac                       >  1:1735836/1‑48 (MQ=255)
                                                                     aGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGt  >  1:2782432/1‑69 (MQ=255)
                                                                     |                                                                    
CCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGT  >  minE/1439788‑1439925

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: