Added iscomplete()

BioJulia · Dec 10, 2019 · 38153a4 · 38153a4
1 parent d70ea7f
commit 38153a4
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 1 deletion.
diff --git a/src/GenomicAnnotations.jl b/src/GenomicAnnotations.jl
@@ -6,7 +6,7 @@ using BioSequences
 using GZip
 
 export Chromosome, Gene, AbstractGene, GeneDataView, Locus
-export readgbk, sequence, iscomplement, addgene!, pushproperty!, printgbk
+export readgbk, sequence, iscomplement, iscomplete, addgene!, pushproperty!, printgbk
 export feature, index, locus
 export @genes
 

diff --git a/src/types.jl b/src/types.jl
@@ -278,6 +278,14 @@ Return `true` if `locus(gene).compliment == '-'`, otherwise return `false`.
 iscomplement(gene::AbstractGene) = locus(gene).strand == '-'
 
 
+"""
+    iscomplete(gene)
+
+Return `true` if `gene` is a complete gene, i.e. not a pseudo gene or partial.
+"""
+iscomplete(gene::AbstractGene) = !any(get(gene, :pseudo, false)) && !any(get(gene, :ribosomal_slippage, false)) && locus(gene).complete_right && locus(gene).complete_left
+
+
 function appendstring(field, v)
     s = ""
     if v isa Bool