Question: Close Bacterial Draft Genome With Other Draft Genome
gravatar for biotech
6.2 years ago by
United States
biotech540 wrote:


I would like to reduce contig number of the bacterial draft genome I'm trying to close. The draft genome is in 47 contigs and has been assembled by me, using Phred, Phrap, Consed.

I'm thinking in using the other draft genome that has been published (same bacterial specie) perform this task.

Do you think I should start from zero, that is, join the not assembled raw data of the two drafts published and see what happens, or use a tool having as input the two assembled drafts?

Thanks for you help, Bernardo

assembly genomics • 2.7k views
ADD COMMENTlink modified 5.7 years ago by mesabioinformatics30 • written 6.2 years ago by biotech540

If you already have a good reference, why not directly use reference based assembly? If you want to identify variations, it might not be a good idea to assembly them together or merge them. It can also be difficult to say which sample your data comes from.

ADD REPLYlink modified 6.2 years ago • written 6.2 years ago by Xingyu Yang260

Do you know any good tool to do reference based assembly?

ADD REPLYlink written 6.2 years ago by biotech540

Sorry for late reply. Here's the pipeline I usually use.

  1. Alignment with bowtie2
  2. samtools mpileup -uf <Fasta of Reference> <Alignment Results>.bam | bcftools view -cg - | vcf2fq > <Assembled Results>.cons.fq
ADD REPLYlink modified 8 months ago by RamRS27k • written 6.2 years ago by Xingyu Yang260

It is unlikely that the other draft genome has any information you can use for your draft genome. I would reckon the assembly will be broken at the same places ie. large repeats such as rRNA locii and insertion sequences that your read length can not span.

ADD REPLYlink written 6.1 years ago by Torst950
gravatar for mesabioinformatics
5.7 years ago by
United States
mesabioinformatics30 wrote:

Try: Ragout—a reference-assisted assembly tool for bacterial genomes


Hope it helps !

ADD COMMENTlink modified 5.7 years ago by Istvan Albert ♦♦ 84k • written 5.7 years ago by mesabioinformatics30
gravatar for Chrispin Chaguza
5.7 years ago by
Wellcome Sanger Institute
Chrispin Chaguza250 wrote:

If you have contigs of your draft genome, you can order and orient them according to the published reference sequence (genome) of the same or closely related species. In order to do this, you can use abacas to contiguate your genome and it can also automatically design primers for you so that you can sequence the 'gapped' regions and close the genome. Another bacterial genome finishing tool that you could use is CONTIGuator.

ADD COMMENTlink modified 7 months ago by RamRS27k • written 5.7 years ago by Chrispin Chaguza250
Please log in to add an answer.


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