A tool to cut off a certain length of nucleotide from reads
2
0
Entering edit mode
4.9 years ago

Hi,

I have a genome of plant provided by a lab for me to assemble. This genome has Pac Bio and Illumina primer sequences that need to be removed beforehand. But I was told that some of the primers do not retain their full sequences, so instead of looking for the full sequences and trimming them out using software, they simply remove 70(?) nucleotide from both ends of every reads. Is this a better approach? And what tool would allow me to do this?

Thanks!

Assembly • 1.4k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
4
Entering edit mode
4.9 years ago
cschu181 ★ 2.6k

In before Brian Bushnell suggests BBDuk. Use BBDuk: http://jgi.doe.gov/data-and-tools/bbtools/bb-tools-user-guide/bbduk-guide/ >:D

ADD COMMENT
0
Entering edit mode

Oh, and I was so close, too!

ADD REPLY
0
Entering edit mode

Thank you! It's time to get back to bbduk then!

ADD REPLY
2
Entering edit mode
4.9 years ago

At least for Illumina reads, that's a bad idea. Especially for libraries under 140 bp, since you'd have nothing left (but assembly would run very fast!)

PacBio libraries don't have read-end adapters like Illumina. Rather, they have "Smart Bell" adapters somewhere in the middle of the read. PacBio's software already removes those for you when creating filtered subreads. Sometimes it doesn't get all of them, so I wrote a tool to remove them removesmartbell.sh), but it's not normally necessary. Therefore for PacBio my recommendation is: Do nothing, though you can try running removesmartbell.sh if you are feeling adventurous or are interested in knowing whether there are any leftover adapter sequences.

For Illumina, you should do adapter-trimming as indicated here. You can also use BBDuk to forcibly trim the first and last 70bp from every read, but I don't recommend that in this case.

ADD COMMENT
0
Entering edit mode

Thank you for your helpful answer! I will follow your suggestion!

ADD REPLY

Login before adding your answer.

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