That's perfect, I can as well use directly CDS length for my purposes, thanks !
I'd be interested in retrieving via Biomart the protein length associated to a transcript, as per screenshot below
I have retrieved the list of attributes available from biomart but can't seem to find the right field. Can someone confirm if this info is accessible programatically ? how ?
Thanks,
2 answers
Hi iatz,
You can't retrieve protein length directly from BioMart. However, you can retrieve the CDS length and divide by 3.
gget seq will return the nucleotide (option "gene") or amino acid (option "transcript") sequence and sequence length for any Ensembl transcript ID:
Alternatively, you can also get the nucleotide sequence length from gget info by calculating the difference between the sequence "start" and "end". (Note: Both gget info and gget seq return information regarding the genomic nucleotide sequence, including all exons and introns.)
All gget tools work from the command line and any Python environment, e.g. JupyterLab.
Log in to answer this question.