Hello,
I tried to get coverage from BAM file in the following format:
chr strand chrom_start chrom_stop count
chr1 + 287 293 3
chr1 + 294 200 8
chr1 + 201 205 1
chr1 - 290 295 5
...
I can use genomeCoverageBed as blow
genomeCoverageBed -bg -ibam *.bam -g chr.size -strand + > pos.txt
genomeCoverageBed -bg -ibam *.bam -g chr.size -strand + > pos.txt
then put them together.
do you know a simpler way?
jian
-
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...07-08-2024, 03:19 PM