What Can We Do For The Reads That Are Not Mapped To The Genome?
2
0
Entering edit mode
10.1 years ago
lanyuhao1991 ▴ 20

I have the reads of the some mammals intestinal microflora. I mapped the reads to the bacteria genomes. And I want to use the rest of reads that are not mapped to do something. I want build some new genomes of bacterium. Can you give me some advice for: 1.how to build some new genomes of bacterium? 2.can I use the rest of the reads to do other things? Thanks!

reads map • 2.2k views
ADD COMMENT
0
Entering edit mode

1.Do some assembly。 2.Map to nt/nt database for classfication.

ADD REPLY
2
Entering edit mode
10.1 years ago

You could use an assembler like IDBA-UD, Spades or Ray (and there are others as well of course) to do de novo assembly in order to "build new genomes" as you write. You may end up with both old and new bacterial contigs, host genome (yak) contigs and perhaps some viruses. Of course you could also try to map your unmapped reads to the host genome to remove those from the analysis, and then map the remaining reads directly to a viral reference.

ADD COMMENT
0
Entering edit mode
10.1 years ago
Haluk ▴ 190

After assembling your unmapped reads, you can BLAST them to reference genome. It increases mapping rate.

ADD COMMENT

Login before adding your answer.

Traffic: 1323 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