awk命指定分隔符输出字符串///使用bgzip遇到的一个报错
awk命指定分隔符输出字符串///使用bgzip遇到的一个报错
用户7010445
发布于 2020-11-13 09:49:01
发布于 2020-11-13 09:49:01
awk指定字符分割字符串、指定分隔符输出字符串
遇到的问题
使用blasr软件将三代测序数据比对到参考序列
blasr longreads.fastq reference.fasta --nproc 16 > blasr.out部分输出结果
m54155_170415_100314/5309390/25118_26816/0_1698 reference 0 1 -3020 75.3097 127858 128847 510317 3 1182
m54155_170415_100314/5243602/0_9742/0_9742 reference 0 1 -1076 79.0576 17916 18284 510317 8946 9296 9742
m54155_170415_100314/5440071/0_9295/0_9295 reference 0 0 -1122 91.2409 470798 471063 510317 0 267 9295 5这个地方不知道为什么 reads 的 ID 多了后面一个部分。如果利用这个ID再来提取比对上的reads时就得不到结果
可以利用awk命令把结尾的部分去掉
参考链接 https://blog.csdn.net/liangbilin/article/details/108593296
cat blasr.out | awk '{print $1}' | awk -F '/' -v OFS="/" '{print $1,$2,$3}' > blasr.out1-F 指定输入文件的的分隔符 -v OFS 指定输出文件的分隔符
bgzip遇到的报错及解决办法
这个服务器上没有bgzip这个命令，我使用conda进行安装
conda install tabix这个安装的是 0.2.6版本
解压fastq文件时
bgzip -d pacbio.sequel.fastq.gz遇到报错
Error: 2然后我卸载，重新安装0.2.5试一下
conda uninstall tabix
conda install tabix=0.2.5再次解压遇到报错
Error: invalid block header以上报错不知道什么原因，搜索一番后看到有人说安装好 htslib后就可以直接使用bgzip了。我试了一下
conda uninstall tabix
conda install htslib果然这次再用bgzip解压就没有报错了
什么原因还是不知道，现在暂时可以用了
本文参与 腾讯云自媒体同步曝光计划，分享自微信公众号。
原始发表：2020-11-08，如有侵权请联系 cloudcommunity@tencent.com 删除
评论
登录后参与评论
推荐阅读
目录
腾讯云开发者
