Trim Galore! (trim_galore)を使うためにはfastqcとcutadaptのインストールが必要となる。cutadaptはpipでインストールできるのだが、そのpipはデフォルトでは入っていない。 pipのインストールには、いくつかのやり方があるが、ここではhomebrewでpython3を入れて、その際に入るpip3を利用する方法をとる。ここでpython3,pip3と末尾に3が付いているのは、pythonのバージョン3のそれ、というこ …
今月(2016年12月)は、AJACSadvanced講習会や共同研究関係で、次世代シークエンサーから得られたデータを自ら直接解析することに関わる事が多 …
ワイルドカードによるファイル指定でそのコマンドに一括で引数で指定できず、数多くのファイルに対して同じ処 …
卒論、修論、追い込みの方も多い今日このごろ、バックアップは重要。大事なデータほどきっちりとバックアップを。もち …
カイコ(Bombyx mori)のEnolaseに関する論文がacceptになりました。研究成果としてすでに公開されている配列情報データベースから候補をドライ解析(in silico)で探 …
Oracle Labs PGX(Parallel Graph Analytics) というグラフに対するクエリや分析のためのフレームワークのユーザー勉強会にお誘いいただいた …
先日のAJACSadvanced講習会で、Trinityで得た結果をアミノ酸配列に翻訳するプログラムをここで言及していないのに気がついたので。Trinityで得られるのは塩基配列レベルの情報で、それをアミノ酸配列に翻訳してくれるのがこのtransdecoder。
brew install -v transdecoder
でインストールされる。
#!/bin/sh
TransDecoder.LongOrfs -t Trinity.fasta …今回の遺伝研研究会では、初の試みとして講習会を連続開催してみた。つまり、研究会をやる前(1日目午前)と後(2日目午後)に講習会を実 …
遺伝研研究会として「次世代モデル生物におけるゲノム情報利用ワークショップ」を2016年12月15,16日にオーガナイズさせていただいた。「ゲノム情報利用ワークショ …
2日目。引き続き、SRAからのRNA-seqデータ回収のパイプライン作り検討。それ以外に、統合TVがflashで再生できなくなっていた問題の対処をお願いするなど。
午後には山口情報芸術センター(YCAM)を見学させて頂く。
今日2016年12月12日から、国内版Biohackathon(BH16.12)、山口県湯田温泉にて。「公共遺伝子発現データ目次AOEを使ったデータ解析事例作成」というテーマで取り組む予定。遺伝子発現デ …
今回が9回目のMishima.syk(ミシマ・エス・ワイ・ケー、三島創薬勉強会)。今回は仮想環境の上のLinuxでTensorFlow（てんさーふろー、と読むらしい)と …
第39回日本分子生物学会年会にて、「いかにして使えるデータベースを維持し続けるか?」というフォーラムをオーガナイズした。本来ワークショップ …
要するに、ArrayExpressのデータをコピった。やった事自体、簡単そうだが、転送量が数十TByteもあると、そう簡単にいかなかった。2016年5月に英国Cambridge郊外のHinxtonにあるEuropean Bioinformatics Institute(EBI)に行って打ち合わせして …
地震で飛び起きたSPARQLthon50の2日目。結局、ArrayExpressのmirrorは終わらず。来月のSPARQLthonへと持ち越し。それ以外のデータ収集手段についてupdateしていただき、より簡単に高速にデータを取得できる …
久しぶりの参加。ArrayExpressのmirrorに関しても目処がついて、次期開発に向けてのちょうど良いタイミング。RDFやSPARQLの利用実例として実用的なものに発展させていければいいな。 50回目のキリ番SPARQLthonのためか、参加者が多く、annexの柏での作業スペースでの参加。仲間と議論するのは楽しいが …
AJACSadvanced三島３でPitagora GalaxyをNGSデータ解析環境として紹介したいと考え、Pitagora Galaxy Meetup 2016-11に参加。仮想環境のダウンロードページにある仮想環境イメージのアップデートのお願いがメ …
普通のftpやrsyncでコピーするとしたら、数ヶ月単位かかるような数十Tbyteのデータを転送する必要があり、先方にAsperaを設定してもらい。ascpというコマンドラインのプログラムをネット上からググって探してダウンロードインストールして、以下のようなコマンドで実行。 [shell] ascp -p -k 1 -QT -l 1G --overwrite=diff+older hoge@xxxx:hoged …
今年(2016年)の３月に急逝したLorenz Poellinger氏の追悼シンポジウム(Professor Lorenz Poellinger Memorial Symposium)で縁の深かった方々によるプレゼンテーションは、泣きそうになる一方、彼らの研究の遍歴がうかがい知れて勉強になっ …
去年は三島でお世話させていただいた、がんとハイポキシア研究会。今年は岐阜で。実は、2年続けての東海地方開催。今年は運営に関わる部分が …
今年度(2016年度)採択していただいた国立遺伝学研究所研究会「次世代モデル生物におけるゲノム情報利用ワークショップ …
cell types in mammals' wordpress_id: 2863 categories:
READとFANTOMDB論文以来、なんと15年ぶりにNucleic Acids Research (NAR) Database issueに名前が入った標記の論文が掲載。
Nucl. Acids Res. (2016) doi: 10.1093/nar/gkw995 First published online: October 27, 2016
Advance accessなのでページ数は付いてないが、DOIのcitationは変わらないので、こちらで是非。 「データ …
去年は参加できなかったJASPUGMだが、今年は参加。今回はとくに発表とかなく。実質アカデミックフリーになって、公式にユーザ …
今週はオープンアクセスウィーク。ということで、連動企画として、2016年のSPARC Japanセミナー2回目はこのタイミングで。
今回もセミナー担当のワーキングメンバーと …
札幌の定山渓にて、10/22-23の一泊二日で。研究会のお題目の通り、オモロイ生き物でさまざまな生命現象を解明しようと取り組んでい …
MacOSXだとSpotlight検索が中身まで検索対象になっていて便利だが、ファイル名でのみ検索したいときがある。しかもコマンドラインで。そういう時に使っているのがlocateコマンド。 [shell] locate hoge [/shell] でhogeを …
今年はトーゴーの日シンポジウム(10/5-6)と1日被り、大会2日目からの参加。企業展示ならびにポスター会場がスカスカでびっくりする。別の …
今年も、10月5日ということで、「トーゴーの日」シンポジウム。今では、この1年の進捗を発表する機会となっているが、統合データベースプロジェクトが始まって早10年。論文で測ること …
Entomology（昆虫学）の国際学会に参加。初参加でアウェーの学会だったが、来てみると全く問題なかった。自分は、いわば系統樹を …
日本生物工学会誌に「バイオインフォマティクスを使い尽くす秘訣教えます!」という連載が開始。2ヶ月に1回、全7回の予定。2016年1月に鹿児島大学医学部で行われた統合データ …
ArrayExpressのdataをrsync開始してから今日でちょうど一ヶ月。前回と変わらず、experimentディレクトリのミラーが終わらない。データ転送量的には、合計1.5Tbyteほど。もう1thread増やしてやっているものの。約44Tbyteと書かれている全体のデータを信じれば、あと2年半。対策を考えねば。
いろいろ書き散らす。2016年の計として元旦に誓ったことでもあるし。まさかこういう形で転ぶとは思っていなかったが。この …
昨日のセミナーで知った言葉だが、どうも気にかかると思ったら、やはり自分たちが果たしてきた役割をいう言葉だ …
前にも書いたとおり、今年度(2016年度)から国際学術情報流通基盤整備事業(SPARC Japan)のセミナー企画ワーキンググループのメン …
ArrayExpressのdataをrsync開始してから3週間。だが、experimentディレクトリのミラーが終わらない。データ転送量的には、合計0.9Tbyteほどで、約44Tbyteと書かれている全体のデータからするとやはりまだまだ。thread数を増やしてみましたが、果たして。
今年度、申請して遺伝研研究会としての開催が決まっている「次世代モデル生物におけるゲノム情報利用ワークショ …
single end readsのTrinityの実行がこけて、hoge.fq.readcountファイルに以下のような記述が。
Error, cannot convert fastq file to fasta since cannot recognize read orientation as /1 or /2 (instead: F)
このエラーメッセージでググると同じ目にあった人のQAが …
かつて関わった厚生労働科学研究費補助金(2009年度〜2011年度)による高齢者のがん研究のデータがNBDCヒトデータベースより公開されたと連絡をいただく。当初は生 …
DRY解析教本ですが、2015年の10月に出版されて早10ヶ月。 その間、「細胞工学」の発刊が最終号を迎えるなどいろいろありましたが、初刷の在庫が無くなってきたとの …
トリミングしたうえでfastqcをかけてくれるソフトがこのtrim_galore。CutadaptとFastQCがすでにインストールされていることが必須(以下の使い方だと)。 singleの場合
[shell] trim_galore --fastqc --trim1 hoge.fq [/shell]
ペアエンドの場合
[shell] trim_galore --fastqc --trim1 --paired hoge_1.fq hoge_2.fq [/shell]
ArrayExpressのdataをrsync開始してから1週間。atlasディレクトリは終わったものの、残りの2つ、arrayとexperimentは終わらない。データ転送量的には、合計0.4Tbyteほどで、約44Tbyteと書かれている全体のデータからすると1/100ほど。ただ、全部ミラーしたらこの容量になるのかどうかも今のとこ …
初めての山の日(8/11)、急逝した研究者仲間のtwittererの初盆に名古屋へ。同じ分野の研究者というわけでもなく、オフラインには学会オフ会に …
ストレージの容量は多くの場合ファイルサイズで問題になることが多いが、たくさんファイルを作るようなプログ …
2016年8月12日日本時間午後2時頃、ArrayExpressのFTPにあるデータをミラー開始。ポイントはFTPではなく、rsyncでやっているというところ。先月のSPARQLthonで調べてわかったことを元に。プロセスがコケていないか、ディスク容量は大丈夫かをちょくちょく見張りつ …
[amazon template=thumbnail&asin=4797382309]
新米探偵本、ようやく通読。出てきたことが一度だけでなく、俵太のまとめと天羽社長の統計学指南で波状 …
先日の沼津でのワークショップで伺った、嬉しいニュース。
アカデミアの方々などに無料でご利用頂けるようになりました。 TIBCO started free of charge Spotfire License. http://spotfire.tibco.com/better-world-donation-program/
本当かと思 …
bam形式なリファレンスゲノムへのアラインメントの情報をwig形式のcoverageファイルへ変換。 実は昔もやったことがあり、ググって出てくるawkのワインライナーでなんとか凌いでいたが、今ネット検索するとbam2wigというコマンド化されたPerlスクリプトをgithubに発見。今回はこ …
東京白金台の北里大学薬学部にて、「ゲノムと遺伝子発現DBの現状ートランスオミクスのための基礎知識」というタイトルで、30分ほど話した。ワークショップの前座。
そのうち統合 …
昨日の続き。直で職場のサーバーにsshできないのはやはり不便だ。 やはり、1日経ってもファイルリストの取得、終わっていない…。定期的にrsyncするように早くしたい。が、昼前に終わったので、実体の取得を開始 …
今回はいつもと違い、DBCLSではなく新潟大学での開催ということで、RDFやSPARQLの入門的な講習会がほぼ1日開催された。こういうスタイルは新しいが、こういう形の地方行脚重要。 こちらはいつものハッカソンで …