Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,639,287 A→G 8.7% intergenic (+57/+53) ypfN → / ← ypfH hypothetical protein/predicted hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,639,2870AG 8.7% 256.1 / 11.4 94intergenic (+57/+53)ypfN/ypfHhypothetical protein/predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base G (2/6);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAA  >  minE/1639219‑1639353
                                                                    |                                                                   
ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                   >  1:4502322/1‑71 (MQ=255)
ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                   <  1:3442558/71‑1 (MQ=255)
ccGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                   <  1:3704294/71‑1 (MQ=255)
 cGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                   <  1:773098/70‑1 (MQ=255)
  gaagaaGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGa                                                                 >  1:511664/1‑71 (MQ=255)
     gaagTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGc                                                             <  1:2470358/71‑1 (MQ=255)
        gTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                           >  1:841144/1‑70 (MQ=255)
          aaGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGAt                                                        >  1:4334209/1‑71 (MQ=255)
            gCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGc                                                      <  1:2063037/71‑1 (MQ=255)
              tGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGa                                                    >  1:217416/1‑71 (MQ=255)
                aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGa                                                    <  1:1472496/69‑1 (MQ=255)
                aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACg                                                  <  1:170100/71‑1 (MQ=255)
                aGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACg                                                  <  1:176796/71‑1 (MQ=255)
                 gTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGc                                                 <  1:754492/71‑1 (MQ=255)
                  tGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCt                                                <  1:3028274/71‑1 (MQ=255)
                       attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            >  1:1694301‑M1/4‑46 (MQ=255)
                       attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            >  1:2786583‑M1/4‑46 (MQ=255)
                       attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            <  1:2626015‑M1/68‑26 (MQ=255)
                       attcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            <  1:4143815‑M1/68‑26 (MQ=255)
                       actcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            <  1:4576253‑M1/68‑26 (MQ=255)
                       actcgtaGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgcta                                                   >  1:800526‑M1/4‑46 (MQ=255)
                       actcgtaGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                            <  1:3539727‑M1/68‑26 (MQ=255)
                       actcgtaGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                            <  1:3142512‑M1/68‑26 (MQ=255)
                        ctcgtAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttat                                            >  1:766880‑M1/3‑45 (MQ=255)
                        ctcgtAGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatc                                           <  1:3266569‑M1/69‑27 (MQ=255)
                         tcgTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatg                                                               >  1:3737166‑M1/2‑44 (MQ=255)
                         tcgTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaa                                                                  >  1:417648‑M1/2‑44 (MQ=255)
                          cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:3560763‑M1/1‑43 (MQ=255)
                          cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                   <  1:2398275‑M1/45‑3 (MQ=255)
                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcagg                                        >  1:3620561‑M1/1‑42 (MQ=255)
                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcag                                         <  1:4664754‑M1/70‑29 (MQ=255)
                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcag                                         <  1:2527874‑M1/70‑29 (MQ=255)
                           gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                          <  1:2292632‑M1/69‑28 (MQ=255)
                           gTAGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatca                                          <  1:2061615‑M1/69‑28 (MQ=255)
                             aGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                     >  1:3526884/1‑71 (MQ=255)
                             aGGCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                            <  1:2228200‑M1/65‑26 (MQ=255)
                              ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                    <  1:3121612/71‑1 (MQ=255)
                              ggCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                    <  1:1179540/71‑1 (MQ=255)
                               gCCGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatcaggccta                                    <  1:2665404‑M1/71‑34 (MQ=255)
                                 cGGATAAGGCGTCCACGCCGCATCCGGCAGTGGTGGcaaatgcctgatgcgctacgcttatcaggccta                                    >  1:3745869‑M1/1‑36 (MQ=255)
                                  ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggccta                                    >  1:1051323‑M1/1‑35 (MQ=255)
                                  ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                          >  1:2540654‑M1/1‑35 (MQ=255)
                                  ggATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatca                                          >  1:1768925‑M1/1‑35 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                                <  1:415688‑M1/71‑38 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                                <  1:3744905‑M1/71‑38 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                   <  1:2947988‑M1/36‑3 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAATGGTGAca                                                                   <  1:860661‑M1/36‑3 (MQ=255)
                                   gATAAGGCGTTAACGCCGCATCCGGCAATGGTGAAc                                                                   <  1:2174536/36‑1 (MQ=255)
                                    aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                            <  1:1256639‑M1/42‑10 (MQ=255)
                                    aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                            <  1:3769719‑M1/42‑10 (MQ=255)
                                    aTAAGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcct                                                            <  1:1726267‑M1/42‑10 (MQ=255)
                                    aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                              <  1:4354912/71‑1 (MQ=255)
                                    aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                              <  1:4210460/71‑1 (MQ=255)
                                    aTAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGcc                              <  1:1902689/71‑1 (MQ=255)
                                     tAAGGCGTTTACGCCGCATCCGGCAACAGTGCCGACTGCCTGATGCGACGCTCGCGCGTCTTATCAGGccg                             >  1:718344/1‑70 (MQ=1)
                                      aaGGCGTTCACGCCGCATCCGGCAATGGTGAcaaatgcctga                                                          >  1:60934‑M1/1‑31 (MQ=255)
                                         gCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTa                         <  1:4425751/71‑1 (MQ=255)
                                          cGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTaa                        >  1:5151097/1‑71 (MQ=255)
                                           gTTCACGCCGCATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctacaga                                <  1:651551‑M1/63‑38 (MQ=255)
                                           gTTAACGCCGCATCCGGCAATGGTGAACGATGCCt                                                           <  1:4731550/35‑1 (MQ=255)
                                            ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTAt                                    <  1:2366829/57‑1 (MQ=255)
                                            ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                               <  1:3638156/62‑1 (MQ=255)
                                            ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                      >  1:1261933/1‑71 (MQ=255)
                                            ttAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATc                      >  1:3411478/1‑71 (MQ=255)
                                             tAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGc                               <  1:841510/61‑1 (MQ=255)
                                              aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCa                                  >  1:2613281/1‑57 (MQ=255)
                                              aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                    >  1:4503588/1‑71 (MQ=255)
                                              aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                    >  1:931845/1‑71 (MQ=255)
                                              aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                    >  1:1180973/1‑71 (MQ=255)
                                              aaCGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTg                    >  1:3506308/1‑71 (MQ=255)
                                                  ccgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTa                         >  1:3891078/1‑62 (MQ=255)
                                                   cgcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCttt               >  1:5414319/1‑71 (MQ=255)
                                                    gcATCCGGCAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctaca                                  >  1:2911165‑M1/1‑17 (MQ=255)
                                                    gcATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                     <  1:5360181/48‑1 (MQ=255)
                                                     cATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGt                          >  1:2983973/1‑58 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:4036882‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:2782819‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:3345929‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:4226855‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:3830605‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAcaaatgcctgatgcgctacgctta                                             >  1:3260030‑M1/1‑13 (MQ=255)
                                                        ccGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTa             >  1:5385728/1‑68 (MQ=255)
                                                         cGGCAGTGGTGGcaaatgcctgatgcgctacgcttat                                            >  1:1302286‑M1/1‑12 (MQ=255)
                                                         cGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcat          <  1:4598506/70‑1 (MQ=255)
                                                         cGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcat          >  1:2368769/1‑70 (MQ=255)
                                                          ggCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTAcatca        <  1:3014156/71‑1 (MQ=255)
                                                            cAATGGTGAcaaatgcctgatgcgctacgcttatcaggcctaca                                  <  1:62743‑M1/44‑36 (MQ=255)
                                                            cAATGGTGAcaaatgcctgatgcgctacgcttatca                                          <  1:3056331‑M1/36‑28 (MQ=255)
                                                            cAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTa                                     >  1:2206842/1‑40 (MQ=255)
                                                            cAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTGTACATCAtt      >  1:3445305/1‑71 (MQ=255)
                                                              aTGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTc    <  1:1408944/71‑1 (MQ=255)
                                                               tGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCa   <  1:2961427/71‑1 (MQ=255)
                                                                 gagAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCaa  >  1:661652/3‑70 (MQ=255)
                                                                 gTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCa                                  >  1:4591110/1‑38 (MQ=255)
                                                                    |                                                                   
CCGAAGAAGTAAGCTGAGTGCTCAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATGGTGAACGATGCCTGATGCGACGCTCACGCGTCTTATCAGGCCAGGTAATCTGTCTTTTACATCATTTCAA  >  minE/1639219‑1639353

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: