Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 137,451 G→C 24.2% V34V (GTC→GTG yadM ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE137,4510GC24.2% 163.0 / 50.2 99V34V (GTC→GTGyadMpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (20/55);  new base C (11/13);  total (31/68)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATC  >  minE/137386‑137517
                                                                 |                                                                  
tAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGAc                                                                  <  1:2739117/68‑1 (MQ=255)
tAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGAc                                                                  <  1:1401065/68‑1 (MQ=255)
           gCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGAc                                                                  <  1:2468453/57‑1 (MQ=255)
              atcatTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAActc                                                   <  1:2229505/69‑1 (MQ=255)
               tcatTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAActc                                                   <  1:1599349/68‑1 (MQ=255)
                cgttCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 <  1:107384/67‑1 (MQ=255)
                catTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 <  1:408743/69‑1 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:128416/1‑68 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:1066546/1‑69 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:64330/1‑69 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:429653/1‑69 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:2698261/1‑69 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:2647122/1‑69 (MQ=255)
                 atTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTAc                                                >  1:1731253/1‑69 (MQ=255)
                    cGTGACGTTATTGCCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:2254913/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:1998354/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:3045384/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:3005284/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:436775/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:3408574/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:2745514/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:84810/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:2376219/1‑65 (MQ=255)
                    cGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                 >  1:2591186/1‑65 (MQ=255)
                             attaCCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGt                                    >  1:2699380/1‑69 (MQ=255)
                             attaCCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGt                                    >  1:1710300/1‑69 (MQ=255)
                              ttaCCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAActc                                                   <  1:2073568/53‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:2701430/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:518944/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:420866/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:2531649/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:3212496/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:114282/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:3408487/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:3023352/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:3354866/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:3201417/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:1345742/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:1278801/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:1756914/69‑1 (MQ=255)
                                   aGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCaaa                              <  1:202356/69‑1 (MQ=255)
                                       aGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCg                                     <  1:1131324/58‑1 (MQ=255)
                                       aGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCg                                     <  1:2498441/58‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:3093264/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:908750/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:35634/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:1995748/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:2065052/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                <  1:1595648/49‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:2505415/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:1683116/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:2962078/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:349147/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:2886700/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:2509429/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:695231/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:129521/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:2118644/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:1972045/69‑1 (MQ=255)
                                                     ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa            <  1:548924/69‑1 (MQ=255)
                                                      gTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataat           <  1:2397726/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1453248/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:2622919/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:937315/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1394993/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1363493/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1615784/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1619772/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1202167/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:412797/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:736502/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1080708/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:1205987/69‑1 (MQ=255)
                                                       ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata          <  1:203971/69‑1 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:36002/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:57643/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:948755/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:1509424/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:154005/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:2722235/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:17339/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:2355059/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:2081026/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:3147849/1‑68 (MQ=255)
                                                        tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag         >  1:1965081/1‑68 (MQ=255)
                                                         gTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        >  1:1193080/1‑69 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:2031412/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:1458436/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:1294528/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:1790956/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:1203936/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:2603347/66‑1 (MQ=255)
                                                            gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt        <  1:2965373/66‑1 (MQ=255)
                                                               ttGACAGTGGCTGTAAATTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:2727718/69‑1 (MQ=255)
                                                               ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:2016070/69‑1 (MQ=255)
                                                               ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:620649/69‑1 (MQ=255)
                                                               ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:1050993/69‑1 (MQ=255)
                                                               ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:2582793/69‑1 (MQ=255)
                                                                tGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc  <  1:2902199/68‑1 (MQ=255)
                                                                 |                                                                  
TAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATC  >  minE/137386‑137517

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: