This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to add new coordinates to the existing .gtf annotation file

I have a few sequences that are not annotated in the reference genome. Using BLAST and STAR I was able to get the coordinates of all those sequences. Using IGV, I could visualize that my transcriptomic data has reads mapped to these coordinates. How can I quantify them?? Do I need to customize my gtf file? If yes, how can I do that?

gtf

1 answer

If it's just a couple regions, you can make a simple annotation file and use featureCounts. If it's anything more, you can make a separate GTF file and cat them together on the command line.

Log in to answer this question.