Question: How to separate sequences of two species from Metagenomics Analysis?
0
gravatar for hafiz.talhamalik
26 days ago by
Pakistan
hafiz.talhamalik40 wrote:

I have sequenced a genomics data. This contains ecoli and naegleria fowleri, and I want to separate both of them so I could assemble naegleria fowleri genome. Any Idea how to separate them from fastq files??

ADD COMMENTlink modified 26 days ago • written 26 days ago by hafiz.talhamalik40
1

You could try profiling with kraken/centrifuge. I believe they will separate out the reads for you too (I could be wrong).

ADD REPLYlink written 26 days ago by jrj.healey10k

OK I will try that and will let you know.

ADD REPLYlink written 26 days ago by hafiz.talhamalik40
1

In this case you can use fastqscreen or bbbbsplit.sh (module available under BBMap package.

Both the tools will map your raw reads on multiple reference genomes.

NOTE: In this case you have to consider unaligned reads (failed to map on any of your reference genomes) for assembly construction.

ADD REPLYlink modified 24 days ago • written 24 days ago by Nitin Narwade380

oh great.! I think these will serve the purpose.

ADD REPLYlink written 23 days ago by hafiz.talhamalik40

Dear hafiz.talhamalik,

If my suggestion somehow helps you to get an idea for resolving your query, could you please hit on upvote button?

Thank you.

ADD REPLYlink written 23 days ago by Nitin Narwade380

Try mapping to E.coli first, check the coverage (they may share some genes), and assemble with unmapped reads. It may be hard to separate If they share some genes.

ADD REPLYlink written 26 days ago by shenwei3564.5k
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: 825 users visited in the last hour