Hello everyone,
I am working on RNA-Seq data analysis, and currently I have generated the bam files (named WT.bam and KO.bam).
Could I ask whether deeptools (https://github.com/deeptools/deepTools/) is available to extract reads that exist in WT.bam but not in KO.bam?
If it can't, is there any other methods (or scripts) to achieve the target?
1 answer
to extract reads that exist in WT.bam but not in KO.bam?
Are you referring to recovering reads from areas of the genome that are covered in WT as opposed to the KO alignment file? If so, you could use bam2bed (https://bedops.readthedocs.io/en/latest/content/reference/file-management/conversion/bam2bed.html ) to first convert the files and then use a set operation like difference (https://bedops.readthedocs.io/en/latest/content/reference/set-operations/bedops.html#difference-d-difference ).
Once you get the read names that are present in one of the files, you can use a tool like filterbyname.sh from BBMap suite to extract the reads from original files.
Log in to answer this question.
-
Performing Differential Analysis with NarrowPeak Files in ATAC-seq: Tumour vs. Control
written by Thom •
-
Seeking Tutorial for Downloading TCGA Data with TCGAbiolinks and Converting to TPM
written by Znow •
-
deeptools k-means clustering
written by biology_inform •
-
Publishing RNA sequencing data analysis scripts on Github before or after Publication?
written by m.sadman.sakib
-
Comparing gene expression data between datasets
written by kpukacz •
-
I need to extract the FASTA reads from a specific region from a chromosome in a BAM file format.
written by screadore •
-
How to operate gene funcational annotation
written by citronxu •
-
Combine microarray and RNA-seq gene expression for differentially expression analysis
written by landscape95
-
Forum: Reuse of publicly available datasets
written by sanathoi •
-
meta analysis of p value data by metaseq package
written by hassan giahi •
