This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Exon falls within a gap in the assembly

Downloaded published gene annotation and after extracting fasta sequences noticed that one exon is only Ns. Intersected annotation with gaps in the assembly and it appears that one exon (short, around 20bp ) falls within a gap.

How to be sure that this exons really falls within a gap before I contact authors?

Edited

Exon starts 10bp from the end of the gap.

gap
• 1,908 views • link
updated by Ram Ram • written by PoGibas PoGibas

1 answer

To make sure that annotations match the genome I would extract actual sequences that correspond to these annotations with say bedtools getfasta and look at not just the Ns but also start/stop codons and also that the annotated genes actually get transcribed to proteins that you expect to see.

• 3 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this