Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 5,624 C→T 5.2% intergenic (+94/+59) yaaX → / ← yaaA hypothetical protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE5,6240CT 5.2% 245.3 / 2.2 77intergenic (+94/+59)yaaX/yaaAhypothetical protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/46);  new base T (2/2);  total (29/48)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTC  >  minE/5555‑5694
                                                                     |                                                                      
ggCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt                                                                       <  1:2009738/71‑1 (MQ=255)
 gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                      >  1:1148074/1‑71 (MQ=255)
     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)
        gcgcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGc                                                                     <  1:2370522/65‑1 (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)
                               tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                      <  1:1395336/41‑1 (MQ=255)
                               tCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGtt                                          <  1:2891438/69‑1 (MQ=255)
                                cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCa                                        <  1:1661665/70‑1 (MQ=255)
                                cTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAc                                       >  1:2892450/1‑71 (MQ=255)
                                 ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                      <  1:3158292/39‑1 (MQ=255)
                                 ttATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg                                                                      <  1:1174162/39‑1 (MQ=255)
                                  tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                     >  1:2386255/1‑71 (MQ=255)
                                  tATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcgc                                     >  1:2381545/1‑71 (MQ=255)
                                   aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                              <  1:3352438/61‑1 (MQ=255)
                                   aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAgg                                              <  1:3764262/61‑1 (MQ=255)
                                   aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg                                      <  1:2396046/69‑1 (MQ=255)
                                   aTCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCAcg                                      <  1:731706/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:545881/48‑1 (MQ=255)
                                              gTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATaa                                                <  1:2316662/48‑1 (MQ=255)
                                               ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                        <  1:2171837/71‑1 (MQ=255)
                                               ttAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAtt                        >  1:2268892/1‑71 (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:1069140/35‑1 (MQ=255)
                                                               ttGAATCTGCATGCTTTTGTAGGCAGGATAAGGCg                                            <  1:3882410/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)
                                                                    tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACt             <  1:583103/61‑1 (MQ=255)
                                                                    tCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctgct   >  1:727503/1‑71 (MQ=255)
                                                                     cTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTc  >  1:2309488/1‑71 (MQ=255)
                                                                     |                                                                      
GGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTC  >  minE/5555‑5694

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