Biostar Beta. Not for public use.
How to blast a paired end fastq file?
0
Entering edit mode
12 months ago

Hello everyone, I have a paired end fastq file and I know that BLAST+ in command line, accepts fasta format. But I don't know how does it work for a paired end fastq file (I mean in two different files R1 and R2). Do I have to interleave them and then convert them to fasta? Or should I concatenate them and convert them to fasta? Or should I treat them as separate files and then convert them to fasta?

Thank you, any help would be appreciated.

ADD COMMENTlink
1
Entering edit mode

Why do you wish to BLAST read sequences? Plus, why do you want to BLAST both pairs and all reads?

ADD REPLYlink
0
Entering edit mode

I already selected these group of reads from all the metagenome sequences by aligning them with bbmap, to the phage_nt db. So before assembling them, I wanted to have a first approach of what I had in the sample by aligning them with BLAST and then export it to MEGAN.

ADD REPLYlink
1
Entering edit mode

You could save some effort and only convert one end to fasta (reformat.sh in=R1.fq.gz out=R1.fa) and then BLAST. PE reads are not going to give you much additional information that MEGAN can use.

ADD REPLYlink
0
Entering edit mode

Thank you! I'm going to try that out :)

ADD REPLYlink
0
Entering edit mode
11 months ago
JC 7.9k
Mexico

Use Magic-Blast

ADD COMMENTlink
0
Entering edit mode

While this would work to do the alignments it would produce SAM format output. I am not sure if MEGAN can use that.

ADD REPLYlink
1
Entering edit mode

magic-blast also reports aligment in tabular format, which is similar to blast outfmt 6

ADD REPLYlink
0
Entering edit mode

Ah I missed that. It may need some manipulation (since it is not quite outfmt 6). I guess OP can let us know.

ADD REPLYlink
0
Entering edit mode
12 months ago
chen ♦ 1.9k
OpenGene

You can use fastp to merge it first, then blast the merged reads

See the intro here: https://github.com/OpenGene/fastp#merge-paired-end-reads

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1