Hello!
I am currently aligning some sequences that contain numerous chimeric reads using bwa (paired-ends). Usually, first third of the R1 is on Chr1 and the rest on Chr18, the R2 is fully mapped to Chr18.
I was wondering how bwa was aligning those reads. If we look at the sam file produced after alignment, what will be the RNAME of the read R1 (chr1 or chr18)? And how is it decided?
Thanks in advance!
•
written
by
VHahaut
1 answer
BWA should create multiple SAM records for the same Read. A first SAM record will be mapped on chr1 and will contains a 'SA:X' attribute containing the alternate alignments, the cigar string will be soft-clipped. The other reads will be the instances of the 'SA:X' tag above = the other alignements, and they will be flagged as 'supplementary alignement'.
Log in to answer this question.
-
Pileup of inner distance for PE reads
written by noodle
-
Number of secondary alignments per read
written by am3
-
Nanopore incomplete fragment sequencing
written by VHahaut
-
plasmid spades R1/R2 files?
written by sarah.goldstein •
-
Obtain Number of PE Reads in BAM file Mapped to Each Chromosome
written by ssm87 •
-
bwa mem fastq order
written by guffrann •
-
How to trim both pair end reads and remove low quality reads
written by Bioinfonext
-
How does Bowtie2 handle chimeric singletons?
written by jnowacki
-
How to obtain insert sequence from paired-end mapping result (SAM/BAM)
written by gundalav
-
Bam To Bed With Read Nucleotide Sequences
written by rhasbunz •
