Generating reference from exome sequence data(bam) to map RNA seq(fastq) from the same individual
Is it possible to generate a reference from exome sequenced/aligned(BWA MEM 0.7.10) to map the RNAseq data from the same individual?The idea is to reduce the mapping bias in RNA seq mapping to the common hg19 reference.
I would appreciate some advice/suggestions on tools at different steps.
Thanks in advance
• 1,837 views
•
link
updated
by
Ram
•
written
by
seqlearner •
•
written
by
seqlearner •
1 answer
If you want to reduce ref bias, I recommend you use a very sensitive aligner such as BBMap rather than mutating the reference in order to align with a less-sensitive aligner. Using a different reference for every experiment can cause major headaches later on when coordinates and base calls no longer match the annotation.
• 0 views
•
link
updated
by
Ram
•
written
by
Brian Bushnell
•
written
by
Brian Bushnell
Log in to answer this question.
More posts like this
-
How to improve the mapping rate of vg giraffe
written by gulin •
-
Substantial reduction in read mapping following VG surject
written by peterdfields •
-
VCF generation of UMI tagged data
written by Saurabh •
-
extending reads before aligning
written by hksk2 •
-
mapping long-reads to a reference library
written by mthm •
-
BWA .alt index file
written by tejaswikoganti
-
Does BWA (BWA-MEM) care about the identity and coverage of the mapped reads?
written by culiu •
-
Identifying contaminants (cow, pig) in oral rinse exome sequencing data
written by AbiB1 •
-
DNA-seq PE150 BWA alignment
written by zhuozhu132 •
-
Paired-End Reads Alignment For Variant Calling ?
written by newDNASeqer
Yes, it is possible. This question has already been discussed several times on this forum. Please use the search feature to find them. At the same time I will keep the question open to get some fresh insights into this topic here.