why bwa-mem only return one read with pair-read end
Hi all,
I have a very basic question. We introduce read1 and read2 into the bwa-mem; why does in the bam file output, they only give 1 sequence in return?
sequencing
bwa-mem
• 3,845 views
•
link
updated
by
rfran010 •
•
written
by
QX •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
bwa mem skip orientation when aligning a small subset of reads
written 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 …
-
nextflow set tuple issue
written by JoeDoasi •Hello, I recently started using Nextflow and I always get this error: > WARN: Input tuple does not match input set cardinality declared by process …
-
can I map multiple (>2) fastq reads against and assembly using bwa?
written by sa91 •Hi I want to map the output of 4 separate illumina smrt cells against a pacbio reference genome using BWA. I've indexed the full genome. …
-
bwa is getting stuck
written by joomi •Hello, I've been running bwa mem to do pair wise alignment to a reference. I have already indexed my reference and it resulted in 5 …
-
How to generate a BAM/SAM file to show a mutation using BWA aligner
written by neranjan •Hi, I am trying to generate Single nucleotide mutations to a sequence, and trying to get SAM files and its related CIGAR sequences. My test …
-
Aligning very short and paired Illumina reads
written by RoxHi everyone, I faced a issue using bwa for some very shorts reads I have (49bp). I have numerous paired Illumina reads that I want …
-
How to find the insertion sequence in .bam file?
written by kimpole1017 •I come up to aligning the fastq file based on the reference and as a result, I got the aligned bam file. For the next …
-
bwa sampe error [fread] Unexpected end of file
written by gsuryawanshi •I am running following command to convert paired end .sai output of bwa aln to one sam file: bwa sampe ref.fa read1.sai read2.sai read1.fastq read2.fastq …
-
bwa mem on a bam file
written by vimartin •Hi, I'm trying to use bwa mem for the first time and I wonder, is it possible to run it on a bam file instead …
-
Questions About Pair-End Reads Which Are Not Match For Bwa
written by Tonyzeng •<p>My question is that my original data or pair end read data are consisted of Read1 and Read2 and they are corresponding to each other. …
Have you tried name sorting the alignment file and/or viewing the alignment in IGV? Paired-end reads are not guaranteed to be next to each other in the file. There may be a discordant alignment for the read pair (because the read captures a chromosomal breakpoint). Or the second read has not mapped.