Coverage of target genes

Question

[gubrins]

Good afternoon, I'm working with target capture sequencing and I would like to know the coverage of my target genes, which I just have in fasta format. I already have my reads aligned into the genome and in bam format and I also already did the SNP calling (with freebayes). So my questions are: - How can I calculate the coverage of my target genes in the bam files I have? - Could be a way to just do the SNP calling of those genes?

I'm aware that with bedtools I could calculate the coverage, but is asking me a bed file that I don't have.

Any type of help will be apreciate it, thanks!

[gubrins]

I'm working with the genome of Podarcis muralis, a lizard species. I'm going to check right now if I can find the bed file, thanks!

[gubrins]

But how can I obtain the genomic coordinates?

[GenoMax]

You can find the genomic coordinates in the GFF file on the genome page for Podarcis here. Look for genome annotation in GFF or tabular format in the box at top of the page.

Answer 1

[GouthamAtla]

What genome are you working with ? For any genome, if you already have genes in fasta, there must be genomic coordinates for those genes in either GTF/GFF or bed format. If you have GTF, convert it to a bed file of exons or gene start end and then calculate coverage using bedtools. Either you can calculate per base pair or per region ( for each gene)