一般情况可以使用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
本站由北京市万商天勤律师事务所王兴未律师提供法律服务