這篇文章給大家分享的是有關linux中fastq文件如何查看指定序列的內容。小編覺得挺實用的,因此分享給大家做個參考,一起跟隨小編過來看看吧。
在處理fastq文件時,可能會遇到序列與質量值長度不一致或格式不正確等問題,這時候軟件會報告哪條序列有問題,并給出序列ID,然后我們就要去查看這條序列是什么問題。但是,fastq文件通常是很大的,使用less命令可能很難搜索到,這時候我們可以用到下面的命令來查看。
zcat /share/nas1/DATA/zx-20180122-15/XM8_H33TNDMXX_L1_2.clean.fq.gz |grep -E "A00164:71:H33TNDMXX:2:2438:30400:20353" -C 4 >xm8-2.txt &
其中zcat可以查看fastq壓縮文件,grep 匹配序列ID,參數 -C 4 可以輸出匹配行的前后四行。
感謝各位的閱讀!關于“linux中fastq文件如何查看指定序列”這篇文章就分享到這里了,希望以上內容可以對大家有一定的幫助,讓大家可以學到更多知識,如果覺得文章不錯,可以把它分享出去讓更多的人看到吧!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。