diff --git a/pyGeno/importation/Genomes.py b/pyGeno/importation/Genomes.py
index <HASH>..<HASH> 100644
--- a/pyGeno/importation/Genomes.py
+++ b/pyGeno/importation/Genomes.py
@@ -337,7 +337,18 @@ def _importGenomeObjects(gtfFilePath, chroSet, genome, batchSize, verbose = 0) :
                     protId = None
                     if verbose > 2 :
                         printf('Warning: no protein_id found in line %s' % gtf[i])
-                
+
+                # Store selenocysteine positions in transcript
+                if regionType == 'Selenocysteine':
+
+                    if store.transcripts[transId].selenocysteine is None:
+                        positions = list()
+                    else:
+                        positions = store.transcripts[transId].selenocysteine
+
+                    positions.append(start)
+                    store.transcripts[transId].set(selenocysteine=positions)
+                        
                 if protId is not None and protId not in store.proteins :
                     if verbose > 1 :
                         printf('\t\tProtein %s...' % (protId))
Stores selenocysteine positions in transcript