Gaius-Augustus / BRAKER

BRAKER is a pipeline for fully automated prediction of protein coding gene structures with GeneMark-ES/ET/EP/ETP and AUGUSTUS in novel eukaryotic genomes
Other
334 stars 80 forks source link

No BUSCO score of BRAKER3 prediction result with protein data only #771

Open phhsieh001 opened 4 months ago

phhsieh001 commented 4 months ago

Dear BRAKER Team,

We ran BRAKER(ver. 3.0.7 docker) with rna-seq only and protein only separately at our termite genome. Here are the BUSCO score at inputs and outputs. <html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns="http://www.w3.org/TR/REC-html40">

  | BUSCO(insecta_odb10) -- | -- Genome | C:99.6%[S:97.3%,D:2.3%] Result of RNA-seq only annotated | C:94.9%[S:74.0%,D:20.9%] Result of protein only annotated | C:0.0%[S:0.0%,D:0.0%] the protein evidenve used in BRAKER | C:99.9%[S:15.4%,D:84.5%]

What are the possible reasons for causing the zero BUSCO in protein only annotation? Thank you!

KatharinaHoff commented 4 months ago

This looks like something went wrong.

Are you sure the file that you applied BUSCO to actually contained any protein sequences at all? (I am currently assuming that you are running BUSCO outside of BRAKER, not using the internal built-in compleasm.)

phhsieh001 commented 4 months ago

Yes, we ran BUSCO outside of BRAKER. We used docker to run BUSCO in transcript mode with BRAKER result (braker.codingseq) as input. There are 34,060 sequences in braker.codingseq file and the BUSCO: 0% (Ran BRAKER with only protein evidence). There are 71,027 sequences in braker.codingseq file and the BUSCO: 94.9% (Ran BRAKER with only RNA-seq evidence). These results are both appeard in BRAKER2 and BRAKER3.