@@ -713,21 +713,21 @@ def __init__(self,prog):
713713 if Transcripts :
714714 output .write ("TRANSCRIPT\t genome\t 1\n " )
715715 #setup EVM run
716- EVM_out = os .path .join (args .out , 'predict_misc' , 'busco.evm.gff3' )
716+ EVM_busco = os .path .join (args .out , 'predict_misc' , 'busco.evm.gff3' )
717717 EVM_script = os .path .join (parentdir , 'bin' , 'funannotate-runEVM.py' )
718718 #get absolute paths for everything
719719 Busco_Weights = os .path .abspath (busco_weights )
720- EVM_out = os .path .abspath (EVM_out )
720+ EVM_busco = os .path .abspath (EVM_out )
721721 Busco_Predictions = os .path .abspath (busco_predictions )
722722 #parse entire EVM command to script, must be absolute paths for everything
723723 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' , 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 ) ]
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' , os .path .abspath (busco_proteins ), '--transcript_alignments' , os .path .abspath (busco_transcripts ), '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_busco ]
725725 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' , os .path .abspath (busco_transcripts ), '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
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' , os .path .abspath (busco_transcripts ), '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_busco ]
727727 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' , os .path .abspath (busco_proteins ), '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), os . path . abspath ( EVM_out ) ]
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' , os .path .abspath (busco_proteins ), '--weights' , Busco_Weights , '--min_intron_length' , str (args .min_intronlen ), EVM_busco ]
729729 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 ), os . path . abspath ( EVM_out ) ]
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_busco ]
731731 #run EVM
732732 if not os .path .isfile (EVM_out ):
733733 subprocess .call (evm_cmd )
0 commit comments