Hello, We have a use-case where we would like to downsample a bam to exactly N reads. One way to do it is to use: bamshuf | view | head -N | sort ... (of course accounting for the length of the header, etc) For this, we would like to ignore read groups and output random read pairs from the original bam. Would it be possible to add a switch to bamshuf to not group by read-groups?
Thank you, Stathis ------------------------------------------------------------------------------ Dive into the World of Parallel Programming The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net/ _______________________________________________ Samtools-help mailing list Samtools-help@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/samtools-help
