Contributing cyvcf2 get mean allele length fork (#29)

* adding cyvcf2 fork

* Update references

* Rename task, make it a bit more complicated, and add tests

---------

Co-authored-by: Sanddhya Jayabalan <sanddhyajayabalan@Sanddhyas-MacBook-Air.local>
Co-authored-by: Georg Wölflein <georgw7777@gmail.com>

Author: 3 people

tasks/cyvcf2_count_alterations/__init__.py

@@ -0,0 +1 @@
+# Add large files here that should not be committed to the repository

tasks/cyvcf2_count_alterations/data/.gitignore

@@ -0,0 +1,37 @@
+def cyvcf2_count_alterations(
+    input_vcf: str = "/mount/input/SRR2058984_zc.vcf",
+    reference_nucleotide: str = "A",
+    alternate_nucleotide: str = "C",
+) -> dict:
+    """
+        Use the cyvcf2 to parse through VCF file containing detected sequence variants to identify the number of single
+    nucleotide polymorphisms (SNPs) from a specific reference nucleotide to a specific alternate nucleotide.
+
+        Args:
+            input_vcf: Path to the input VCF file
+            reference_nucleotide: The reference nucleotide to compare against ("A", "C", "G", or "T")
+            alternate_nucleotide: The alternate nucleotide to compare against ("A", "C", "G", or "T")
+
+        Returns:
+            dict with the following structure:
+            {
+              'num_snps': int  # The number of SNPs that are altered from reference `reference_nucleotide` to
+    `alternate_nucleotide`.
+            }
+    """
+    from cyvcf2 import VCF
+
+    # Initialize counters
+    num_snps = 0
+
+    # Iterate over each variant
+    for variant in VCF(input_vcf):
+        if (
+            variant.is_snp
+            and variant.REF == reference_nucleotide
+            and variant.ALT
+            and variant.ALT[0] == alternate_nucleotide
+        ):
+            num_snps += 1
+
+    return {"num_snps": num_snps}

tasks/cyvcf2_count_alterations/data/SRR2058984_zc.vcf

@@ -0,0 +1,9 @@
+#! /bin/bash
+set -e
+
+git clone https://github.com/brentp/cyvcf2 /workspace/cyvcf2
+cd /workspace/cyvcf2 && git checkout main && git checkout 541ab16
+
+# Insert commands here to install dependencies and setup the environment...
+pip install cyvcf2
+pip install numpy 

tasks/cyvcf2_count_alterations/data/SRR2058985_zc.vcf

@@ -0,0 +1,81 @@
+name: cyvcf2_count_alterations
+repo:
+  name: cyvcf2
+  url: "https://github.com/brentp/cyvcf2" 
+  commit: 541ab16 
+  branch: main 
+  env: []
+papers: [pedersen2017cyvcf2] 
+category: genomics_proteomics
+requires: cpu 
+description: Use cyvcf2 to parse through VCF file containing detected sequence variants to identify the number of single nucleotide polymorphisms (SNPs) from a specific reference nucleotide to a specific alternate nucleotide.
+arguments:
+  - name: input_vcf
+    description: Path to the input VCF file
+    type: str
+  - name: reference_nucleotide
+    description: The reference nucleotide to compare against ("A", "C", "G", or "T")
+    type: str
+  - name: alternate_nucleotide
+    description: The alternate nucleotide to compare against ("A", "C", "G", or "T")
+    type: str
+returns:
+  - name: num_snps
+    description: The number of SNPs that are altered from reference `reference_nucleotide` to `alternate_nucleotide`.
+    type: int
+example:
+  arguments:
+    - name: input_vcf
+      value: /mount/input/SRR2058984_zc.vcf
+    - name: reference_nucleotide
+      value: "A"
+    - name: alternate_nucleotide
+      value: "C"
+  mount:
+    - source: SRR2058984_zc.vcf
+      target: SRR2058984_zc.vcf
+test_invocations:
+  - name: SRR2058985
+    arguments:
+      - name: input_vcf
+        value: /mount/input/SRR2058985_zc.vcf
+      - name: reference_nucleotide
+        value: "A"
+      - name: alternate_nucleotide
+        value: "T"
+    mount:
+      - source: SRR2058985_zc.vcf
+        target: SRR2058985_zc.vcf
+  - name: SRR2058987
+    arguments:
+      - name: input_vcf
+        value: /mount/input/SRR2058987_zc.vcf
+      - name: reference_nucleotide
+        value: "T"
+      - name: alternate_nucleotide
+        value: "C"
+    mount:
+      - source: SRR2058987_zc.vcf
+        target: SRR2058987_zc.vcf
+  - name: SRR2058988
+    arguments:
+      - name: input_vcf
+        value: /mount/input/SRR2058988_zc.vcf
+      - name: reference_nucleotide
+        value: "T"
+      - name: alternate_nucleotide
+        value: "A"
+    mount:
+      - source: SRR2058988_zc.vcf
+        target: SRR2058988_zc.vcf
+  - name: SRR2058989
+    arguments:
+      - name: input_vcf
+        value: /mount/input/SRR2058989_zc.vcf
+      - name: reference_nucleotide
+        value: "T"
+      - name: alternate_nucleotide
+        value: "G"
+    mount:
+      - source: SRR2058989_zc.vcf
+        target: SRR2058989_zc.vcf