Question: Variant and Genotype calling
gravatar for dorokhin.alex.v
16 months ago by
dorokhin.alex.v0 wrote:


1) I have a bam-file (sorted, with duplicates removed using samtools). It is important that the reads do not completely cover the reference genome.

2) I also have a vcf-file obtained from this bam-file (using bcftools tools) with options (consensus calling)

3) For the reference genome there is a gtf-file with annotation (from ensembl).

I would like to find out some obvious things:

a) What genes have been sequenced in this case (number and name of genes)?

b) I want to analyze the variants for these genes. (How different are these versions from the reference ones, are there any problems - stop codons, deletions, inversions, etc.).

Tell me, please examples of specific commands, what tools is better to achieve this (bedtools? bedops?)

It is worth noting that my data is not on the human genome, so I want to know the solution in general, and not with the help of specific for humans.


snp gene • 381 views
ADD COMMENTlink modified 16 months ago • written 16 months ago by dorokhin.alex.v0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1857 users visited in the last hour