bacterial-annotator 0.1.2 → 0.1.3
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/lib/bacterial-annotator/genbank-manip.rb +1 -0
- data/lib/bacterial-annotator.rb +3 -1
- metadata +1 -1
checksums.yaml
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
---
|
2
2
|
SHA1:
|
3
|
-
metadata.gz:
|
4
|
-
data.tar.gz:
|
3
|
+
metadata.gz: 746fa68fef3d2d9c1c6df6fef95c354a818e55e1
|
4
|
+
data.tar.gz: e0bbb9f28e87ef0daaf745e1eba9c1bf215d1f63
|
5
5
|
SHA512:
|
6
|
-
metadata.gz:
|
7
|
-
data.tar.gz:
|
6
|
+
metadata.gz: 0cd5f9a866b7f5bff37cd176ba941a11bf59fb51cba37f626da458b38ae9cb1acc03693848f0b1d7bbf276580f1b252c9ea3040cd4098ff61e3a6c4086a0ca1c
|
7
|
+
data.tar.gz: 35fe32dcca6c1a5f5316718fb48b5fa38fbb69069a06230595b1c024e654ba3319122183260dcc6beb6c93083c9c75a1784e845a0b2c6f94e18aab61c93ff918
|
@@ -146,6 +146,7 @@ class GenbankManip
|
|
146
146
|
ftArray.push(qProd)
|
147
147
|
end
|
148
148
|
|
149
|
+
# check if there is a reference genome.. reference_locus shouldn't be nil in that case
|
149
150
|
if locus != nil
|
150
151
|
qNote = Bio::Feature::Qualifier.new('note', "correspond to #{locus} locus (#{pId}% identity) from #{reference_locus.entry_id}")
|
151
152
|
ftArray.push(qNote)
|
data/lib/bacterial-annotator.rb
CHANGED
@@ -228,7 +228,9 @@ class BacterialAnnotator
|
|
228
228
|
@contig_annotations.each do |contig, contig_prot_annotations|
|
229
229
|
gbk_path = @fasta.prodigal_files[:gbk_path]
|
230
230
|
gbk_to_annotate = GenbankManip.new("#{gbk_path}/#{contig}.gbk", "#{gbk_path}")
|
231
|
-
|
231
|
+
reference_locus = nil
|
232
|
+
reference_locus = @refgenome.gbk.locus if @with_refence_genome
|
233
|
+
gbk_to_annotate.add_annotation contig_prot_annotations, gbk_path, 0, reference_locus
|
232
234
|
end
|
233
235
|
|
234
236
|
end # end of method
|