awk命指定分隔符输出字符串///使用bgzip遇到的报错是怎样的

awk命指定分隔符输出字符串///使用bgzip遇到的报错是怎样的,相信很多没有经验的人对此束手无策,为此本文总结了问题出现的原因和解决方法,通过这篇文章希望你能解决这个问题。

成都创新互联专注于企业营销型网站、网站重做改版、威信网站定制设计、自适应品牌网站建设、H5场景定制成都做商城网站、集团公司官网建设、成都外贸网站建设公司、高端网站制作、响应式网页设计等建站业务,价格优惠性价比高,为威信等各大城市提供网站开发制作服务。

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解压就没有报错了

看完上述内容,你们掌握awk命指定分隔符输出字符串///使用bgzip遇到的报错是怎样的的方法了吗?如果还想学到更多技能或想了解更多相关内容,欢迎关注创新互联行业资讯频道,感谢各位的阅读!


当前标题:awk命指定分隔符输出字符串///使用bgzip遇到的报错是怎样的
转载来源:http://pwwzsj.com/article/pechdd.html