Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,067,040 T→C 29.0% K55K (AAA→AAG yqjH ← predicted siderophore interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,067,0400TC29.0% 84.2 / 38.4 86K55K (AAA→AAGyqjHpredicted siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base C (19/6);  total (54/32)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAA  >  minE/2066980‑2067101
                                                            |                                                             
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:3061887/69‑1 (MQ=255)
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:49423/69‑1 (MQ=255)
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:173203/69‑1 (MQ=255)
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:2987085/69‑1 (MQ=255)
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:741750/69‑1 (MQ=255)
gATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:2355654/69‑1 (MQ=255)
gATGCCCTCTTCAGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:586333/69‑1 (MQ=255)
 aTGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:169165/68‑1 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:235317/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:825597/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:790484/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:768125/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:1904596/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:3324863/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:576271/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:3025988/1‑68 (MQ=255)
      ctctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGt                                                  >  1:3109872/1‑68 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1902313/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1810514/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:373616/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:277129/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1753308/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1744987/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:531835/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1587899/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:597710/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:91041/67‑1 (MQ=255)
             gTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGc                                            <  1:1011321/67‑1 (MQ=255)
                   gTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGa                                     >  1:1357891/1‑68 (MQ=255)
                   gTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGa                                     >  1:282568/1‑68 (MQ=255)
                      ggcgtcACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCtgtg                                                       <  1:2854951/47‑1 (MQ=255)
                                agtgagtgTCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:2590153/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCtgt                                                        >  1:1069864/1‑36 (MQ=16)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTTATCGTCAAAGCCGCGCGATGTaa                                >  1:1806214/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:3259281/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:551254/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:85450/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:886790/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:889580/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:936224/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:1109421/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:2006998/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:220288/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:1802781/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:1703368/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:2350046/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:240370/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:1232864/1‑60 (MQ=255)
                                agtgagtgCCTGGTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTaa                                >  1:1203107/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:1348374/69‑1 (MQ=255)
                                           ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg            <  1:770758/69‑1 (MQ=255)
                                           ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg            <  1:867961/69‑1 (MQ=255)
                                           ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg            <  1:1453730/69‑1 (MQ=255)
                                           ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg            <  1:386997/69‑1 (MQ=255)
                                           ggTTGAGGAAAGAAGAGCTTGCTGTGATCGTCAAAGCCGCGCGATGTAAAACCGTCCAGCGCCTcgccg            <  1:2216858/69‑1 (MQ=255)
                                                  gAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAgg      >  1:1876973/1‑68 (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:63133/69‑1 (MQ=255)
                                                   aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc    <  1:1561267/69‑1 (MQ=255)
                                                   aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc    <  1:519573/69‑1 (MQ=255)
                                                   aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc    <  1:1765239/69‑1 (MQ=255)
                                                   aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc    <  1:411595/69‑1 (MQ=255)
                                                   aaaGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGAc    <  1:20667/69‑1 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTcgccgc           >  1:433389/1‑61 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2720657/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2740915/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2575059/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:238486/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:1156844/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2464048/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:1560990/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:1578203/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:1597664/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2884252/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2971613/1‑69 (MQ=255)
                                                    aagaagAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:3382185/1‑69 (MQ=255)
                                                    aagaagAGTATGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACa   >  1:2453807/1‑69 (MQ=255)
                                                     agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa  >  1:15992/1‑69 (MQ=255)
                                                     agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa  >  1:2470781/1‑69 (MQ=255)
                                                     agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa  >  1:3377203/1‑69 (MQ=255)
                                                     agaagaGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACaa  >  1:321234/1‑69 (MQ=255)
                                                            |                                                             
GATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAA  >  minE/2066980‑2067101

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: