import cyvcf2 # Open VCF file vcf = cyvcf2.VCF('filename.vcf') # Print number of samples in VCF file print("Number of samples:", vcf.num_samples) # Iterate through each variant in the VCF file for variant in vcf: print(variant.CHROM, variant.POS, variant.REF, variant.ALT)
import cyvcf2 # Open VCF file vcf = cyvcf2.VCF('filename.vcf') # Iterate through each variant in the VCF file for variant in vcf: # Iterate through each sample for sample in variant.samples: # Print sample ID and genotype print(sample.sample, sample.gt_bases)Package library - cyvcf2 is a Python package library.