This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Reformatting Fastq to Fasta for Metagenomic Analysis

Hi there, I'm trying to use VirFinder with shotgun metagenomic sequencing reads. My reads came in FASTQ format, and VirFinder requires a Fasta input format. I've read you can use tools like BBMap reformat.sh, but these don't map the reads. If I'm using metagenomic data, is this an appropriate way to reformat into a fasta? Do I need the reads to be mapped?

Additionally, what is the difference between unmapped and mapped reads. Thanks!

metagenomics fasta fastq
• 1,432 views • link
updated by colindaven colindaven • written by echolley • echolley

2 answers

Fasta and fastq are distinct formats, with the latter having information about the basecall Q-score values for each base.

My reads came in FASTQ format, and VirFinder requires a Fasta input format.

Then you will have to convert the reads into fasta format.

If I'm using metagenomic data, is this an appropriate way to reformat into a fasta?

There is only one way to convert a fastq format sequence read to fasta and reformat.sh will do that as you already note.

but these don't map the reads

VirFinder will likely do that.

Do I need the reads to be mapped?

Depends on what you are trying to do. If a software needs mapped reads then you will have to do that up front.

what is the difference between unmapped and mapped reads

Exactly what those designations say. Reads are either mapped or will not be mapped to a particular set of reference sequences used for the alignments.

• 0 views • link
written by GenoMax GenoMax

Another option to convert fastq to fasta formats is seqtk

seqtk seq -A 

Usage:   seqtk seq [options] <in.fq>|<in.fa>

         -A        force FASTA output (discard quality)
• 0 views • link
written by colindaven colindaven

Log in to answer this question.

More posts like this