Count the number of exons in a trasncript

Question

[speycast]

Hello,

I'd like to know how to count the number of exons in a given transcript using NCBI eutils? for example NM_020742.3.

thanks very much!

Answer 1

[GenoMax]

Using Entrezdirect:

$ esearch -db nuccore -query "NM_020742.3" | efetch -format ft | grep - w "exon" | wc -l
       6

If you need starts and stops:

$ esearch -db nuccore -query "NM_020742.3" | efetch -format ft | grep -w "exon"
1   159 exon
160 936 exon
937 1089    exon
1090    1275    exon
1276    2065    exon
2066    5706    exon

[speycast]

Thanks so much genomax!!!!