bin/homokaryot_primers.rb in bio-polyploid-tools-0.1.0 vs bin/homokaryot_primers.rb in bio-polyploid-tools-0.2.3
- old
+ new
@@ -80,10 +80,10 @@
snps = Array.new
#0. Load the fasta index
fasta_reference_db = nil
if reference_file
- fasta_reference_db = Bio::DB::Fasta::FastaFile.new(reference_file)
+ fasta_reference_db = Bio::DB::Fasta::FastaFile.new({:fasta=>reference_file})
fasta_reference_db.load_fai_entries
p "Fasta reference: #{reference_file}"
end
#1. Read all the SNP files
#All the SNPs should be on the same chromosome as the first SNP.
\ No newline at end of file