山内のサイト

岸本 DNA-Reseq(Breseq) 追加データ 手順をメモっておく 2019-10-26

受け取ったデータをコピー

~/src/KishimotoRNA4/下

gzip strain1_1.fastq.gz
gzip strain1_2.fastq.gz
gzip strain2_2.fastq.gz
gzip strain2_1.fastq.gz

ungzipの後、md5チェック

fastqcによる品質チェック

~/src/KishimotoRNA4/下

#!/usr/bin/bash
for f in *.fastq
do
  fastqc --nogroup -o qc $f &
done

トリミング

Trimmomaticでトリミングを行う。「ここ」と同じ処理で、Pair EndでPEを指定する。入力ファイルは2つ(forwardとbackward)で、出力ファイルは4つ(pairedとunpaired × forwardとbackword)。またリード長が150なのでMINLENは120にセット。入力はgzipしてから食わすことにする。

java -jar -Xmx256g /usr/local/Trimmomatic/trimmomatic-0.38.jar PE\
 -threads 16 -phred33 -trimlog strain1.trimlog\
 strain1_1.fastq.gz strain1_2.fastq.gz\
 strain1_1.paired.fastq.gz strain1_2.paired.fastq.gz\
 strain1_1.unpaired.fastq.gz strain1_2.unpaired.fastq.gz\
 ILLUMINACLIP:/usr/local/Trimmomatic/adapters/TruSeq3-PE.fa:2:30:10\
 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:120
java -jar -Xmx256g /usr/local/Trimmomatic/trimmomatic-0.38.jar PE\
 -threads 16 -phred33 -trimlog strain2.trimlog\
 strain2_1.fastq.gz strain2_2.fastq.gz\
 strain2_1.paired.fastq.gz strain2_2.paired.fastq.gz\
 strain2_1.unpaired.fastq.gz strain2_2.unpaired.fastq.gz\
 ILLUMINACLIP:/usr/local/Trimmomatic/adapters/TruSeq3-PE.fa:2:30:10\
 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:120


Breseq処理

breseq -j 32 -o strain1 -r ./AP012030.gb strain1_1.fastq.gz strain1_2.fastq.gz >& strain1.out &
breseq -j 32 -o strain2 -r ./AP012030.gb strain2_1.fastq.gz strain2_2.fastq.gz >& strain2.out &

  出力


トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2019-10-27 (日) 09:28:45 (22d)