Coverage of target genes
1
0
Entering edit mode
12 months ago
gubrins ▴ 70

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!

coverage target capture target capturing BAM VCF • 263 views
ADD COMMENT
0
Entering edit mode

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!

ADD REPLY
0
Entering edit mode

But how can I obtain the genomic coordinates?

ADD REPLY
1
Entering edit mode

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.

ADD REPLY
1
Entering edit mode
12 months ago

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)

ADD COMMENT

Login before adding your answer.

Traffic: 1286 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