Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 137,499 T→C 20.9% S18S (TCA→TCG yadM ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE137,4990TC20.9% 163.9 / 33.7 120S18S (TCA→TCGyadMpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (23/71);  new base C (12/13);  total (35/85)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-02

Only 100 of 121 total aligned reads displayed.

GGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGT  >  minE/137439‑137566
                                                            |                                                                   
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:2118644/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:1972045/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:2505415/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:2886700/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:2962078/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:695231/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:1683116/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:129521/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:349147/69‑1 (MQ=255)
ggTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataa                                                             <  1:548924/69‑1 (MQ=255)
 gTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataat                                                            <  1:2397726/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:937315/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1619772/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:412797/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1394993/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1453248/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1080708/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1363493/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:736502/69‑1 (MQ=255)
  ttgttgCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaata                                                           <  1:1205987/69‑1 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:2081026/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:2355059/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:1965081/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:2722235/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:17339/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:3147849/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:57643/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:36002/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:154005/1‑68 (MQ=255)
   tGTTGCTTTCACGGTAGCGGTAAAATCGACATTAACATCCGTTGCAAATACCGAAGGCGataacaatag                                                          >  1:1509424/1‑68 (MQ=255)
    gTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         >  1:1193080/1‑69 (MQ=255)
       gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         <  1:1203936/66‑1 (MQ=255)
       gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         <  1:1790956/66‑1 (MQ=255)
       gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         <  1:2965373/66‑1 (MQ=255)
       gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         <  1:2603347/66‑1 (MQ=255)
       gCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATAtt                                                         <  1:2031412/66‑1 (MQ=255)
          ttGACAGTGGCTGTAAATTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                                                   <  1:2727718/69‑1 (MQ=255)
          ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                                                   <  1:620649/69‑1 (MQ=255)
          ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                                                   <  1:2582793/69‑1 (MQ=255)
          ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                                                   <  1:1050993/69‑1 (MQ=255)
          ttGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATc                                                   <  1:2016070/69‑1 (MQ=255)
             aCAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAg                                                 <  1:3274234/68‑1 (MQ=255)
             aCAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAg                                                 <  1:1770254/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTTCATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:1249801/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:2252145/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:354407/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:1853288/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:2478274/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:762881/68‑1 (MQ=255)
              cAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGt                                                <  1:3258783/68‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:543667/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:1261579/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:2809841/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:2979760/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:1810045/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:3021218/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:777085/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:2443371/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:686867/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:2136840/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCAc                                           <  1:2312126/69‑1 (MQ=255)
                  ggCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGGGATAATAATATTCCCATCAGTATCAc                                           <  1:2746434/69‑1 (MQ=255)
                         aaCTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTa                                    <  1:2301480/69‑1 (MQ=255)
                             ctACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg                                >  1:1637191/1‑69 (MQ=255)
                             ctACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg                                >  1:2109422/1‑69 (MQ=255)
                             ctACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg                                >  1:2835535/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:2816048/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:613231/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:1012159/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:2322438/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:2525884/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:3096564/1‑69 (MQ=255)
                                  ttAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgtt                           >  1:1676095/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)
                                    aaTATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg                                >  1:1581070/1‑62 (MQ=255)
                                    aaTATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCgttgt                         >  1:1872682/1‑69 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:3007369/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:2538534/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:115802/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:1901442/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:1593612/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:677885/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:2587091/43‑1 (MQ=255)
                                              gCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                         <  1:3035693/43‑1 (MQ=255)
                                                   tACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGAc           >  1:322608/1‑68 (MQ=255)
                                                   tACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACg          >  1:2793502/1‑69 (MQ=255)
                                                   tACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACg          >  1:3058567/1‑69 (MQ=255)
                                                          ggTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCtt    <  1:2031807/68‑1 (MQ=255)
                                                          ggTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCt     <  1:2661074/67‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCt     <  1:1494852/66‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:1276964/69‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:649513/69‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:3420789/69‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:3056368/69‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:773305/69‑1 (MQ=255)
                                                           gTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:189200/69‑1 (MQ=255)
                                                           gTGATAATAATATTACCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGt  <  1:2293335/69‑1 (MQ=255)
                                                            tGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCTTTGTTTTTATCATTTGACGTCCCTTGt  <  1:2979951/68‑1 (MQ=255)
                                                            |                                                                   
GGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGT  >  minE/137439‑137566

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: