BWA, fastq reads mapped to multi ref fasta. ALL vs ALL
I have a fasta with two or more references. I want to map ALL my reads (fq) to both references.
What happens now, bwa chooses which is the best reference for each paired read, and it's discarded for the other reference.
I want exactly the same behavior than mapping separately using two different fasta files (with one ref each).
Is it possible?
Cheers
• 984 views
•
link
updated
by
colindaven
•
written
by
erneheay •
•
written
by
erneheay •
1 answer
Just map to each fasta ref completely separately , i.e. split the fasta before alignment ? BTW, alignment is NOT assembly, see your tags.
• 0 views
•
link
written
by
colindaven
Log in to answer this question.
More posts like this
-
Bowtie - multiple alignments and seed length vs read length in miRNA mapping
written by MH85 •
-
How to not align reverse complementary reads?
written by dcheng1 •
-
Merging BAM vs concatenate FASTQ
written by nhaus •
-
BWA mapping to multiple reference genome
written by dsklsdoiwld •
-
Running STAR aligner with paired-end and single-end reads simultaneously
written by nanoide
-
beginner question BWA program
written by reza
-
get a consensus out of a Sorted bam file
written by lait
-
BWA barcode trimming and labeling
written by igor
-
some problem with my RRBS data when dealt with bismark
written by hua.peng1314
-
Using Bwa Sampe Using Illumna Pair End Reads
written by Juliofdiaz
