This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Close Bacterial Draft Genome With Other Draft Genome

Hi,

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

genomics assembly
• 5,067 views • link
updated by Ram Ram • written by biotech biotech

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.

• 0 views • link
written by Xingyu Yang Xingyu Yang

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

• 0 views • link
written by biotech biotech

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 - | vcfutils.pl vcf2fq > <Assembled Results>.cons.fq
• 0 views • link
updated by Ram Ram • written by Xingyu Yang Xingyu Yang

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.

• 0 views • link
written by Torst Torst

2 answers

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

https://github.com/fenderglass/Ragout

Hope it helps!

• 0 views • link
updated by Ram Ram • written by mesabioinformatics • mesabioinformatics

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.

• 0 views • link
updated by Ram Ram • written by Chrispin Chaguza Chrispin Chaguza

Log in to answer this question.

More posts like this