Base call from BAM file
0
0
Entering edit mode
4.7 years ago
Simo ▴ 50

Hi, I have a BAM file and a reference (let's say 1000 genomes data), I would like to extract from the BAM file all the possible overlapping bases with my reference dataset (not only the variants. Basically I would like to perform a base calling and not a variant calling). This extraction should be done according to a random selection of one read per base. Is there any way to do it?

Thanks

BAM Vcf • 1.6k views
ADD COMMENT
0
Entering edit mode

You could pipe the output of samtools mpileup to a script that randomly selects a single base (no clue why you'd want to do that).

ADD REPLY
0
Entering edit mode

all the possible overlapping bases with my reference dataset

How do you define indel overlap?

I would like to perform a base calling and not a variant calling

How does this differ from setting the flag on your favorite variant caller to also output homozygous reference calls?

This extraction should be done according to a random selection of one read per base.

Why are you downsampling? Why aren't you extracting all read bases that align to any given position and processing that?

ADD REPLY
0
Entering edit mode
  • I do not take into account indels.

  • I don't have any experience on this, that's why I was asking for help.

  • I'm working on a low coverage aDNA.

Thanks

ADD REPLY

Login before adding your answer.

Traffic: 2814 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6