Skip to content

Commit a398940

Browse files
khetherintcezard
authored andcommitted
handling for if no reference provided and reference_seq not in GVF
1 parent aa8eea1 commit a398940

File tree

1 file changed

+5
-1
lines changed

1 file changed

+5
-1
lines changed

convert_gvf_to_vcf/convertGVFtoVCF.py

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -516,7 +516,11 @@ def get_ref(self):
516516
if "Reference_seq" in self.vcf_value.keys():
517517
return self.vcf_value["Reference_seq"] # attributes:reference_seq
518518
else:
519-
reference_allele = extract_reference_allele(self.assembly, self.chrom, self.pos)
519+
if self.assembly:
520+
reference_allele = extract_reference_allele(self.assembly, self.chrom, self.pos)
521+
else:
522+
print("WARNING: No reference provided. Placeholder inserted for Reference allele.")
523+
reference_allele = "."
520524
return reference_allele
521525

522526
def __str__(self):

0 commit comments

Comments
 (0)