lib/bio/PolyploidTools/ExonContainer.rb in bio-polyploid-tools-0.10.1 vs lib/bio/PolyploidTools/ExonContainer.rb in bio-polyploid-tools-1.0.0
- old
+ new
@@ -74,11 +74,10 @@
path = opts[:alig_path]
chromosomes[name] = Bio::DB::Fasta::FastaFile.new({:fasta=>path})
end
def add_snp(snp)
- #TODO: add to the snp the maximum number of hits?
snp.max_hits = self.max_hits
@snp_map[snp.gene] = Array.new unless @snp_map[snp.gene]
@snp_map[snp.gene] << snp
end
@@ -139,10 +138,11 @@
#puts snp.inspect
begin
file.puts snp.aligned_sequences_fasta
rescue Exception=>e
+ #puts snp.inspect
@missing_exons << snp.to_s
$stderr.puts "print_fasta_snp_exones:" + snp.to_s + ":" + e.to_s
$stderr.puts "Local position: #{snp.local_position}"
$stderr.puts "Local position: #{snp.parental_sequences.to_s}"
$stderr.puts e.backtrace
@@ -158,11 +158,11 @@
snp_array.each do |snp|
string = ""
begin
primer_3_min_seq_length
string = snp.primer_3_string( snp.chromosome, parental )
- #TODO: add tan error to the SNP this snp has more than max_hits. Or maybe inside the SNP file.
- #puts "print_primer_3_exons: #{string.size}"
+ #TODO: add tan error to the SNP this snp has more than max_hits.
+ #Or maybe inside the SNP file.
if string.size > 0
file.puts string
added += 1
end
rescue Exception=>e
\ No newline at end of file