Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 357,299 T→C 16.4% F5L (TTT→CTT)  sfmA → predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE357,2990TC16.4% 179.3 / 28.6 73F5L (TTT→CTT) sfmApredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (5/7);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

AAAGTAATCGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATG  >  minE/357231‑357365
                                                                    |                                                                  
aaaGTAATCGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGatt                                                                   >  1:33361/1‑70 (MQ=255)
  aGTAATCGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGatttat                                                                <  1:2892439/71‑1 (MQ=255)
        cGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTc                                                          >  1:2706860/1‑71 (MQ=255)
        cGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTc                                                          >  1:1568075/1‑71 (MQ=255)
          gTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTg                                                        >  1:1179309/1‑71 (MQ=255)
           ttATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTgc                                                       <  1:2695004/71‑1 (MQ=255)
           ttATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTgc                                                       <  1:997192/71‑1 (MQ=255)
           ttATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTgc                                                       <  1:1542555/71‑1 (MQ=255)
              tataCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGat                                                                    <  1:1957937/55‑1 (MQ=255)
              tataCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgct                                                    <  1:1033397/71‑1 (MQ=255)
              tataCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgct                                                    <  1:1284589/71‑1 (MQ=255)
                 aCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGatt                                                                   <  1:502996/53‑1 (MQ=255)
                 aCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGatt                                                                   <  1:1361049/53‑1 (MQ=255)
                 aCGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGatt                                                                   <  1:2209614/53‑1 (MQ=255)
                  cGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgctggct                                                >  1:1467821/1‑71 (MQ=255)
                  cGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCgctggct                                                >  1:1071774/1‑71 (MQ=255)
                        tgGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGc                                           <  1:1978191/70‑1 (MQ=255)
                        tgGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGc                                           <  1:1971107/70‑1 (MQ=255)
                          gAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCAc                                         <  1:2747028/70‑1 (MQ=255)
                           aGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTa                                       <  1:122978/71‑1 (MQ=255)
                           aGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTa                                       <  1:894979/71‑1 (MQ=255)
                           aGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTa                                       <  1:1494916/71‑1 (MQ=255)
                           aGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTa                                       <  1:1181402/71‑1 (MQ=255)
                             tCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgctg                                                   <  1:79753/57‑1 (MQ=255)
                             tCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTa                                       >  1:2141905/1‑69 (MQ=255)
                             tCGATAAATGAGATTGAAGGAATACATATGAAATTAAGATTTATTTCGTCTGCgctg                                                   <  1:869746/57‑1 (MQ=255)
                              cGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCt                                                         >  1:2081392/1‑50 (MQ=255)
                               gATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgctggc                                                 >  1:3210833/1‑57 (MQ=255)
                               gATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTc                                    <  1:2847632/70‑1 (MQ=255)
                                  aaaTGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccg                                >  1:3068961/1‑71 (MQ=255)
                                   aaTGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccgc                               <  1:2959634/71‑1 (MQ=255)
                                   aaTGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccgc                               <  1:129386/71‑1 (MQ=255)
                                   aaTGAGATTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccgc                               <  1:900621/71‑1 (MQ=255)
                                    aTGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccgc                               >  1:3097874/1‑70 (MQ=255)
                                    aTGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccgc                               >  1:78122/1‑70 (MQ=255)
                                    aTGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCAGCGCTGGCTGCCGCACTATTcgccg                                >  1:3048186/1‑69 (MQ=255)
                                       agaTTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgctggc                                                 >  1:805751/1‑49 (MQ=255)
                                        gaTTGAAGGAATATATATGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACgg                          >  1:416488/1‑71 (MQ=255)
                                          ttGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACggg                         >  1:972500/1‑70 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCgctggc                                                 >  1:2101308/1‑44 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccg                                >  1:532332/1‑61 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTa                             >  1:2285994/1‑64 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTa                             >  1:2333659/1‑64 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTa                       <  1:887176/70‑1 (MQ=255)
                                            gAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAg                      >  1:992410/1‑71 (MQ=255)
                                             aaGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTa                       <  1:348331/69‑1 (MQ=255)
                                             aaGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGt                     <  1:2558864/71‑1 (MQ=255)
                                              aGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACggg                         <  1:439716/66‑1 (MQ=255)
                                              aGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACggg                         <  1:435649/66‑1 (MQ=255)
                                              aGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGtt                    <  1:483656/71‑1 (MQ=255)
                                              aGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGtt                    <  1:2059406/71‑1 (MQ=255)
                                               ggAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgccg                                >  1:3164207/1‑58 (MQ=255)
                                               ggAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACggg                         >  1:1182695/1‑65 (MQ=255)
                                               ggAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACggg                         >  1:572919/1‑65 (MQ=255)
                                                gAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTAt                  >  1:2579184/1‑71 (MQ=255)
                                                gAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTAt                  >  1:1623424/1‑71 (MQ=255)
                                                 aaTATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTAtt                                     >  1:1337837/1‑51 (MQ=255)
                                                 aaTATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTAtg                 <  1:496300/71‑1 (MQ=255)
                                                  atatatatGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTAtgc                >  1:1950927/1‑71 (MQ=255)
                                                   tatatatGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGt                     <  1:81526/65‑1 (MQ=255)
                                                   tatatatGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTcgcc                                 <  1:2850285/53‑1 (MQ=255)
                                                    atatatGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGc                                           <  1:3057525/42‑1 (MQ=255)
                                                    atatatGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGc                                           <  1:2682291/42‑1 (MQ=255)
                                                    atatatGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGc                                           <  1:2544209/42‑1 (MQ=255)
                                                     tatatGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTAtgctgc             >  1:1970658/1‑71 (MQ=255)
                                                        atGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGt          >  1:2392904/1‑71 (MQ=255)
                                                         tGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTa                       <  1:1038598/57‑1 (MQ=255)
                                                         tGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACCGGTAGTTATGCTGCCGtt         >  1:3039752/1‑71 (MQ=255)
                                                         tGAAATTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGt          <  1:1022605/70‑1 (MQ=255)
                                                          gAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACg                           >  1:472008/1‑52 (MQ=255)
                                                           aaaTTAAGATTTATTTCGTCTGCGCTGGCTGCCGCGCTATTCGCCGCTACGGGTAGTTATGCTGCCGTTg        <  1:2956785/70‑1 (MQ=255)
                                                           aaaTTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTg        <  1:1036977/70‑1 (MQ=255)
                                                           aaaTTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTg        >  1:2076720/1‑70 (MQ=255)
                                                            aaTTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTa      <  1:3320451/71‑1 (MQ=255)
                                                             aTTAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGATGCCGTTg        >  1:1191361/1‑68 (MQ=255)
                                                              ttAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTg        >  1:3353369/1‑67 (MQ=255)
                                                              ttAAGACTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTg        <  1:1592032/67‑1 (MQ=255)
                                                                aaGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAg     >  1:298483/1‑68 (MQ=255)
                                                                aaGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATg  >  1:1801756/1‑71 (MQ=255)
                                                                    |                                                                  
AAAGTAATCGGTTATATACGGATGTGGAGTCGATAAATGAGATTGAAGGAATATATATGAAATTAAGATTTATTTCGTCTGCGCTGGCTGCCGCACTATTCGCCGCTACGGGTAGTTATGCTGCCGTTGTAGATG  >  minE/357231‑357365

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: