[khmer] Sweep-reads3.py

C. Titus Brown ctbrown at ucdavis.edu
Wed Dec 16 06:16:20 PST 2015


Hi Raf,

yes, the sandbox/ scripts aren't subjected to our continued and relentless
automated testing regime :).

What is the inital format of your read names?

best,
--titus

On Wed, Dec 16, 2015 at 07:54:08AM +0000, Winand, Raf wrote:
> Dear Khmer mailing list
> 
> After running the sweep-reads3.py script from the sandbox I ran into a problem trying to extract the paired reads after reinflating. During the reinflation the /1 and /2 at the end of the read names get removed and the reads are separated in the output file. Trying to extract the paired reads from the files as mentioned in the tutorial (http://khmer-protocols.readthedocs.org/en/v0.8.4/metagenomics/3-partition.html) therefore always results in an error that no paired reads could be found. As this script is found in the sandbox I'm not sure if it is still supposed to work or of it is deprecated?
> 
> Best regards
> Raf
> 
> ###########################################################
> Checked as being free of known viruses.
> 
> Scientific Institute of Public Health, Brussels, Belgium
> Wetenschappelijk Instituut Volksgezondheid, Brussel, Belgi?
> Institut scientifique de Sant? publique, Bruxelles, Belgique
> Visit our website: http://www.wiv-isp.be 
> ########################################################################
> DISCLAIMER: Please see https://www.wiv-isp.be/Pages/EmailDisclaimer.aspx
> ########################################################################

> _______________________________________________
> khmer mailing list
> khmer at lists.idyll.org
> http://lists.idyll.org/listinfo/khmer


-- 
C. Titus Brown, ctbrown at ucdavis.edu



More information about the khmer mailing list