Linux查看基因数据压缩包.vcf.gz文件内容

已于 2022-03-26 10:14:27 修改
阅读量5.4k 收藏 7
点赞数 3
分类专栏： Timo学习笔记 文章标签： shell
于 2022-03-26 10:13:21 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/timo_csdn/article/details/123751243
版权

得到超大的基因数据压缩包.vcf.gz，我们想知道他的id是怎样的，里面都包含了什么信息，所以想查看一下他的前几行。
但是用一般的看压缩包的方式是看不到的，只会看到乱码。

less -S yourfile.vcf.gz

即可显示文件前几行的内容。
图为1000 genome的数据
在这里插入图片描述
按 f或 空格键，继续print到屏幕上。
按q退出，返回命令行。

野柚子__
关注
  • 3
    点赞
  • 踩
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
从 Excel 导出 联系人名片文件（.vcf）安卓
08-10
从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓从 Excel 导出 联系人名片文件（.vcf）安卓
Linux 三剑客：grep、awk、sed
bushengxinxiaobai的博客
08-28 207
Linux 三剑客：grep、awk、sed (介绍并不仔细，只是分享曾经用过的命令行，也作为自己的一个笔记，后续会不断填充命令) 输入文件：testchr01.vcf.gz 输出文件：chr01.vcf 1.grep 优势：查找功能 举例： 1.1 输出 testchr01.vcf.gz中带##（表头部分） less -S testchr01.vcf.gz | grep -v ’ ^## ’ > chr01.vcf 1.2 输出已用命令行中关于grep的命令 history | grep ‘gre
linux 查看vcf文件,Linux生信练习4-vcf
weixin_28887297的博客
04-29 772
原始数据准备cd biosoft/bowtie2/bowtie2-2.4.1-linux-x86_64/example/reads/../../bowtie2 -x ../index/lambda_virus -1 reads_1.fq -2 reads_2.fq | samtools sort -@ 5 -o tmp.bam -bcftools mpileup -f ../reference/l...
Linux压缩/解压命令
Jason的博客
03-10 508
1、tar 压缩（归档） tar 用来压缩和解压文件, tar 命令本身没有压缩、解压的功能，是调用其他压缩功能实 现的。tar 压缩后的文件也叫归档文件。 语法: tar 参数 要压缩或解压的文件或目录 常用参数： z : 使用 gzip 解压缩程序，生成的文件名是 xxx.tar.gz 这是 linux 中常用的压缩格式。 c : 创建压缩文档 v : 显示压缩，解压过程中处理的...
php简单读取.vcf格式文件的方法示例
10-19
主要介绍了php简单读取.vcf格式文件的方法,结合具体实例形式分析了php自定义函数读取vcf格式文件的具体实现方法与相关注意事项,需要的朋友可以参考下
如何打开vcf.gz文件
weixin_53682198的博客
04-22 1370
解压后的文件会保存在当前目录中，文件名为原始文件名去掉"gz"的后缀。如果要指定解压后的文件名和保存路径，可以使用"-c"选项和重定向符”＞"来实现，例如：gzip -dc file.gz > /path/to/output/file。
VCF变异文件读取和详细
热门推荐
weixin_42594968的博客
08-02 1万+
下载了千人基因组variants数据，是.vcf.gz和.vcf.gz.tbi文件格式，需要在linux上打开，此时我还是一个生信小菜鸡，记录一下vcf文件读取和详细解释。 在服务器先解压再打开， gunzip 文件.vcf.gz cat 文件.vcf 结果如下，60G+的文件太大了，只好读取一部分： VCF各列意义说明 各列之间用tab空白隔开；前面9列为固定列，第10列开始为样品信息列，可以无限多个； #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT 后面
Linux环境下利用perl脚本批量筛选VCF文件指定样本
青笋的博客
11-18 938
今天分享一个实用的数据提取思路：从VCF文件中提取指定列，比如有1000个样本，需要从中提取出200个制定样本，形成新的vcf文件。需要在Linux环境下进行操作，主要用到perl脚本和shell脚本。 原始数据 原始文件vcf ：开始为若干行注释信息以##开头，然后是表头信息由#开始（该行包括样本ID） 待提取的样本ID：txt文件，每个样本ID一行（list.txt） 什么是VCF文件？ 伴随着大规模的基因分型及测序数据的产生，迫切需要一种新的格式来记录高效的记录这些变异信息。VCF(Vari
linux 查看vcf文件,转载-VCF格式详解
weixin_35619710的博客
04-29 1144
CHROM(chromosome):染色体POS - position:参考基因组variant碱基位置，如果是INDEL(插入缺失)，位置是INDEL的第一个碱基位置ID - identifier:variant的ID。比如在dbSNP中有该SNP的id，则会在此行给出；若没有，则用’.'表示其为一个novel variant。REF - reference base(s):参考碱基，染色体上面...
vcf-0.9.8.tar.gz
01-16
vcf-0.9.8.tar.gz 牛逼的C++ ui库vcf-0.9.8.tar.gz 源代码
vcf2DIYABC.py:将vcf文件转换为基因型数据矩阵以输入DIYABC程序的脚本
05-22
vcf2DIYABC.py 目的 转换DIYABC输入文件中的VCF文件以获取SNP数据的工具 要求 python 2.x 分支python3.x上提供了与python 3.x兼容的版本。 Windows用户：如果您的计算机上未安装python，则此脚本无法运行。 我强烈建议您通过以下来完成此步骤，因为它可能对许多其他生物信息学目的很有用。 待办事项：我目前正在尝试生成一个自包含的可执行文件来运行脚本，而无需使用进行python安装。 下载/安装 如果git用户使用以下命令在其系统上安装了git，则可以克隆该存储库： git clone https://github.com/loire/vcf2DIYABC.snp.git 此存储库的zip归档文件的链接在此页面的右下角。 要安装，只需将存档解压缩到本地目录中即可。 用法 Unix和Mac用户可以使用终端和类型 python
基因型文件转换为vcf文件.rar_vcf中./. 基因型_基因型文件转换为vcf文件
09-15
将基因型文件转换为vcf文件，可以直接运行，压缩文件中有例子
vcf.iobio.io
05-01
vcf.iobio.io
Contacts-2022-05-18 (1).vcf
05-19
Contacts-2022-05-18 (1).vcf
压缩算法的原理丨基因型vcf文件为什么压缩后发生了什么？
最新发布
青笋的博客
07-28 388
本文中提到的vcf文件是存储基因型变异数据的常用格式，VCF（Variant Call Format）通常用于存储个体或群体的基因组数据，其中包含了多个基因突变位点的信息。在VCF文件中，可能存在大量的重复信息，比如在多个位点上的质量分数或者过滤标记可能会重复出现，gzip能够识别这些重复模式，并使用更短的编码方式来表示它们，从而减少存储空间。接下来的两列分别是两个样本的基因型数据。因为gzip是一种有效的压缩算法，它可以利用重复出现的模式和冗余信息，将数据进行编码和压缩，从而减少数据的存储空间。
实操 | 合并VCF文件的几种方法及注意事项
悟道西方
08-15 9502
背 景 在基因组分析领域的很多不同场景中，需要合并VCF文件。VCF文件。简单来说，就是记录样本基因型的文件。但多数VCF文件不只记录了基因型，也包含有关该基因型的来源的细节。其它文件。VCF文件的上游是BAM文件，主要记录Reads与参考基因组的比对信息；更上游的，就是FASTQ测序数据，以及物种的参考基因组。不同类型的VCF文件。VCF文件有单样本的、多样本的；也有普通VCF文件...
linux 查看vcf文件,VCF格式文件的shell小练习
weixin_42488542的博客
04-29 969
首先使用bowtie2软件自带的测试数据生成sam/bam文件，还有vcf文件代码如下：mkdir -p ~/biosoftcd ~/biosoftwget https://mirrors.tuna.tsinghua.edu.cn/anaconda/miniconda/Miniconda3-latest-Linux-x86_64.shbash Miniconda3-latest-Linux-x86...
linux下查看vcf文件格式,plink格式文件与vcf格式文件相互转化
weixin_31243809的博客
05-12 954
有的时候在做数据分析的时候，会用到各种各样的分析软件，而每个分析软件又有与之相匹配的输入文件格式。格式不对，会很麻烦，有时候会出现错误的结果，甚至就不能运行程序。而许多软件之间的文件格式是可以相互转化的。直接转化就可以方便很多。今天就来讲一讲plink的ped和map文件与vcf文件的相互转化。因为在许多情况下这两种文件格式是需要转化的。首先介绍一下由plink转vcfplink --file t...
gatk GenotypeGVCFs -R ref.fa -V chrx.merged.g.vcf.gz -O chrx.genotype.vcf.gz各个参数的意义，并举例
06-10
1. `-R ref.fa`：指定参考基因组文件，`ref.fa`是参考基因组文件的路径和文件名。 2. `-V chrx.merged.g.vcf.gz`：指定输入的GVCF文件，`chrx.merged.g.vcf.gz`是输入文件的路径和文件名。 3. `-O chrx.genotype.vcf.gz`：指定输出文件的路径和文件名，`chrx.genotype.vcf.gz`是输出文件的路径和文件名。 举例来说，假设我们有一个参考基因组文件`hg38.fa`，一个包含多个样本的GVCF文件`sample1.g.vcf.gz, sample2.g.vcf.gz, sample3.g.vcf.gz`，我们想要对这些样本进行基因型调用并输出到一个VCF文件中，我们可以使用以下命令： ``` gatk GenotypeGVCFs -R hg38.fa -V sample1.g.vcf.gz -V sample2.g.vcf.gz -V sample3.g.vcf.gz -O genotypes.vcf.gz ``` 其中，`-R`指定参考基因组文件，`-V`指定输入的GVCF文件，可以在命令中使用多个`-V`参数来指定多个输入文件，最终的基因型调用结果将会输出到`genotypes.vcf.gz`文件中。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值