Removing contaminations from PacBio reads
0
2
Entering edit mode
7.6 years ago
Ric ▴ 440

Hi I have PacBio RS II and Sequal reads. I am looking for a tool which removes common contaminations and allow also to provide fasta files with chloroplast and mitochondrial genomes and it would remove these genomes from the reads too.

Thank you in advance.

pacbio sequencing Assembly fastq • 3.1k views
ADD COMMENT
0
Entering edit mode

If you know what contaminant genomes are present then you could see if bbsplit.sh from BBMap suite may work. Otherwise use mapPacBio.sh (again from BBMap) to get the reads you want (assuming you have a reference) and leave the rest behind.

ADD REPLY
0
Entering edit mode

Hi Ric, did you managed to solve your issue using Genomax answer ? I'm still stuck with this question, because I don't now yet all organism that contaminated my assembly. Any others suggestions ? I was thinking about blasting maybe a sequence highly specific to bacteria, which is shared among bacteria (maybe following the busco criterias) and then remove all contigs that contains theses sequences... Does that sound good, or is it a bit harsh ?

ADD REPLY

Login before adding your answer.

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