投稿

初めて年末ジャンボ宝くじあたった

 気づいたときにたまに買っている宝くじが今回はじめてあたった。といっても下一桁があたった7等で、300円なんだけど。 宝くじが割に合わないのは知っているけど、買う楽しみがあるのも確か。だからいつも一枚だけ買っている。当たらなくても慈善事業に使われるのなら惜しくはない。 それが今年は初めてあたった。7等は下一桁だから1/10で当たるんだけど、それでもなかなか当たらないんだよね。 今回は初めてあたったということで、今年の運は良かったか来年の運が良くなるのかよくわからないけど、バンザイ。

ログイン・パスワード地獄

 ネットでなにかすると、すぐにログインしてくださいってなる。それが本当にストレス。なんで店の予約するときに生年月日とか住所とか登録しないといけなんだよ。一方電話予約とか直接予約だと名前だけでいいとかあるし。色々最近腹立つことが多いな。

ニトリに珪藻土バスマットを交換しに行く

 ニトリの珪藻土バスマットがうちにあったので交換しに行った。ネット上のフローで回収商品とわかり、割れて使っていなかった物があった。割れているのでちょっと怖いけど。自転車で行くので結構疲れた。 専用のカウンターがあり、5−10人ほど待機していた。店員も大変だな。

アンリのアップルパイ

イメージ
近くのパン屋のアップルパイ  

仕事納め・・・ならず

 12月25日が仕事納め。朝から出張で、超特急の仕事をした。そこから夜八時に終わろうかと思ったが、納得いかず結局次の日の土曜日も休日出勤する羽目に。さらに、日曜日も家でできる仕事をすることに。 仕事と趣味が同じような人ならいいけど、少なくとも今やってる作業はそういうわけでもないからなー。

fastqを横一列に並べる

Before @SEQ_ID GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT + !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 After @SEQ_ID        GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTT        +        !''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65 上だとよくわからないけど、改行がタブになってます。 コマンド for i in 1 2; do for j in 1 2 3 4; do zcat FASTQ_S1_L001_R${i}_001.fastq.gz | awk '{if(NR%4==$j){print}}' > R${i}-${i}; done; done paste R1-1 R1-2 R1-3 R1-4 > R1 paste R2-1 R2-2 R2-3 R2-4 > R2 これが賢いやり方だとは思えないけど・・・

fastqファイルからTotalseq抗体のFeature barcodeの配列だけを抜き出す

イメージ
zcat FASTQ_S1_L001_R2_001.fastq.gz | awk '{if(NR%4==2){print}}' | cut -c 11-25