Hello,
I need to extract reads from a bam file - for each genomic position, I need to extract a different number of reads.
Eg.
CHR START END #READS
2 12345 12350 10
2 14567 14568 7
I have gone through samtools, and cannot see a way to sample a specified number of reads for each position. I have also ran into some elegant perl scripts that will do this by position and %, but not variable for each position.
Can anyone help?
Thanks in advance!
I need to extract reads from a bam file - for each genomic position, I need to extract a different number of reads.
Eg.
CHR START END #READS
2 12345 12350 10
2 14567 14568 7
I have gone through samtools, and cannot see a way to sample a specified number of reads for each position. I have also ran into some elegant perl scripts that will do this by position and %, but not variable for each position.
Can anyone help?
Thanks in advance!
Comment