Low number of reads mapped, looking for a relaxed paired-end Illumina mapper

Hi,

I mapped 2x75 Illumina to a bacterial genome using BWA with default settings. Very low number of reads were mapped. I was used to TMAP tool when worked with IonTorrent data, with high number of mapping reads always.

I'm looking for a relaxed mapping tool. Maybe could be useful to relax BWA thresholds before?

Thanks for your help,

Bernardo

snp rna-seq illumina

• 3,427 views • link

updated 4.8 years ago by Ram

• written 12.1 years ago by biotech

1 answer

You put rna-seq among your tags. Are you aligning RNA-Seq data with bwa? If so, this could partially explain your low mapping rate and I would try a splice-aware aligner like tophat. Otherwise, bwa mem might be an option to consider if you expect gaps.

• 0 views • link

written 12.1 years ago by dariober

I'm doing some checks and will update the post soon.

• 0 views • link

written 12.1 years ago by biotech

Seems the problem is not the aligner. I have very low number of RNA-seq data from the organism I'm mapping to, only 10%.

• 0 views • link

written 12.1 years ago by biotech

You are studying bacteria so splice aware aligner may not be really necessary.

• 0 views • link

updated 4.8 years ago by Ram

• written 12.1 years ago by Ashutosh Pandey

More posts like this

Low % reads mapped and high % genome coverage
written 3.3 years ago by natural_decline •

Bacterial isolate sequenced on Illumina platform, fastq's are decontaminated prior to mapping and alignment. Coverage analysis run with samtools and then mapping stats calculated with …
How to reduce cross mapping alignments with paired-end reads ?
written 8.9 years ago by David

Hi, I´m trying to align paired-end reads (illumina) coming from a metagenomics experiment to a custom database of bacterial species. I´m using bwa-mem (althought could …
Prokaryotic RNA-seq / Mixing IonTorrent and Illumina data
written 10.8 years ago by caramba •

I have sequenced RNA from a bacterial culture grown under different conditions. The control sample was sequenced with IonTorrent, the other three other treatments were …
bwa misalignment for TruSeq Amplicon reads
written 11.0 years ago by idedios •

<img alt="" src="http://i.imgur.com/3BCFpgu.png" style="height:300px; width:533px" /> I'm having the following issues using bwa mem on paired-end TruSeq Amplicon reads. I used the default settings with …
Removing low complexity reads from RNA-seq
written 11.1 years ago by Asaf

I'm mapping bacterial RNA-seq to the genome and found something very weird. A very abundant RNA has a AATAATAAT repeat somewhere in the middle and …
RNA-seq reads mapping to paralogs with high MAPQ. How to ensure this?
written 11.7 years ago by biotech

I've mapped 75 bp paired-end Illumina RNA-seq data to a bacterial draft genome. We are interested in a family of genes belonging to trimeric autotransporter …
BWA Small Indels vs Roche Newbler Assembler
written 12.1 years ago by pld

I've been experimenting with using BWA to assemble 454 (1kb reads), comparing previous assemblies performed with the proprietary Roche assembler and I've noticed that BWA …
R1 R2 orientation (RNA-seq)
written 12.1 years ago by biotech

Hi, Will always R2 file contain sense transcript and R1 reverse complement of transcript? If R2 file contains sense transcript, I will have to run …
Counting sense reads in bacterial paired-end RNA-seq data
written 12.1 years ago by biotech

Hi, I'm trying to count reads mapping to sense strand. I have doubts which counts file I should chose from this pipeline. I think is …
Too Many Deg With EdgeR Output?
written 12.8 years ago by biotech

Hi there, I'm dealing with bacterial RNA-seq analysis. My experiment is very simple. Two samples to compare and no replicates. Reads were generated with Ion …