Answers to: trans-ABySSのエラー（おそらくBlatのところ）

Comment by かどた on かどた's 質問

かどた — Mon, 28 Feb 2011 15:49:22 +0900

おかげさまで無事実行できました。 pysamのre-installでうまく動きました。まだstage1でkの値が二つだけで試しただけですが、とりあえずはエラーを吐かずにいけたことがなによりです。 pysam問題ですが、おそらく、他のエラーが出ていた際に、pythonのre-installをしたので、pysamの以前のインストールが無効になってしまっていたのでしょう。

Comment by かどた on かどた's 質問

かどた — Tue, 22 Feb 2011 19:52:31 +0900

そうなんですよ。その後「Error in stage1, perhaps merge.pl」上で開発者に上記の-fastMapのところを外せばいい、というアドバイスを受けてやったところ、blatのところでは（相当時間がかかりましたが）エラーがでなかったようですが、そのあとのところでエラーがでました。。。「ImportError: No module named pysam」みたいなエラーが出て。えー、pysamいれてるのにー、なんでよ、みたいな。。。今は意地でtrans-ABySSを動かそうとしていますが、別のプログラムに乗り換えたほうがいいのか、、、難しいところです。どなたか私このプログラム使って無事動かせましたよ、っていうのがありましたら、コメントよろしくお願いしますm(_ _)m

Comment by mya_ on かどた's 質問

mya_ — Tue, 22 Feb 2011 15:10:06 +0900

fastMapの言及しているのが、configs/projects.cfg の上２行だけなんで

 
[default]
blat_1_2: FA2 FA1 OUT_DIR/KOUT1-KOUT2.psl -minIdentity=100 -maxGap=0 -fastMap
blat_2_1: FA1 FA2 OUT_DIR/KOUT2-KOUT1.psl -minIdentity=100 -maxGap=0 -fastMap

のオプション削ればfastMap無しで受け付けそうなんですが、だいたいこういうパイプラインってデフォルト設定が前提のコーディングしていてmergerがおかしくなるとか出そうなんですよね。

Comment by かどた on かどた's 質問

かどた — Thu, 17 Feb 2011 15:06:07 +0900

Blat with -fastMap option is used as an alignment tool in trans-ABySS とtrans-ABySSの開発者が書いてあるので、このオプションをつけなければ25,000bpまで大丈夫だが、このオプションを使ったときに5000bpという長さ制限を加えることで早くしているとか、、、そういうこともあるかもしれませんね。今日はまだ自分のデータを眺められてはいませんが、もしコンティグ長が全て25,000bp以下ならtrans-ABySSのプログラムをいじって、-fastMapオプションをつけないようにすることも一つの手ですね。ありがとうございます！ひょっとして、このようなプログラムを使いこなすための会合がBiohackathon(つづりあやしい)だったりして。。。

Answer by mn3

mn3 — Thu, 17 Feb 2011 12:09:08 +0900

blat は 5000bp より長いクエリを受け付けています。たとえば、http://genome.ucsc.edu/cgi-bin/hgBlat?command=startでは、

Only DNA sequences of 25,000 or fewer bases and protein or translated sequence of 10000 or fewer letters will be processed. Up to 25 sequences can be submitted at the same time. The total limit for multiple sequence submissions is 50,000 bases or 25,000 letters.

といのことで、すくなくとも 25,000 塩基までは受け付けています。ですので、クエリの長さ制限以外のところにエラーがあるようですね。BLAT Program Specifications にもクエリ長さ制限については記載がありませんでした。

Comment by かどた on かどた's 質問

かどた — Thu, 17 Feb 2011 09:19:05 +0900

ひょっとしたらBlatでqueryできる最大の配列長が5000bpで自分のデータは1106370番目のコンティグ長が6471bpでリミットを超えちゃっているからではないかと思い始めました。ちょっと調べてみます。ちなみに私はblatをこれまで一度も使ったことがないのですが、blatってちょっと短めの配列用？！だから5000bp程度までしかできないような設定になってるっていうようなことがあるんですかね。。。ほとんど独り言ですが。あ、mn3さま編集どうもありがとうございますm(_ _)m

Comment by かどた on かどた's 質問

かどた — Wed, 16 Feb 2011 18:58:37 +0900

さらにやらかしてますねm(_ _)m あのURLだけだとどのスレッド（という表現でいいのか？！）かわかりません。「Error in stage1, perhaps merge.pl」というタイトルのやりとりです。

Comment by かどた on かどた's 質問

かどた — Wed, 16 Feb 2011 18:23:04 +0900

あ、すいません。 GoogleディスカッショングループのURLをつけ忘れました。 http://groups.google.com/group/trans-abyss です。