Skip to content

Commit 6de6920

Browse files
committed
handling for if no reference provided and reference_seq not in GVF
1 parent eca888e commit 6de6920

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)