Read a VCF file

public function readVcf(VcfReaderInterface $vcfReader)
{
    // The example file of the VCF 4.3 specification
    $oResult = $vcfReader->read(file($this->getParameter('kernel.project_dir') . '/public/data/sample.vcf'));
    $aVariants = $oResult->getVariants();
 
    // Only the variants that passed every filter
    $aPassing = array_filter($aVariants, fn(VcfVariant $oVariant) => $oVariant->isPass());
 
    // Malformed lines are reported and skipped, the others are read
    $oFaulty = $vcfReader->read([
        "#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO",
        "20\t14370\trs6054257\tG\tA\t29\tPASS\tDP=14",
        "20\tabc\t.\tT\tA\t3\tq10\tDP=11",
        "20\t17330\t.\tT\tA\t3\tq10",
    ]);
 
    return $this->render('default/variantsReadVcf.html.twig', [
        'variants'  => $this->describe($aVariants), // type of each alternate allele (SNV, insertion...)
        'passing'   => count($aPassing),
        'faulty'    => $oFaulty,
    ]);
}