Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,792 G→A 6.7% Y170Y (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,7920GA 6.7% 292.6 / 7.5 89Y170Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (51/32);  new base A (2/4);  total (53/36)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

CGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGT  >  minE/1439725‑1439862
                                                                   |                                                                      
cGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                      <  1:3797834/70‑1 (MQ=255)
cGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTAc                                                                     >  1:4069528/1‑71 (MQ=255)
cGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTAc                                                                     >  1:1155823/1‑71 (MQ=255)
cGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTAc                                                                     >  1:3203480/1‑71 (MQ=255)
cGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTAc                                                                     >  1:433690/1‑71 (MQ=255)
  tCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                   <  1:910524/71‑1 (MQ=255)
  tCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                   <  1:117296/71‑1 (MQ=255)
   cACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                   >  1:1626694/1‑70 (MQ=255)
    aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                    >  1:3738870/1‑68 (MQ=255)
    aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                 >  1:1955896/1‑71 (MQ=255)
    aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                 >  1:3684169/1‑71 (MQ=255)
    aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                 >  1:426148/1‑71 (MQ=255)
    aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGCGTTGCCGTTTTTACCCTGGTACTGaa                                                                 >  1:456000/1‑71 (MQ=255)
       aCCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                                                               <  1:1147877/70‑1 (MQ=255)
         cGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                            <  1:2041681/71‑1 (MQ=255)
           ttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                            >  1:3323924/1‑69 (MQ=255)
              ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                       <  1:205604/71‑1 (MQ=255)
              ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt                                                        <  1:2678959/70‑1 (MQ=255)
               tagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                       >  1:297701/1‑70 (MQ=255)
                agttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTg                                                                   <  1:291908/57‑1 (MQ=255)
                    aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                 >  1:2159222/1‑71 (MQ=255)
                    aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                 >  1:2028957/1‑71 (MQ=255)
                     cGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                 <  1:2429539/70‑1 (MQ=255)
                       ccACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                      <  1:2652228/47‑1 (MQ=255)
                        cACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                    >  1:1060351/1‑48 (MQ=255)
                         aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTc                                              <  1:3982698/69‑1 (MQ=255)
                         aCTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                            <  1:2590613/71‑1 (MQ=255)
                          cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                           >  1:511295/1‑71 (MQ=255)
                          cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                           <  1:1435957/71‑1 (MQ=255)
                           tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                                                               <  1:145155/50‑1 (MQ=255)
                           tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                           >  1:3863923/1‑70 (MQ=255)
                           tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                           >  1:3321781/1‑70 (MQ=255)
                              tAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaa                                                           <  1:2296268/51‑1 (MQ=255)
                              tAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                       <  1:4068518/55‑1 (MQ=255)
                                aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                      >  1:1481632/1‑38 (MQ=255)
                                aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                     >  1:2133078/1‑55 (MQ=255)
                                aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGa                                       >  1:814612/1‑69 (MQ=255)
                                aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGAcc                                     >  1:2521302/1‑71 (MQ=255)
                                   ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                          >  1:2042108/1‑63 (MQ=255)
                                   ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAcc                                          >  1:3573862/1‑63 (MQ=255)
                                   ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                  <  1:4038060/71‑1 (MQ=255)
                                   ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                  <  1:939048/71‑1 (MQ=255)
                                       cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTc                                                      >  1:895348/1‑47 (MQ=255)
                                       cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                            >  1:3797908/1‑57 (MQ=255)
                                       cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                            >  1:1057791/1‑57 (MQ=255)
                                       cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                              >  1:4096180/1‑71 (MQ=255)
                                        aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                             >  1:2956021/1‑71 (MQ=255)
                                        aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                             <  1:2022812/71‑1 (MQ=255)
                                        aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                             >  1:2897085/1‑71 (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)
                                               gggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAgggc                                                 >  1:1191602/1‑42 (MQ=255)
                                               gggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                    >  1:2163995/1‑41 (MQ=255)
                                               gggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                       >  1:1224135/1‑70 (MQ=255)
                                                ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                      <  1:1277748/70‑1 (MQ=255)
                                                 gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                          <  1:3481406/65‑1 (MQ=255)
                                                 gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                    >  1:2736074/1‑71 (MQ=255)
                                                  ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                         <  1:956409/49‑1 (MQ=255)
                                                   cttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACg                   <  1:56719/67‑1 (MQ=2)
                                                    gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  <  1:3717165/70‑1 (MQ=255)
                                                    gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  <  1:1424771/70‑1 (MQ=255)
                                                    gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  >  1:829369/1‑70 (MQ=255)
                                                    gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  >  1:3802273/1‑70 (MQ=255)
                                                       gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                               >  1:3198704/1‑54 (MQ=255)
                                                       gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaagaa                               >  1:3146695/1‑54 (MQ=255)
                                                        tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                           <  1:1887122/41‑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:1804202/71‑1 (MQ=255)
                                                           ttACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc          <  1:3110226/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:1374860/1‑55 (MQ=255)
                                                                cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                     >  1:2445327/1‑55 (MQ=255)
                                                                cTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc     <  1:2598011/71‑1 (MQ=255)
                                                                 tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                             >  1:3092619/1‑46 (MQ=255)
                                                                 tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                             >  1:3023380/1‑46 (MQ=255)
                                                                  ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGaa                                      >  1:110833/1‑35 (MQ=255)
                                                                   gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  >  1:3478811/1‑55 (MQ=255)
                                                                   gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                  >  1:2716471/1‑55 (MQ=255)
                                                                   gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGGCGCGAAGcc    <  1:2959033/69‑1 (MQ=255)
                                                                   gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGAGTTACGGTAGGTCGCGAAGCCGt  >  1:1755088/1‑71 (MQ=255)
                                                                   |                                                                      
CGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGT  >  minE/1439725‑1439862

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: