BWA option for FR sens only

Hi,

I am looking to bwa mem to map my reads.

However, I would like to keep only the real Forward-reverse paired end and to discard any weird sens (FF, RR, RF).

Is there any option in bwa to do this ?

Thanks

bwa mem

• 2,615 views • link

written 9.7 years ago by Picasa

maybe align them first then extract them using samtools view -hf 0x2 alignments.bam

• 0 views • link

written 9.7 years ago by Medhat

Filtering by flag is a good idea. In case you need the flag concept explained, check the SAM format specifications and also this website, which can decode a flag:

Bitwise Flag explained

• 0 views • link

written 9.7 years ago by ATpoint

You could also pipe bwa directly into samtools view, avoiding intermediate files.

• 0 views • link

written 9.7 years ago by WouterDeCoster

0 answers

No answers yet.

More posts like this

bwa mem skip orientation when aligning a small subset of reads
written 3.1 years ago by Yihang •

For pair-end read files, `read1.fastq`, `read2.fastq`, I use the following command `cat read1.fastq | head -n 4000 > read1_first1000.fastq` (and `cat read2.fastq | head -n …
BWA-MEM : skip orientation FR as there are not enough pairs ?
written 3.6 years ago by Saran •

Hello, I have Paired End Data that I am trying to align to a very short reference(222 bp); I expect that the paired-end data have …
skip orientation FF/RF/RR by bwa mem
written 5.9 years ago by 14301010010 •

Hi everyone I extracted DNA fragments from human somatic tissue, sequenced them, and then try to align data to the ref genome (hg38) using bwa …
Question: BWA mem skip orientation
written 6.8 years ago by BellaK •

Hello everyone, I am very new to genomic analysis and alignment, and I am posting this with the hope that someone with more knowledge can …
BWA mem skip orientation
written 7.3 years ago by williamsbrian5064

Hi, So I am not entirely sure if this is an issue or not. So to give a little background, I used bwa, picard, samtools, …
Strand specific pair end sequencing
written 9.4 years ago by Bioinfonext

I have a strand-specific pair-end Illumina sequencing data. I am not able to understand which option to use in the different software. For each library …
bowtie --fr/--rf/--ff reporting
written 10.0 years ago by ThePresident

Hello, I am mapping Illumina generated paired-end reads with bowtie. By default, bowtie reports --fr, i.e. paired mates are in Forward/Reverse position (which is normally …
bwa -mem log messages
written 10.4 years ago by aleka

Hello, I am using bwa mem to aligh paired ends and I get the following errors: [M::process] read 1057202 sequences (80000122 bp)... [M::mem_pestat] # candidate …
BWA MEM 0.7.12 Error : [gzread] <fd:4>: invalid distance code
written 10.7 years ago by gsr9999

Dear BioStars Leaders, I have been running BWA MEM 0.7.12 on a handful of paired-end sample fastq files and it worked out pretty well. I …
BWA MEM indexing
written 10.9 years ago by david.4.ray •

I'm a relative newbie and trying t work my way through using bwa mem. I've successfully aligned PE reads to a reference genome previously using …