Exceptionally high peak/reads in ATACseq bam and peak files
Hi there, could anyone tell me how to remove the exceptionally high peak in the ATACseq? I got those high reads in both bam files and peak files? Thanks a lot!


sequencing
• 1,765 views
•
link
updated
by
i.sudbery
•
written
by
bright602 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Sorting BAM files never fully finished
written by dolphin6906 •I have some large BAM files that were created based on large indexes. I trying to sort and index them and threading the sort with …
-
SILAC protein data analysis
written by yueli7Hello I am new in analysis SILAC protein data. I have three patients samples which collected from three different points of time, and no replicates. …
-
vlm.adjust_totS_totU -- Input contains NaN, infinity or a value too large for dtype('float64')
written by ericrkofman •I am getting close to getting the velocyto pipeline working. However, I am running into some issues I think having to do with sparse data. …
-
Bioinformatics: tSNE Perplexity parameter
written by tom5 •Hi, I hope you're well. I am going through the 'Orchestrating Single-Cell Analysis with Bioconductor' tutorial found here [tutorial][1] and am unsure what the perplexity …
-
Unbalanced design and large difference in number of input reads for DESeq2
written by JW •Hi all, My understanding is that DESeq2 should by default normalize for differences in the number of reads. However, I am concerned that the magnitude …
-
GSVA R packages
written by Morris_ChairHello everyone, I'm trying to do a gene set varian analysis using R to detect a specific gene set signature of a specific pathway from …
-
How would you trim this RNA-seq data
written by williamsbrian5064Hi, I am working on an RNA-Seq project. We just got back some sequencing and I had a couple questions about it. I figure someone …
-
How to make plots with gene clusters based on the gene length and orientation?
written by xupbuy •I have >30 gene clusters in fasta files downloaded from NCBI. I want to make plots like this: ![enter image description here][1] [1]: https://i.postimg.cc/zf9N0p9C/Screen-Shot-2018-11-05-at-10-23-48-AM.png How …
-
Making a counts
written by Nick •Hi everyone, I'm a newbie in R and RNA seq analysis, so don't judge me so much) I'm trying to perform RNA-seq analysis as it …
-
Forum: Review of the CIGAR string format
written by Juke34**--EDIT-- For a corrected review see next answer below** I had difficulties to retrieve the meaning of **R** in a cigar string from the annotation …
Can you confirm that you have deduplicated the data and excluded the ENCODE blacklist regions?
I removed the PCR duplicates through Picard. For the species I am working on, ENCODE does not have the blacklist regions yet.
Did you set
REMOVE_DUPLICATES parameter to true?yes, I did set the REMOVE_DUPLICATES=true
Removed mitochondrial reads? Even if not I would ignore it and proceed with downstream analysis and only come back if you see odd things downstream.
Thanks, I have removed the mitochondrial reads too.