Proper Bowtie2 settings for snp calling applications

Hi,

I'd like to map illumina paired-end reads to the reference genome by Bowtie2 for downstream snp calling with Freebayes. The reads are already adapters and quality trimmed. I plan to specify the following Bowtie2 options:

--end-to-end
--very-sensitive
--no-mixed
--no-discordant

I'm wondering which options i have to include/change in my Bowtie2 command line to produce the most reliable results.

Regards,
Denis

illumina snp freebayes bowtie2

• 1,305 views • link

updated 3.0 years ago by Ram

• written 3.0 years ago by Denis

1 answer

My recommendation is to use bwa mem with default settings which has proven a gold standard over the years. Bwa is expected and tested by almost any variant calling pipeline. No need to use anything else.

• 0 views • link

written 3.0 years ago by ATpoint

More posts like this

Drastic drop in RNA-seq read mapping rate when disabling gaps in Bowtie2
written 1.2 years ago by triplee0305 •

Hi all, I have some RNA-seq reads from a grasshopper species. I mapped them to a transcriptome reference (from the same species) using Bowtie2. When …
Block size of SSD file storage
written 1.4 years ago by Denis

Hi there, I'm creating a `virtual machine` for different bioinformatic applications (`de novo assembly` and `snp calling` primarily) and have to specify a block size …
It is weird there are no reads less than 100 bp in qualimap.
written 4.1 years ago by Dan •

For ATAC-seq, I trimmed the reads using ``` trim_galore --illumina --paired --fastqc \ fastq/${FILE}/*R1_001.fastq.gz fastq/${FILE}/*R2_001.fastq.gz \ -o trim/${FILE} ``` When I run Bowtie2 using ``` …
rsem-calculate-expression error
written 5.1 years ago by Bhoomika •

**I have prepared reference using the command:** rsem-prepare-reference -p 5 --bowtie2 cd_hit_est.fasta ref_output_12 **And the files i obtained are:** ref_12.1.bt2 ref_12.2.bt2 ref_12.3.bt2 ref_12.4.bt2 ref_12.grp ref_12.idx.fa …
Low mapping percent (whole metagenome shotgun assembly)
written 5.4 years ago by mewgia •

Hello! I have 4 samples of whole metagenome shotgun sequencing results (soil communities). Illumina, single end, read length ~76 bp. The raw reads were preprocessed …
Plotting range/type of discordance in alignment data
written 6.8 years ago by rbronste

I was wondering if there was a way, in a Bowtie2 output, to plot the range and type of discordance? The really obvious kind when …
RSEM ignores quality scores?
written 9.0 years ago by AW

Dear all, I want to quantify gene expression for my Trinity transcriptome assembly using RSEM. I am getting confusing results depending on how I specify …
Bowtie2/HISAT using single-end to generate FPKM???
written 9.2 years ago by jamieson.pierce •

Hi all, I am familiar with the fact that paired-end reads are used to generate FPKM, and single-end reads are used to make RPKM. I …
Bowtie2 error with single end mapping
written 9.4 years ago by H.Hasani

Dear all, I'm facing a tedious error and not sure what causing it. I have single-end reads and simply trying to align it using Bowtie2 …
How to obtain insert sequence from paired-end mapping result (SAM/BAM)
written 9.9 years ago by gundalav

I have paired-end reads R1 and R2. |-----R1-----|----------------------------|-----R2-----| A B C Total insert here is A + B + C. How can I get this …