How to read vcf.bgz or vcf.gz in python using PyVcf?

When I try to do simply like this:

vcf_reader = vcf.Reader(open("input.vcf.gz", 'r')) vcf_writer = vcf.Writer(open('input.vcf.gz', 'w'), vcf_reader)

there is an error:

UnicodeDecodeError: 'utf-8' codec can't decode byte 0x8b in position 1: invalid start byte

pyvcf vcf python

• 6,285 views • link

updated 6.7 years ago by dariober

• written 6.7 years ago by ja4123 •

1 answer

I haven't tried it by I think you should use the filename parameter:

vcf_reader = vcf.Reader(filename='input.vcf.gz')

• 0 views • link

written 6.7 years ago by dariober

More posts like this

Filter BAM by chromosome using PySam
written 3.8 years ago by pablosolar.r •

Hi all! I am trying to filter a BAM by chromosome using PySam. Using **Samtools**: samtools view -h mybam.bam 22 > mybam.bam.chr22.sam samtools view -bS …
Error when opening CEL files in python
written 4.9 years ago by Paula •

Hi! I am trying to open a CEL file using: with open("pathtofile/file.CEL", "r") as handle: c = CelFile.read(handle) and it is returning the error: --------------------------------------------------------------------------- …
Issue reading fasta file with Biopython
written 5.0 years ago by Rox

Hello everyone, This should be very easy and I know it, but I am stuck with it and I cannot pinpoint my mistake. I wanted …
Ensembl encoding problem
written 5.1 years ago by giammafer •

Hi everybody, I have a problem with two files downloaded from Ensembl FTP Homo_sapiens.GRCh38.pep.all.fa.gz - http://ftp.ensembl.org/pub/release-104/fasta/homo_sapiens/pep/ Homo_sapiens.GRCh38.104.chr.gtf.gz - http://ftp.ensembl.org/pub/release-104/gtf/homo_sapiens/ I need to open it on …
GTFtools with Homo_sapiens.GRCh38.103.gtf encoding error in python
written 5.3 years ago by Doreen •

I adapted [GTF tools][1] to work with python3. The demo.gtf runs fine, but if I try using it with Homo_sapiens.GRCh38.103.gtf from [ensembl][2] I get Decoding …
Demultiplexing fastq files
written 6.1 years ago by sidrah.maryam •

I am doing demultiplexing of fastq files using demultiplex. The command is: system2("demultiplex demux /home/sidrah19220/Thesis/new/fastqfile/atac_v1_pbmc_10k_fastqs/index_fastq.gz /home/sidrah19220/Thesis/new/fastqfile/atac_v1_pbmc_10k_fastqs/pe1_fastq.gz /home/sidrah19220/Thesis/new/fastqfile/atac_v1_pbmc_10k_fastqs/pe2_fastq.gz", ) But it is showing error": **demultiplex: error: …
VCF filtering by allele frequency
written 6.6 years ago by Jimpix •

Hi, I am trying to filtering a vcf file by allele frequency. I'm using a file downloaded from the gnomad v2 liftover. I am using …
Reading vcf file using python gives UnicodeDecodeError
written 7.5 years ago by Medhat

python version `Python 3.5.4` Trying to process vcf file as shown in the code below: with open(myfile, 'r') as datain: for myline in datain: if …
Failed to liftover from hg38 to hg19
written 7.6 years ago by BAGeno

Hi, I am converting my file from hg38 to hg19. For this purpose I used `crossmap` and `picard LiftoverVcf`. But I am getting this error …
IDR on blacklisted region-filtered narrowPeak
written 8.3 years ago by endicytosis •

Hi guys, I am having a problem performing IDR on my blacklist region-filtered file. After I called my narrow peaks using MACS2, I used bedtools …