安装

pip install PyVCF

使用

#如果有多个样本,可以这样访问最后两列信息,例如双样本,0为normal,1为tumor

record.samples[0]

对vcf进行筛选,并重新生成vcf

切记用完一定要close()

import vcffilaname1='xxx.vcf'
filename2='xxx.vcf'
tumor_idx=1
vcf_reader = vcf.Reader(open(filename1, 'r'))
vcf_writer = vcf.Writer(open(filename2,'w'),vcf_reader)
for record in vcf_reader:if len(record.FILTER) == 0: #passaf = record.samples[tumor_idx]['AF']if af >= 0.01:vcf_writer.write_record(record)
vcf_writer.close()

也可以读取vcf.gz

import vcf vcf_name = "xx.vcf.gz"
vcf_reader = vcf.Reader(filename = vcf_name)
for record in vcf_reader:print(record.CHROM)  ##strprint(record.POS)    ##intprint(record.REF)    ##strprint(record.ALT)    ##listprint(record.INFO)   ##dictprint(record.FORMAT)   ##str

python操作VCF文件相关推荐

  1. Python操作Zip文件

    Python操作Zip文件 需要使用到zipfile模块 读取Zip文件 随便一个zip文件,我这里用了bb.zip,就是一个文件夹bb,里面有个文件aa.txt.

  2. python处理excel表格实例-使用Python操作excel文件的实例代码

    使用的类库 pip install openpyxl 操作实现 •工作簿操作 # coding: utf-8 from openpyxl import Workbook # 创建一个excel工作簿 ...

  3. python入门教程2word-入门干货:Python操作Word文件经验分享

    原标题:入门干货:Python操作Word文件经验分享 导读:Microsoft Word在当前使用中是占有巨大优势的文字处理器,这使得Word专用的档案格式Word 文件(.docx)成为事实上最通 ...

  4. 手机使用python操作图片文件

    手机使用python操作图片文件 起因 前几天去国图拍了一本书,一本心理学方面的书,也许你问我为什么不去买一本,或者去网上找pdf. 其实吧,关于心理学方面的书可以说在市面上一抓就是一堆,至于拍这本书 ...

  5. 详解Python操作Excel文件

    前言 本篇文章主要总结了一下利用python操作Excel文件的第三方库和方法. 常见库简介 1.xlrd xlrd是一个从Excel文件读取数据和格式化信息的库,支持.xls以及.xlsx文件. 地 ...

  6. python实现xmind_Python 使用Python操作xmind文件

    使用Python操作xmind文件 by:授客QQ:1033553122 测试环境 Win10 Python 3.5.4 XMind-1.2.0.tar.gz 下载地址: 创建及更新xmind文件 # ...

  7. 如何使用 Python 操作 .npy 文件?详细教程分享

    如何使用 Python 操作 .npy 文件?详细教程分享 .npy 文件是指 NumPy 存储数组数据的文件格式,它具有快速.高效和易存储等优点.在科学计算和机器学习领域中,.npy 文件的使用非常 ...

  8. Python操作Excel文件:插入一列数据

    问题描述:已有Excel文件,其中包含5列数据,要求在第3列前插入一列数据,保存为新文件. 参考代码: 运行结果: --------董付国老师Python系列教材-------- 1)<Pyth ...

  9. Python 操作Excel文件

    Python 操作Excel文件 文章目录 Python 操作Excel文件 1. 用到的三方库 2. 读excel文件 3. 写excel 文件 简单的总结 Excel 是我们在日常办公中可谓是无处 ...

  10. python 操作PDF文件 之 A3页面转A4

    python 操作PDF文件 A3页面转A4页面 文章目录 1. 需求概述 2. 代码实现 1. 需求概述 最近接到一份PDF资料需要打印,奈何页面是如图所示的A3格式的,奈何目前条件只支持打印A4. ...

最新文章

  1. ActiveMQ死信队列使用
  2. Python网络爬虫笔记:下载、分析京东P20销售数据
  3. 记录工作中常用的CSS3
  4. Android开发日记(六)
  5. 中如何构造有参和无惨_导数中构造函数的神力你应该借助一下
  6. LeetCode 简单等级
  7. Js字符串与十六进制的相互转换 【转】
  8. 关系型数据库学习手记——初见倾心PostgreSQL、MySQL、SQLite、MongoDB
  9. 基于阿里云服务器使用宝塔面板进行Hexo建站
  10. 独立安装 lamp 兄弟连
  11. IntelliJ IDEA设置豆沙绿背景
  12. 如何用仅用C语言判断编译器的大小端
  13. 该怎么复习安徽省考计算机专业课
  14. opencv中calcOpticalFlowPyrLK实现的光流法(Lucas-Kanade Method for Sparse Optical Flow)原理解析 (摘要翻译)
  15. 通过镜像安装torch,torchvision
  16. 苹果展开新显示器带动高阶需求:Mini LED背光技术
  17. 冥想的重大功能——人类21在21世纪的伟大发现
  18. 【渝粤题库】陕西师范大学152205 西方行政学说史 作业(专升本)
  19. Smba服务器的配置
  20. VS中项目属性win32 application和win32 console application的切换

热门文章

  1. 【vue学习笔记】vue-cli安装问题
  2. 饭谈:失眠,还有梦魇,第二天要上班应该怎么办?
  3. Python:实现monte carlo蒙特卡罗算法(附完整源码)
  4. 日紫白飞星算法_年、月、日、时紫白飞星法!
  5. Android 讯飞语音合成、语音播报(详细步骤+源码)
  6. 网页实现数据离线永久保存localStorage、storage
  7. 简单实现网页登录框中的“记住我”
  8. 【网络安全】GitHub项目监控,teemo子域名查询
  9. 失业的程序员(十二):潜意识的智商
  10. 关于DSP28335CCS6编译时出现error #10099-D: program will not fit into available memory.