More posts like this
-
VCFfiltering
written by Peter •Hi, A) I have a large set of samples diploid, tetraploid, and hexaploid mixed. I have used FreeBayes to build a VCF file. I want …
-
What causes very low read mapping? CUT&RUN with Bowtie2
written by smurph50 •Generally, what causes low read mapping efficiency? Some samples have 70% and others are 8-20%. I am experienced with bowtie2 and tested with/without trimming. These …
-
Get rid of multimapped and unmapped reads for snp calling improvement
written by DenisI'm doing `snp` calling for several projects in which either `BWA` or `Bowtie2` was used to map reads to the reference genome. I use only …
-
VG Toolkit for human WGS analysis
written by venkat •Hello! I have been developing a computational pipeline with VG Toolkit for human WGS analysis. I have built a genome graph with GRCH38 and variants …
-
bowtie2 -all ... multiple alignment Chip-seq
written by poltora4enko •Hello! I'm trying to use bowtie2 with options `bowtie2 -p 8 -a -x bt2_file -U file.fq -S file.sam`, because I need multiple alignment for my …
-
Calculate ChIP-seq correlation in promoter regions
written by dcheng1 •I read a few papers that calculate Pearson correlation coefficient of ChIP-seq replicates **in promoter regions**. However, even for active marks like H3K4me3, not 100% …
-
Filter reads by calculating the distribution of phred score in each individual read
written by Nagesh •Dear All, I would like to know if there is any tool or script to calculate the distribution of PHRED score for each individual read …
-
bowtie2 filter for maximum alignment score
written by rishi •Is it possible to set a filter/argument in which bowtie2 reports just one alignment for each seed, the one which is of the highest alignment …
-
High confidence junction events in complete genomics data
written by MAPK<p>I have been trying to understand translocation data from complete genomics. They have high confidence and alljunction data(includes low confidence events). I want to understand …
-
Practice Of Filtering Vcf Files (From Gatk)
written by newDNASeqer<p>I use GATK to make variants calling on exome sequencing data from human tumor samples, and have been using GATK for a few months now. …
I personally use MAPQ 20 for standard tasks like ATAC-seq or ChIP-seq. But don't ask why 20 and not 25, 15 or 10. I think I stumbled over 20 which was used by a reputable tool when I got started with NGS so I took it over. Unless you are investigating difficult genomic regions specifically, choice should not matter too much as long as it is > 0 to avoid multimappers.
Thanks for your help!
The ENCODE pipeline uses a MAPQ 30 cut-off