Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,591,000 T→A 36.0% V17V (GTT→GTA 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,0000TA36.0% 93.2 / 93.0 111V17V (GTT→GTAhdeAstress response protein acid‑resistance protein
Reads supporting (aligned to +/- strand):  ref base T (28/43);  new base A (21/19);  total (49/62)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 111 total aligned reads displayed.

GAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAC  >  minE/2590938‑2591054
                                                              |                                                      
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:244894/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3192685/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3226304/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3144030/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3076710/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3244466/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:472894/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:473425/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:3060001/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:2407770/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:2988217/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:2578241/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:136742/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:1275258/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:1102020/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:878390/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:1188474/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:1196130/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:864053/69‑1 (MQ=255)
gAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:812233/69‑1 (MQ=255)
   gATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCaa                                                <  1:1532185/68‑1 (MQ=255)
   gATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCaa                                                <  1:3054221/68‑1 (MQ=255)
      atTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCAAGTTGTGAGCAATGCa                                            >  1:1270660/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:1543785/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:1474259/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:506189/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:2174039/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:29474/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:781144/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:1079973/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:951160/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:2098156/1‑69 (MQ=255)
       ttACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAg                                           >  1:2667646/1‑69 (MQ=255)
         aCGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCg                                         <  1:966778/69‑1 (MQ=255)
         aCGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCg                                         <  1:3243709/69‑1 (MQ=255)
         aCGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCg                                         <  1:1603992/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTGTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:673046/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:2060144/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:1987403/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:575666/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:2903176/69‑1 (MQ=255)
                   aaGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCaaa                               <  1:1026971/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCTCaaaa                              <  1:3017012/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:629653/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:672743/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:736031/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:766539/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:2844800/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:2841692/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:2823466/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1652952/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1693330/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1447452/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1376486/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1284357/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:1254205/69‑1 (MQ=255)
                    aGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTAGTGAGCAATGCAGCGGATGCACaaaa                              <  1:2461095/69‑1 (MQ=255)
                        ttAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAagca                          >  1:284852/1‑69 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAg                                                   >  1:2955137/1‑35 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:1137903/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:1902596/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:712010/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:940293/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:3098746/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:2165765/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  >  1:303476/1‑36 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                  >  1:2216069/1‑68 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                  >  1:1211019/1‑68 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                  >  1:1411392/1‑68 (MQ=255)
                                 aTTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATaa                  >  1:289044/1‑68 (MQ=255)
                                  ttCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGc                                                  <  1:316036/35‑1 (MQ=255)
                                  ttCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCg                                         <  1:3204345/44‑1 (MQ=255)
                                             cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc     <  1:495866/69‑1 (MQ=255)
                                             cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc     <  1:952910/69‑1 (MQ=255)
                                             cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc     <  1:3228198/69‑1 (MQ=255)
                                             cTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc     <  1:247092/69‑1 (MQ=255)
                                               gCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTc     <  1:491841/67‑1 (MQ=255)
                                                cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc  >  1:668183/1‑69 (MQ=255)
                                                cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc  >  1:1042619/1‑69 (MQ=255)
                                                cttcttCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAc  >  1:1556102/1‑69 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:1617078/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:458694/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:3271568/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2896481/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2057454/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2608252/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:924728/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2060678/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2136088/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:603186/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:843540/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2383094/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:2431618/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGa                     >  1:920692/1‑47 (MQ=255)
                                                   cttctGCCAGTAGTGAGCAATGCAGCGGAAGCACAAAAAGCAGCTGa                     >  1:2140641/1‑47 (MQ=255)
                                                    ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa    >  1:723856/1‑63 (MQ=255)
                                                    ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa    >  1:1072789/1‑63 (MQ=255)
                                                    ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa    >  1:2038352/1‑63 (MQ=255)
                                                    ttctGCCAGTAGTGAGCAATGCAGCGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa    >  1:2730912/1‑63 (MQ=255)
                                                    ttctGCCAGTAGTGAGCAATGCAGAGGATGCACAAAAAGCAGCTGATAACAAAAAACCGGTCa    >  1:1389426/1‑63 (MQ=255)
                                                              |                                                      
GAGGATATTACGATGAAAAAAGTATTAGGCGTTATTCTTGGTGGTCTGCTTCTTCTGCCAGTTGTGAGCAATGCAGCGGATGCGCAAAAAGCAGCTGATAACAAAAAACCGGTCAAC  >  minE/2590938‑2591054

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: