Bowtie Assembly Parameters
1
0
Entering edit mode
11.0 years ago
biolab ★ 1.4k

Hi, everyone, I am trying Bowtie for reference based assembly. My draft genome dataset is not from a model organism, but a related species. This data contains ~20 million 88 bp illumina reads. Because it's not resequencing of reference genome, I need to use a relatively loose parameters. Actually I only focus on genes without caring about intergenic regions. A CDS dataset of model organism can be used as a reference in this assembly. What Bowtie parameters you suggest to use? Thank you very much.

assembly • 4.2k views
ADD COMMENT
1
Entering edit mode

Bowtie is an aligner, not an assembler. Why don't you just use a splice-aware aligner, like tophat, or just try de novo assembly (e.g., with trinity)?

ADD REPLY
0
Entering edit mode

Can you actually use trinity for genome assembly? I thought it's specifically for transcript assembly.

ADD REPLY
0
Entering edit mode

Since biolab mentioned aligning to CDS, I'm assuming that he/she is actually doing RNAseq (though, upon rereading, this may well be incorrect!). Otherwise, yes, I think you're correct.

ADD REPLY
1
Entering edit mode
11.0 years ago
Adrian Pelin ★ 2.6k

I concur with the first comment, specifically I suggest an assembler. Your data is 88bp illumina, is it paired end? That will improve your assembly statistics. Here are my suggestions:

1) Try SPAdes or Velvet, read here about how they differ. 2) You can map your resulting contigs to your close species. This is better, since contigs are much much longer than reads, and you can virtually see which parts of the genome remain homologous.

ADD COMMENT
0
Entering edit mode

thanks for ur suggestions.

ADD REPLY
0
Entering edit mode

It depends on how similar your reference is. bwa is quite good for losely mapping, or "last" given you have enough computational resource. Then just parse the outputs incorporating snps/indels.

ADD REPLY
0
Entering edit mode

thanks for your suggestions, all helpful!

ADD REPLY

Login before adding your answer.

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