Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,628 A→G 5.5% intergenic (+98/+55) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,6280AG 5.5% 209.1 / 2.6 73intergenic (+98/+55)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/41);  new base G (2/2);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTA  >  minE/5560‑5697
                                                                    |                                                                     
tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATg                                                                     >  1:2321021/1‑71 (MQ=255)
 cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGAATCTGCATGc                                                                    >  1:772863/1‑71 (MQ=255)
 cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGc                                                                    >  1:1961496/1‑71 (MQ=255)
  aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt                                                                   >  1:3159980/1‑71 (MQ=255)
     gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt                                                                >  1:1413923/1‑71 (MQ=255)
      ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt                                                                <  1:1291823/70‑1 (MQ=255)
          aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                           <  1:2561009/71‑1 (MQ=255)
          aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg                                                            <  1:376284/70‑1 (MQ=255)
           tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                           >  1:3251649/1‑70 (MQ=255)
              gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg                                                           <  1:2569631/67‑1 (MQ=255)
                cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa                                                      >  1:3668316/1‑70 (MQ=255)
                cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                                     <  1:2723348/71‑1 (MQ=255)
                 gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt                                                     >  1:3482613/1‑70 (MQ=255)
                 gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa                                                    <  1:1389208/71‑1 (MQ=255)
                    ggCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa                                                      <  1:841959/66‑1 (MQ=255)
                     gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc                                                >  1:3199098/1‑71 (MQ=255)
                       gcgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGt                                              >  1:3160887/1‑71 (MQ=255)
                        cgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                             <  1:1265795/71‑1 (MQ=255)
                          tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                             <  1:2891438/69‑1 (MQ=255)
                           cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                           <  1:1661665/70‑1 (MQ=255)
                           cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc                                          >  1:2892450/1‑71 (MQ=255)
                             tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                        >  1:2386255/1‑71 (MQ=255)
                             tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                        >  1:2381545/1‑71 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                                 <  1:3764262/61‑1 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                                 <  1:3352438/61‑1 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg                                         <  1:731706/69‑1 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg                                         <  1:2396046/69‑1 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg                                         <  1:1736853/69‑1 (MQ=255)
                              aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc                                          <  1:3514893/68‑1 (MQ=255)
                                cAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg                                               >  1:479535/1‑61 (MQ=255)
                                  ggCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCAt                                   <  1:1797126/71‑1 (MQ=255)
                                   gCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                  >  1:569159/1‑71 (MQ=255)
                                     cTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTa                                                             >  1:2274993/1‑42 (MQ=255)
                                      tACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                               >  1:2444461‑M1/1‑71 (MQ=255)
                                        cgctaatTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCA                             <  1:556490‑M1/68‑1 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCtttt                                                                <  1:2291463/35‑1 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCtttt                                                                <  1:1661963/35‑1 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                   <  1:2316662/48‑1 (MQ=255)
                                         gTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                   <  1:545881/48‑1 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                           >  1:2268892/1‑71 (MQ=255)
                                          ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                           <  1:2171837/71‑1 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                           <  1:1609033/69‑1 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg                          <  1:3008581/70‑1 (MQ=255)
                                            aaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTg                          <  1:325016/70‑1 (MQ=255)
                                              ttCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACt                       <  1:1625082/71‑1 (MQ=255)
                                               tCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTg                      >  1:443700/1‑71 (MQ=255)
                                                 tGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc                     >  1:816602/1‑70 (MQ=255)
                                                 tGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc                     >  1:1338375/1‑70 (MQ=255)
                                                  gCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                           <  1:2632500/63‑1 (MQ=255)
                                                    aaTATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc                                                <  1:3248275/40‑1 (MQ=255)
                                                    aaTATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAAc                 <  1:2682855/71‑1 (MQ=255)
                                                     atatatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                               <  1:1018048‑M1/56‑1 (MQ=255)
                                                      tatatTGATTCTGCATTCTTTTGTAGGCAGTATAAGGCGtt                                             <  1:3766705/41‑1 (MQ=255)
                                                      tatatTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcac               >  1:2025025‑M1/1‑58 (MQ=255)
                                                       atatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAT                            >  1:1835342‑M1/1‑57 (MQ=255)
                                                        tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa                    <  1:274357/64‑1 (MQ=255)
                                                        tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTa              <  1:745164/70‑1 (MQ=255)
                                                        tatTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGG                               <  1:2988205‑M1/53‑1 (MQ=255)
                                                          ttGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactca            <  1:3227803‑M1/70‑17 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg                                               <  1:3882410/35‑1 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg                                               <  1:1069140/35‑1 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                             >  1:320899/1‑37 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgccgc                                     >  1:2448972/1‑45 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACg           >  1:2275236/1‑71 (MQ=255)
                                                          ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACg           <  1:976455/71‑1 (MQ=255)
                                                            gAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcact              >  1:1493863‑M1/1‑52 (MQ=255)
                                                            gAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGAc                        <  1:727554/56‑1 (MQ=255)
                                                            gAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgc          <  1:710464/70‑1 (MQ=255)
                                                             aaTTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaaca         <  1:1098384‑M1/70‑20 (MQ=255)
                                                              atttgcgTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATgaataaagcgcactcaacaat       <  1:258357‑M1/68‑22 (MQ=255)
                                                               tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACt                <  1:583103/61‑1 (MQ=255)
                                                               tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctgct      >  1:727503/1‑71 (MQ=255)
                                                                cTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTc     >  1:2309488/1‑71 (MQ=255)
                                                                 tGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCa                    >  1:1755539/1‑55 (MQ=255)
                                                                    aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGAc                        >  1:2661154/1‑48 (MQ=255)
                                                                    aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCg    >  1:2256088/1‑68 (MQ=255)
                                                                    aTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTa  >  1:372398/1‑70 (MQ=255)
                                                                    |                                                                     
TCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTA  >  minE/5560‑5697

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: