Merge pull request #21 from PacificBiosciences/cg_sites_dict-fix

bug fix in cg_sites_from_fasta()
PacificBiosciences · Jun 21, 2022 · afd63e8 · afd63e8
2 parents 8b51033 + e22c576
commit afd63e8
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/aligned_bam_to_cpg_scores.py b/aligned_bam_to_cpg_scores.py
@@ -163,6 +163,8 @@ def cg_sites_from_fasta(input_fasta, ref):
     :return cg_sites_ref_set: Set with all CG ref positions. (set)
     """
     # open fasta with BioPython and iterated over records
+    cg_sites_ref_set = set()
+
     with open(input_fasta) as fh:
         for record in SeqIO.parse(fh, "fasta"):
             # if record name matches this particular ref,