把sra資料移動到我們工作目錄後,我們開始sra轉faq。
正式執行**之前,必須先拿乙個樣品測試下**能否執行成功,這點很關鍵,因為這步就算成功執行也特別慢,要是**再出錯了就更浪費時間了。
ls srr5315196.sra |fastq-dump -gzip --split-3 -o ./ srr5315196.srasra-tools 裡的 fastq-dump工具可以將srr檔案轉換為fastq格式,–split-3引數表示如果是雙端測序就自動拆分,如果是單端不受影響。也就是說,–split-3引數可以將pe的sra檔案解壓後的fastq檔案拆分成_1.fastq和_2.fastq,如果示例資料集是se測序,不會進行拆分。–gzip轉換fastq為壓縮檔案,節省空間。
單個測試成功,那我們就寫迴圈進行批量轉換格式。
cat >sra.sh #寫指令碼 ls srr* |while read id; do (fastq-dump -gzipd表示–split-3 -o ./ id)
1>./
) 1>./
id)1
>./
.sra.log 2>&1;done #多個迴圈 nohup bash sra.sh & #掛後台執行
輸出(o
utpu
t)到當
前資料夾
下,檔案
名字首不
變(這樣
的話輸出
張這樣s
rr5315196.fa
stq.
gz)1
>./
表示輸出(output)到當前資料夾下,檔名字首不變(這樣的話輸出張這樣srr5315196.fastq.gz) 1>./
id表示輸出
(out
put)
到當前文
件夾下,
檔名字首不變(
這樣的話
輸出張這
樣srr
5315
196.
fast
q.gz
)1>./
.sra.log 2>&1 (0.標準輸入;1.標準輸出;2.標準錯誤)這個命令是說,重定向標準輸出到當前檔案下srr5315197.sra.sra.log檔案,且標準輸出、標準錯誤到乙個檔案中(2>&1)
製作config檔案
C語言實戰課 將FASTQ轉成FASTA
我這次學習c語言的目的非常單純,就是嘗試將c語言應用在日常的分析任務重,解決實際問題。既然如此,那麼第一課就不是列印 hello world 了,畢竟說了那麼多次世界你好,依舊寫不好 我們第一課直接就來處理乙個實際的小需求,讀取fastq,將其轉成乙個fasta。處理這個問題和把大象放進冰箱裡一樣,...
fastq資料質控過濾軟體 soapnuke 使用
安裝2.0version 參照readme.md for 2.x make soapnuke 使用測試 soapnuke filter 1 home extend u r exam newbgiseq500 1.fq.gz 2 home extend u r exam newbgiseq500 2....
將扁平資料轉成樹狀資料
將扁平資料轉成樹狀資料,轉化中間可以新增自己想要的字段 序號等等 export const flattranstree ids,list else brancharr.length 0?father.children brancharr return father.pid 0 return newa...