Following my other post, I was wondering how do the short reads mappers (bwa & bowtie2) handle the "-" strand of the reference at the step to build index. As short reads are always in 5'-->3' direction for either strand, but from the mapping result SAM file there are reads mapped to both strands.
1) I assumed the "-" strand of the reference is concatenated to the "+" so that the final length/index is simply doubled;
2) Or, it is needed to reverse complemented each read at mapping step so that mapping was done twice;
Mapper STAR uses the first approach as the author told me, but I could not figure this out from the source codes of both bwa or bowtie2. Can somebody confirm this for me? Thanks.
0 answers
No answers yet.
Log in to answer this question.
-
GSNAP insufficient memory error
written by Hood •
-
fastq reads cleaning
written by blueskypie •
-
Bowtie2 orientation mapping question
written by Andrés Ribone •
-
How is the FM-index built for multiple entries of a fasta file?
written by yifangt86 •
-
How do mappers behave at end of reference sequence?
written by nickp60 •
-
Identifying contaminants (cow, pig) in oral rinse exome sequencing data
written by AbiB1 •
-
Mapping mate-pairs reads with bowtie2
written by Biomonika (Noolean)
-
RNA-seq short-read mapping with wildcards (IUPAC codes) using GEM mapper
written by vitor.aguiar •
-
Should Rna-Seq Reads Map To One Strand Of Cdna Reference?
written by chjiao3456 •
-
Tool: Aligning Sequence Reads, Clone Sequences And Assembly Contigs With Bwa-Mem
written by lh3
I think bwa and bowtie index both the forward and reverse reference, and map once.
BBMap reverse-complements the read and does mapping twice.
Thanks Brian!
Is it possible for you to elaborate the indexing procedure for both the forward and reverse strand, technically? What I can guess is the "-" strand is concatenated to the "+" strand to have a single string. Then, how the offsets for each chromosome are distinguished from each other (2 strands and multiple chromosomes)? I saw the outputs of bwa index step for mouse genome, there are
Are those files with "r" (*.rbwt, *rpac, *.rsa) for the reverse strand? However, if I index the small genome as lambda_virus, I did not see the pattern.
Does this mean if the genome is big (Mouse) then there will be separate index for the "-" strand of the genome?
From bowtie2-build, I have:
Are the .rev.1.bt2, *.rev.2.bt2 for the reverse strand? But, how came there is no .rev.3.bt2 or *.rev.4.bt2?
Thanks a lot!
I have not looked at the code, so my assumption of how they operate was indeed based on the file extension generated when indexing :) I'm not sure why sometimes the "reverse" ones are not present.