Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,067,061 A→G 38.2% R48R (CGT→CGC yqjH ← predicted siderophore interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,067,0610AG38.2% 59.8 / 77.1 115R48R (CGT→CGCyqjHpredicted siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base G (38/6);  total (73/42)
Fisher's exact test for biased strand distribution p-value = 5.50e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 117 total aligned reads displayed.

GTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTC  >  minE/2066999‑2067129
                                                              |                                                                    
gTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGa                                                                 >  1:282568/1‑68 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTTATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1806214/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1802781/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1109421/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1703368/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1203107/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:1232864/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:889580/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:2006998/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:2350046/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:886790/1‑60 (MQ=255)
             agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:240370/1‑60 (MQ=255)
             aGTGAGAGCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                                            >  1:439305/1‑60 (MQ=255)
                     tCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTcg                                           >  1:236402/1‑69 (MQ=255)
                       aGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTcgcc                                         >  1:1968067/1‑69 (MQ=255)
                        ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg                                        <  1:867961/69‑1 (MQ=255)
                        ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg                                        <  1:1453730/69‑1 (MQ=255)
                        ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg                                        <  1:770758/69‑1 (MQ=255)
                        ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg                                        <  1:1348374/69‑1 (MQ=255)
                        ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg                                        <  1:2216858/69‑1 (MQ=255)
                               gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAgg                                  >  1:1399632/1‑68 (MQ=255)
                               gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAgg                                  >  1:1275533/1‑68 (MQ=255)
                               gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGACGCCGAgg                                  >  1:2883148/1‑68 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:1561267/69‑1 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:20667/69‑1 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:63133/69‑1 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:1765239/69‑1 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:519573/69‑1 (MQ=255)
                                aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc                                <  1:411595/69‑1 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTcgccgc                                       >  1:433389/1‑61 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:2720657/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:3382185/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:2575059/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:2884252/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:1597664/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:1578203/1‑69 (MQ=255)
                                 aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:2464048/1‑69 (MQ=255)
                                 aagaagAGTATGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa                               >  1:2453807/1‑69 (MQ=255)
                                  agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa                              >  1:321234/1‑69 (MQ=255)
                                  agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa                              >  1:2470781/1‑69 (MQ=255)
                                  agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa                              >  1:15992/1‑69 (MQ=255)
                                  agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa                              >  1:3377203/1‑69 (MQ=255)
                                                tgATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCg                >  1:3066796/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGcc                                              >  1:730865/1‑37 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGa                                 >  1:3093800/1‑50 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCgg               >  1:671073/1‑68 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:698383/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:3359354/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:1293349/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:303298/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:624367/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:755215/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:856016/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:2726691/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:1012230/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:2640051/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:2239593/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:895952/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:787924/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:2031646/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:1664281/1‑69 (MQ=255)
                                                  aTCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCAAGGACAATGCGCTGAAAACCGgc              >  1:180641/1‑69 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2854668/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3389561/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:490821/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:376017/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1030672/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1086404/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3121054/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3041233/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1144652/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:632123/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2929388/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1897533/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1548398/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2374722/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:844737/69‑1 (MQ=255)
                                                     gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2131518/69‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:78021/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2704261/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3131173/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3075496/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1247169/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:3008563/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1630763/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:1631401/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:2339553/67‑1 (MQ=255)
                                                       cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt           <  1:415335/67‑1 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGa         >  1:774860/1‑68 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:2472862/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:1015605/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:1368421/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:2833538/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:2974504/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:2978340/1‑69 (MQ=255)
                                                        aaaGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGAt        >  1:370735/1‑69 (MQ=255)
                                                          aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATg       <  1:266080/68‑1 (MQ=255)
                                                          aGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATg       <  1:1482697/68‑1 (MQ=255)
                                                              aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTc  <  1:1539182/69‑1 (MQ=255)
                                                              aCGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTc  <  1:3014326/69‑1 (MQ=255)
                                                              |                                                                    
GTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTC  >  minE/2066999‑2067129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: