|
今ある菌のde novo assemblを行っています。 先日454で読んだ配列をABySSでアセンブルして、コンタミのDNAを除きました。 その後コンソーシアムのある454で読んだESTやSangerで読んだDNA配列と一緒にアセンブルしたいのですがどうしたらよいですか？長さは40bpのもあれば、長くて100kのもあります。平均700bpぐらいだと思います。 試しにSOAP denovoを使ったところ Floating point exceptionとかsegmentation faultのエラーで止まってしまいます。 ABySS使っても Colour-space assembly Assertion failed: (isdigit(seq[0])), function loadSequences, file AssemblyAlgorithms.cpp, line 79. Abort trap で止まってしまいます。 違う原因で止まってしまっているのかも知れませんが、そこも踏まえてご助言いただけると助かります。 |
|
454の配列があると言う事は、お近くに454純正アセンブラのnewbler(gsAssembler)が動く環境があるかと思います。 newblerは一般的なfasta形式の配列を扱えますので、sffと混ぜてアセンブルされるのが ソフトウェア導入でご苦労されるより手っ取り早いのではないでしょうか。 |
確証はありませんが、SOAP denovoの"segmentation fault"はコンパイルが失敗していたり、バイナリが環境に合っていない可能性もあるような気がするので、実行環境と、バージョンも書いた方が良いかもしれません。 ABySSの方は入力ファイル依存の問題のような気がするので、ここに入力ファイルの一部でもさらすのが無理なら、ご自分でエラー原因を特定するしかないような気がします。