Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,846 G→A 9.7% Y152Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8460GA 9.7% 327.1 / 18.6 113Y152Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (55/47);  new base A (5/6);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 117 total aligned reads displayed.

TGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATT  >  minE/1439776‑1439913
                                                                      |                                                                   
tGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                                     >  1:1771873/1‑71 (MQ=255)
   cGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                                    <  1:3047435/69‑1 (MQ=255)
   cGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAg                                                                  >  1:629630/1‑71 (MQ=255)
    gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                                                                 >  1:2364191/1‑71 (MQ=255)
    gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                                                                 >  1:2133154/1‑71 (MQ=255)
     tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                                    <  1:2949539/67‑1 (MQ=255)
     tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                                >  1:1662705/1‑71 (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)
           cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                            <  1:2765595/69‑1 (MQ=255)
           cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                            >  1:202365/1‑69 (MQ=255)
           cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                            >  1:1247060/1‑69 (MQ=255)
           cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGa                                                           <  1:1979738/70‑1 (MQ=255)
           cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaac                                                          <  1:1605507/71‑4 (MQ=255)
           cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaac                                                          <  1:2843971/71‑4 (MQ=255)
           cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa                                                           <  1:212156/70‑3 (MQ=255)
           cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa                                                           <  1:3432106/70‑3 (MQ=255)
            ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTATGGTAGGTCGCGAAg                                                         <  1:541593/71‑1 (MQ=255)
            ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                                                                    >  1:1132752/1‑60 (MQ=255)
            ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                         <  1:3639299/71‑1 (MQ=255)
            ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                         <  1:2413742/71‑1 (MQ=255)
            ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg                                                         <  1:3569742/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:2767405/71‑1 (MQ=255)
                gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                     >  1:2406348/1‑71 (MQ=255)
                gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                     >  1:3323222/1‑71 (MQ=255)
                gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                                                     >  1:277325/1‑71 (MQ=255)
                 tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                                    <  1:294004/71‑1 (MQ=255)
                 tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                                                    <  1:165070/71‑1 (MQ=255)
                  aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                                   >  1:1549251/1‑70 (MQ=255)
                  aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg                                                   >  1:1124359/1‑70 (MQ=255)
                   cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                                  <  1:1444877/71‑1 (MQ=255)
                      aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                                               >  1:3771385/1‑71 (MQ=255)
                        cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg                                                      >  1:3010504/1‑62 (MQ=255)
                         aGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTTGCAACACCATTAGCACgac                                            >  1:1383453/1‑69 (MQ=255)
                            aaaGTTCAGGCCGTCATCCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                                <  1:2793821/48‑1 (MQ=255)
                            aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                                                <  1:3204583/48‑1 (MQ=255)
                            aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                         >  1:459440/1‑71 (MQ=255)
                            aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                                         >  1:288432/1‑71 (MQ=255)
                             aaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                                                     >  1:1433708/1‑42 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTa                                                   >  1:1529600/1‑58 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc                                                  >  1:2840793/1‑59 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCtgctgc                                       >  1:3128980/1‑70 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                      >  1:3775932/1‑71 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                      >  1:1738353/1‑71 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                      <  1:3667624/71‑1 (MQ=255)
                               gTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCa                                      <  1:1545053/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: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)
                                   aGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                                              >  1:2600455/1‑43 (MQ=255)
                                    ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                 >  1:649939/1‑71 (MQ=255)
                                    ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg                                 >  1:1162079/1‑71 (MQ=255)
                                     gCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                                                            <  1:1387196/43‑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:3459522/71‑1 (MQ=255)
                                       cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg                              >  1:2465204/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)
                                              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:2064845/71‑1 (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)
                                                                      |                                                                   
TGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATT  >  minE/1439776‑1439913

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: