snp-search 0.2.0 → 0.3.0

data/Gemfile

@@ -3,6 +3,10 @@ source "http://rubygems.org"
 # Example:
 #   gem "activesupport", ">= 2.3.5"
 
+ gem "activerecord"
+ gem "bio"
+ gem "slop"
+
 # Add dependencies to develop your gem here.
 # Include everything needed to run rake, tests, features, etc.
 group :development do
@@ -10,4 +14,5 @@ group :development do
   gem "bundler", "~> 1.0.0"
   gem "jeweler", "~> 1.6.4"
   gem "rcov", ">= 0"
-end
+ 
+  end

data/Gemfile.lock

@@ -1,12 +1,30 @@
 GEM
   remote: http://rubygems.org/
   specs:
+    activemodel (3.1.0)
+      activesupport (= 3.1.0)
+      bcrypt-ruby (~> 3.0.0)
+      builder (~> 3.0.0)
+      i18n (~> 0.6)
+    activerecord (3.1.0)
+      activemodel (= 3.1.0)
+      activesupport (= 3.1.0)
+      arel (~> 2.2.1)
+      tzinfo (~> 0.3.29)
+    activesupport (3.1.0)
+      multi_json (~> 1.0)
+    arel (2.2.1)
+    bcrypt-ruby (3.0.0)
+    bio (1.4.2)
+    builder (3.0.0)
     diff-lcs (1.1.3)
     git (1.2.5)
+    i18n (0.6.0)
     jeweler (1.6.4)
       bundler (~> 1.0)
       git (>= 1.2.5)
       rake
+    multi_json (1.0.3)
     rake (0.9.2.2)
     rcov (0.9.11)
     rspec (2.3.0)
@@ -17,12 +35,17 @@ GEM
     rspec-expectations (2.3.0)
       diff-lcs (~> 1.1.2)
     rspec-mocks (2.3.0)
+    slop (2.4.0)
+    tzinfo (0.3.29)
 
 PLATFORMS
   ruby
 
 DEPENDENCIES
+  activerecord
+  bio
   bundler (~> 1.0.0)
   jeweler (~> 1.6.4)
   rcov
   rspec (~> 2.3.0)
+  slop

data/Rakefile

@@ -21,6 +21,7 @@ Jeweler::Tasks.new do |gem|
   gem.description = %Q{Use the snp-search toolset to query the SNP database}
   gem.email = "ali.al-shahib@hpa.org.uk"
   gem.authors = ["Ali Al-Shahib", "Anthony Underwood"]
+  gem.executables = ["snp-search"]
   # dependencies defined in Gemfile
 end
 Jeweler::RubygemsDotOrgTasks.new
@@ -38,10 +39,9 @@ end
 
 task :default => :spec
 
-require 'rake/rdoctask'
-Rake::RDocTask.new do |rdoc|
+require 'rdoc/task'
+RDoc::Task.new do |rdoc|
   version = File.exist?('VERSION') ? File.read('VERSION') : ""
-
   rdoc.rdoc_dir = 'rdoc'
   rdoc.title = "snp-search #{version}"
   rdoc.rdoc_files.include('README*')

data/VERSION

@@ -1 +1 @@
-0.2.0
+0.3.0

data/bin/snp-search

