How to blast a paired end fastq file?
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.
• 7,278 views
•
link
updated
by
chen
•
written
by
chillchik29 •
•
written
by
chillchik29 •
2 answers
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
Use Magic-Blast
Log in to answer this question.
More posts like this
-
Best way to handle samples with both single end and paired end?
written by ichbinlynn93 •
-
MACS for a paired-end sequencing
written by KT •
-
Should I merge the Paired-end RNA-seq reads for BLAT for intron identification
written by praasu •
-
Paired end data with R1 demultiplexed but not R2
written by singlegenome •
-
Kallisto | bustools workflow for fastq with R1 and R2 (forward and reverse)?
written by Farah
-
FASTQ --> FASTA Conversion
written by Khosh.0 •
-
Trim Paired-end Fastq Files
written by yuabrahamliu •
-
hisat2 parameters for paired end and single end reads from same samples
written by bioinfo17 •
-
which tool is used for separating forward and reverse from paired end fastq
written by ranjini.shakthi •
-
BWA MEM alignment output of splited fastq files differ from the original(unsplit) fastq file
written by manojkumar_bhosale
Why do you wish to BLAST read sequences? Plus, why do you want to BLAST both pairs and all reads?
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.
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.Thank you! I'm going to try that out :)