|
とある生物のゲノムを解読しました。 遺伝子を同定するためにRNA-seqを行いましたが、発現が低い遺伝子や特異的に発現する遺伝子は見逃していると思います。 そこで、以下のサイトの通りに、RNA-seqで得られたイントロン情報をAUGUSTUSに読み込ませて、RNA-seqでは同定が難しい遺伝子を予測しようとしました。 http://bioinf.uni-greifswald.de/bioinf/wiki/pmwiki.php?n=IncorporatingRNAseq.Tophat 5.のイントロン情報の作成まで行い、6.パラメーターの設定は入っていたファイルをコピーし、実際AUGUSTUSを実行してみると、途中まではきちんと予測してくれるのですが、ある程度進むと「Segmentation fault」で中断されてしまいます。 上記の方法は（１）のように実行していますが、諦めてイントロン情報を入れない（２）の方法ならば最後まで実行されるので、入れたイントロン情報か設定が悪いのかと思っていますが、何が原因なのか分かりません。アドバイス等いただけましたら幸甚です。 （１）augustus --species=...... --extrinsicCfgFile=extrinsic.cfg --alternatives-from-evidence=true --hintsfile=hints.gff --allow_hinted_splicesites=atac --introns=on --genemodel=complete genome.fa > aug1.out （２）augustus --species=....... --genemodel=complete genome.fa > aug1.out また、CufflinksやAUGUSTUS以外に遺伝子予測するのに、オススメのプログラム等ありましたら教えて下さい。 よろしくお願いします。 |
