Dear all,
I want to align my ChIP-Seq samples but with two strategies:
- Align the
fastq.gzfiles and remove the reads that match to more than one site - Align the
fastq.gzfiles and, for those reads that are duplicated that match more than one site, assign a random site for them.
Can someone tell me hoe to do this?
As far as I know, I can use bowtie2: For second situation (random assign multi-matching reads):
- Use the
-aargument forbowtie2to look for all the "best" sites for a read and selects the binding randomly from the "best" sites.
But I don't know how to perform the first situation (delete multi-matching reads).
•
written
by
kuragari.ch •
2 answers
Hi, I can recommend Bowtie1 for this purpose. It has an option for exactly what you want:
-m <int> suppress all alignments if > <int> exist (def: no limit)
So in your case, you want to set -m 1
I don't have experience with Bowtie2, (I've looked in the help but couldn't find a similar option), maybe this options is hidden in a different parameter, or it is not possible.
For you 2nd strategy, in Bowtie1 you can use:
-M <int> like -m, but reports 1 random hit (MAPQ=0); requires --best
Again, for your purpose -M 1
For both cases you can use bowtie2 as it will map your reads using random mapping but you can filter them afterwards. Be sure to use an assembly that contains all unassigned regions (those _random). Otherwise some regions will appear unique when they are not.
To remove multi-reads, filter the mapping file file by mapping quality. Multi-map reads usually have a low quality score of 1. This is as far as I understand the best way to remove multi-reads.
I will discourage the use of bowtie1 if your reads are longer than 36bp.
Log in to answer this question.
-
How to count the number of multimapped reads from a bam file?
written by O.rka
-
Dealing with multimapped reads with Chip-seq and ATAC-seq data when studying Transposable Elements
written by ahnje770 •
-
Bowtie2 error :: (ERR): bowtie2-align died with signal 6 (ABRT) (core dumped)
written by shrutidabral •
-
inconsistent results of rsem quantification
written by gnilihzeux •
-
How to remove duplicate sites and inconsistent sites from a VCF?
written by Sharon
-
featureCounts allowMultiOverlap - fractional or not?
written by Lucy •
-
Finding alignment score for multi-read frombam/sam file
written by EVR
-
How to account for multiple mapping reads while calling variants?
written by kirannbishwa01
-
Bowtie2 or RSEM (with Bowtie2) alignment of PE reads with reference transcriptome?
written by deanpett •
-
Count multi-maped reads independent of mapper
written by Tomislav Ilicic •
