複数の候補アプタマー配列から選択された選択アプタマー配列を含む被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測装置であって、入力手段と、自由エネルギーの算出手段と、目的アプタマー配列の濃縮の有無を予測する予測手段とを備え、前記入力手段は、被対照プールに含まれる被対照アプタマー配列群の配列情報、および対照アプタマー配列群の配列情報を入力する入力工程を実行し、前記算出手段は、前記被対照アプタマー配列群の配列情報および前記対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程を実行し、前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを前記目的アプタマー配列が濃縮されたプールと予測する予測工程を実行し、前記対照アプタマー配列群が、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする予測装置。前記対照アプタマー配列群が、前記1候補アプタマー配列群であり、前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記1候補アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記1候補アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する、請求項1記載の予測装置。前記対照アプタマー配列群が、前記2仮想アプタマー配列群であり、前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記2仮想アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記2仮想アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する、請求項1記載の予測装置。さらに、前記仮想アプタマー配列の生成手段を備え、前記生成手段は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成する生成工程を実行する、請求項3記載の予測装置。前記対照アプタマー配列群が、前記1候補アプタマー配列群、前記2仮想アプタマー配列群、および下記3仮想アプタマー配列群であり、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群3前記候補アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記候補アプタマー配列群由来の仮想アプタマー配列群前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記1候補アプタマー配列群の自由エネルギーと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記被対照アプタマー配列群の自由エネルギーの減少程度が、前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記1候補アプタマー配列群の自由エネルギーの減少程度よりも大きい場合、前記被対照プールを前記目的アプタマー配列が濃縮されたプールと予測する前記予測工程を実行する、請求項1から4のいずれか一項に記載の予測装置。さらに、前記仮想アプタマー配列の生成手段を備え、前記生成手段は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成し、前記候補アプタマー配列群の配列情報から、前記候補アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記3仮想アプタマー配列群を生成する生成工程を実行する、請求項5記載の予測装置。前記仮想アプタマー配列群の配列数が、100個~1000万個である、請求項1から6のいずれか一項に記載の予測装置。前記被対照アプタマー配列群および前記候補アプタマー配列群の配列数が、それぞれ、5個~1億個である、請求項1から7のいずれか一項に記載の予測装置。前記自由エネルギーが、自由エネルギーの平均および分散の少なくとも一方である、請求項1から8のいずれか一項に記載の予測装置。さらに、出力手段を備え、前記出力手段は、前記予測手段で予測された予測結果を出力する、請求項1から9のいずれか一項に記載の予測装置。複数候補アプタマー配列から選択された選択アプタマー配列を有する被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測方法であって、被検対照プールに含まれる被対照アプタマー配列群の配列情報および対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程と、前記被対照アプタマー配列群の自由エネルギーと、前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが、前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを、目的アプタマー配列が濃縮されたプールと予測する予測工程とを有し、前記対照アプタマー配列群が、1複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする予測方法。前記対照アプタマー配列群が、前記1の候補アプタマー配列群であり、前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記1候補アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記1候補アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する、請求項11記載の予測方法。前記対照アプタマー配列群が、前記2被対照アプタマー配列群由来の仮想アプタマー配列群であり、前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記2仮想アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記2仮想アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する、請求項11記載の予測方法。さらに、仮想アプタマー配列の生成工程を有し、前記生成工程は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成する工程である、請求項13記載の予測方法。前記対照アプタマー配列群が、前記1候補アプタマー配列群、前記2仮想アプタマー配列群、および下記3の仮想アプタマー配列群であり、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群3前記候補アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記候補アプタマー配列群由来の仮想アプタマー配列群前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記1候補アプタマー配列群の自由エネルギーと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記被対照アプタマー配列群の自由エネルギーの減少程度が、前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記1候補アプタマー配列群の自由エネルギーの減少程度よりも大きい場合、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する、請求項11から14のいずれか一項に記載の予測方法。さらに、生成工程を有し、前記生成工程は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成し、前記候補アプタマー配列群の配列情報から、前記候補アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記3仮想アプタマー配列群を生成する工程である、請求項15記載の予測方法。前記仮想アプタマー配列群の配列数が、100個~1000万個である、請求項11から16のいずれか一項に記載の予測方法。前記被対照アプタマー配列群および前記候補アプタマー配列群の配列数が、それぞれ、5個~1億個である、請求項11から17のいずれか一項に記載の予測方法。前記自由エネルギーが、自由エネルギーの平均および分散の少なくとも一方である、請求項11から18のいずれか一項に記載の予測方法。請求項11から19のいずれか一項に記載の予測方法を、コンピュータ上で実行可能なことを特徴とするプログラム。請求項20記載のプログラムを記録していることを特徴とする記録媒体。そこで、本発明は、ライブラリーから目的アプタマー配列を選択する際に、容易に、前記目的アプタマー配列の濃縮の有無を予測できる予測装置、予測方法、プログラムおよび記録媒体を提供することを目的とする。前記目的を達成するために、本発明の予測装置は、複数の候補アプタマー配列から選択された選択アプタマー配列を含む被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測装置であって、入力手段と、自由エネルギーの算出手段と、目的アプタマー配列の濃縮の有無を予測する予測手段とを備え、前記入力手段は、被対照プールに含まれる被対照アプタマー配列群の配列情報、および対照アプタマー配列群の配列情報を入力する入力工程を実行し、前記算出手段は、前記被対照アプタマー配列群の配列情報および前記対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程を実行し、前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを前記目的アプタマー配列が濃縮されたプールと予測する予測工程を実行し、前記対照アプタマー配列群が、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする。また、本発明の予測方法は、複数候補アプタマー配列から選択された選択アプタマー配列を有する被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測方法であって、被検対照プールに含まれる被対照アプタマー配列群の配列情報および対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程と、前記被対照アプタマー配列群の自由エネルギーと、前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが、前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを、目的アプタマー配列が濃縮されたプールと予測する予測工程とを有し、前記対照アプタマー配列群が、1複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする。また、本発明のプログラムは、前記本発明の予測方法を、コンピュータ上で実行可能なことを特徴とする。また、本発明の記録媒体は、前記本発明のプログラムを記録していることを特徴とする。本発明によれば、前記被対照アプタマー配列群と前記対照アプタマー配列群との自由エネルギーを比較するのみで、容易に前記目的アプタマー配列の濃縮の有無を予測できる。このため、例えば、目的アプタマーを濃縮するための選択処理を、必要以上のラウンド数繰り返したり、濃縮のためのラウンド数が不足すること等を防止できる。また、本発明により濃縮の有無を予測すれば、例えば、必要と考えられるプールについてのみ、ウェットでの確認試験を行えば足りるため、労力やコストを軽減できる。本発明者らは、鋭意研究の結果、目的アプタマー配列が濃縮されたプールは、目的アプタマー配列が濃縮されていないプールと比較して、自由エネルギーが低い値を示すことを見出した。その結果、任意の2種類のプールについて、それぞれの自由エネルギーを比較し、一方のプールが、他方のプールよりも、低い自由エネルギーを示していれば、前者のプールは、後者のプールよりも、目的アプタマー配列が濃縮されていると予測することが可能となった。また、さらに、プールに含まれるアプタマー配列群と同じ塩基組成となる仮想アプタマー配列を含む仮想アプタマー配列群を生成した場合、目的アプタマー配列が濃縮されたプールは、前記プールに含まれるアプタマー配列群の自由エネルギーが、前記仮想アプタマー配列群の自由エネルギーよりも低い値を示すことを見出した。その結果、任意のプールについて、前記プールに含まれるアプタマー配列群の自由エネルギーと、その仮想アプタマー配列群の自由エネルギーとを比較し、前記プールに含まれるアプタマー配列群が、前記仮想アプタマー配列群よりも低い自由エネルギーを示していれば、前記プールは、目的アプタマー配列が濃縮されていると予測することが可能となった。本発明において、目的アプタマー配列とは、任意のターゲットに結合可能な核酸の配列であり、候補アプタマー配列とは、ターゲットに対する結合性が不明な核酸の配列であって、前記目的アプタマー配列の候補となる配列であり、選択アプタマー配列とは、候補アプタマー配列のうち、任意の条件を満たすアプタマー配列である。本発明において、目的アプタマー配列の濃縮とは、例えば、複数の候補アプタマー配列を含むプールPnから、任意の条件を満たすアプタマー配列を選択して、前記選択アプタマー配列を含むプールPmを得た際、前記選択前の前記プールPnにおける目的アプタマー配列の割合よりも、前記選択後の前記プールPmにおける目的アプタマー配列の割合重複数ともいうが高くなることを意味する。また、目的アプタマー配列の濃縮は、例えば、プールに含まれる核酸の配列が、目的アプタマー配列に収束するということもできる。本発明において、被対照プールは、目的アプタマー配列の濃縮の有無を予測する、予測対象のプールであり、複数の候補アプタマー配列から選択された選択アプタマー配列を含む。本発明において、被対照アプタマー配列群は、前記被対照プールにおける選択アプタマー配列を含むアプタマー配列群であり、例えば、以下、選択アプタマー配列群ともいう。前記被対照アプタマー配列群は、例えば、前記被対照プールに含まれる全ての選択アプタマー配列からなるアプタマー配列群でもよいし、一部の選択アプタマー配列からなるアプタマー配列群でもよい。本発明において、対照アプタマー配列群は、前記被対照プールについて前記予測を行うにあたり、前記被対照アプタマー配列群の比較対象となるアプタマー配列群である。前記対照アプタマー配列群の詳細は、後述する。本発明は、複数の候補アプタマー配列からなる対照プールPnに、任意の回数m回、任意のスクリーニングを施し、前記候補アプタマー配列の選択を行い、選択された候補アプタマー配列選択アプタマー配列からなる被対照プールPnmについて、目的アプタマー配列が、前記対照プールPnよりも濃縮されているか否かを予測するために使用できる。前記対照プールPnのnは、そのプールPnが、任意のスクリーニングによって選択をうけたプールであるか否か、および、スクリーニングの回数ラウンドを示し、0および正の整数を示す。すなわち、例えば、P0n=0は、任意のスクリーニングをうけていない未選択の候補アプタマー配列からなるプールである。また、P1n=1は、例えば、前記P0に、1回の任意のスクリーニングを施すことにより選択された候補アプタマー配列からなるプールである。また、P2n=2は、前記P1に、さらに1回の任意のスクリーニングを施すことにより選択された候補アプタマー配列からなるプールである。他方、前記プールPmは、そのプールPmが、前記プールPnから、m回の任意のスクリーニングによって選択を受けたことを示し、mは、正の整数を示す。すなわち、nが0、mが1の場合、Pnm=P1は、任意のスクリーニングをうけていない未選択の候補アプタマー配列からなるプールP0に、1回の任意のスクリーニングを施すことにより選択された選択アプタマー配列からなるプールである。前記任意のスクリーニングは、制限されず、例えば、SELEX法における選択処理があげられる。本発明は、SELEX法において、各ラウンドで得られるプールについて、目的アプタマー配列の濃縮の有無を予測する際に、特に有用である。本発明の予測装置は、前述のように、複数の候補アプタマー配列から選択された選択アプタマー配列を含む被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測装置であって、入力手段と、自由エネルギーの算出手段と、目的アプタマー配列の濃縮の有無を予測する予測手段とを備え、前記入力手段は、被対照プールに含まれる被対照アプタマー配列群の配列情報、および対照アプタマー配列群の配列情報を入力する入力工程を実行し、前記算出手段は、前記被対照アプタマー配列群の配列情報および前記対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程を実行し、前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを前記目的アプタマー配列が濃縮されたプールと予測する予測工程を実行し、前記対照アプタマー配列群が、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする。本発明の予測装置は、例えば、前記対照アプタマー配列群が、前記1候補アプタマー配列群であり、前記予測手段は、算出された前記被対照アプタマー配列群の自由エネルギーと算出された前記候補アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記1候補アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。本発明の予測装置は、例えば、前記対照アプタマー配列群が、前記2被対照アプタマー配列群由来の仮想アプタマー配列群であり、前記予測手段は、算出された前記被対照アプタマー配列群の自由エネルギーと算出された前記2仮想アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記2仮想アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。本発明の予測装置は、例えば、さらに、前記仮想アプタマー配列の生成手段を備え、前記生成手段は、前記入力された前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成する生成工程を実行する。本発明の予測装置は、例えば、前記対照アプタマー配列群が、前記1候補アプタマー配列群、前記2仮想アプタマー配列群、および下記3仮想アプタマー配列群であり、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群3前記候補アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記候補アプタマー配列群由来の仮想アプタマー配列群前記予測手段は、前記被対照アプタマー配列群の自由エネルギーと前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記1候補アプタマー配列群の自由エネルギーと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記被対照アプタマー配列群の自由エネルギーの減少程度が、前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記1候補アプタマー配列群の自由エネルギーの減少程度よりも大きい場合、前記被対照プールを前記目的アプタマー配列が濃縮されたプールと予測する前記予測工程を実行する。本発明の予測装置は、例えば、さらに、前記仮想アプタマー配列の生成手段を備え、前記生成手段は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成し、前記候補アプタマー配列群の配列情報から、前記候補アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記3仮想アプタマー配列群を生成する生成工程を実行する。本発明の予測装置は、例えば、前記仮想アプタマー配列群の配列数が、100個~1000万個である。本発明の予測装置は、例えば、前記被対照アプタマー配列群および前記候補アプタマー配列群の配列数が、それぞれ、5個~1億個である。本発明の予測装置は、例えば、前記自由エネルギーが、自由エネルギーの平均および分散の少なくとも一方である。本発明の予測装置は、例えば、さらに、出力手段を備え、前記出力手段は、前記予測手段で予測された予測結果を出力する。本発明の予測方法は、前述のように、複数候補アプタマー配列から選択された選択アプタマー配列を有する被対照プールについて、目的アプタマー配列の濃縮の有無を予測する予測方法であって、被検対照プールに含まれる被対照アプタマー配列群の配列情報および対照アプタマー配列群の配列情報から、それぞれ、前記被対照アプタマー配列群の自由エネルギーおよび前記対照アプタマー配列群の自由エネルギーを算出する算出工程と、前記被対照アプタマー配列群の自由エネルギーと、前記対照アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが、前記対照アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを、目的アプタマー配列が濃縮されたプールと予測する予測工程とを有し、前記対照アプタマー配列群が、1複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群の少なくとも一方であることを特徴とする。本発明の予測方法は、前記各工程が、全てコンピュータ上で実行される。本発明の予測方法は、例えば、前記対照アプタマー配列群が、前記1の候補アプタマー配列群であり、前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記1候補アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記1候補アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。本発明の予測方法は、前記対照アプタマー配列群が、前記2被対照アプタマー配列群由来の仮想アプタマー配列群であり、前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記2仮想アプタマー配列群の自由エネルギーとを比較して、前記被対照アプタマー配列群の自由エネルギーが前記2仮想アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。本発明の予測方法は、さらに、仮想アプタマー配列の生成工程を有し、前記生成工程は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成する工程である。本発明の予測方法は、前記対照アプタマー配列群が、前記1候補アプタマー配列群、前記2仮想アプタマー配列群、および下記3仮想アプタマー配列群であり、1前記複数の候補アプタマー配列を含む候補アプタマー配列群、および、2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群3前記候補アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記候補アプタマー配列群由来の仮想アプタマー配列群前記予測工程は、前記被対照アプタマー配列群の自由エネルギーと前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記1候補アプタマー配列群の自由エネルギーと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較し、前記2被対照アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記被対照アプタマー配列群の自由エネルギーの減少程度が、前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーに対する前記1候補アプタマー配列群の自由エネルギーの減少程度よりも大きい場合、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。本発明の予測方法は、さらに、生成工程を有し、前記生成工程は、前記被対照アプタマー配列群の配列情報から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記2仮想アプタマー配列群を生成し、前記候補アプタマー配列群の配列情報から、前記候補アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成することにより、前記3仮想アプタマー配列群を生成する工程である。本発明の予測方法は、例えば、前記仮想アプタマー配列群の配列数が、100個~1000万個である。本発明の予測方法は、例えば、前記被対照アプタマー配列群および前記候補アプタマー配列群の配列数が、それぞれ、5個~1億個である。本発明の予測方法は、例えば、前記自由エネルギーが、自由エネルギーの平均および分散の少なくとも一方である。本発明のプログラムは、前記本発明の予測方法を、コンピュータ上で実行可能なことを特徴とする。本発明の記録媒体は、本発明のプログラムを記録していることを特徴とする。以下、図面を参照しながら本発明のさらに具体的な実施形態について説明する。ただし、本発明は、以下の実施形態に限定されない。[実施形態1]本実施形態は、前記対照アプタマー配列群を、前記1候補アプタマー配列群とする一例である。前記候補アプタマー配列群は、複数の候補アプタマー配列を含む。本実施形態によれば、例えば、前記複数の候補アプタマー配列からなるプールPnから、選択により前記選択アプタマー配列からなるプールPnmを得た際に、前記目的アプタマー配列が、前記プールPnmにおいて、前記プールPnよりも濃縮されているか否かを予測できる。本実施形態において、前記複数の候補アプタマー配列からなるプールPnを、前記非対照プールの比較対象となる対照プールまたは候補プールという。そして、前記対照アプタマー配列群、すなわち候補アプタマー配列群は、前記候補プールにおける前記候補アプタマー配列のうち、任意の前記候補アプタマー配列からなる群である。前記候補アプタマー配列群は、前記候補プールにおける全ての前記候補アプタマー配列からなる群でもよいし、一部の前記候補アプタマー配列からなる群でもよい。前記被対照プールPnmは、前記複数の候補アプタマー配列から選択された選択アプタマー配列からなるプールであり、選択プールともいう。そして、前記被対照アプタマー配列群は、前記被対照プールにおける前記選択アプタマー配列のうち、任意の前記選択アプタマー配列からなる群であり、選択アプタマー配列群ともいう。前記被対照アプタマー配列群は、前記被対照プールにおける全ての前記選択アプタマー配列からなる群でもよいし、一部の前記選択アプタマー配列からなる群でもよい。以下、本実施形態において、前記対照アプタマー配列は候補アプタマー配列、前記対照アプタマー配列群は候補アプタマー配列群、前記被対照アプタマー配列は選択アプタマー配列、前記被対照アプタマー配列群は選択アプタマー配列群ともいう。図1Aに、本実施形態の予測装置の一例の構成を示す。図1Aに示すとおり、予測装置10は、入力手段11と自由エネルギー算出手段12と濃縮予測手段13と出力手段14とを備える。自由エネルギー算出手段12と濃縮予測手段13とは、例えば、同図に示すように、ハードウェアであるデータ処理装置15に組み込まれていてもよい。データ処理装置15は、例えば、CPU等を備えていてもよい。入力手段11は、例えば、自由エネルギー算出手段12に電気的に接続されている。自由エネルギー算出手段12は、例えば、濃縮予測手段13に電気的に接続されている。濃縮予測手段13は、例えば、出力手段14に電気的に接続されている。入力手段11は、自由エネルギー算出手段12に、前記候補アプタマー配列群対照アプタマー配列群および前記選択アプタマー配列群被対照アプタマー配列群の配列情報を入力する手段である。前記入力手段11は、特に制限されず、例えば、キーボード、マウス等のコンピュータに備わる通常の入力装置、入力ファイルおよび他のコンピュータ等を用いることができる。入力手段11は、例えば、図1Bに示すように、データベースに格納された、前記候補アプタマー配列群の配列情報および前記選択アプタマー配列群の配列情報を読み出す手段でもよい。この場合、例えば、同図に示すように、予めサーバ16に格納された配列情報が、回線網17を通じて、入力手段11に呼び出される。また、入力手段11は、例えば、図1Bに示すように、通信インターフェース111を含んでもよい。図1Bにおいて、図1Aと同一部分には同一符号を付している。入力手段11により配列情報を入力するアプタマー配列の数は、特に制限されない。前記候補アプタマー配列群の配列情報として、例えば、前記対照プールに含まれる全ての候補アプタマー配列の配列情報を入力してもよいし、一部の候補アプタマー配列の配列情報を入力してもよい。また、前記選択アプタマー配列群の配列情報として、例えば、前記被対照プールに含まれる全ての選択アプタマー配列の配列情報を入力してもよいし、一部の選択アプタマー配列の配列情報を入力してもよい。具体例として、入力手段11により、配列情報を入力する前記候補アプタマー配列群の配列数および前記選択アプタマー配列群の配列数は、それぞれ、例えば、5個~1億個であり、好ましくは、10個~10万個であり、より好ましくは、20個~200個である。前記目的アプタマー配列の濃縮は、通常、アプタマー配列のライブラリーについて、複数ラウンドのスクリーニングを行うことで実施できる。このため、濃縮にあたって、例えば、スクリーニング開始前の初期プールラウンド0および各ラウンドのプールが得られる。本発明においては、例えば、スクリーニングを行っていない初期プールおよびスクリーニングを行った各ラウンドのプールから、前記対照プールおよび前記被対照プールを設定できる。前記対照プールPnは、いずれのラウンドのプールを設定してもよい。前記被対照プールPnmは、前述のように、選択された前記選択アプタマー配列を含むプールであることから、例えば、前記対照プールのラウンド以降のプールを設定できる。前記対照プールPnおよび前記被対照プールPnmにおいて、nおよびmは、スクリーニングの回数ラウンドを示し、nは、0および正の整数を示し、mは、正の整数を示す。具体的には、前記対照プールPnは、例えば、スクリーニングを行っていないラウンドn=0のプールP0または1回以上のスクリーニングを行ったラウンドn≧1のプールPn≧1を設定できる。一方、前記被対照プールPnmは、例えば、さらにm回のスクリーニングを行った、ラウンドn以降のラウンドnmのプールPnmnは、0および正の整数、mは、正の整数を設定できる。前記被対照プールPnmは、例えば、前記対照プールがラウンド0のプールP0の場合、ラウンド0以降のラウンドのプールPm≧1を設定でき、具体的には、ラウンド1のプールP1および/またはラウンド1以降のラウンドのプールPm≧2を設定できる。本実施形態において、例えば、前記対照プールおよび前記被対照プールとして設定した、異なる任意のラウンドの2種類の各プールについて、アプタマー配列の配列情報を入力するだけでなく、同様に、その他の異なる任意のラウンドの1種類以上のプールについて、アプタマー配列の配列情報を入力してもよい。また、例えば、全てのラウンドのプールについて、アプタマー配列の配列情報を入力してもよい。このような場合、例えば、入力後に、前記対照プールおよび前記被対照プールを設定して、後述するように、それぞれの自由エネルギーを算出してもよい。また、例えば、入力した各ラウンドのプールについて、自由エネルギーを算出した後、前記対照プールおよび前記被対照プールを設定して、後述するように、自由エネルギーの比較を行ってもよい。なお、前述のように、各プールについて、入力する配列情報は、例えば、前記プールに含まれる全てのアプタマー配列の配列情報でもよいし、一部のアプタマー配列の配列情報でもよい。前記候補アプタマー配列および前記選択アプタマー配列の配列情報は、例えば、ウェット系で実際に作製したアプタマーの配列情報でもよいし、ドライ系insilicoで設計したアプタマーの配列情報でもよい。具体例として、前記候補アプタマー配列がウェット系で作製したアプタマーの配列の場合、例えば、ウェット系で、前記候補アプタマーから任意の条件により選択したアプタマーを、選択アプタマーとし、その配列の情報を、前記選択アプタマー配列の配列情報とすることができる。ウェット系でアプタマーを選択する方法は、特に制限されず、例えば、SELEX法があげられる。他方、前記候補アプタマー配列がドライ系で設計したアプタマーの配列情報の場合、例えば、ドライ系で、前記候補アプタマー配列から任意の条件により選択したアプタマー配列を、選択アプタマー配列とすることができる。前記ドライ系でアプタマー配列を選択する方法は、特に制限されず、例えば、前記候補アプタマー配列の配列情報から、任意の特徴を有するアプタマー配列を選択することにより行える。前記特徴は、特に制限されず、適宜設定できる。前記特徴は、例えば、塩基の並び方、各塩基の出現頻度、推定二次構造、各塩基の含有率、連続する塩基の出願頻度、モチーフ配列の有無、融解温度等があげられる。自由エネルギー算出手段12は、入力された前記対照プールの候補アプタマー配列群および入力された前記被対照プールの選択アプタマー配列群の配列情報から、それぞれ、前記候補アプタマー配列群の自由エネルギーおよび前記選択アプタマー配列群の自由エネルギーを算出する。前記候補アプタマー配列群の自由エネルギーは、例えば、前記対照プールの自由エネルギーということもでき、前記選択アプタマー配列群の自由エネルギーは、例えば、前記被対照プールの自由エネルギーということもできる。前記自由エネルギーは、例えば、入力した配列情報から算出できる。前記自由エネルギーは、例えば、前記配列に関する、塩基の並び、推定二次構造、各塩基の含有率等の配列情報、および温度、塩濃度等の溶媒情報等に基づいて算出できる。前記算出方法は、特に制限されず、従来公知の方法により行うことができ、例えば、MfoldMZukerMfoldwebserverfornucleicacidfoldingandhybridizationpredictionNucleicAcidsRes313406152003、RNAfoldW70W74Jul2008を用いて、自由エネルギー最小化法等を使用する方法等があげられる。前記自由エネルギーは、例えば、最小自由エネルギーMinimumFreeEnergy:MFEが好ましい。前記候補アプタマー配列群の自由エネルギーは、例えば、前記候補アプタマー配列群に含まれる各候補アプタマー配列の自由エネルギーの平均でもよいし、分散でもよい。具体的には、入力された複数の候補アプタマー配列の配列情報について、それぞれ、自由エネルギーを算出し、さらに、それらの平均値またはそれらの分散値を算出して、前記候補アプタマー配列群の自由エネルギーの平均または分散としてもよい。また、前記選択アプタマー配列の自由エネルギーも、同様に、例えば、前記選択アプタマー配列群に含まれる各選択アプタマー配列の自由エネルギーの平均でもよいし、分散でもよい。具体的には、入力された複数の選択アプタマー配列の配列情報について、それぞれ、自由エネルギーを算出し、さらに、それらの平均値またはそれらの分散値を算出して、前記選択アプタマー配列群の自由エネルギーの平均または分散としてもよい。自由エネルギーを算出する前記候補アプタマー配列の数および前記選択アプタマー配列の数は、特に制限されず、それぞれ、例えば、5個~1億個であり、好ましくは、10個~10万個であり、より好ましくは、20個~200個である。濃縮予測手段13は、自由エネルギー算出手段12から、前記自由エネルギーの情報を取り出し、前記対照プールの候補アプタマー配列群の自由エネルギーと、前記被対照プールの選択アプタマー配列群の自由エネルギーとを比較する。そして、濃縮予測手段13は、前記選択アプタマー配列群の自由エネルギーが前記候補アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。つまり、前記対照プールよりも、前記被対照プールは、目的アプタマーが濃縮されたプールであると予測できる。また、予想結果の信頼性を向上できることから、例えば、前記被対照プールの選択アプタマー配列群の自由エネルギーが、前記対照プールの候補アプタマー配列群の自由エネルギーよりも有意に低い場合に、前記被対照プールを、目的アプタマー配列が濃縮されたプールと予測することが好ましい。有意か否かの判断方法は、特に制限されず、例えば、t検定、KS検定等があげられる。例えば、前記t検定の場合、有意水準5P=0.05で有意と判断できる。出力手段14は、濃縮予測手段13の予測結果を出力する手段であればよい。前記出力手段14は、特に制限されず、例えば、ディスプレイ装置、印刷装置等のコンピュータに備わる通常の出力装置、出力ファイル、および、他のコンピュータ等を使用できる。つぎに、図2および図3のフローチャートを参照し、本実施形態の予測方法を説明する。本実施形態の予測方法は、自由エネルギー算出工程S21と、濃縮予測工程S22とを含む。[自由エネルギー算出工程]前記対照プールの候補アプタマー配列群の配列情報および前記被対照プールの選択アプタマー配列群の配列情報から、それぞれ、前記候補アプタマー配列群の自由エネルギーCRおよび前記選択アプタマー配列群の自由エネルギーSRを算出するS301。前記自由エネルギー算出工程は、例えば、自由エネルギー算出手段12によって実行できる。以下、CRを、前記対照プールの自由エネルギー、SRを、前記被対照プールの自由エネルギーともいう以下、同様。[濃縮予測工程]前記対照プールの候補アプタマー配列群の自由エネルギーCRと前記被対照プールの選択アプタマー配列群の自由エネルギーSRとを比較し、前記選択アプタマー配列群の自由エネルギーSRが前記候補アプタマー配列群の自由エネルギーCRよりも低い値であるか否か、すなわち、SR<CRであるか否かを判断するS302。SR;選択アプタマー配列群の自由エネルギーCR:候補アプタマー配列群の自由エネルギーそして、前記選択アプタマー配列群の自由エネルギーSRが、前記候補アプタマー配列群の自由エネルギーCRよりも低い値であれば、前記被対照プールは、前記対照プールよりも目的アプタマーが濃縮されていると予測しS303、前記候補アプタマー配列群の自由エネルギーCRと同じまたは高い値であれば、前記被対照プールは、目的アプタマーが濃縮されていないと予測するS304。以上の工程により、濃縮予測が完了する。また、前述のように、例えば、2種類以上または全てラウンドのプールについて、アプタマー配列群の配列情報を入力した場合、例えば、複数のプールのアプタマー配列群から、任意のプールのアプタマー配列群を、対照プールの候補アプタマー配列群および被対照プールの選択アプタマー配列群として選択する工程、選択した前記候補アプタマー配列群および前記選択アプタマー配列群についての自由エネルギー算出工程、および濃縮予測工程を、繰り返し行ってもよい。また、例えば、前記複数のプールの各アプタマー配列群について、自由エネルギーの算出工程を行った後、任意のプールのアプタマー配列群を、対照プールの候補アプタマー配列群および被対照プールの選択アプタマー配列群として選択する工程、および濃縮予測工程を、繰り返し行ってもよい。この際、前記対照プールは、例えば、各濃縮予測工程において同じプールでもよいし、異なるプールであってもよい。これによって、前記複数のプールから、目的アプタマー配列が濃縮されたプールを予測可能である。予測結果は、例えば、ユーザーが指定した条件に従って、出力手段14で出力してもよい。さらに、出力手段14は、例えば、前記予測結果の他に、前記候補アプタマー配列の配列情報、前記選択アプタマー配列の配列情報、自由エネルギー等を出力する機能を備えていてもよい。[実施形態2]図4に、本実施形態の予測装置40の一例の構成を示す。同図において、図1と同一部分には同一符号を付している。特に示さない限り、本実施形態は、前記実施形態1と同様である。図4に示すとおり、予測装置40は、さらに、配列記憶手段41および自由エネルギー記憶手段42を備えること以外は、実施形態1の予測装置と同様の構成を有する。すなわち、予測装置40は、入力手段11と自由エネルギー算出手段12と濃縮予測手段13と出力手段14と配列記憶手段41と自由エネルギー記憶手段42とを備える。自由エネルギー算出手段12と濃縮予測手段13とは、例えば、図4に示すように、ハードウェアであるデータ処理装置15に組み込まれてもよい。自由エネルギー算出手段12と濃縮予測手段13とは、例えば、ソフトウェアそのものでも良いし、ハードウェアにソフトウェアが組み込まれたものでもよい。データ処理装置15は、例えば、CPU等を備えてもよい。配列記憶手段41と自由エネルギー記憶手段42とは、例えば、図4に示すように、ハードウェアである記憶装置43に組み込まれてもよい。図4に示すように、配列記憶手段41は、例えば、入力手段11および自由エネルギー算出手段12に電気的に接続されている。また、自由エネルギー記憶手段42は、例えば、自由エネルギー算出手段12および濃縮予測手段13に電気的に接続されている。配列記憶手段41は、例えば、入力手段11から入力された、アプタマー配列の情報を記憶し、自由エネルギー算出手段12に前記情報を出力する。自由エネルギー記憶手段42は、自由エネルギー算出手段12で算出された自由エネルギーを、それぞれのプール、すなわち、対照プールおよび被対照プールに関連付けて記憶し、濃縮予測手段13に前記自由エネルギーの情報を出力する。また、前述のように、入力手段11により、複数のラウンドのプールについて、アプタマー配列群の配列情報を入力した場合、配列記憶手段41は、例えば、全てのラウンドのプールまたは2種類以上のラウンドのプールについて、アプタマー配列の配列情報を記憶し、自由エネルギー算出手段12に前記情報を出力してもよい。そして、自由エネルギー記憶手段42は、自由エネルギー算出手段12で算出された自由エネルギーを、それぞれのプールに関連付けて記憶し、その中から、前記対照プールおよび前記被対照プールとして選択した任意の2種類のプールについて、濃縮予測手段13に前記自由エネルギーの情報を出力してもよい。[実施形態3]本実施形態は、前記対照アプタマー配列群を、前記2仮想アプタマー配列群とする例である。前記2仮想アプタマー配列群は、前記仮想アプタマー配列群は、前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群である。以下、本実施形態において、前記対照アプタマー配列は仮想アプタマー配列、前記対照アプタマー配列群は仮想アプタマー配列群、前記選択アプタマー配列を含む被対照アプタマー配列群は選択アプタマー配列群ともいう。図5Aに、本実施形態の予測装置の一例の構成を示す。図5Aにおいて、図1Aと同一部分には同一符号を付している。特に示さない限り、本実施形態は、前記実施形態1と同様である。図5Aに示すとおり、予測装置50は、仮想アプタマー生成手段54を有する以外は、前記図1Aに示す実施形態1の予測装置10と同様である。具体的には、予測装置50は、入力手段11と仮想アプタマー生成手段54と自由エネルギー算出手段12と濃縮予測手段13と出力手段14とを備える。仮想アプタマー生成手段54と自由エネルギー算出手段12と濃縮予測手段13とは、例えば、図5Aに示すように、ハードウェアであるデータ処理装置15に組み込まれてもよい。仮想アプタマー生成手段54と自由エネルギー算出手段12と濃縮予測手段13とは、例えば、ソフトウェアそのものでも良いし、ハードウェアにソフトウェアが組み込まれたものでもよい。入力手段11は、例えば、図5Bに示すように、データベースに格納された、前記候補アプタマー配列群の配列情報および前記選択アプタマー配列群の配列情報を読み出す手段であってもよい。この場合、例えば、同図に示すように、予めサーバ16に格納された配列情報が、回線網17を通じて、入力手段11に呼び出される。また、入力手段11は、例えば、図5Bに示すように、通信インターフェース111を含んでもよい。図5Bにおいて、図5Aと同一部分には同一符号を付している。図5Aに示すように、仮想アプタマー生成手段54は、例えば、入力手段11および自由エネルギー算出手段12に電気的に接続されている。仮想アプタマー生成手段54は、入力された前記被対照プールの選択アプタマー配列群の配列情報から、前記選択アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を生成し、前記仮想アプタマー配列を含む前記仮想アプタマー配列群を生成する。前記仮想アプタマー配列の生成は、例えば、前記被対照プールの選択アプタマー配列群の塩基組成に基づいて行うことができる。具体的には、例えば、前記選択アプタマー配列群に含まれる各選択アプタマー配列の塩基配列の情報から、前記選択アプタマー配列群における各塩基の出現頻度を求め、これを塩基組成として設定する。そして、この塩基組成となり得る塩基配列をランダムに生成することで、同じ塩基組成である複数の仮想アプタマー配列を生成して、前記仮想アプタマー配列からなる前記仮想アプタマー配列群を得ることができる。また、本実施形態において、入力手段11は、例えば、前記被対照プールの選択アプタマー配列群の配列情報として、塩基組成を入力してもよい。この場合、例えば、入力された塩基組成に基づいて、仮想アプタマー配列を生成できる。生成する仮想アプタマーの数は、特に制限されないが、例えば、100個~1000万個であり、好ましくは、1000個~100万個であり、より好ましくは、1万個~10万個である。また、前記仮想アプタマー配列群の配列情報は、例えば、入力手段11により入力してもよい。自由エネルギー算出手段12は、入力された前記被対照アプタマー配列群由来の仮想アプタマー配列群の配列情報から、前記仮想アプタマー配列群の自由エネルギーを算出する。前記自由エネルギーの算出は、例えば、前述の通りである。濃縮予測手段13は、自由エネルギー算出手段12から、前記自由エネルギーの情報を取り出し、前記選択アプタマー配列群の自由エネルギーと、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較する。そして、濃縮予測手段13は、前記選択アプタマー配列群の自由エネルギーが前記仮想アプタマー配列群の自由エネルギーよりも低い場合に、前記被対照プールを目的アプタマー配列が濃縮されたプールと予測する。また、予想結果の信頼性を向上できることから、例えば、前記選択アプタマー配列群の自由エネルギーが、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーよりも有意に低い場合に、前記被対照プールを、目的アプタマー配列が濃縮されたプールと予測することが好ましい。有意か否かの判断方法は、特に制限されず、前述と同様である。また、本実施形態においては、例えば、前記実施形態1に示すように、前記被対照プールの選択アプタマー配列群の自由アプタマーと、前記対照プールの候補アプタマー配列群の自由エネルギーとの比較を、あわせて行ってもよい。本実施形態の予測装置は、例えば、前記実施形態2における図4に示す予測装置40と同様に、さらに、配列記憶手段41および自由エネルギー記憶手段42を備えてもよい。配列記憶手段41は、例えば、前記被対照プールの選択アプタマー配列群の配列情報の他に、仮想アプタマー生成手段で生成された前記仮想アプタマー配列群の配列情報を記憶し、自由エネルギー算出手段12に前記情報を出力してもよい。また、配列記憶手段41は、さらに、前記対照プールの候補アプタマー配列群の配列情報を記憶し、自由エネルギー算出手段12に前記配列情報を出力してもよい。自由エネルギー記憶手段42は、例えば、自由エネルギー算出手段12で算出された、前記被対照プールの選択アプタマー配列群の自由エネルギーの他に、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーを、前記被対照プールに関連付けて記憶し、濃縮予測手段13に前記自由エネルギーの情報を出力してもよい。また、自由エネルギー記憶手段42は、さらに、前記対照プールの候補アプタマー配列群の自由エネルギーを、前記対照プールに関連付けて記憶し、濃縮予測手段13に出力してもよい。つぎに、図6および図7のフローチャートを参照し、本実施形態の予測方法を説明する。本実施形態の予測方法は、仮想アプタマー生成工程S60、自由エネルギー算出工程S61と、濃縮予測工程S62とを含む。[仮想アプタマー生成工程]前記被対照プールの選択アプタマー配列群から、前記被対照アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を含む前記仮想アプタマー配列群を生成するS700。[自由エネルギー算出工程]そして、前記被対照プールの選択アプタマー配列群の配列情報および前記選択アプタマー配列群由来の仮想アプタマー配列群の配列情報から、それぞれ、前記選択アプタマー配列群の自由エネルギーSRおよび前記仮想アプタマー配列群の自由エネルギーSVを算出するS701。この自由エネルギー算出工程は、例えば、自由エネルギー算出手段12によって実行できる。以下、SRを、前記被対照プールの自由エネルギー、SVを、前記対照プールの自由エネルギーともいう以下、同様。[濃縮予測工程]算出した前記被対照プールの選択アプタマー配列群の自由エネルギーSRと前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVとを比較し、前記選択アプタマー配列群の自由エネルギーSRが、前記仮想アプタマー配列群の自由エネルギーSVよりも低い値であるか否か、すなわち、SR<SVであるか否かを判断するS702。SR:選択アプタマー配列群の自由エネルギーSV:仮想アプタマー配列群の自由エネルギーそして、前記被対照プールの選択アプタマー配列群の自由エネルギーSRが、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVよりも低い値であれば、前記被対照プールは、目的アプタマーが濃縮されていると予測しS703、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVと同じまたは高い値であれば、前記被対照プールは、目的アプタマーが濃縮されていないと予測するS704。以上の工程により、濃縮予測が完了する。本実施形態では、例えば、前記被対照プールの選択アプタマー配列群の自由エネルギーと、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較するのみで、濃縮を予測できる。また、本実施形態では、例えば、前記実施形態1と同様に、あわせて、前記自由エネルギー算出工程において、前記対照プールの候補アプタマー配列群の自由エネルギーCRを算出してもよい。そして、前記濃縮予測工程において、前記被対照プールの選択アプタマー配列群の自由エネルギーSRと前記対照プールの候補アプタマー配列群の自由エネルギーCRとの比較により、濃縮予測を行うこともできる。[実施形態4]本実施形態は、前記対照アプタマー配列群を、前記1候補アプタマー配列群、前記2仮想アプタマー配列群、および前記3仮想アプタマー配列群とする例である。なお、特に示さない限り、前記実施形態1、実施形態2および実施形態3と同様にして行うことができる。1前記複数の候補アプタマー配列を含む候補アプタマー配列群2前記被対照アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記被対照アプタマー配列群由来の仮想アプタマー配列群3前記候補アプタマー配列群と同じ塩基組成の仮想アプタマー配列を含む、前記候補アプタマー配列群由来の仮想アプタマー配列群本実施形態の予測方法は、仮想アプタマー生成工程、自由エネルギー算出工程と、濃縮予測工程とを含む。[仮想アプタマー生成工程]前記被対照プールの選択アプタマー配列群から、前記選択アプタマー配列群と同じ塩基組成となる仮想アプタマー配列を含む、前記2選択アプタマー配列群由来の仮想アプタマー配列群を生成する。他方、前記1対照プールの候補アプタマー配列群から、前記候補アプタマー配列群と同じ組成となる仮想アプタマー配列を含む、前記3候補アプタマー配列群由来の仮想アプタマー配列群を生成する。[自由エネルギー算出工程]そして、前記被対照プールの選択アプタマー配列群および前記2選択アプタマー配列群由来の仮想アプタマー配列群から、それぞれ、前記選択アプタマー配列群の自由エネルギーSRおよび前記2選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVを算出する。他方、前記1対照プールの候補アプタマー配列群および前記3候補アプタマー配列群由来の仮想アプタマー配列群から、それぞれ、前記1候補アプタマー配列群の自由エネルギーCRおよび前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーCVを算出する。これらの自由エネルギー算出工程は、例えば、自由エネルギー算出手段12によって実行できる。[濃縮予測工程]前記被対照プールの選択アプタマー配列群の自由エネルギーSRと前記2選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVとを比較し、前記選択アプタマー配列群の自由エネルギーSRが前記2仮想アプタマー配列群の自由エネルギーSVよりも低い値であるか否か、すなわち、SR<SVであるか否かを判断する。SR:選択アプタマー配列群の自由エネルギーSV:仮想アプタマー配列群の自由エネルギー他方、算出した前記1対照プールの候補アプタマー配列群の自由エネルギーCRと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーCVとを比較し、前記1候補アプタマー配列群の自由エネルギーCRが前記3仮想アプタマー配列群の自由エネルギーCVよりも低い値であるか否か、すなわち、CR<CVであるか否かを判断する。CR:候補アプタマー配列群の自由エネルギーCV:仮想アプタマー配列群の自由エネルギーつぎに、前記被対照プールについて、前記2選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーSVに対する前記選択アプタマー配列群の自由エネルギーSRの減少程度SV-SRを求める。また、前記対照プールについて、前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーCVに対する前記1候補アプタマー配列群の自由エネルギーCRの減少程度CV-CRを求める。そして、前記被対照プールの減少程度SV-SRが、前記対照プールの減少程度CV-CRよりも大きい場合、前記被対照プールを、前記対照プールよりも、目的アプタマー配列が濃縮されたプールと予測できる。[実施形態5]前記被対照プールについて、前記選択アプタマー配列群の自由エネルギーと前記2選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較する場合、例えば、以下の方法によって、両者の自由エネルギーの差が有意であるかを判断できる。これらは一例であって、本発明を制限するものではない。また、前記対照プールについて、前記1候補アプタマー配列群の自由エネルギーと前記3候補アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーとを比較する場合も、同様である。1二群の比較前記被対照プールに含まれる選択アプタマー配列群Sの塩基A、G、C、ならびに、Tおよび/またはUの頻度を基準に、同じ頻度となるランダムな配列を生成し、これを仮想アプタマー配列群Srとする。選択アプタマー配列群Sから計算される自由エネルギーをgとし、仮想アプタマー配列群Srから計算される自由エネルギーをgrとする。目的アプタマー配列の濃縮が起こったプールでは、前述のように、選択アプタマー配列群の自由エネルギーの分布が、仮想アプタマー配列群の自由エネルギーの分布よりも有意に低い。したがって、選択アプタマー配列群の自由エネルギーgの分布が、仮想アプタマー配列群の自由エネルギーgrよりも低い分布であれば、プールにおいて目的アプタマー配列が濃縮されているといえる。ここで、選択アプタマー配列群の自由エネルギーgの分布が、仮想アプタマー配列群の自由エネルギーgrの分布よりも低いことを示すため、下記式1に示す仮説をたてる。下記式において、μおよびμrは、それぞれ、選択アプタマー配列群の自由エネルギーgの平均値および仮想アプタマー配列群の自由エネルギーgrの平均値である。この仮説が棄却されれば、選択アプタマー配列群の自由エネルギーgの平均は、仮想アプタマー配列群の自由エネルギーgrの平均よりも、有意に低いことが示される。ここでは、平均値の差の検定として、Welchのt検定を行う。t値は、下記式2で表わされ、自由度vのt分布に従う。この際、nは、選択アプタマー配列群Sの選択アプタマー配列数、nrは、仮想アプタマー配列群の仮想アプタマー配列数とする。また、vは、分散の推定値を用いて、下記式3とする。2母集団との比較前記1と同様に、被対照プールに含まれる選択アプタマー配列群Sから計算される自由エネルギーをgとする。ここで、gが正規分布に従うと仮定すると、下記式4となる。つぎに、前記1と同様に、選択アプタマー配列群Sの塩基頻度を基準に、同じ頻度となるランダムな配列を生成し、これを仮想アプタマー配列群Srとする。そして、仮想アプタマー配列群Srから計算される自由エネルギーをgrとし、その平均値をμ0とする。ここで、下記式5に示す仮説をたてる。この時、t値は、下記式6となり、自由度n-1のt分布に従う。この仮説が棄却されれば、選択アプタマー配列群の自由エネルギーgの平均μは、仮想アプタマー配列群の自由エネルギーgrの平均μ0よりも、有意に低いことが示される。3分布の平行移動前記被対照プールに含まれる選択アプタマー配列群の自由エネルギーが、前記選択アプタマー配列群由来の仮想アプタマー配列群の自由エネルギーよりも有意に低い場合、その差が十分であるかは、例えば、以下の方法により、判断可能である。ここで、下記式7に示す仮説をたてる。下記式において、μは、選択アプタマー配列群Sから計算される自由エネルギーgの平均値であり、μ0は、仮想アプタマー配列群Srから計算される自由エネルギーgrの平均値であり、δは、どの程度の平均値に差があるかを示すパラメータである。この時、t値は、下記式8となり、自由度n-1のt分布に従う。この仮説が棄却されれば、選択アプタマー配列群の自由エネルギーgの平均μは、仮想アプタマー配列群の自由エネルギーgrの平均μ0よりも、十分に有意に低いことが示される。[実施形態6]前記対照プールの候補アプタマー配列群の自由エネルギーと、前記被対照プールの選択アプタマー配列群の自由エネルギーとを比較する場合、例えば、以下の方法によって、両者の自由エネルギーの差が有意であるかを判断できる。1ラウンド間の比較i番目のラウンドi-thのプールを対照プールとし、それに含まれる候補アプタマー配列群Siから計算される自由エネルギーの分布を、μiとする。他方、前記ith以降の任意のj番目のラウンドjthのプールを被対照プールとし、それに含まれる選択アプタマー配列群をSjとし、選択アプタマー配列群Sjから計算される自由エネルギーの分布を、μjとする。ただし、j>i。ここで、下記式9に示す仮説をたて、Welchのt検定を行う。検定の結果、有意であれば、j-thラウンドの自由エネルギーの分布μjの平均は、i-thラウンドの自由エネルギーの分布μiの平均よりも低い。つまり、目的アプタマー配列の濃縮が進んでいると判断できる。本発明のプログラムは、前述の本発明の予測装置の使用方法を、同装置に実行させるプログラムである。本発明のプログラムは、例えば、記録媒体に記録されていてもよい。前記記録媒体は、特に制限されず、例えば、HDD、CD-ROMCD-R、CD-RW、DVD、メモリーカード等があげられる。また、本発明のプログラムは、例えば、予め、前述の本発明の予測装置にインストールされていてもよいし、前記記録媒体またはインターネット等の回線網を介して、インストールされてもよい。なお、本発明のプログラムは、必ずしも、前述の本発明の予測装置にインストールされている必要はなく、例えば、サーバに格納されている本発明のプログラムにより、前述の使用方法を、前述の装置に実行させてもよい。つぎに、本発明の実施例について説明する。ただし、本発明は、下記の実施例により制限されない。[実施例1]本実施例では、複数のアプタマーを含む初期プールを用いて、SELEX法により、ヒトHMGB1Highmobilitygroupbox1;商品名H4652HMG-1、シグマ社製に対して結合可能な目的プライマーの濃縮を行った。前記ラウンド0の初期プールは、複数のアプタマーを含むプールを使用した。前記初期プールにおけるアプタマーは、3’末端領域および5’末端領域にプライマー領域を有し、前記3’末端領域および5’末端領域の間にランダム領域を有し、各アプタマーにおいて、前記プライマー領域の配列は、完全に一致し、前記ランダム領域の配列は、異なる配列とした。SELEX法では、ターゲットに結合したアプタマーの選択および結合したアプタマーの増幅を1セットとして、前記初期プールR0から8回繰り返し行い、各ラウンドR1~R8でプールを得た。また、各ラウンドのプールについて、BiacoreX商品名;GEHealthcareUKLtd.を用いた表面プラズモン共鳴分子間相互作用解析により、HMGB1との結合能Rmax:単位RUを測定した。各ラウンドのプールの情報を下記表1に示す。下記表1において、個数は、各ラウンドのプールに含まれる全アプタマーの個数配列数であり、塩基組成は、各ラウンドのプールについて、各プールに含まれる全アプタマーのランダム領域における塩基組成各塩基の出現頻度であり、Rmaxは、各ラウンドのプールとHMGB1との結合能を示し、MFE平均およびMFE分散は、各ラウンドのプールの最小自由エネルギーMFEの平均kcal/molおよび分散kcal/mol2を示す。また、後述する自由エネルギーによる解析の信頼性を確認するため、前記各ラウンドのプールについて、実際に、表面プラズモン共鳴分子間相互作用解析を行い、いずれのラウンドでHMGB1に結合する目的アプタマーが濃縮されているかを確認した。具体的には、各ラウンドのプールについて、重複配列プール内に有意に複数個存在するアプタマーの配列の有無と、各ラウンドのプールのHMGB1への結合を確認した。また、表面プラズモン共鳴分子間相互作用解析により、R8のプールから、HMGB1と実際に結合するアプタマーHMGB1結合アプタマーを選択し、前記HMGB1結合アプタマーの塩基配列決定を行った。さらに、前記HMGB1結合アプタマーの配列に基づいて、各プールに含まれるアプタマーについて、モチーフ解析を行った。その結果、R1~R3のプールにおいては、有意に存在する重複配列は見られなかった。また、R4のプールでは、前記HMGB1結合アプタマー配列のモチーフを有する重複配列が検出されたが、HMGB1には結合しなかったことから、前記HMGB1結合アプタマー配列の濃縮は、起こっていなかった。R5のプールでは、実際に、複数種類の重複配列が検出され、前記プールにおける頻度が低い重複配列に、HMGB1に結合するアプタマーが確認された。このことから、R5のプールは、前記HMGB1結合アプタマーの濃縮が起こっていることがわかった。また、R6およびR7のプールにおいても、複数種類の重複配列が確認され、前記プールにおける頻度が上位である重複配列最上位の重複配列を除くに、HMGB1に結合するアプタマーが確認された。このことから、R6および7のプールは、R5のプールよりも、さらに前記HMGB1結合アプタマーの濃縮が起こっていることがわかった。また、R8のプールにおいても、複数種類の重複配列が確認され、前記プールにおける頻度が最上位を含む上位である複数種類の重複配列が、HMGB1に結合するアプタマーであることが確認された。このように、最上位を含む上位の重複配列がHMGB1に結合したことから、R8のプールは、R6およびR7のプールよりも一層、前記HMGB1結合アプタマーの濃縮が起こっていることがわかった。以上の結果から、後述する自由エネルギーに基づく解析において、R1~R4の検定が有意となった場合は、偽陽性となり、R5~R8の検定が有意とならなかった場合は、偽陰性となる。1-1仮想アプタマー配列の自由エネルギー前記表1に示す各ラウンドの塩基組成に基づいて、同じ組成となる仮想アプタマー配列を、コンピュータ上で、1000個、1万個、10万個生成した。そして、各仮想アプタマー配列から自由エネルギーを求め、各ラウンドのプールについて、自由エネルギーの平均値および分散値を算出した。この仮想アプタマー配列の生成、ならびに、自由エネルギーの平均および分散の算出を、合計10回n=10行った。図8に、10回の試行における、各ラウンドの自由エネルギーの平均を示す。図8Aは、ラウンドあたり1万個の仮想アプタマー配列を生成した結果であり、図8Bは、ラウンドあたり1000個の仮想アプタマー配列を生成した結果である。図8において、縦軸が、平均自由エネルギー単位kcal/molを示し、横軸が、試行回n=1~10を示す。同図に示すように、ラウンドあたり1000個または1万個の仮想アプタマー配列を生成した際、前記自由エネルギーの平均は、10回の試行によっても、ほぼ同等の値を示した。中でも、1万個の仮想アプタマー配列の生成より、極めて安定に自由エネルギーの平均を算出できることがわかった。図9に、10回の試行における、各ラウンドの自由エネルギーの分散を示す。図9は、ラウンドあたり10万個の仮想アプタマー配列を生成した結果である。図9において、縦軸が、自由エネルギーの分散単位kcal/mol2を示し、横軸が、試行回n=1~10を示す。同図に示すように、ラウンドあたり10万個の仮想アプタマー配列を生成した際、前記仮想アプタマー配列の自由エネルギーは、10回の試行によっても、ほぼ同等の自由エネルギーを示した。1-2サンプリングアプタマー配列と仮想アプタマー配列との比較前記表1に示す各ラウンドのプールから、所定数のアプタマー配列をサンプリングした。前記サンプリングしたアプタマー配列数は、20、40、60、80、100、120、140、160、180、200個とした。前記各ラウンドからサンプリングした所定数の全サンプリングアプタマー配列のランダム領域の塩基組成を解析し、各ラウンドの塩基組成に基づいて、同じ組成となる仮想アプタマー配列を、コンピュータ上で、10000個生成した。続いて、前記各サンプリングアプタマー配列および前記各仮想アプタマー配列から、RNAfoldW70W74Jul2008を用いて自由エネルギーを求め、各ラウンドのプールについて、前記サンプリングアプタマー配列の自由エネルギーの平均、および、前記仮想アプタマー配列の自由エネルギーの平均を算出した。自由エネルギーの算出の条件は、温度摂氏37°C、最小ステム長2塩基対とした以下、同様。この仮想アプタマー配列の生成および自由エネルギーの平均の算出を、合計100回n=100試行した。そして、100回の試行のうち、前記サンプリングアプタマー配列の平均自由エネルギーSRが、前記仮想アプタマー配列の平均自由エネルギーSVよりも低かった回数をカウントした。下記表2に、これらの結果を示す。下記表2は、各ラウンドおよび各サンプリング数において、前記サンプリングアプタマー配列の平均自由エネルギーSRが、前記仮想アプタマー配列の平均自由エネルギーSVよりも低かった回数を示す。また、下記表2に、各ラウンドのプールの結合能Rmaxを合わせて示す。自由エネルギーの解析の結果、前記表2に示すように、サンプリング数にかかわらず、SR<SVの回数が、R5で増加して100回に近づき、R6以降で100回となった。これらの結果は、実際にHMGB1との結合を確認した前述の表面プラズモン共鳴分子間相互作用解析の結果に対応していた。このことから、自由エネルギーの解析によって、標的物質と結合可能なアプタマー配列の濃縮の有無を、優れた感度で判断可能であるといえる。さらに、前述の100回の試行のうち、t検定により、有意水準5で、前記サンプリングアプタマー配列の平均自由エネルギーSRが、前記仮想アプタマー配列の平均自由エネルギーSVよりも低かった回数をカウントした。下記表3に、これらの結果を示す。下記表3は、各ラウンドおよび各サンプリング数において、有意水準5で、前記サンプリングアプタマー配列の平均自由エネルギーSRが、前記仮想アプタマー配列の平均自由エネルギーSVよりも低かった回数を示す。また、下記表3において、各ラウンドのプールの結合能Rmaxを合わせて示す。この結果、前記表3に示すように、サンプリング数にかかわらず、SR<SV有意水準5の回数は、R4までほぼ0となった。そして、R5で増加が見られ、R6以降で、ほぼ100回となり、サンプリング数40回以上の場合、完全に100回の結果が得られた。これらの結果は、実際にHMGB1との結合を確認した前述の表面プラズモン共鳴分子間相互作用解析の結果に対応していた。特に、t検定を併用することによって、R4の結果もほぼ0回となったため、R4に関する、偽陽性の判断を回避でき、優れた信頼性で、濃縮の有無を判断可能であるといえる。また、サンプリング数を40以上とすることで、結果が安定することから、偽陰性の判断を回避でき、さらに信頼性を向上できるといえる。さらに、サンプリング数20、40、80、180の各ラウンドのプールについて、自由エネルギーの平均の標準偏差を、下記表4に示す。前記表4に示すように、サンプリング数の増加に伴い、標準偏差値の減少が見られた。この結果から、例えば、セレクションのどの段階においても、サンプリング数を増加させることによって、より信頼度の高い判定が可能るといえる。1-3サンプリングアプタマー配列のラウンド間の比較およびサンプリングアプタマー配列と仮想アプタマー配列との比較各ラウンドのプールについて、前記プールに含まれる全アプタマー配列をサンプリングアプタマー配列とした以外は、前記1-2と同様にして、10000個の仮想アプタマー配列の生成、前記サンプリングアプタマー配列および前記仮想アプタマー配列の自由エネルギーの平均の算出を行った。これらの結果を、図10の箱ひげ図に示す。図10において、縦軸は、自由エネルギー単位kcal/molを示し、横軸が、各ラウンドR1~8を示す。図10において、左半分が、サンプリングアプタマー配列の結果であり、右半分が、仮想アプタマー配列の結果である。まず、各ラウンドのサンプリングアプタマー配列の自由エネルギーの中央値を比較した。この結果、R1~R5と比較して、R6以降において、著しい自由エネルギーの低下が見られた。R6は、前述のように、HMGB1結合アプタマーの濃縮が顕著に確認されたラウンドである。このことから、自由エネルギーの低下を確認することで、目的アプタマーの濃縮の有無を確認できることがわかった。つぎに、各ラウンドについて、サンプリングアプタマー配列の自由エネルギーの中央値と、仮想アプタマー配列の自由エネルギーの中央値とを比較した。この結果、R5において、サンプリングアプタマー配列の中央値と比較して仮想アプタマー配列の中央値の低下が確認され、R6以降において、サンプリングアプタマー配列の中央値と比較して、より大きな仮想アプタマー配列の中央値の低下が確認された。このことから、サンプリングアプタマー配列の自由エネルギーと、仮想アプタマー配列の自由エネルギーとを比較することで、目的アプタマーの濃縮の有無を確認できることがわかった。以上、実施形態を参照して本願発明を説明したが、本願発明は、上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解しうる様々な変更をすることができる。この出願は、2010年7月2日に出願された日本出願特願2010-152020を基礎とする優先権を主張し、その開示の全てをここに取り込む。