Hello All,
Say I have 2 BAM files from 2 diploid reference genome(simulated maternal and paternal one). Now I need to extract the best reads alignments from both the BAMs and merged into oneBAM file. I have sorted both Bam files.
I would like to know the best ways to proceed, at the end I need to get the coverage of each variants from the BAM file.
Any suggestion would be greatly appreciated.
•
written
by
a.james
May be it would be easier to map the reads on a reference made of both parental genome, so they compete as template for the alignment. From there, selecting the best alignment is straightforward.
Now if you really have to start from the two bam files, here is a possible strategy:
@Carlo Thanks for the detailed answer, I have doubt. For example, the 3rd point comparing based on the MAPQ value and in that case which threshold is good, (sy, MAPQ 30 or 20 etc)? For the 4th point is there a tool to save or extract the reads from the BAM file based on the given list of read_ids?
Thank you
which threshold is good for the MAPQ value ?
You don't really need a treshold, just that one score is bigger than the other, right ?
is there a tool to extract reads from the BAM file based on read_ids?
Yes, for instance picard-tools FilterSamReads with the option FILTER=includeReadList
Log in to answer this question.