@@ -0,0 +1,63 @@
+require 'snp-search'
+require 'snp_db_connection'
+require 'snp_db_models'
+require 'snp_db_schema'
+require 'slop'
+
+
+
+opts = Slop.new :help do
+  banner "ruby snp-search [OPTIONS]"
+
+  on :V, :verbose, 'Enable verbose mode'
+  on :n, :name=, 'Name of database', true 
+  on :r, :reference_file=, 'Path for the reference database, in gbk or embl file format'
+  on :v, :vcf_file=, 'Path for the .vcf file', true
+  on :s, :strain, 'Path for the list of strains text file', true
+  on :c, :cuttoff_snp=, 'cuttoff for SNP quality' 
+  on :t, :cuttoff_genotype=, 'cuttoff for genotype quality'
+  on_empty do
+    puts help
+  end
+end
+
+opts.parse
+
+strains = []
+File.read(opts[:strain]).each_line do |line|
+  strains << line.chop
+end
+
+# Enter the name of your database 
+establish_connection(opts[:name])
+
+# Schema will run here
+db_schema
+
+ref = opts[:reference_file]
+
+sequence_format = guess_sequence_format(ref)
+
+      case sequence_format
+      when :genbank
+        sequence_flatfile = Bio::FlatFile.open(Bio::GenBank,opts[:reference_file]).next_entry
+      when :embl
+        sequence_flatfile = Bio::FlatFile.open(Bio::EMBL,opts[:reference_file]).next_entry
+      else
+        puts "All sequence files should be of genbank or embl format"
+        exit
+      end
+
+
+# path for vcf file here
+vcf_mpileup_file = opts[:vcf_file]
+
+
+# The populate_strains method populates the strains in the db. It uses the strain names in array.
+populate_strains(strains)
+
+# The populate_features_and_annotations method populates the features and annotations.  It uses the embl/gbk file.
+populate_features_and_annotations(sequence_flatfile)
+
+#The populate_snps_alleles_genotypes method populates the snps, alleles and genotypes.  It uses the strain names (array) and vcf file.
+populate_snps_alleles_genotypes(strains, vcf_mpileup_file, opts[:cuttoff_snp].to_i, opts[:cuttoff_genotype].to_i)

data/lib/{snp-search.rb → snp-search/snp-search.rb}

@@ -1,7 +1,19 @@
 require 'rubygems'
 require 'bio'
 require  'snp_db_models'
-establish_connection
+#establish_connection
+
+def guess_sequence_format(reference_genome)
+  file_extension = File.extname(reference_genome).downcase
+  file_format = nil
+  case file_extension
+  when ".gbk", ".genbank", ".gb"
+    file_format = :genbank
+  when ".embl", ".emb"
+    file_format = :embl
+  end
+  return file_format
+end
 
 # A method to populate the strain names in the Strain table.  strain_names is an array of strain names.
 def populate_strains(strain_names)
@@ -16,8 +28,8 @@ end
 # We include all features that are not 'source' or 'gene' as they are repetitive info.  'CDS' is the gene.
 # The annotation table includes also the start and end coordinates of the CDS.  The strand is also included.  the 'locations' method is defined in bioruby under genbank.  It must be required at the top (bio).
 # Also, the qualifier and value are extracted from the embl file and added to the database.
-def populate_features_and_annotations(embl_ncbi_file)
-	embl_ncbi_file.features.each do |feature|
+def populate_features_and_annotations(sequence_file)
+	sequence_file.features.each do |feature|
 		unless feature.feature == "source" || feature.feature == "gene"
 			db_feature = Feature.new
 			db_feature.start = feature.locations.first.from
@@ -43,7 +55,7 @@ end
 #This method populates the rest of the information, i.e. SNP information, Alleles and Genotypes.
 # It requires the strain_names as array and the output (vcf file) from mpileup-snp identification algorithm.
 
-def populate_snps_alleles_genotypes(strain_names, vcf_file)
+def populate_snps_alleles_genotypes(strain_names, vcf_file, cuttoff_snp, cuttoff_genotype)
 	strains = Array.new
  	strain_names.each do |strain_name|
   		strain = Strain.find_by_name(strain_name) # equivalent to Strain.find.where("strains.name=?", strain_name).first
@@ -84,13 +96,13 @@ def populate_snps_alleles_genotypes(strain_names, vcf_file)
 		    genotypes.each_with_index do |gt, index|
 		    	if gt == "1/1"
 			        variant_genotypes << index
-			        if genotypes_qualities[index].to_i >= 30
+			        if genotypes_qualities[index].to_i >= cuttoff_genotype
 			          high_quality_variant_genotypes << index
 			        end
 		    	end
 			end
 		    
