Question: paired end metagenomics analysis
0
gravatar for saraswati.awasthi22
21 days ago by
saraswati.awasthi220 wrote:

Hello, I am Saraswati and I am new in the field of metagenomics. I have to do taxonomic classification(Archea, Bacteria, Eukaryotes and Viruses) of the shotgun sequence which is 3.5 GB in size by using tools DIAMOND and MEGAN6.

I downloaded the .sra file from ncbi and splitted into two files, forward and reverse .fastq format using sra toolkit. Now I have to do blastx of the sequence using DIAMOND but I want to know if I should merge forward and reverse reads into a single file or I should do blastx with forward read only? If I should merge then which tool I should use for that purpose?

I want to do blastx against nr database I tried doing it with the forward end .fasta file but it is taking a long time(more than 7 days and still continuing) so how should I make it faster? Also if there are some other tools or strategy to the same for taxonomy in short time and better way then please suggest me.

ADD COMMENTlink modified 13 days ago by Biostar ♦♦ 20 • written 21 days ago by saraswati.awasthi220

Yes, you should assemble paired-end reads (assuming they overlap) with any paired-end assembler tool, such as PEAR (this is just one among many tools available - paper) before trying to annotate them. There are software written to annotate taxonomically metagenomes such as kraken2 (site) and centrifuge (site) among others. You can, of course depending on your objective, try to assemble the metagenomes into genomes (the so called MAGs) annotating this after.

If you're blasting locally, you can provide -num_threads (depending on the threads that your computer/server has) to parallelize the work. If you're blasting remotely this is not possible as far as I know. Blasting entire fastq files without performing any kind of clustering etc, it will take some time, assuming that each file just have a few million of reads.

ADD REPLYlink written 21 days ago by antonioggsousa1.9k

IDseq is also a good option for taxonomic assignment.

ADD REPLYlink written 13 days ago by Shaminur100
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1927 users visited in the last hour
_