Question: Transposons and Pacbio
0
gravatar for burnsro
4.2 years ago by
burnsro20
Austria
burnsro20 wrote:

What tools are available to call TEs using Pacbio data? My understanding is that most work using paired-end reads from short read Illumina data.

pacbio • 1.2k views
ADD COMMENTlink modified 4.2 years ago by nterhoeven110 • written 4.2 years ago by burnsro20

What are you trying to do exactly? Map to a genome or infer repeat abundance straight from the reads? You are correct though, Pacbio is not as pervasive and most tools have been developed with Illumina or 454 data.

ADD REPLYlink written 4.2 years ago by SES8.2k

I was planning to do de novo assembly, and infer repeats from the reads

ADD REPLYlink written 4.2 years ago by burnsro20

I would not normally try to assemble the reads, though this may be a reasonable approach with this data (depending on the species). Try Transposome, it works very well with long single-end 454 reads so I'm sure it will do fine with Pacbio. I'm the author of that tool so you can ask me any questions about it if you like.

ADD REPLYlink written 4.2 years ago by SES8.2k
1
gravatar for nterhoeven
4.2 years ago by
nterhoeven110
nterhoeven110 wrote:

The advantage of PacBio reads is that they are quite long, so it may be possible to use methods designed for finding TEs in assemblies. I tested a workflow based on this protocol with Arabidopsis data and it worked quite well.

ADD COMMENTlink written 4.2 years ago by nterhoeven110
1

Did you use raw or corrected PacBio reads?

ADD REPLYlink written 4.2 years ago by thackl2.7k

corrected with prooveread :-)

ADD REPLYlink written 4.2 years ago by nterhoeven110

That protocol is for identifying repeats in a genome assembly. OP wants to infer repeats from the reads.

ADD REPLYlink written 4.2 years ago by SES8.2k

yes, I know.  Read the second sentence of my answer ;-)

ADD REPLYlink written 4.2 years ago by nterhoeven110
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1968 users visited in the last hour