# HG changeset patch # User jjohnson # Date 1511988959 18000 # Node ID a8218b11216f5e9b524769fc295d443dd28a3293 Uploaded diff -r 000000000000 -r a8218b11216f ._ensembl_cdna_translate.py Binary file ._ensembl_cdna_translate.py has changed diff -r 000000000000 -r a8218b11216f ._ensembl_cdna_translate.xml Binary file ._ensembl_cdna_translate.xml has changed diff -r 000000000000 -r a8218b11216f digest.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/digest.py Wed Nov 29 15:55:59 2017 -0500 @@ -0,0 +1,158 @@ +# Copyright 2012 Anton Goloborodko, Lev Levitsky +# +# Licensed under the Apache License, Version 2.0 (the "License"); +# you may not use this file except in compliance with the License. +# You may obtain a copy of the License at +# +# http://www.apache.org/licenses/LICENSE-2.0 +# +# Unless required by applicable law or agreed to in writing, software +# distributed under the License is distributed on an "AS IS" BASIS, +# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +# See the License for the specific language governing permissions and +# limitations under the License. + +import re +from collections import deque +import itertools as it + +def cleave(sequence, rule, missed_cleavages=0, min_length=None): + """Cleaves a polypeptide sequence using a given rule. + + Parameters + ---------- + sequence : str + The sequence of a polypeptide. + + .. note:: + The sequence is expected to be in one-letter uppercase notation. + Otherwise, some of the cleavage rules in :py:data:`expasy_rules` + will not work as expected. + + rule : str or compiled regex + A regular expression describing the site of cleavage. It is recommended + to design the regex so that it matches only the residue whose C-terminal + bond is to be cleaved. All additional requirements should be specified + using `lookaround assertions + `_. + :py:data:`expasy_rules` contains cleavage rules for popular cleavage agents. + missed_cleavages : int, optional + Maximum number of allowed missed cleavages. Defaults to 0. + min_length : int or None, optional + Minimum peptide length. Defaults to :py:const:`None`. + + ..note :: + This checks for string length, which is only correct for one-letter + notation and not for full *modX*. Use :py:func:`length` manually if + you know what you are doing and apply :py:func:`cleave` to *modX* + sequences. + + Returns + ------- + out : set + A set of unique (!) peptides. + + Examples + -------- + >>> cleave('AKAKBK', expasy_rules['trypsin'], 0) == {'AK', 'BK'} + True + >>> cleave('GKGKYKCK', expasy_rules['trypsin'], 2) == \ + {'CK', 'GKYK', 'YKCK', 'GKGK', 'GKYKCK', 'GK', 'GKGKYK', 'YK'} + True + + """ + return set(_cleave(sequence, rule, missed_cleavages, min_length)) + +def _cleave(sequence, rule, missed_cleavages=0, min_length=None): + """Like :py:func:`cleave`, but the result is a list. Refer to + :py:func:`cleave` for explanation of parameters. + """ + peptides = [] + ml = missed_cleavages+2 + trange = range(ml) + cleavage_sites = deque([0], maxlen=ml) + cl = 1 + for i in it.chain([x.end() for x in re.finditer(rule, sequence)], + [None]): + cleavage_sites.append(i) + if cl < ml: + cl += 1 + for j in trange[:cl-1]: + seq = sequence[cleavage_sites[j]:cleavage_sites[-1]] + if seq: + if min_length is None or len(seq) >= min_length: + peptides.append(seq) + return peptides + +def num_sites(sequence, rule, **kwargs): + """Count the number of sites where `sequence` can be cleaved using + the given `rule` (e.g. number of miscleavages for a peptide). + + Parameters + ---------- + sequence : str + The sequence of a polypeptide. + rule : str or compiled regex + A regular expression describing the site of cleavage. It is recommended + to design the regex so that it matches only the residue whose C-terminal + bond is to be cleaved. All additional requirements should be specified + using `lookaround assertions + `_. + labels : list, optional + A list of allowed labels for amino acids and terminal modifications. + + Returns + ------- + out : int + Number of cleavage sites. + """ + return len(_cleave(sequence, rule, **kwargs)) - 1 + +expasy_rules = { + 'arg-c': r'R', + 'asp-n': r'\w(?=D)', + 'bnps-skatole' : r'W', + 'caspase 1': r'(?<=[FWYL]\w[HAT])D(?=[^PEDQKR])', + 'caspase 2': r'(?<=DVA)D(?=[^PEDQKR])', + 'caspase 3': r'(?<=DMQ)D(?=[^PEDQKR])', + 'caspase 4': r'(?<=LEV)D(?=[^PEDQKR])', + 'caspase 5': r'(?<=[LW]EH)D', + 'caspase 6': r'(?<=VE[HI])D(?=[^PEDQKR])', + 'caspase 7': r'(?<=DEV)D(?=[^PEDQKR])', + 'caspase 8': r'(?<=[IL]ET)D(?=[^PEDQKR])', + 'caspase 9': r'(?<=LEH)D', + 'caspase 10': r'(?<=IEA)D', + 'chymotrypsin high specificity' : r'([FY](?=[^P]))|(W(?=[^MP]))', + 'chymotrypsin low specificity': + r'([FLY](?=[^P]))|(W(?=[^MP]))|(M(?=[^PY]))|(H(?=[^DMPW]))', + 'clostripain': r'R', + 'cnbr': r'M', + 'enterokinase': r'(?<=[DE]{3})K', + 'factor xa': r'(?<=[AFGILTVM][DE]G)R', + 'formic acid': r'D', + 'glutamyl endopeptidase': r'E', + 'granzyme b': r'(?<=IEP)D', + 'hydroxylamine': r'N(?=G)', + 'iodosobenzoic acid': r'W', + 'lysc': r'K', + 'ntcb': r'\w(?=C)', + 'pepsin ph1.3': r'((?<=[^HKR][^P])[^R](?=[FLWY][^P]))|' + r'((?<=[^HKR][^P])[FLWY](?=\w[^P]))', + 'pepsin ph2.0': r'((?<=[^HKR][^P])[^R](?=[FL][^P]))|' + r'((?<=[^HKR][^P])[FL](?=\w[^P]))', + 'proline endopeptidase': r'(?<=[HKR])P(?=[^P])', + 'proteinase k': r'[AEFILTVWY]', + 'staphylococcal peptidase i': r'(?<=[^E])E', + 'thermolysin': r'[^DE](?=[AFILMV])', + 'thrombin': r'((?<=G)R(?=G))|' + r'((?<=[AFGILTVM][AFGILTVWA]P)R(?=[^DE][^DE]))', + 'trypsin': r'([KR](?=[^P]))|((?<=W)K(?=P))|((?<=M)R(?=P))' + } +""" +This dict contains regular expressions for cleavage rules of the most +popular proteolytic enzymes. The rules were taken from the +`PeptideCutter tool +`_ +at Expasy. +""" + diff -r 000000000000 -r a8218b11216f ensembl_cdna_translate.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/ensembl_cdna_translate.py Wed Nov 29 15:55:59 2017 -0500 @@ -0,0 +1,426 @@ +#!/usr/bin/env python +""" +# +#------------------------------------------------------------------------------ +# University of Minnesota +# Copyright 2017, Regents of the University of Minnesota +#------------------------------------------------------------------------------ +# Author: +# +# James E Johnson +# +#------------------------------------------------------------------------------ +""" + +import argparse +import sys +from time import sleep + +from Bio.Seq import translate + +import requests + +import digest + + +server = "https://rest.ensembl.org" +ext = "/info/assembly/homo_sapiens?" +max_region = 5000000 + + +def ensembl_rest(ext, headers): + if True: print >> sys.stderr, "%s" % ext + r = requests.get(server+ext, headers=headers) + if r.status_code == 429: + print >> sys.stderr, "response headers: %s\n" % r.headers + if 'Retry-After' in r.headers: + sleep(r.headers['Retry-After']) + r = requests.get(server+ext, headers=headers) + if not r.ok: + r.raise_for_status() + return r + + +def get_species(): + results = dict() + ext = "/info/species" + req_header = {"Content-Type": "application/json"} + r = ensembl_rest(ext, req_header) + for species in r.json()['species']: + results[species['name']] = species + print >> sys.stdout,\ + "%s\t%s\t%s\t%s\t%s"\ + % (species['name'], species['common_name'], + species['display_name'], + species['strain'], + species['taxon_id']) + return results + + +def get_biotypes(species): + biotypes = [] + ext = "/info/biotypes/%s?" % species + req_header = {"Content-Type": "application/json"} + r = ensembl_rest(ext, req_header) + for entry in r.json(): + if 'biotype' in entry: + biotypes.append(entry['biotype']) + return biotypes + + +def get_toplevel(species): + coord_systems = dict() + ext = "/info/assembly/%s?" % species + req_header = {"Content-Type": "application/json"} + r = ensembl_rest(ext, req_header) + toplevel = r.json() + for seq in toplevel['top_level_region']: + if seq['coord_system'] not in coord_systems: + coord_systems[seq['coord_system']] = dict() + coord_system = coord_systems[seq['coord_system']] + coord_system[seq['name']] = int(seq['length']) + return coord_systems + + +def get_transcripts_bed(species, refseq, start, length,params=None): + bed = [] + param = params if params else '' + req_header = {"Content-Type": "text/x-bed"} + regions = range(start, length, max_region) + if not regions or regions[-1] < length: + regions.append(length) + for end in regions[1:]: + ext = "/overlap/region/%s/%s:%d-%d?feature=transcript;%s"\ + % (species, refseq, start, end, param) + start = end + 1 + r = ensembl_rest(ext, req_header) + if r.text: + bed += r.text.splitlines() + return bed + + +def get_seq(id, seqtype,params=None): + param = params if params else '' + ext = "/sequence/id/%s?type=%s;%s" % (id, seqtype,param) + req_header = {"Content-Type": "text/plain"} + r = ensembl_rest(ext, req_header) + return r.text + + +def get_cdna(id,params=None): + return get_seq(id, 'cdna',params=params) + + +def get_cds(id,params=None): + return get_seq(id, 'cds',params=params) + + +def bed_from_line(line): + fields = line.rstrip('\r\n').split('\t') + (chrom, chromStart, chromEnd, name, score, strand, + thickStart, thickEnd, itemRgb, + blockCount, blockSizes, blockStarts) = fields[0:12] + bed_entry = BedEntry(chrom=chrom, chromStart=chromStart, chromEnd=chromEnd, + name=name, score=score, strand=strand, + thickStart=thickStart, thickEnd=thickEnd, + itemRgb=itemRgb, + blockCount=blockCount, + blockSizes=blockSizes.rstrip(','), + blockStarts=blockStarts.rstrip(',')) + return bed_entry + + +class BedEntry(object): + def __init__(self, chrom=None, chromStart=None, chromEnd=None, + name=None, score=None, strand=None, + thickStart=None, thickEnd=None, itemRgb=None, + blockCount=None, blockSizes=None, blockStarts=None): + self.chrom = chrom + self.chromStart = int(chromStart) + self.chromEnd = int(chromEnd) + self.name = name + self.score = int(score) if score is not None else 0 + self.strand = '-' if str(strand).startswith('-') else '+' + self.thickStart = int(thickStart) if thickStart else self.chromStart + self.thickEnd = int(thickEnd) if thickEnd else self.chromEnd + self.itemRgb = str(itemRgb) if itemRgb is not None else r'100,100,100' + self.blockCount = int(blockCount) + if isinstance(blockSizes, str) or isinstance(blockSizes, unicode): + self.blockSizes = [int(x) for x in blockSizes.split(',')] + elif isinstance(blockSizes, list): + self.blockSizes = [int(x) for x in blockSizes] + else: + self.blockSizes = blockSizes + if isinstance(blockStarts, str) or isinstance(blockSizes, unicode): + self.blockStarts = [int(x) for x in blockStarts.split(',')] + elif isinstance(blockStarts, list): + self.blockStarts = [int(x) for x in blockStarts] + else: + self.blockStarts = blockStarts + self.seq = None + self.pep = None + + def __str__(self): + return '%s\t%d\t%d\t%s\t%d\t%s\t%d\t%d\t%s\t%d\t%s\t%s' % ( + self.chrom, self.chromStart, self.chromEnd, + self.name, self.score, self.strand, + self.thickStart, self.thickEnd, str(self.itemRgb), self.blockCount, + ','.join([str(x) for x in self.blockSizes]), + ','.join([str(x) for x in self.blockStarts])) + + # (start, end) + def get_subrange(self, tstart, tstop, debug=False): + chromStart = self.chromStart + chromEnd = self.chromEnd + if debug: + print >> sys.stderr, "%s" % (str(self)) + r = range(self.blockCount) + if self.strand == '-': + r.reverse() + bStart = 0 + bEnd = 0 + for x in r: + bEnd = bStart + self.blockSizes[x] + if bStart <= tstart < bEnd: + if self.strand == '+': + chromStart = self.chromStart + self.blockStarts[x] +\ + (tstart - bStart) + else: + chromEnd = self.chromStart + self.blockStarts[x] +\ + self.blockSizes[x] - (tstart - bStart) + if bStart <= tstop < bEnd: + if self.strand == '+': + chromEnd = self.chromStart + self.blockStarts[x] +\ + (tstop - bStart) + else: + chromStart = self.chromStart + self.blockStarts[x] +\ + self.blockSizes[x] - (tstop - bStart) + if debug: + print >> sys.stderr,\ + "%3d %s\t%d\t%d\t%d\t%d\t%d\t%d"\ + % (x, self.strand, bStart, bEnd, + tstart, tstop, chromStart, chromEnd) + bStart += self.blockSizes[x] + return(chromStart, chromEnd) + + # get the blocks for sub range + def get_blocks(self, chromStart, chromEnd): + tblockCount = 0 + tblockSizes = [] + tblockStarts = [] + for x in range(self.blockCount): + bStart = self.chromStart + self.blockStarts[x] + bEnd = bStart + self.blockSizes[x] + if bStart > chromEnd: + break + if bEnd < chromStart: + continue + cStart = max(chromStart, bStart) + tblockStarts.append(cStart - chromStart) + tblockSizes.append(min(chromEnd, bEnd) - cStart) + tblockCount += 1 + return (tblockCount, tblockSizes, tblockStarts) + + def trim(self, tstart, tstop, debug=False): + (tchromStart, tchromEnd) =\ + self.get_subrange(tstart, tstop, debug=debug) + (tblockCount, tblockSizes, tblockStarts) =\ + self.get_blocks(tchromStart, tchromEnd) + tbed = BedEntry( + chrom=self.chrom, chromStart=tchromStart, chromEnd=tchromEnd, + name=self.name, score=self.score, strand=self.strand, + thickStart=tchromStart, thickEnd=tchromEnd, itemRgb=self.itemRgb, + blockCount=tblockCount, + blockSizes=tblockSizes, blockStarts=tblockStarts) + if self.seq: + ts = tchromStart-self.chromStart + te = tchromEnd - tchromStart + ts + tbed.seq = self.seq[ts:te] + return tbed + + +def __main__(): + parser = argparse.ArgumentParser( + description='Retrieve Ensembl cDNAs and three frame translate') + parser.add_argument( + '-s', '--species', default='human', + help='Ensembl Species to retrieve') + parser.add_argument( + '-B', '--biotypes', action='append', default=[], + help='Restrict Ensembl biotypes to retrieve') + parser.add_argument( + '-i', '--input', default=None, + help='Use BED instead of retrieving cDNA from ensembl (-) for stdin') + parser.add_argument( + '-t', '--transcripts', default=None, + help='Path to output cDNA transcripts.bed (-) for stdout') + parser.add_argument( + '-r', '--raw', action='store_true', + help='Report transcript exacty as returned from Ensembl') + parser.add_argument( + '-f', '--fasta', default=None, + help='Path to output translations.fasta') + parser.add_argument( + '-b', '--bed', default=None, + help='Path to output translations.bed') + parser.add_argument( + '-m', '--min_length', type=int, default=7, + help='Minimum length of protein translation to report') + parser.add_argument( + '-e', '--enzyme', default=None, + help='Digest translation with enzyme') + parser.add_argument( + '-a', '--all', action='store_true', + help='Include reference protein translations') + parser.add_argument('-v', '--verbose', action='store_true', help='Verbose') + parser.add_argument('-d', '--debug', action='store_true', help='Debug') + args = parser.parse_args() + # print >> sys.stderr, "args: %s" % args + species = args.species + input_rdr = None + if args.input is not None: + input_rdr = open(args.input, 'r') if args.input != '-' else sys.stdin + tx_wtr = None + if args.transcripts is not None: + tx_wtr = open(args.transcripts, 'w')\ + if args.transcripts != '-' else sys.stdout + fa_wtr = open(args.fasta, 'w') if args.fasta is not None else None + bed_wtr = open(args.bed, 'w') if args.bed is not None else None + + enzyme = digest.expasy_rules.get(args.enzyme,args.enzyme) + + # print >> sys.stderr, "args biotypes: %s" % args.biotypes + biotypea = ['biotype=%s' % bt.strip() for biotype in args.biotypes for bt in biotype.split(',')] + # print >> sys.stderr, "args biotypes: %s" % biotypea + biotypes = ';'.join(['biotype=%s' % bt.strip() for biotype in args.biotypes for bt in biotype.split(',') if bt.strip()]) + # print >> sys.stderr, "biotypes: %s" % biotypes + + translations = dict() # start : end : seq + + def unique_prot(tbed, seq): + if tbed.chromStart not in translations: + translations[tbed.chromStart] = dict() + translations[tbed.chromStart][tbed.chromEnd] = [] + translations[tbed.chromStart][tbed.chromEnd].append(seq) + elif tbed.chromEnd not in translations[tbed.chromStart]: + translations[tbed.chromStart][tbed.chromEnd] = [] + translations[tbed.chromStart][tbed.chromEnd].append(seq) + elif seq not in translations[tbed.chromStart][tbed.chromEnd]: + translations[tbed.chromStart][tbed.chromEnd].append(seq) + else: + return False + return True + + def translate_bed(bed): + translate_count = 0 + if any([fa_wtr, bed_wtr]): + transcript_id = bed.name + refprot = None + if not args.all: + try: + cds = get_cds(transcript_id) + if len(cds) % 3 != 0: + cds = cds[:-(len(cds) % 3)] + refprot = translate(cds) if cds else None + except: + refprot = None + cdna = get_cdna(transcript_id) + cdna_len = len(cdna) + for offset in range(3): + seqend = cdna_len - (cdna_len - offset) % 3 + aaseq = translate(cdna[offset:seqend]) + aa_start = 0 + while aa_start < len(aaseq): + aa_end = aaseq.find('*', aa_start) + if aa_end < 0: + aa_end = len(aaseq) + prot = aaseq[aa_start:aa_end] + if enzyme and refprot: + frags = digest._cleave(prot,enzyme) + for frag in reversed(frags): + if frag in refprot: + prot = prot[:prot.rfind(frag)] + else: + break + if len(prot) < args.min_length: + pass + elif refprot and prot in refprot: + pass + else: + tstart = aa_start*3+offset + tend = aa_end*3+offset + prot_acc = "%s_%d_%d" % (transcript_id, tstart, tend) + tbed = bed.trim(tstart, tend) + if args.all or unique_prot(tbed, prot): + translate_count += 1 + tbed.name = prot_acc + bed_wtr.write("%s\t%s\n" % (str(tbed), prot)) + bed_wtr.flush() + fa_id = ">%s\n" % (prot_acc) + fa_wtr.write(fa_id) + fa_wtr.write(prot) + fa_wtr.write("\n") + fa_wtr.flush() + aa_start = aa_end + 1 + return translate_count + + if input_rdr: + translation_count = 0 + for i, bedline in enumerate(input_rdr): + try: + bed = bed_from_line(bedline) + translation_count += translate_bed(bed) + except: + print >> sys.stderr, "BED format error: %s\n" % bedline + if args.debug or (args.verbose and any([fa_wtr, bed_wtr])): + print >> sys.stderr,\ + "%s\tcDNA translations:%d" % (species, translation_count) + else: + coord_systems = get_toplevel(species) + if 'chromosome' in coord_systems: + for ref in sorted(coord_systems['chromosome'].keys()): + length = coord_systems['chromosome'][ref] + if not any([tx_wtr, fa_wtr, bed_wtr]): + print >> sys.stderr,\ + "%s\t%s\tlength: %d" % (species, ref, length) + continue + if args.debug: + print >> sys.stderr,\ + "Retrieving transcripts: %s\t%s\tlength: %d"\ + % (species, ref, length) + translation_count = 0 + start = 0 + regions = range(start, length, max_region) + if not regions or regions[-1] < length: + regions.append(length) + for end in regions[1:]: + bedlines = get_transcripts_bed(species, ref, start, end, params=biotypes) + if args.verbose or args.debug: + print >> sys.stderr,\ + "%s\t%s\tstart: %d\tend: %d\tcDNA transcripts:%d"\ + % (species, ref, start, end, len(bedlines)) + # start, end, seq + for i, bedline in enumerate(bedlines): + try: + bed = bed_from_line(bedline)\ + if any([not args.raw, fa_wtr, bed_wtr])\ + else None + if tx_wtr: + tx_wtr.write(bedline if args.raw else str(bed)) + tx_wtr.write("\n") + tx_wtr.flush() + if bed: + translation_count += translate_bed(bed) + except Exception as e: + print >> sys.stderr,\ + "BED error (%s) : %s\n" % (e, bedline) + start = end + 1 + + if args.debug or (args.verbose and any([fa_wtr, bed_wtr])): + print >> sys.stderr,\ + "%s\t%s\tlength: %d\tcDNA translations:%d"\ + % (species, ref, length, translation_count) + + +if __name__ == "__main__": + __main__() diff -r 000000000000 -r a8218b11216f ensembl_cdna_translate.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/ensembl_cdna_translate.xml Wed Nov 29 15:55:59 2017 -0500 @@ -0,0 +1,248 @@ + + using Ensembl REST API + + requests-cache + biopython + + + + + = 0: + --transcripts + #if str($output_choice).find('translation') >= 0: + '-' | tee '$transcript_bed' + #else + '$transcript_bed' + #end if + #elif str($output_choice).find('translation') >= 0: + --transcripts '-' + #end if + #end if + #if str($output_choice).find('translation') >= 0: + | python '$__tool_directory__/ensembl_cdna_translate.py' -i '-' + #if $input and str($output_choice).find('transcript_bed') >= 0: + --transcripts '$transcript_bed + #end if + #if str($output_choice).find('translation_bed') >= 0: + --bed '$translation_bed' + #end if + #if str($output_choice).find('translation_fasta') >= 0: + --fasta '$translation_fasta' + #end if + #if $enzyme: + --enzyme '$enzyme' + #end if + #end if + + ]]> + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 'transcript_bed' in output_choice + + + 'translation_bed' in output_choice + + + 'translation_fasta' in output_choice + + + + + + + + + + + + + + + + + + + + + + 10.1093/bioinformatics/btu613 + 10.1093/nar/gku1010 + + diff -r 000000000000 -r a8218b11216f test-data/human_transcripts.bed --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/human_transcripts.bed Wed Nov 29 15:55:59 2017 -0500 @@ -0,0 +1,19 @@ +chr1 14403 29570 ENST00000488147 1000 - 14402 14402 0,0,0 11 98,34,152,159,198,136,137,147,99,154,37, 0,601,1392,2203,2454,2829,3202,3511,3864,10334,15130, WASH7P-201 none none -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1, unprocessed_pseudogene ENSG00000227232 WASH7P unprocessed_pseudogene +chr1 29553 31097 ENST00000473358 1000 + 29552 29552 0,0,0 3 486,104,122, 0,1010,1422, MIR1302-2HG-202 none none -1,-1,-1, lincRNA ENSG00000243485 MIR1302-2HG lincRNA +chr1 30266 31109 ENST00000469289 1000 + 30265 30265 0,0,0 2 401,134, 0,709, MIR1302-2HG-201 none none -1,-1, lincRNA ENSG00000243485 MIR1302-2HG lincRNA +chr1 30365 30503 ENST00000607096 1000 + 30364 30364 0,0,0 1 138, 0, MIR1302-2-201 none none -1, miRNA ENSG00000284332 MIR1302-2 miRNA +chr1 34553 36081 ENST00000417324 1000 - 34552 34552 0,0,0 3 621,205,361, 0,723,1167, FAM138A-201 none none -1,-1,-1, lincRNA ENSG00000237613 FAM138A lincRNA +chr1 35244 36073 ENST00000461467 1000 - 35243 35243 0,0,0 2 237,353, 0,476, FAM138A-202 none none -1,-1, lincRNA ENSG00000237613 FAM138A lincRNA +chr1 52472 53312 ENST00000606857 1000 + 52471 52471 0,0,0 1 840, 0, AL627309.6-201 none none -1, unprocessed_pseudogene ENSG00000268020 AL627309.6 unprocessed_pseudogene +chr1 57597 64116 ENST00000642116 1000 + 57596 57596 0,0,0 3 56,157,1201, 0,1102,5318, OR4G11P-202 none none -1,-1,-1, processed_transcript ENSG00000240361 OR4G11P transcribed_unprocessed_pseudogene +chr1 62948 63887 ENST00000492842 1000 + 62947 62947 0,0,0 1 939, 0, OR4G11P-201 none none -1, transcribed_unprocessed_pseudogene ENSG00000240361 OR4G11P transcribed_unprocessed_pseudogene +chr1 65418 71585 ENST00000641515 1000 + 69090 70008 0,0,0 3 15,54,2549, 0,101,3618, OR4F5-202 cmpl cmpl -1,-1,0, protein_coding ENSG00000186092 OR4F5 protein_coding +chr1 69054 70108 ENST00000335137 1000 + 69090 70008 0,0,0 1 1054, 0, OR4F5-201 cmpl cmpl 0, protein_coding ENSG00000186092 OR4F5 protein_coding +chr1 131024 134836 ENST00000442987 1000 + 131023 131023 0,0,0 1 3812, 0, CICP27-201 none none -1, processed_pseudogene ENSG00000233750 CICP27 processed_pseudogene +chr1 139789 140339 ENST00000493797 1000 - 139788 139788 0,0,0 2 58,265, 0,285, AL627309.2-201 none none -1,-1, antisense_RNA ENSG00000239906 AL627309.2 antisense_RNA +chr1 157783 157887 ENST00000410691 1000 - 157782 157782 0,0,0 1 104, 0, RNU6-1100P-201 none none -1, snRNA ENSG00000222623 RNU6-1100P snRNA +chr1 187890 187958 ENST00000612080 1000 - 187889 187889 0,0,0 1 68, 0, MIR6859-2-201 none none -1, miRNA ENSG00000273874 MIR6859-2 miRNA +chr1 263014 297502 ENST00000424587 1000 - 263013 263013 0,0,0 4 5190,150,105,158, 0,5652,26251,34330, AP006222.1-206 none none -1,-1,-1,-1, processed_transcript ENSG00000228463 AP006222.1 transcribed_processed_pseudogene +chr1 347981 348366 ENST00000458203 1000 - 347980 347980 0,0,0 1 385, 0, RPL23AP24-201 none none -1, processed_pseudogene ENSG00000236679 RPL23AP24 processed_pseudogene +chr1 439869 440232 ENST00000437905 1000 + 439868 439868 0,0,0 1 363, 0, WBP1LP7-201 none none -1, processed_pseudogene ENSG00000269732 WBP1LP7 processed_pseudogene +chr1 450702 451697 ENST00000426406 1000 - 450739 451678 0,0,0 1 995, 0, OR4F29-201 cmpl cmpl 0, protein_coding ENSG00000284733 OR4F29 protein_coding