bin/bfr.rb in bio-polyploid-tools-0.1.0 vs bin/bfr.rb in bio-polyploid-tools-0.2.3
- old
+ new
@@ -1,5 +1,6 @@
+#!/usr/bin/env ruby
require 'rubygems'
#require 'extensions/all'
require 'bio-samtools'
require 'optparse'
@@ -68,29 +69,23 @@
chunk = options[:chunk]
chunk_size = options[:chunk_size]
output_filename = options[:output_filename]
stats_file = options[:stats_file]
-
-#reference = ARGV[6]
-
-
min = chunk * chunk_size
max = min + chunk_size
-#AvocetS
parental_1=options[:parent_1]
-#AvocetS (Yr15)
parental_2=options[:parent_2]
bulk_1 = options[:bulk_1]
bulk_2 = options[:bulk_2]
-fasta_db = Bio::DB::Fasta::FastaFile.new(reference)
+fasta_db = Bio::DB::Fasta::FastaFile.new({:fasta=>reference})
fasta_db.load_fai_entries
if chunk_size == 0
min = 0