Hello everyone, I am confused about a strategy hired to mapping paired-end sequencing data which was employed in the paper Accurate identification of human Alu and non-Alu RNA editing sites. In brief, they mapped each of the paired-end reads separately using the commands “bwa aln fastqfile” and “bwa samse -n4”, and then merged these reads into a single file. But bwa has the command bwa sampe, why they just not use that command?
•
written
by
Proxy Agent •
1 answer
Previously BWA didn't have support for direct alignment of paired-end sequence data. So, first, you had to align the files in .sai(suffix array index) format separately then merge the matching pairs.
Now, with bwa aln -1 <read1.fq> -2 <read2.fq> you can directly align paired end reads.
Log in to answer this question.
-
Alignment tool for very short reads (25bp)
written by Anisur Rahman •
-
Working with VEP in anaconda
written by luzglongoria •
-
BWA mapping parameters for DNA capture sequencing protocol
written by graeme.thorn
-
Forum: BWA command guide
written by Chen Sun
-
RNA editing mapping step with BWA
written by AleF •
-
Help, I'm having problem with bwa mapping
written by davidcopen •
-
PE alignment with bwa aln and conversions
written by raplayer •
-
Can Anyone Suggest Me A Script Based Pipeline For Exome Sequencing With Paired End Reads Generated …
written by ivivek_ngs
-
Tutorial: Reference Assembly - Mapping Reads To A Reference Genome
written by Joseph Hughes
-
Tutorial: Mapping Reads With Bwa And Bowtie
written by MSU course 2011 •
It looks weird, not to make use of the paired information. Probably it has to do something with the fact that they are studying editing, and therefore there can be a difference between reads and template.