Chiseq files for deseq2
Hello everyone, I have 4 peaks bed files and i want to create count matrix for analysis in deseq2 I want to see the difference acetylation between them, I want in the column the 4 samples, and in the rows peaks I understand i can make it by using this function:
featureCounts -p -O -T 8 -a $gtf -o $outfn --extraAttributes gene_name -Q 10 $*
But i am not sure how to do that, I am really confused and i would like to get any details about that
Thanks!!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
htseq, stringtie, featurecounts all not working on nanopore direct RNA
written by CL •
-
featureCounts problems with GTF generated by myself
written by soda •
-
use ROSE to identify super enhancer
written by CL •
-
counting EXOM reads using subread featureCounts
written by adR
-
bedmap to annotate bed file with GTF that was converted to BED produces way too many results
written by Tawny
-
featureCounts strand specify parameter
written by ek699 •
-
How to connect htseq-count to DESeq2?
written by takoyaki
-
Should I use gene_id or gene_name in featureCounts for differential expression gene analysis?
written by biock •
-
featuresCount order sam file?
written by noeD
-
Is there any differences between tophat, cufflinks command with and without GTF file?
written by bioinforesearchquestions
Check the SAF format in the manual: https://subread.sourceforge.net/SubreadUsersGuide.pdf You will need to convert the BED files into a format compatible with featureCounts. SAF is simpler than GTF but some scripting will still be required. Last but not least: some peaks in your four BED files may be identical or overlap with peaks from other files. You will need to decide how to treat these overlapping peaks. BTW do you also have input/control BAMs?