More posts like this
-
How to iterate through paired-end reads in pysam?
written by Tao •
-
Which side do I need to trim when having adapter residues in my reads?
written by gogeni5529 •
-
STAR split reads into two strands
written by Alexander •
-
Trim bam file reads to the same length
written by rrapopor •
-
Sampling random reads from BAM files with allowing same reads
written by SSK •
-
Trimmomatic; How to get only q20 < q < q30 reads
written by SSK •
-
pysam: Use read ID to access location
written by bd2810 •
-
Iterating through paired-end reads in samtools / pysam
written by olavur
-
Trim Reads From BAM files
written by fusion.slope
-
cutadapt for trimming reads to specific length
written by Varun Gupta
Your question lacks any details. What do you want to trim? What is
not whole reads but each read? Please give representative examples.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.
Ex.
[pre bam]
read1 150b
read2 150b
read3 150b
[trimmed bam]
read1 100b (trimmed 50b)
read2 150b (trimmed 0b)
read3 20b (trimmed 130b)
Do you understand what I mean?