@@ -671,6 +671,7 @@ def __init__(self,prog):
671
671
end = int (cols [4 ]) + 100
672
672
if end > ContigSizes .get (cols [2 ]): #check that it doesn't go over length
673
673
end = ContigSizes .get (cols [2 ])
674
+ print end
674
675
bedfile .write ('%s\t %i\t %i\t %s\n ' % (cols [2 ],start ,end ,cols [0 ]))
675
676
#now get BUSCO GFF models
676
677
busco_augustus_tmp = os .path .join (args .out , 'predict_misc' , 'busco_augustus.tmp' )
@@ -719,15 +720,15 @@ def __init__(self,prog):
719
720
Busco_Weights = os .path .abspath (busco_weights )
720
721
EVM_out = os .path .abspath (EVM_out )
721
722
Busco_Predictions = os .path .abspath (busco_predictions )
722
- #parse entire EVM command to script
723
+ #parse entire EVM command to script, must be absolute paths for everything
723
724
if Exonerate and Transcripts :
724
- evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--protein_alignments' , busco_proteins , '--transcript_alignments' , busco_transcripts , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_out ]
725
+ evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--protein_alignments' , os . path . abspath ( busco_proteins ) , '--transcript_alignments' , os . path . abspath ( busco_transcripts ) , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
725
726
elif not Exonerate and Transcripts :
726
- evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ),str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--transcript_alignments' , busco_transcripts , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_out ]
727
+ evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ),str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--transcript_alignments' , os . path . abspath ( busco_transcripts ) , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
727
728
elif not Transcripts and Exonerate :
728
- evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--protein_alignments' , busco_proteins , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_out ]
729
+ evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--protein_alignments' , os . path . abspath ( busco_proteins ) , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
729
730
elif not any ([Transcripts ,Exonerate ]):
730
- evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_out ]
731
+ evm_cmd = [sys .executable , EVM_script , os .path .join (args .out , 'logfiles' , 'funannotate-EVM_busco.log' ), str (args .cpus ), '--genome' , MaskGenome , '--gene_predictions' , Busco_Predictions , '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
731
732
#run EVM
732
733
if not os .path .isfile (EVM_out ):
733
734
subprocess .call (evm_cmd )
0 commit comments