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.
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
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.
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.
Log in to answer this question.