What are the methods and recources of filtering contigs after de novo genome assembly ?
3
0
Entering edit mode
7.8 years ago

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 • 2.3k views
ADD COMMENT
1
Entering edit mode
7.8 years ago
k.kathirvel93 ▴ 310

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.

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode
7.8 years ago

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.

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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