One read of the pair-ended reads is not able to map to the region when the reference is shorter using BWA
Hi All, I am not sure if anyone has the same experience as mine. When I perform BWA mapping using whole genome as reference, I could get both of a specific pair-ended reads mapped to the desired regions. However, when I extracted the regions from the target site (extended 500 bp front and back) and performed BWA mapping, one of the pair-ended reads is unable to map to it. The only difference is for whole genome as reference, I used -a bwtsw option and for short sequence as reference, I used -a is option. I wonder what the reason is behind it. Appreciate if anyone may know the answer. Thanks.
• 2,187 views
•
link
written
by
Louis Kok •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Aligning whole genome sequencing pair-end reads with a previous made STACKS catalog from ddRADseq p…
written by rujiporn.sun •
-
Paired-end amplicon sequencing data
written by bioplanet •
-
How to know which reads were not aligned to a reference genome?
written by Charles Yin
-
samtools flagstat result is different from samtools view result
written by Louis Kok •
-
Low quality bases at the starts of HiSeq 2 x 100bp reads
written by Louis Kok •
-
Single end or paired end for variant calling
written by memory_donk
-
Fastq reads mapped to the desired position on reference genome
written by Louis Kok •
-
Merging pe and contig alignments with Samtools
written by zggayk •
-
How Are Human Sex Chromosomes Handled When Mapping Reads?
written by robert •
-
Almost Nothing Mapped Using Bwa Or Bowtie. But A Lot Of Mapping In Blast
written by lwc628
