Skip to content

Commit 5235cc5

Browse files
committed
tidy up comments
1 parent d024351 commit 5235cc5

File tree

1 file changed

+1
-3
lines changed

1 file changed

+1
-3
lines changed

convert_gvf_to_vcf/convertGVFtoVCF.py

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -232,9 +232,7 @@ def extract_reference_allele(fasta_file, chromosome_name, position):
232232
:param position: position
233233
:return: reference_allele: base found at this chromosome_name at this position within this fasta_file
234234
"""
235-
# commented out the below due to low memory efficiency: https://biopython.org/docs/1.76/api/Bio.SeqIO.html#input-multiple-records
236-
# with open(fasta_file) as assembly:
237-
# records_dictionary = SeqIO.to_dict(SeqIO.parse(assembly, "fasta"))
235+
# using .index instead of .todict for memory efficiency: https://biopython.org/docs/1.76/api/Bio.SeqIO.html#input-multiple-records
238236
records_dictionary = SeqIO.index(fasta_file, "fasta")
239237
zero_indexed_position = position - 1 # minus one because zero indexed
240238
reference_allele = records_dictionary[chromosome_name].seq[zero_indexed_position]

0 commit comments

Comments
 (0)