Features of gene on different strands in pool-seq analysis
Entering edit mode
3.8 years ago
Dracaena ▴ 50

Hi everyone,

Currently I am working with pool-seq data and I am specifically interested in SNPs and underlying genes with significant differentiation. I was able to find several thousands of SNPs that show significant differentiation after correcting for multiple testing. Now I want to look at the genes underlying these SNPs. Since I am working with pool-seq data I am working with the Popoolation(2) software package. This package contains a script to create a genewise-sync file which enables you to check for genes that are particularly enriched with high differentiation SNPs. This is the point where I am stuck. When trying to create the genewise sync file I get the following error: "all features of a gene have to be on the same strand; Problem for BG8803808". When I check the gene, I see that a CDS and an exon are indeed annotated to the opposite strand from the rest of the features. So here's my question: is it common for an annotated (draft) genome to have features of the same gene on both strands? And how can I solve this issue (maybe write a costum script to find the genes using my set of SNPs)?

Thanks in advance,

popoolation pool-seq annotation genomics • 731 views

Login before adding your answer.

Traffic: 1366 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6