您好,欢迎来到欧得旅游网。
搜索
您的当前位置:首页使用linux命令 根据序列ID快速提取fastq序列

使用linux命令 根据序列ID快速提取fastq序列

来源:欧得旅游网

一般情况可以使用grep命令,就可以实现

grep -A 3 seq_id fastq.file
# 匹配seq_id 并向下取3行

但是,如果需要的序列很多,在使用grep就会很慢了,所以这里给出了 awk 的命令,速度简直快的飞起。

awk -F ' ' 'NR==FNR {a[$1]=1; next} { if (a[substr($1,2)]) {print $0; getline b; print b ;getline b; print b ;  getline b; print b ;}}'  list.txt  fastq.file

list.txt:中是要提取的序列名

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- ovod.cn 版权所有 湘ICP备2023023988号-4

违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务