MAKER2 no FASTAs generated · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

MAKER2 no FASTAs generated

I have recently used MAKER2 for genome annotation. The run appears to have been successful, but when I came to use fasta_merge -d there was no output! gff3_merge -d worked fine, producing a new file AGND01.1.all.gff. I decided to check if any FASTAs had been generated using find . -type f -name "*.fa" in the AGND01.1_datastore directory generated by MAKER, but there were no results. Does this mean that MAKER did not produce and FASTAs for the contigs it annotated? If so, how do I fix this? Is there an option I must choose in the .ctl files?

Thanks,

Jamie

maker fasta fasta_merge genome annotation

• 2,200 views • link

updated 2.2 years ago by Kinoppy •

• written 6.2 years ago by jamie.pike

is your find cmdline accurate enough? I mean it's not the fasta files are called .FA or .fasta or such ('cus then your find will not report them)

• 0 views • link

written 6.2 years ago by lieven.sterck

Ah, that seems obvious now you say it!! I have just tried find . -type f -name "*.fasta" which has produced results. It still does not explain why fasta_merge -d AGND01.1_master_datastore_index.log is not working though. Would you perhaps know why?

• 0 views • link

written 6.2 years ago by jamie.pike

no it does not indeed :)

no errors/warnings or such, I assume?

• 0 views • link

written 6.2 years ago by lieven.sterck

Nope, no errors or warning at all. I have had a look at where the fastas are and they are in a subdirectory e.g. AGND01.1.maker.output/AGND01.1_datastore/1D/3B/AGND01000482.1/theVoid.AGND01000482.1/query.fasta

not where the gff outputs are e.g.

AGND01.1.maker.output/AGND01.1_datastore/1D/3B/AGND01000482.1/AGND01000482%2E1.gff

So I am still wondering if MAKER generated the correct fastas?

• 0 views • link

written 6.2 years ago by jamie.pike

1 answer

Hi Jamie, I have also encountered the same issue. Were you by any chance able to find a solution or understand why MAKER was giving that error?

Thank you

• 0 views • link

written 2.2 years ago by Kinoppy •

Log in to answer this question.

More posts like this

vg index -j hangs for weeks (minigraph-cactus pangenome, dist file incomplete or missing distance i…
written 6 months ago by RianeZ •

Hi everyone, I'm having trouble building the distance index (.dist) for a large plant pangenome graph generated by Minigraph-Cactus, and would appreciate any advice. I'm …
Sub-sampling a BAM to a fixed number of reads
written 3.2 years ago by matt81rd •

Hi all, I have recently been trying to write a simple python script that can produce contaminated synthetic paired end reads for a tool my …
Methods of annotating a genome?
written 6.3 years ago by jamie.pike

Hello, I have some fungal genomes that I downloaded from NCBI, it appears that 5 of the genomes have not been annotated. I have a …
Matching Query ID and Subject ID in a new file
written 6.3 years ago by jamie.pike

Hi, I am trying to write a script that takes a list of genomes, performs BLASTN, creates a file containing the sequence of the BLASTN …
Extracting upstream and downstream bases from a Blastn hit.
written 6.4 years ago by jamie.pike

I have the output from BLASTN searches and want to extract 2500 bases upstream and downstream of each BLASTN hit from an assembled genome. I …
News: Discovering the world of the Long-read nanopore Sequencing
written 6.7 years ago by Physalia-courses

![enter image description here][1] Nanopore sequencing technology is an emerging and a promising sequencing technology, which promises high sequencing throughput, low cost and longer read …
MAKER AED score
written 9.0 years ago by Tye Kahn •

Hi everyone, I recently have performed a **new genome annotation** with EVM, thanks to *ab initio* predictions, RNA-Seq data and homology protein alignments. As there …
bwa aln. Reads align to 'N' characters in reference file.
written 9.4 years ago by izaak.coleman1 •

Hi All, This question seems quite simple, but, it seems to be very difficult to get an answer. I have a file `cns_pairs.fastq` containing reads. …
Method for generating a reference-based genome assembly in fasta with no reference genome bases
written 10.8 years ago by memory_donk

Hi Biostars, I'm trying to produce a genome assembly of a microbat species which is fairly closely related to the previously published *Myotis lucifugus.* The …
Forum: Isaac Asimov Accurately Predicts Woes Of Bioinformatician
written 13.2 years ago by Obi Griffith

I recently re-read Isaac Asimov's excellent Foundation series. I decided to read them in the order they were released rather than chronological order of the …