MAPQ mapping quality

I would like to filter only mapped reads with hight MAPQ. I use bowtie2 What value should i give to MAPQ? i would like to make variant calling after filtering good reads from bam.

snp mpileup call

• 3,684 views • link

written 6.0 years ago by sami •

1 answer

Useful links:

• 0 views • link

written 6.0 years ago by Ram

More posts like this

Why Is My Uniquely Mapped Read Count Higher in Samtools Compared to Bowtie2?
written 1.5 years ago by daffodil •

Hi, I aligned my paired-end sequencing data using Bowtie2 and processed the resulting BAM files with Picard MarkDuplicates. However, I noticed a discrepancy between the …
Which threshold to be put in the MAPQ when doing mpileup of BCF tools on a STAR-generated BAM file ?
written 3.7 years ago by mohsamir2016 •

Dear All Sorry for bothering. I am doing variant calling using BCFtools mpileup then call and I would like to add a filteration step in …
bbmap, mapq and uniquelly mapped reads filtering
written 5.3 years ago by boczniak767

Hi, I'd like to know, how do postalignment filtering to get only *uniquelly mapped* reads. I know that it is possible at the `bbmap` level, …
FAIRE for non-standard cultivar, mapping to reference, MAPQ
written 5.4 years ago by boczniak767

Hi, I have FAIRE-seq data from two cultivars: reference B73 and less known S68. I assumed I can align reads from both cultivars to B73 …
bowtie2 maximum value of MAPQ
written 6.0 years ago by sami •

I would like to know the maximum value of MAPQ generated by bowtie2 ? Thank you
Remove duplicate reads before variant calling?
written 6.6 years ago by O.rka

I have a very large fastq file that I want to do variant calling on. Would it make sense to remove duplicate reads before mapping …
Samtools: How to extract reads with MAPQ 0?
written 7.3 years ago by stefano.meucci •

After performing BWA MEM alignment, I want to extract only the reads with MAPQ 0, how to do it? The "-q" parameter only filter out …
Bowtie2 - How can MAPQ be used to obtain uniquely mapped reads?
written 8.4 years ago by ej

Hi, I am interested in extracting uniquely mapped reads from a bam file by filtering based on the MAPQ field. I ran some tests to …
Filtering high quality reads from bowtie2
written 11.0 years ago by yulia.gordon •

Hi everyone! I've encountered an interesting discrepancy when trying to filter the results of my paired end CHIP-SEQ. When using bowtie2 in Galaxy interface it …
Population Rate Calculation Programme
written 13.9 years ago by jackuser1979

I have RNA seq of illumina paired reads mapped to reference genome (bam file) and VCF file from variant calling software. I would like to …