Reformatting Fastq to Fasta for Metagenomic Analysis
2
0
Entering edit mode
5 weeks ago
echolley ▴ 20

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 • 332 views
ADD COMMENT
2
Entering edit mode
5 weeks ago
GenoMax 144k

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.

ADD COMMENT
0
Entering edit mode
26 days ago

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)
ADD COMMENT

Login before adding your answer.

Traffic: 2689 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6