Error: qual length () differs from seq length () for fastq record !

CAGE-seqのマッピングをするときに、
Cutadapt, rRNAdust, removeNをかませるとたまに長さがあわない配列や、クオリティスコアのない配列が混ざってきて、結果タイトルのエラーに出くわす。

しかもなったりならなかったりして大変困るので、長さがあっている配列（とスコア）だけ抜き出すワンライナーとスクリプトを書いたのでおいておきます。





check_read.pyを保存して、ターミナルにcheck_qs.shをぺたっとしてくだし。 @SRRから始まる4行を抽出して、長さが違うものはpythonスクリプトで捨てます。 SRRXXXXXX_extracted.fastqが訂正されたfastqファイルになる予定です。

コメント

コメントを投稿

このブログの人気の投稿

旧姓を失った研究者が海外でも活躍するためには

イメージ
突然ですが結婚しました！ そして苗字の珍しさで負けたので、夫の姓にすることに話し合いで決まりました。 今まで名前と生年月日という情報をもとに個人が識別されていたのに、名前が変わるって大丈夫なの！？と思ったらやっぱり大変でした。 ということで、それに伴う名前変更作業などについて記録しておきます。 （これを参考にして何か問題が起きても責任はとれませんので悪しからず） 　 　 そもそも仕事での名前をどうするか？ 研究者にとって、今までの名前で検索して業績がひっかからないというのは致命的です。 私は一応生物系（？？？）の研究者なので基本的にPubmedを参照することが多いですが、 知り合いの研究者の名前をPubmedで検索してみると、論文誌によって保持される著者名情報が違うためミドルネームの情報ではヒットしない場合が多く、苗字の変更はクリティカルになります。 ↑ 姓名などについて詳しく条件づける検索方法は特にないようです 例えばHetta HaraさんがHetta Okaneさんになったとき こんなとき、名前をどうするかについては知り合いの（主に女性の）研究者の人たちをみるといくつかのパターンがありました。 旧姓のまま（Hetta Hara） メリット　当然楽 であり、今までの活動をまるっとそのまま引き継げる。 デメリット　戸籍名を強制する表記があったときに、それらを統合する術がない。また、その本人であることを証明することも難しくなる。税金関係や未確認だが卒業論文、証明などは戸籍名になるという説がある。パスポートも原則戸籍名です。 新姓に変更する（Hetta Okane） メリット　戸籍名と一致するため本人確認で身分証を必要するときにも問題なし デメリット　もしそれまでに研究活動をしていたらそれらを統合する術がない。医師免許や博士号などの名前と異なるので、証明が難しい。 ミドルネームに旧姓を入れる（Hetta Hara Okane） メリット　一応公式の名前は戸籍名と一致しており、旧姓の情報もわかるので親切な人であれば昔の論文も検索してくれるかもしれない。 デメリット　Pubmedの検索では名前とミドルネーム情報はイニシャルでしか含まれないことが多い（Pubmedの仕様なのか元の...
続きを読む

再現性なんてないさ（？）

イメージ
※この投稿は 「今年読んだ一番好きな論文2016」 Advent calendar の12月21日の記事です。 （わざわざ手描きにしたのに結局著作権が危うい絵を描いている） どうもこんにちは、 博士論文（仮）を提出したばかりのぴちぴちのD3、@calicolicali です。学生最後ということで、院生しか登録できないAdvent calendarに登録させていただきました！ ʕ◔ϖ◔ʔ 　「統計について」という事前コメントをしましたが、 私の専門はアルゴリズム寄りの Computational Biology で、これまではRNA二次構造予測ソフトウェア ParasoR の開発や、配列アラインメントソフトウェア LAST の開発に関わらせていただいたりしてきました。 （ご興味あればぜひ！） しかしこの分野では「網羅的な遺伝子データや配列情報を扱うためには統計や機械学習が必要！」ということで、日々これらの分野の勉強をしております。 　そんな私が紹介する、今年読んだ一番好きな論文はこちら(｀・ω・´)つ "Estimating the reproducibility of phychological science" Open Science Collaboration. Science 349.6251 (2015): aac4716. 去年の８月にサイエンス誌に投稿された論文ながら、Google scholarではすでに652本のサイテーションがなされています。 ざっくり解説  　この論文は心理学の分野の3つのジャーナルに2008年に出た 約100の実験的な相関解析に対して再実験を行い 、どのような実験・結果であれば再現性が得られるのかということを調べた論文です。元論文の97%で有意なp値が得られていたものに対して、 再現実験では有意とされたものは37% 、95%信頼区間に入ったものが47%、元論文のデータと組み合わせても 最終的に68%までしか有意性を確かめられなかった 、という報告がなされました。 　この論文は発表後様々な研究者によって正当性を検証するコメントや論文がでており、そして、最近Natureなど大手の雑誌でも心理学に加えて生物学など様々な分野での Reprod...
続きを読む

WLSのssh

WLS上でのsshが以下のように突然拒否される病にかかった。 ssh: connect to host xxx port 22: Resource temporarily unavailable おそらくWindows updateでFirewall設定が変更になるかしたからだと思われるが、 ~/.ssh/configに Host *              ProxyCommand nc %h %p %r と入れて、 chmod 600 ~/.ssh/configにしておいたところうまくいった。 回答をよく読むとFirewall設定ではなく、内部でのProxyを利用したjumpがうまくいっていない模様？ サーバーの方でRAMを追加するかなにかしたらしいので、内部でどこかに飛ばされる仕様になったのだろうか。 わたしはふいんきでsshをしている。 参考： https://superuser.com/questions/1098526/windows-10-linux-subsystem-ssh-client-resource-temporarily-unavailable
続きを読む