Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 137,477 T→C 20.7% I26V (ATT→GTT)  yadM ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE137,4770TC20.7% 161.4 / 32.1 115I26V (ATT→GTT) yadMpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (20/71);  new base C (11/13);  total (31/84)
Fisher's exact test for biased strand distribution p-value = 3.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 115 total aligned reads displayed.

ATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGT  >  minE/137415‑137543
                                                              |                                                                  
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:420866/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:3212496/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:3354866/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:3408487/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:1345742/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:2701430/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:2531649/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:1278801/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:518944/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:1756914/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:3201417/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:114282/69‑1 (MQ=255)
      aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                                                        <  1:202356/69‑1 (MQ=255)
          aGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCg                                                               <  1:2498441/58‑1 (MQ=255)
          aGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCg                                                               <  1:1131324/58‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:1595648/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:908750/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:2065052/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:1995748/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:3093264/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                          <  1:35634/49‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:2886700/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:129521/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:349147/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:1683116/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:2509429/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:2962078/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:548924/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:2118644/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:695231/69‑1 (MQ=255)
                        ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                      <  1:1972045/69‑1 (MQ=255)
                         gTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataat                                     <  1:2397726/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:2622919/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:203971/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:736502/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:937315/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1080708/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1394993/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1619772/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1205987/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1615784/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1202167/69‑1 (MQ=255)
                          ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                    <  1:1453248/69‑1 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:36002/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:17339/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:2355059/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:2722235/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:154005/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:2081026/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:3147849/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:1509424/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:57643/1‑68 (MQ=255)
                           tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                   >  1:1965081/1‑68 (MQ=255)
                            gTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  >  1:1193080/1‑69 (MQ=255)
                               gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  <  1:1458436/66‑1 (MQ=255)
                               gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  <  1:2965373/66‑1 (MQ=255)
                               gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  <  1:1790956/66‑1 (MQ=255)
                               gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  <  1:2603347/66‑1 (MQ=255)
                               gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                  <  1:1294528/66‑1 (MQ=255)
                                  ttGACAGTGGCTGTAAATTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:2727718/69‑1 (MQ=255)
                                  ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:1050993/69‑1 (MQ=255)
                                  ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:620649/69‑1 (MQ=255)
                                  ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:2582793/69‑1 (MQ=255)
                                  ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:2016070/69‑1 (MQ=255)
                                   tGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                            <  1:2902199/68‑1 (MQ=255)
                                     aCAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAg                          <  1:1770254/68‑1 (MQ=255)
                                     aCAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAg                          <  1:3274234/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTTCATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:1249801/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:1853288/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:3139786/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:301913/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:762881/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:2252145/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:354407/68‑1 (MQ=255)
                                      cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                         <  1:2478274/68‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:1261579/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:543667/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:3021218/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:2979760/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:2809841/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:2443371/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:777085/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:2312126/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:2136840/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                    <  1:1810045/69‑1 (MQ=255)
                                          ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGGGATAATAATATTCCCATCAGTATCAc                    <  1:2746434/69‑1 (MQ=255)
                                                     ctACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg         >  1:1637191/1‑69 (MQ=255)
                                                     ctACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg         >  1:2109422/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:613231/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:1012159/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:3096564/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:2816048/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:2775755/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:1676095/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:2525884/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:2322438/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGTGATAATAATAATCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:1441640/1‑69 (MQ=255)
                                                          ttAATATCCGTTGCAAATACTGAGGGCGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt    >  1:1703811/1‑69 (MQ=255)
                                                            aaTATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTAt                       >  1:1921932/1‑48 (MQ=255)
                                                            aaTATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgttgt  >  1:1821416/1‑69 (MQ=255)
                                                              |                                                                  
ATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGT  >  minE/137415‑137543

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: