Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,356,166 G→C 5.0% L154V (CTA→GTA)  yehA ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,356,1660GC 5.0% 262.8 / 3.2 79L154V (CTA→GTA) yehApredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base C (2/2);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCAT  >  minE/1356100‑1356235
                                                                  |                                                                     
gtgCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGaaa                                                                    <  1:1289140/70‑1 (MQ=255)
   ccAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTg                                                                 >  1:20913/1‑70 (MQ=255)
   ccAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTg                                                                 <  1:2228833/70‑1 (MQ=255)
     aaTGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGaaa                                                                    <  1:2140971/65‑1 (MQ=255)
     aaTGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTAc                                                              <  1:2351861/71‑1 (MQ=255)
      aTGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGaaa                                                                    <  1:233423/64‑1 (MQ=255)
        gCATTTCTGTAGCTATTTACACATGCATCCCATGAGTTGCCACCTAATAAACACAGTAGAAAt                                                                   >  1:655185/1‑63 (MQ=255)
        gCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAAt                                                                   >  1:1958561/1‑63 (MQ=255)
         cATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGaaa                                                                    >  1:810472/1‑61 (MQ=255)
           tttCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGt                                                                <  1:2944320/63‑1 (MQ=255)
            ttCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCa                                                       >  1:533468/1‑71 (MQ=255)
             tCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCa                                                       >  1:2384146/1‑70 (MQ=255)
             tCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCAc                                                      >  1:488702/1‑71 (MQ=255)
             tCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCAc                                                      >  1:1491498/1‑71 (MQ=255)
             tCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCAc                                                      >  1:1958222/1‑71 (MQ=255)
             tCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCAc                                                      >  1:1030773/1‑71 (MQ=255)
              cTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCAc                                                      <  1:674544/70‑1 (MQ=255)
              cTGTAGCTATTTACACATGCATCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTc                                                    >  1:297341/1‑71 (MQ=255)
                 tAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACt                                                     <  1:1163974/68‑1 (MQ=255)
                   gCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCcgc                                                >  1:1067019/1‑71 (MQ=255)
                    cTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCcgcg                                               >  1:2455474/1‑71 (MQ=255)
                      aTTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCgtg                                             >  1:2929546/1‑71 (MQ=255)
                      aTTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCgtg                                             >  1:1116607/1‑71 (MQ=255)
                       tttACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAAt                                                                   <  1:2187096/48‑1 (MQ=255)
                       tttACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCgt                                              <  1:1846071/69‑1 (MQ=255)
                         tACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTc                                           >  1:2587130/1‑70 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                          >  1:2259948/1‑54 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                          >  1:934749/1‑54 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                          >  1:756941/1‑54 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCg                                                          >  1:760459/1‑54 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCg                                          <  1:1313823/70‑1 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCg                                          <  1:846364/70‑1 (MQ=255)
                          acacATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCg                                          <  1:1256894/70‑1 (MQ=255)
                           cacaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCgt                                              <  1:1134271/65‑1 (MQ=255)
                            acaTGCATCCCATGACTTGCCACCTAATAAACGCAGTACaatt                                                                   <  1:468601/43‑3 (MQ=38)
                            acaTGCATCCCATGACTTGCCACCTAATAAACGCAGTACaatt                                                                   <  1:123414/43‑3 (MQ=38)
                            acaTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCgatg                                       <  1:1967017/71‑1 (MQ=255)
                              gtGCATCCCATGACTTGCCACCTAATAAACACAGTACAAtttc                                                                 >  1:2869799/2‑42 (MQ=255)
                              aTGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCgatgat                                     >  1:1050815/1‑71 (MQ=255)
                              aTGCATCCCATGACTTGCCACCTAATAAACACAGTACAAtttc                                                                 >  1:1508034/1‑42 (MQ=255)
                               tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                    <  1:3024132/71‑1 (MQ=255)
                               tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                    <  1:2935425/71‑1 (MQ=255)
                               tGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                    <  1:1171239/71‑1 (MQ=255)
                                 cATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                    >  1:2148252/1‑69 (MQ=255)
                                  aTCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTg                                                            <  1:842924/44‑1 (MQ=255)
                                    cccATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGatgatg                                    <  1:2878202/66‑1 (MQ=255)
                                     ccATGACTTGCCACCTAATAAACACAGTAGAAATTGTa                                                               >  1:658419/1‑38 (MQ=255)
                                     ccATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCgg                                                         <  1:1854079/44‑1 (MQ=255)
                                     ccATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCAc                              >  1:590283/1‑71 (MQ=255)
                                      cATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCgatg                                       <  1:2875805/61‑1 (MQ=255)
                                        tGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGaa                                  >  1:1551197/1‑64 (MQ=255)
                                          aCTTGCCACCTAATAAACACATTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGAt                         <  1:1238988/71‑1 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCgtgt                                            >  1:2171089/1‑52 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCga                                         >  1:2600324/1‑55 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCga                                         >  1:1080423/1‑55 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCga                                         >  1:2467030/1‑55 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTg                            >  1:2488579/1‑68 (MQ=255)
                                          aCTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGa                          >  1:857630/1‑70 (MQ=255)
                                            ttGCCACCTAATAAACACAGTAGAAATTGTACTGc                                                           <  1:291521/35‑1 (MQ=255)
                                               ccACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGt                    <  1:2399214/71‑1 (MQ=255)
                                                cACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGAt                         <  1:382908/65‑1 (MQ=255)
                                                  ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCg                  >  1:2174184/1‑70 (MQ=255)
                                                  ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCg                  >  1:3032580/1‑70 (MQ=255)
                                                  ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCg                  >  1:1879206/1‑70 (MQ=255)
                                                  ccTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCg                  >  1:468241/1‑70 (MQ=255)
                                                   cTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCa                <  1:904004/71‑1 (MQ=255)
                                                     aataaACACAGTAGAAATTGTACTGCGGCACTCCCGCgtgt                                            >  1:3051317/1‑41 (MQ=255)
                                                       taaACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAacac            <  1:264785/71‑1 (MQ=255)
                                                        aaaCACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCgatga                                      >  1:2653440/1‑44 (MQ=255)
                                                         aaCACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCa          >  1:1442937/1‑71 (MQ=255)
                                                          acacAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGaa                                  >  1:1512177/1‑46 (MQ=255)
                                                              aGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCa                               >  1:810197/1‑45 (MQ=255)
                                                               gTAGAAATTGTACTGCGGCACTCCCGCGTGTCgatgat                                     <  1:2220567/38‑1 (MQ=255)
                                                               gTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAAt                      >  1:1053423/1‑53 (MQ=255)
                                                                tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATaa                       <  1:2706929/51‑1 (MQ=255)
                                                                tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATaa                       <  1:1468108/51‑1 (MQ=255)
                                                                tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAatca   >  1:2376121/1‑71 (MQ=255)
                                                                tAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAatca   >  1:2853540/1‑71 (MQ=255)
                                                                 aGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAatcat  >  1:1174606/1‑71 (MQ=255)
                                                                  |                                                                     
GTGCCAATGCATTTCTGTAGCTATTTACACATGCATCCCATGACTTGCCACCTAATAAACACAGTAGAAATTGTACTGCGGCACTCCCGCGTGTCGATGATGAAGCACTGGATAATGTCGCAACACCATTAATCAT  >  minE/1356100‑1356235

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: