Overlapping paired end - tophat

[ACTGangster]

Running 2.7 and this doesn't happen. Do you have setuptools installed?

Try typing import setuptools in the python console.

If not, download setup_tools and install it using python2.7

[gpcr]

Thanks a tonn.
It worked after importing setuptools.

[ACTGangster]

no problem. email me if you need anything else: adavisr at ufl dot edu.

Also, please cite the website if it works. https://github.com/audy/stitch

[jstjohn]

I have a working draft C program to do adapter stripping/ paired-end merging similar to stitch and SHERA for large gzipped datasets. It looks like it processes somewhere around 20M 100x2 pairs per hour. Its available here if anyone is interested:

https://github.com/jstjohn/SeqPrep

I don't have correctness statistics available, but the program can copy a subset of the merged reads into a human-readable aligned format so you can sanity check the settings.

[christophe]

Originally Posted by Cole Trapnell

As of TopHat 1.0.13, you should be able to specify a negative inner distance of -30. TopHat does map the reads independently, and has a different algorithm from Bowtie for handling the ends. The coverage.wig file display depth of read coverage, not depth of physical coverage, so those bases will be double counted, as you suggest. However, Cufflinks operates at the fragment level, not the read level, and so should do the right thing here.

Hi Cole,
First, According to the quote above, we may use negative inner distance. If I set it up at -30 and the stddev at 10, this mean tophat will look for overlapping of size 20 to 40.
Am I correct?


Second
if I have a mix of pair-ended reads such as some read pairs are overlapping some are not;
---> May I set up the inner distance to 10 and then the std deviation for the inner distance at 50 for example?
---> In that case, will Tophat look for inner distance from -40 to +60? (-40 meaning an overlap of 40)

Thanks a lot in advance for your help.
Chris***

[luc]

I just wanted to say thanks for both "Stitch" and "SeqPrep". They seem to be very useful !!