-			if snp_qual.to_i >= 90 && genotypes.include?("1/1") &&  ! high_quality_variant_genotypes.empty? && high_quality_variant_genotypes.size == variant_genotypes.size # first condition checks the overall quality of the SNP is >=90, second checks that at least one genome has the 'homozygous' 1/1 variant type with quality >= 30 and informative SNP
+			if snp_qual.to_i >= cuttoff_snp && genotypes.include?("1/1") &&  ! high_quality_variant_genotypes.empty? && high_quality_variant_genotypes.size == variant_genotypes.size # first condition checks the overall quality of the SNP is >=90, second checks that at least one genome has the 'homozygous' 1/1 variant type with quality >= 30 and informative SNP
 		    	if  genotypes.include?("0/0") && !genotypes.include?("0/1") # exclude SNPs which are all 1/1 i.e something strange about ref and those which have confusing heterozygote 0/1s
 			        good_snps +=1
 			        # puts good_snps

data/lib/{snp_db_connection.rb → snp-search/snp_db_connection.rb}

@@ -1,5 +1,5 @@
 require 'active_record'
-def establish_connection(db_location= "snp.db.sqlite3")
+def establish_connection(db_location)
   ActiveRecord::Base.establish_connection(
     :adapter => "sqlite3",
     :database => db_location,

data/lib/{snp_db_schema.rb → snp-search/snp_db_schema.rb}

@@ -1,6 +1,4 @@
-require 'snp_db_connection'
-establish_connection
-
+def db_schema
 ActiveRecord::Schema.define do
   unless table_exists? :strains
     create_table :strains do |t|
@@ -83,4 +81,5 @@ ActiveRecord::Schema.define do
   unless index_exists? :annotations, :feature_id
     add_index :annotations, :feature_id
   end
+end
 end

data/snp-search.gemspec

@@ -5,13 +5,14 @@
 
 Gem::Specification.new do |s|
   s.name = "snp-search"
-  s.version = "0.2.0"
+  s.version = "0.3.0"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.authors = ["Ali Al-Shahib", "Anthony Underwood"]
-  s.date = "2011-11-25"
+  s.date = "2011-11-30"
   s.description = "Use the snp-search toolset to query the SNP database"
   s.email = "ali.al-shahib@hpa.org.uk"
+  s.executables = ["snp-search"]
   s.extra_rdoc_files = [
     "LICENSE.txt",
     "README",
@@ -26,12 +27,13 @@ Gem::Specification.new do |s|
     "README.rdoc",
     "Rakefile",
     "VERSION",
-    "lib/snp-search.rb",
-    "lib/snp_db_connection.rb",
-    "lib/snp_db_models.rb",
-    "lib/snp_db_schema.rb",
-    "lib/user_entry_file.rb",
+    "bin/snp-search",
+    "lib/snp-search/snp-search.rb",
+    "lib/snp-search/snp_db_connection.rb",
+    "lib/snp-search/snp_db_models.rb",
+    "lib/snp-search/snp_db_schema.rb",
     "snp-search.gemspec",
+    "snp-search_test.rb",
     "spec/snp-search_spec.rb",
     "spec/spec_helper.rb"
   ]
@@ -45,17 +47,26 @@ Gem::Specification.new do |s|
     s.specification_version = 3
 
     if Gem::Version.new(Gem::VERSION) >= Gem::Version.new('1.2.0') then
+      s.add_runtime_dependency(%q<activerecord>, [">= 0"])
+      s.add_runtime_dependency(%q<bio>, [">= 0"])
+      s.add_runtime_dependency(%q<slop>, [">= 0"])
       s.add_development_dependency(%q<rspec>, ["~> 2.3.0"])
       s.add_development_dependency(%q<bundler>, ["~> 1.0.0"])
       s.add_development_dependency(%q<jeweler>, ["~> 1.6.4"])
       s.add_development_dependency(%q<rcov>, [">= 0"])
     else
+      s.add_dependency(%q<activerecord>, [">= 0"])
+      s.add_dependency(%q<bio>, [">= 0"])
+      s.add_dependency(%q<slop>, [">= 0"])
       s.add_dependency(%q<rspec>, ["~> 2.3.0"])
       s.add_dependency(%q<bundler>, ["~> 1.0.0"])
       s.add_dependency(%q<jeweler>, ["~> 1.6.4"])
       s.add_dependency(%q<rcov>, [">= 0"])
     end
   else
+    s.add_dependency(%q<activerecord>, [">= 0"])
+    s.add_dependency(%q<bio>, [">= 0"])
+    s.add_dependency(%q<slop>, [">= 0"])
     s.add_dependency(%q<rspec>, ["~> 2.3.0"])
     s.add_dependency(%q<bundler>, ["~> 1.0.0"])
     s.add_dependency(%q<jeweler>, ["~> 1.6.4"])

data/snp-search_test.rb

@@ -0,0 +1,207 @@
+require 'rubygems'
+require 'bio'
+require  'snp_db_models'
+require 'snp_db_connection'
+require 'snp_db_models'
+require 'snp_db_schema'
+require 'highline/import'
+require 'pp'
+require 'slop'
+
+
+opts = Slop.new :help do
+  banner "ruby snp-search.rb [OPTIONS]"
+
+  on :v, :version, 'Display the version of App' do
+    puts 'Version 1.5!'
+    exit
+  end
+
+  on :V, :verbose, 'Enable verbose mode'
+  on :n, :name=, 'Name of database', true 
+  on :r, :reference=, 'Path for the reference database', true 
+  on :f, :vcf=, 'Path for the .vcf file', true
+  on :s, :strain, 'Path for the list of strains text file', true
+  on :c, :cuttoff, 'cuttoff for SNP quality', true 
+  on_empty do
+    puts help
+  end
+end
+
+opts.parse
+
+strains = []
+File.read(opts[:strain]).each_line do |line|
+  strains << line.chop
+end
+
+
+# A method to populate the strain names in the Strain table.  strain_names is an array of strain names.
+def populate_strains(strain_names)
+	strain_names.each do |strain|
+	s = Strain.new
+	s.name = strain
+	s.save
+	end
+end
+
+# A method to populate the database with the features (genes etc) and the annotations from the embl file.  
+# We include all features that are not 'source' or 'gene' as they are repetitive info.  'CDS' is the gene.
+# The annotation table includes also the start and end coordinates of the CDS.  The strand is also included.  the 'locations' method is defined in bioruby under genbank.  It must be required at the top (bio).
+# Also, the qualifier and value are extracted from the embl file and added to the database.
+def populate_features_and_annotations(embl_ncbi_file)
+	embl_ncbi_file.features.each do |feature|
+		unless feature.feature == "source" || feature.feature == "gene"
+			db_feature = Feature.new
+			db_feature.start = feature.locations.first.from
+			db_feature.end = feature.locations.first.to
+			db_feature.strand = feature.locations.first.strand
+			db_feature.name = feature.feature
+			db_feature.save
+			puts "populated #{db_feature.name}, start: #{db_feature.start}, end: #{db_feature.end}, strand: #{db_feature.strand} for feature: #{db_feature.id}"
+			# Populate the Annotation table with qualifier information from the genbank file
+			feature.qualifiers.each do |qualifier|
+				a = Annotation.new
+				a.qualifier = qualifier.qualifier
+				a.value = qualifier.value
+				a.save
+				db_feature.annotations << a
+				puts "populated #{a.qualifier} for feature: #{db_feature.id}"
+			end
+		end
+	end
+end
+
+
+#This method populates the rest of the information, i.e. SNP information, Alleles and Genotypes.
+# It requires the strain_names as array and the output (vcf file) from mpileup-snp identification algorithm.
+
+def populate_snps_alleles_genotypes(strain_names, vcf_file)
+	strains = Array.new
+ 	strain_names.each do |strain_name|
+  		strain = Strain.find_by_name(strain_name) # equivalent to Strain.find.where("strains.name=?", strain_name).first
+  		strains << strain
+	end
+
+# open vcf file and parse each line
+	File.open(vcf_file) do |f|
+	  # header names
+	  header = f.gets
+	  header2 = f.gets.chomp
+	  column_headings = header2.split("\t")
+	  sample_names = column_headings[9..-1]
+	 
+	  good_snps = 0
+	  # start parsing snps
+		while line = f.gets
+		    details = line.split("\t")
+		    ref = details[0]
+		    ref_pos = details[1]
+		    ref_base = details[3]
+		    snp_base = details[4]
+		    snp_qual = details [5]
+		    samples = details[9..-1]
+		     
+		    genotypes = samples.map do |s| 
+		      pl, gt, gq = s.chomp.split(":")
+		      gt
+		    end
+
+		    genotypes_qualities = samples.map do |s| 
+		      pl, gt, gq = s.chomp.split(":")
+		      gq
+		    end
+
+		    high_quality_variant_genotypes = Array.new # this will be filled with the indicies of genotypes that are "1/1" and have a quality >= 30
+		    variant_genotypes = Array.new
+		    genotypes.each_with_index do |gt, index|
+		    	if gt == "1/1"
+			        variant_genotypes << index
+			        if genotypes_qualities[index].to_i >= 30
+			          high_quality_variant_genotypes << index
+			        end
+		    	end
+			end
+		    
+			if snp_qual.to_i >= opts[:cuttoff] && genotypes.include?("1/1") &&  ! high_quality_variant_genotypes.empty? && high_quality_variant_genotypes.size == variant_genotypes.size # first condition checks the overall quality of the SNP is >=90, second checks that at least one genome has the 'homozygous' 1/1 variant type with quality >= 30 and informative SNP
+		    	if  genotypes.include?("0/0") && !genotypes.include?("0/1") # exclude SNPs which are all 1/1 i.e something strange about ref and those which have confusing heterozygote 0/1s
+			        good_snps +=1
+			        # puts good_snps
+			        #create snp
+			        s = Snp.new
+			        s.ref_pos = ref_pos
+			        s.save
+			        puts "Adding Reference SNP position: #{ref_pos}"
+
+			        # create ref allele
+			        ref_allele = Allele.new
+			        ref_allele.base = ref_base
+			        ref_allele.snp = s
+			        ref_allele.save
+
+			        puts "Adding Reference SNP base: #{ref_base}"
+
+			        s.reference_allele = ref_allele
+			        s.save
+
+			        # create snp allele
+			        snp_allele = Allele.new
+			        snp_allele.base = snp_base
+			        snp_allele.snp = s
+			        snp_allele.save
+
+			        puts "Adding SNP base: #{snp_base}"
+
+			        
+
+			        genotypes.each_with_index do |gt, index|
+			          genotype = Genotype.new
+			          genotype.strain = strains[index]
+			          puts index if strains[index].nil?
+			          # print "#{gt}(#{genotypes_qualities[index]}) "
+			          if gt == "0/0" # wild type
+			            genotype.allele = ref_allele
+			          elsif gt == "1/1" # snp type
+			            genotype.allele = snp_allele
+			          else
+			            puts "Strange SNP #{gt}"
+			          end
+			          genotype.save
+			        end
+		      	end
+		    end
+
+		end
+	end
+	#Here we link the features to snps.
+	Snp.all.each do |snp|
+		x = Feature.where("features.start <= ? AND features.end >= ?", snp.ref_pos, snp.ref_pos).first
+		snp.feature = x
+		snp.save
+	end
+end
+
+
+#puts opts[:name]
+# Enter the name of your database 
+establish_connection(opts[:name])
+
+# # Schema will run here
+db_schema
+
+# path for embl file here
+#path_for_embl_file = ask("Please enter the full path for the embl reference file")
+genome_sequence = Bio::FlatFile.open(Bio::EMBL,opts[:reference]).next_entry
+
+# # path for vcf file here
+vcf_mpileup_file = opts[:vcf]
+
+
+# # The populate_strains method populates the strains in the db. It uses the strain names in array.
+populate_strains(strains)
+
+# # The populate_features_and_annotations method populates the features and annotations.  It uses the embl/gbk file.
+populate_features_and_annotations(genome_sequence)
+
+# #The populate_snps_alleles_genotypes method populates the snps, alleles and genotypes.  It uses the strain names (array) and vcf file.
+populate_snps_alleles_genotypes(strains, vcf_mpileup_file)

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: snp-search
 version: !ruby/object:Gem::Version
-  version: 0.2.0
+  version: 0.3.0
   prerelease: 
 platform: ruby
 authors:
@@ -10,11 +10,44 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2011-11-25 00:00:00.000000000Z
+date: 2011-11-30 00:00:00.000000000Z
 dependencies:
+- !ruby/object:Gem::Dependency
+  name: activerecord
+  requirement: &2159603360 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: *2159603360
+- !ruby/object:Gem::Dependency
+  name: bio
+  requirement: &2159602700 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: *2159602700
+- !ruby/object:Gem::Dependency
+  name: slop
+  requirement: &2159602100 !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: *2159602100
 - !ruby/object:Gem::Dependency
   name: rspec
-  requirement: &2161702220 !ruby/object:Gem::Requirement
+  requirement: &2159601440 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -22,10 +55,10 @@ dependencies:
         version: 2.3.0
   type: :development
   prerelease: false
-  version_requirements: *2161702220
+  version_requirements: *2159601440
 - !ruby/object:Gem::Dependency
   name: bundler
-  requirement: &2161701620 !ruby/object:Gem::Requirement
+  requirement: &2159600840 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -33,10 +66,10 @@ dependencies:
         version: 1.0.0
   type: :development
   prerelease: false
-  version_requirements: *2161701620
+  version_requirements: *2159600840
 - !ruby/object:Gem::Dependency
   name: jeweler
-  requirement: &2161701080 !ruby/object:Gem::Requirement
+  requirement: &2159600240 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ~>
@@ -44,10 +77,10 @@ dependencies:
         version: 1.6.4
   type: :development
   prerelease: false
-  version_requirements: *2161701080
+  version_requirements: *2159600240
 - !ruby/object:Gem::Dependency
   name: rcov
-  requirement: &2161699820 !ruby/object:Gem::Requirement
+  requirement: &2159599640 !ruby/object:Gem::Requirement
     none: false
     requirements:
     - - ! '>='
@@ -55,10 +88,11 @@ dependencies:
         version: '0'
   type: :development
   prerelease: false
-  version_requirements: *2161699820
+  version_requirements: *2159599640
 description: Use the snp-search toolset to query the SNP database
 email: ali.al-shahib@hpa.org.uk
-executables: []
+executables:
+- snp-search
 extensions: []
 extra_rdoc_files:
 - LICENSE.txt
@@ -73,12 +107,13 @@ files:
 - README.rdoc
 - Rakefile
 - VERSION
-- lib/snp-search.rb
-- lib/snp_db_connection.rb
-- lib/snp_db_models.rb
-- lib/snp_db_schema.rb
-- lib/user_entry_file.rb
+- bin/snp-search
+- lib/snp-search/snp-search.rb
+- lib/snp-search/snp_db_connection.rb
+- lib/snp-search/snp_db_models.rb
+- lib/snp-search/snp_db_schema.rb
 - snp-search.gemspec
+- snp-search_test.rb
 - spec/snp-search_spec.rb
 - spec/spec_helper.rb
 homepage: http://github.com/hpa-bioinformatics/snp-search
@@ -96,7 +131,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: 3531722145808918413
+      hash: -4259420816018168147
 required_rubygems_version: !ruby/object:Gem::Requirement
   none: false
   requirements:

data/lib/user_entry_file.rb

@@ -1,21 +0,0 @@
-require 'snp-search'
-
-#path for embl file here
-genome_sequence = Bio::FlatFile.open(Bio::EMBL, "path_for_embl_file_here").next_entry
-
-#path for vcf file here
-vcf_mpileup_file = "path_for_vcf_file_here"
-
-#array of strain names here
-strains = ["STRAIN_NAME_1", "STRAIN_NAME_2"]
-
-# Thats it, you job is done here.
-
-# The populate_strains method populates the strains in the db. It uses the strain names in array.
-populate_strains(strains)
-
-# The populate_features_and_annotations method populates the features and annotations.  It uses the embl/gbk file.
-populate_features_and_annotations(genome_sequence)
-
-# The populate_snps_alleles_genotypes method populates the snps, alleles and genotypes.  It uses the strain names (array) and vcf file.
-populate_snps_alleles_genotypes(strains, vcf_mpileup_file)