@@ -338,7 +338,7 @@ def __init__(self,prog):
338
338
blat2hints = os .path .join (AUGUSTUS_BASE , 'scripts' , 'blat2hints.pl' )
339
339
subprocess .call ([blat2hints , b2h_input , b2h_output , '--minintronlen=20' , '--trunkSS' ], stdout = FNULL , stderr = FNULL )
340
340
else :
341
- lib .log .error ("No transcripts availble to generate Augustus hints, provide --transcript_evidence" )
341
+ lib .log .error ("No transcripts available to generate Augustus hints, provide --transcript_evidence" )
342
342
343
343
#check for protein evidence/format as needed
344
344
p2g_out = os .path .join (args .out , 'predict_misc' , 'exonerate.out' )
@@ -452,7 +452,7 @@ def __init__(self,prog):
452
452
shutil .rmtree (os .path .join (args .out , 'predict_misc' , 'braker' ))
453
453
os .rename ('braker' , os .path .join (args .out , 'predict_misc' , 'braker' ))
454
454
#okay, now need to fetch the Augustus GFF and Genemark GTF files
455
- aug_out = os .path .join (args .out , 'predict_misc' , 'braker' , aug_species , 'augustus.gff ' )
455
+ aug_out = os .path .join (args .out , 'predict_misc' , 'braker' , aug_species , 'augustus.gff3 ' )
456
456
gene_out = os .path .join (args .out , 'predict_misc' , 'braker' , aug_species , 'GeneMark-ET' , 'genemark.gtf' )
457
457
#now convert to EVM format
458
458
Augustus = os .path .join (args .out , 'predict_misc' , 'augustus.evm.gff3' )
@@ -759,7 +759,7 @@ def __init__(self,prog):
759
759
lib .log .error ("GeneMark predictions failed, proceeding with only Augustus" )
760
760
761
761
#if hints used for Augustus, get high quality models > 80% coverage to pass to EVM
762
- if os .path .isfile (hints_all ):
762
+ if os .path .isfile (hints_all ) and not args . rna_bam :
763
763
lib .log .info ("Pulling out high quality Augustus predictions" )
764
764
hiQ_models = []
765
765
with open (aug_out , 'rU' ) as augustus :
0 commit comments