sra 資料轉成 fastq並改名

2021-09-25 07:33:17 字數 1301 閱讀 2046

把sra資料移動到我們工作目錄後,我們開始sra轉faq。

正式執行**之前,必須先拿乙個樣品測試下**能否執行成功,這點很關鍵,因為這步就算成功執行也特別慢,要是**再出錯了就更浪費時間了。

ls srr5315196.sra |fastq-dump -gzip --split-3 -o ./ srr5315196.sra

sra-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 -gzip

–split-3 -o ./ id)

1>./

) 1>./

id)1

>./

.sra.log 2>&1;done #多個迴圈 nohup bash sra.sh & #掛後台執行

d表示

輸出(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...