This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What are the methods and recources of filtering contigs after de novo genome assembly ?

Hello,

I have a plant chloroplast genome that I de novo assembled using velvet. Now I want to filter the contigs.

Kindly suggest me which methods will be the best for this purpose and what tools should I use that are not much complicated.

Thankyou.

Best Regards, Faizan Saleem

assembly genome

3 answers

Using BLASTn against nt is very helpful. Download the nt DataBase from SRA. and do Blast against with your query contig seq file. Here is the comment : blastn -query /home/assembly.fa -db nt -max_target_seqs 1 -outfmt '6 qseqid pident evalue staxids sscinames scomnames sskingdoms stitle' -out outblast.txt. All the Best.

Thanks alot. Someone told me that I have to filter my contigs according to the DNA threshold value or something like that and I dont understand how to do it. can you help me please.

Typically, I'd suggest BLASTing them against nt or similar, and removing the ones that hit suspicious things. RefSeq also has a plastid dataset; you could align your contigs to that also.

Furthermore, the chloroplast contigs should have similar coverage. Once you know the chloroplast coverage, you can usually just throw away contigs with very different coverage. For example, if the chloroplast is 500x on average, and you get some contigs with 100x coverage, those are probably something else (like the plant main genome). You can determine coverage by mapping the reads to the assembly (e.g. with BBMap: bbmap.sh in=reads.fq ref=assembly.fa covstats=covstats.txt) or usually by looking at the contig names, though that's less accurate.

Thanks alot. Someone told me that I have to filter my contigs according to the DNA threshold value or something like that and I dont understand how to do it. can you help me please.

There is no recipe for decontamination. You need to use your judgement.

Thanks alot. Someone told me that I have to filter my contigs according to the DNA threshold value or something like that and I dont understand how to do it. can you help me please.

Log in to answer this question.