I'm trying to trim each reads in bam.
I want to crop not whole reads but each read, so I don't think trimBam is good way.
How can I do this ?
Your question lacks any details. What do you want to trim? What is not whole reads but each read? Please give representative examples.
not whole reads but each read
Sorry, my words was not enough.
I want to do bootstrap method. Thus I sampled randomly from bam file, and made subsampled.bam. On the top of that, I want to crop read randomly from subsampled.bam.
read1 100b (trimmed 50b)
read2 150b (trimmed 0b)
read3 20b (trimmed 130b)
Do you understand what I mean?
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy