BRAKER miss some genes despite of RNA-Seq evidence
Hi,
I ran braker.pl --UTR=on --genome=$1 --bam=$2 --species=$3 --cores 20 --softmasking --gff3 but in few cases, BRAKER did not predict any genes where there is RNA-Seq alignment.
Is there a way to get still gene predictions for those missing areas?
Thank you in advance,
rna-seq
genome
assembly
braker
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
StringTie fails with “Fatal error: Exit code 1” when using GFF3 annotation for a non-model halophyte
written by Chamara Lakshitha •Hello everyone, I am analyzing RNA-seq data from a non-model halophyte plant. Because no reference genome is available, I aligned the reads to an available …
-
Low busco scores after annotating genome using same model as another species in genus
written by Wilber0x •I am annotating the genome of a non model species of monocot. I have successfully annotated the sister species using BRAKER3 with a lot of …
-
Falsely high busco score after genome annotation with braker
written by Wilber0x •I have a genome that I am annotating, and there are no reference assemblies or closely related species with genomes sequenced. I am very pleased …
-
Need help with providing RNA Seq data to run braker.pl in ETP mode.
written by Vijith •I have some queries before running `Braker` in `ETP mode`, where I will provide both RNA Seq and protein data. I am looking forward to …
-
Error in file braker.pl at line 5404. Failed to execute gmes_petap.pl
written by natt1621 •Hi, I'm trying to run Braker to predict genes, but I'm getting this error with gmes_petap.pl script. When I try to look out on 'errors' …
-
Centrally installed Braker2 error: AUGUSTUS_CONFIG_PATH/species is not writeable.
written by Timotheus •Hi, I'm trying to run Braker2 to annotate my genome with the help of RNA-seq data: ``` braker.pl \ --genome=$ASSEMBLY \ --bam=$RNA_BAM \ --gff3 \ …
-
Issues with BRAKER2 (v.2.1.6)--not producing "genemark_hintsfile.gff.prot.tmp"
written by vizzerraa •Hi, BRAKER2 (v.2.1.6) is not producing the "genemark_hintsfile.gff.prot.tmp" file after ProtHint finishes, what could be causing this issue? The full command was /data/selaginella/anaconda3/envs/braker2_conda_env/BRAKER-2.1.6/scripts/braker.pl --genome=/data/selaginella/DATA/annotations/braker_wallacei_test/a.lines.fa --bam=/data/selaginella/DATA/annotations/braker_wallacei_test/wallacei.sorted.bam …
-
braker2 output: gff file is missing and species folder is empty
written by praasu •I have to annotate novel diplonema genome. So that I chose to use braker. I have run command prothint to generate and prothint_augustus.gff used to …
-
gff3_sp_complement_annotations.pl parameter settings
written by RicI ran `gff3_sp_complement_annotations.pl --ref stringtie.gff3 --add augustus.hints_utr.gff3 --out gaas.gff3` and got the below new IDs: NbV1Ch01 transdecoder five_prime_UTR 99705 99705 . - . ID=nbis_NEW-five_prime_utr-4254;Parent=STRG2.t1 NbV1Ch01 …
-
Genome annotation with BRAKER: how to interpret the results?
written by svitlana.lukicheva •Hello, I am new to the genome annotation and I'm lost with the interpretation of results produced by BRAKER. - I have a de novo …
Rna-seq are not necessarily from protein coding genes. Have you seen any ORF that make you think that a gene is missing? Otherwise you should check. If there is no ORF Augustus will not be able to predict anything.
apart from what Juke-34 already mentioned, you have to keep in mind that all these tools (approaches) have the goal to correctly predict as much genes as possible. This of course does not mean that all genes will be correctly predicted (if at all) but that the majority of the genes should be correct (or as good as). A consequence of this can be that some genes will be perhaps missed if by that eg. a 100 others might be correct