Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,591,021 G→A 47.0% A24A (GCG→GCA hdeA → stress response protein acid‑resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,591,0210GA47.0% 9.8 / 69.2 84A24A (GCG→GCAhdeAstress response protein acid‑resistance protein
Reads supporting (aligned to +/- strand):  ref base G (19/25);  new base A (21/18);  total (40/44)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

AAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTGGCTGTGGACG  >  minE/2590957‑2591088
                                                                |                                                                   
aaGTATTAGGCGTTATTCTTGGTGTTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:353442/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTGTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:673046/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:1026971/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:2060144/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:1987403/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:2903176/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:313558/69‑1 (MQ=255)
aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                                                                 <  1:575666/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCTCaaaa                                                                <  1:3017012/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1693330/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:2461095/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:2654817/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:2823466/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:2841692/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:2844800/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:821940/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:766539/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1376486/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1652952/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:736031/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1254205/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1284357/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:629653/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:1447452/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:639297/69‑1 (MQ=255)
 aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:672743/69‑1 (MQ=255)
  gTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                                                                <  1:984325/68‑1 (MQ=255)
     ttAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAagca                                                            >  1:284852/1‑69 (MQ=255)
              aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                                                    >  1:1211019/1‑68 (MQ=255)
              aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                                                    >  1:2216069/1‑68 (MQ=255)
              aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                                                    >  1:1411392/1‑68 (MQ=255)
              aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                                                    >  1:289044/1‑68 (MQ=255)
                          cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc                                       <  1:952910/69‑1 (MQ=255)
                          cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc                                       <  1:495866/69‑1 (MQ=255)
                          cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc                                       <  1:247092/69‑1 (MQ=255)
                          cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc                                       <  1:3228198/69‑1 (MQ=255)
                            gCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc                                       <  1:491841/67‑1 (MQ=255)
                             cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc                                    >  1:1042619/1‑69 (MQ=255)
                             cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc                                    >  1:668183/1‑69 (MQ=255)
                             cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc                                    >  1:1556102/1‑69 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2896481/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2608252/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:3271568/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:1617078/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2383094/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2431618/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:458694/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:843540/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2057454/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:920692/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:924728/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2031931/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:603186/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2136088/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                                                       >  1:2060678/1‑47 (MQ=255)
                                cttctGCCAGTAGTGAGCAATGCAGCGGAAGCACAAAAAGCAGCTGa                                                       >  1:2140641/1‑47 (MQ=255)
                                 ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa                                      >  1:2730912/1‑63 (MQ=255)
                                 ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa                                      >  1:723856/1‑63 (MQ=255)
                                 ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa                                      >  1:1072789/1‑63 (MQ=255)
                                 ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa                                      >  1:2038352/1‑63 (MQ=255)
                                 ttctGCCAGTAGTGAGCAATGCAGAGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa                                      >  1:1389426/1‑63 (MQ=255)
                                            gtgAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGaa                     <  1:10809/69‑1 (MQ=255)
                                            gtgAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGaa                     <  1:3039422/69‑1 (MQ=255)
                                            gtgAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGaa                     <  1:2908252/69‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:69790/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:2735714/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:454453/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:3430298/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:3272051/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:1882285/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:3080646/68‑1 (MQ=255)
                                              gAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:1846202/68‑1 (MQ=255)
                                               aGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAg                    <  1:433781/67‑1 (MQ=255)
                                                 cAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGAttt                >  1:2092634/1‑69 (MQ=255)
                                                  aaTGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTc               >  1:687258/1‑69 (MQ=255)
                                                  aaTGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTc               >  1:2658297/1‑69 (MQ=255)
                                                  aaTGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTc               >  1:1564223/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAGCCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:114008/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAGCCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:2739694/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAGCCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:1995971/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAGCCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:1649713/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAGCCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:2188594/1‑69 (MQ=255)
                                                     gCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTg            >  1:781299/1‑69 (MQ=255)
                                                           gATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTGGCtgtg      >  1:3171141/1‑69 (MQ=255)
                                                               cgcAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTGGCTGTGGACg  >  1:762039/1‑69 (MQ=255)
                                                               cgcAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTGGCTGTGGACg  >  1:2581780/1‑69 (MQ=255)
                                                                |                                                                   
AAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAACTCCTGGACCTGTGAAGATTTCCTGGCTGTGGACG  >  minE/2590957‑2591088

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: