Entering edit mode
8 months ago
qudrat.nii
▴
10
Hello,
I have ATAC-Seq bam files and I want to know the quality of files, so I run featureCounts. I found the following result...
Successfully assigned alignments : 54255620 (40.3%)
Is it OK or what should be the percentage of assigned alignment?
If this is the percentage of reads overlapping callable peaks then this actually pretty good. Some context on commands and what you use as annotations here to count would be helpful though.
I use following command:
If you are using
-p
for paired end reads then you will want to use following as well to count fragments.As said, if this is reads in peaks (FRiPs) then this is fine.
Sort of impossible to answer without more specific information. What organism? What annotation? Has anyone ever done ATAC Seq in this organism before? Do you have a positive control experiment (even in another organism to control for your analysis pipeline) that you can compare the number to?