Mercurial > repos > devteam > cuffcompare
changeset 4:806c27c97df7 draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tool_collections/cufflinks/cuffcompare commit a0b0845a9d1b3e7ecdeacd1e606133617e3918bd"
author | iuc |
---|---|
date | Tue, 16 Jun 2020 16:56:55 +0000 |
parents | 5aac9b9d6f2a |
children | a739ea0fef9d |
files | cuff_macros.xml cuffcompare_wrapper.py cuffcompare_wrapper.xml test-data/cuffcompare_out7.txt test-data/output.combined.gtf test-data/output.cuffcompare_in1.gtf.refmap test-data/output.cuffcompare_in1.gtf.tmap test-data/output.cuffcompare_in2.gtf.refmap test-data/output.cuffcompare_in2.gtf.tmap test-data/output.loci test-data/output.stats test-data/output.tracking |
diffstat | 12 files changed, 459 insertions(+), 160 deletions(-) [+] |
line wrap: on
line diff
--- a/cuff_macros.xml Sun Feb 19 12:11:14 2017 -0500 +++ b/cuff_macros.xml Tue Jun 16 16:56:55 2020 +0000 @@ -66,20 +66,26 @@ </token> <xml name="cufflinks_gtf_inputs"> <param format="gtf" name="inputs" type="data" label="GTF file(s) produced by Cufflinks" help="" multiple="true" /> - <repeat name="additional_inputs" title="Additional GTF Inputs (Lists)"> - <param format="gtf" name="additional_inputs" type="data_collection" label="GTF file(s) produced by Cufflinks" help="" /> - </repeat> </xml> + <token name="@CUFFLINKS_LINK_GTF_INPUTS@"><![CDATA[ + ## Inputs. + #for $i, $input_file in enumerate($inputs): + ln -s '${input_file}' input_$i && + #end for + ]]></token> <token name="@CUFFLINKS_GTF_INPUTS@"> ## Inputs. - #for $input_file in $inputs: - '${input_file}' - #end for - #for $additional_input in $additional_inputs: - #for $input_file in $additional_input.additional_inputs: - '${input_file}' - #end for + #for $i, $input_file in enumerate($inputs): + 'input_$i' #end for </token> <token name="@HAS_MULTIPLE_INPUTS@">getattr(inputs, "__len__", [].__len__)() >= 2</token> + + <xml name="citations"> + <citations> + <citation type="doi">10.1038/nbt.1621</citation> + <yield/> + </citations> + </xml> + </macros>
--- a/cuffcompare_wrapper.py Sun Feb 19 12:11:14 2017 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,117 +0,0 @@ -#!/usr/bin/env python - -# Supports Cuffcompare versions v1.3.0 and newer. - -import optparse -import os -import shutil -import subprocess -import sys -import tempfile - - -def stop_err(msg): - sys.stderr.write('%s\n' % msg) - sys.exit() - - -def __main__(): - parser = optparse.OptionParser() - parser.add_option('-r', dest='ref_annotation', help='An optional "reference" annotation GTF. Each sample is matched against this file, and sample isoforms are tagged as overlapping, matching, or novel where appropriate. See the refmap and tmap output file descriptions below.') - parser.add_option('-R', action="store_true", dest='ignore_nonoverlap_reference', help='If -r was specified, this option causes cuffcompare to ignore reference transcripts that are not overlapped by any transcript in one of cuff1.gtf,...,cuffN.gtf. Useful for ignoring annotated transcripts that are not present in your RNA-Seq samples and thus adjusting the "sensitivity" calculation in the accuracy report written in the transcripts accuracy file') - parser.add_option('-Q', action="store_true", dest='ignore_nonoverlap_transfrag', help='If -r was specified, this option causes cuffcompare to consider only the input transcripts that overlap any of the reference transcripts (Sp correction); Warning: this will discard all "novel" loci!)') - - parser.add_option('-s', dest='use_seq_data', action="store_true", help='Causes cuffcompare to look into for fasta files with the underlying genomic sequences (one file per contig) against which your reads were aligned for some optional classification functions. For example, Cufflinks transcripts consisting mostly of lower-case bases are classified as repeats. Note that <seq_dir> must contain one fasta file per reference chromosome, and each file must be named after the chromosome, and have a .fa or .fasta extension.') - - parser.add_option('-M', action="store_true", dest='discard_single_exon_all', help='discard (ignore) single-exon transfrags and reference transcript') - parser.add_option('-N', action="store_true", dest='discard_single_exon_ref', help='discard (ignore) single-exon reference transcripts') - parser.add_option('-e', dest='max_dist_exon', help='Max. Distance for assessing exon accuracy" help="max. distance (range) allowed from free ends of terminal exons of reference transcripts when assessing exon accuracy. Default: 100') - parser.add_option('-d', dest='max_dist_group', help='Max.Distance for transcript grouping" help="max. distance (range) for grouping transcript start sites. Default: 100') - parser.add_option('-F', action="store_true", dest='discard_redundant_intron_transfrags', help='Discard intron-redundant transfrags if they share the 5-prime end (if they differ only at the 3-prime end)') - - # Wrapper / Galaxy options. - parser.add_option('', '--index', dest='index', help='The path of the reference genome') - parser.add_option('', '--ref_file', dest='ref_file', help='The reference dataset from the history') - - # Outputs. - parser.add_option('', '--combined-transcripts', dest='combined_transcripts') - - (options, args) = parser.parse_args() - - # Set/link to sequence file. - if options.use_seq_data: - if options.ref_file: - # Sequence data from history. - # Create symbolic link to ref_file so that index will be created in working directory. - seq_path = "ref.fa" - os.symlink(options.ref_file, seq_path) - else: - if not os.path.exists(options.index): - stop_err('Reference genome %s not present, request it by reporting this error.' % options.index) - seq_path = options.index - - # Build command. - - # Base. - cmd = "cuffcompare -o cc_output " - - # Add options. - if options.ref_annotation: - cmd += " -r %s " % options.ref_annotation - if options.ignore_nonoverlap_reference: - cmd += " -R " - if options.ignore_nonoverlap_transfrag: - cmd += " -Q " - if options.use_seq_data: - cmd += " -s %s " % seq_path - if options.discard_single_exon_all: - cmd += " -M " - if options.discard_single_exon_ref: - cmd += " -N " - if options.max_dist_exon: - cmd += " -e %i " % int(options.max_dist_exon) - if options.max_dist_group: - cmd += " -d %i " % int(options.max_dist_group) - if options.discard_redundant_intron_transfrags: - cmd += " -F " - # Add input files. - - # Need to symlink inputs so that output files are written to the current directory. - for i, arg in enumerate(args): - input_file_name = "./input%i" % (i + 1) - os.symlink(arg, input_file_name) - cmd += "%s " % input_file_name - - # Run command. - try: - with tempfile.NamedTemporaryFile(dir=".") as tmp_stderr: - returncode = subprocess.call(args=cmd, stderr=tmp_stderr, shell=True) - - # Error checking. - if returncode != 0: - # Get stderr, allowing for case where it's very large. - buffsize = 1048576 - stderr = '' - with open(tmp_stderr.name) as tmp_stderr2: - try: - while True: - stderr += tmp_stderr2.read(buffsize) - if not stderr or len(stderr) % buffsize != 0: - break - except OverflowError: - pass - raise Exception(stderr) - - # Copy outputs. - shutil.move("cc_output.combined.gtf", options.combined_transcripts) - - # check that there are results in the output file - cc_output_fname = "cc_output.stats" - if len(open(cc_output_fname).read().strip()) == 0: - raise Exception('The main output file is empty, there may be an error with your input file or settings.') - except Exception as e: - stop_err('Error running cuffcompare: %s' % e) - - -if __name__ == "__main__": - __main__()
--- a/cuffcompare_wrapper.xml Sun Feb 19 12:11:14 2017 -0500 +++ b/cuffcompare_wrapper.xml Tue Jun 16 16:56:55 2020 +0000 @@ -1,31 +1,34 @@ -<tool id="cuffcompare" name="Cuffcompare" version="@VERSION@.1"> +<tool id="cuffcompare" name="Cuffcompare" version="@VERSION@.2"> <description>compare assembled transcripts to a reference annotation and track Cufflinks transcripts across multiple experiments</description> <macros> <import>cuff_macros.xml</import> </macros> <expand macro="requirements" /> - <version_command>cuffcompare 2>&1 | head -n 1</version_command> + <version_command><![CDATA[ +cuffcompare 2>&1 | head -n 1 + ]]></version_command> <command detect_errors="aggressive"><![CDATA[ - python '$__tool_directory__/cuffcompare_wrapper.py' + @CUFFLINKS_LINK_GTF_INPUTS@ + cuffcompare + -o output ## Use annotation reference? #if $annotation.use_ref_annotation == "Yes": - -r '$annotation.reference_annotation' + -r '$annotation.reference_annotation' #if $annotation.ignore_nonoverlapping_reference: -R #end if #if $annotation.ignore_nonoverlapping_transfrags: -Q #end if - #end if ## Use sequence data? #if $seq_data.use_seq_data == "Yes": -s #if $seq_data.seq_source.index_source == "history": - --ref_file '$seq_data.seq_source.ref_file' + '$seq_data.seq_source.ref_file' #else: - --index '${seq_data.seq_source.index.fields.path}' + '${seq_data.seq_source.index.fields.path}' #end if #end if @@ -34,15 +37,10 @@ -e $max_dist_exon -d $max_dist_group - #if $discard_intron_redundant_transfrags: - -F - #end if - - ## Outputs. - --combined-transcripts '${transcripts_combined}' + $discard_intron_redundant_transfrags @CUFFLINKS_GTF_INPUTS@ - ]]></command> +]]></command> <inputs> <expand macro="cufflinks_gtf_inputs" /> <conditional name="annotation"> @@ -86,47 +84,48 @@ </when> </conditional> <param type="select" name="discard_single_exon" label="discard (ignore) single-exon transcripts"> - <option value="" selected="True">No</option> - <option value="-M">Discard single-exon transfrags and reference transcripts</option> - <option value="-N">Discard single-exon reference transcripts</option> + <option value="" selected="True">No</option> + <option value="-M">Discard single-exon transfrags and reference transcripts</option> + <option value="-N">Discard single-exon reference transcripts</option> </param> <param type="integer" name="max_dist_exon" value="100" label="Max. Distance for assessing exon accuracy" help="max. distance (range) allowed from free ends of terminal exons of reference transcripts when assessing exon accuracy. Default: 100" /> <param type="integer" name="max_dist_group" value="100" label="Max.Distance for transcript grouping" help="max. distance (range) for grouping transcript start sites. Default: 100" /> <param type="boolean" name="discard_intron_redundant_transfrags" label="discard intron-redundant transfrags sharing 5'" + value="false" truevalue="-F" falsevalue="" help="Discard intron-redundant transfrags if they share the 5' end (if they differ only at the 3' end)" /> </inputs> <outputs> <data format="txt" name="transcripts_accuracy" label="${tool.name} on ${on_string}: transcript accuracy" - from_work_dir="cc_output.stats" /> + from_work_dir="output.stats" /> <data format="tabular" name="input1_tmap" label="${tool.name} on ${on_string}: data ${inputs[0].hid} tmap file" - from_work_dir="cc_output.input1.tmap" /> + from_work_dir="output.input_0.tmap" /> <data format="tabular" name="input1_refmap" label="${tool.name} on ${on_string}: data ${inputs[0].hid} refmap file" - from_work_dir="cc_output.input1.refmap"> + from_work_dir="output.input_0.refmap"> <filter>annotation['use_ref_annotation'] == 'Yes'</filter> </data> - <data format="tabular" name="input2_tmap" label="${tool.name} on ${on_string}: data ${inputs[1].hid} tmap file" from_work_dir="cc_output.input2.tmap"> + <data format="tabular" name="input2_tmap" label="${tool.name} on ${on_string}: data ${inputs[1].hid} tmap file" from_work_dir="output.input_1.tmap"> <filter>@HAS_MULTIPLE_INPUTS@</filter> </data> <data format="tabular" name="input2_refmap" label="${tool.name} on ${on_string}: data ${inputs[1].hid} refmap file" - from_work_dir="cc_output.input2.refmap"> + from_work_dir="output.input_1.refmap"> <filter>annotation['use_ref_annotation'] == 'Yes' and @HAS_MULTIPLE_INPUTS@</filter> </data> - <data format="tabular" name="transcripts_tracking" label="${tool.name} on ${on_string}: transcript tracking" from_work_dir="cc_output.tracking"> + <data format="tabular" name="transcripts_tracking" label="${tool.name} on ${on_string}: transcript tracking" from_work_dir="output.tracking"> <filter>@HAS_MULTIPLE_INPUTS@</filter> </data> - <data format="gtf" name="transcripts_combined" label="${tool.name} on ${on_string}: combined transcripts"/> + <data format="gtf" name="transcripts_combined" label="${tool.name} on ${on_string}: combined transcripts" from_work_dir="output.combined.gtf"/> </outputs> <tests> <!-- cuffcompare -r cuffcompare_in3.gtf -R cuffcompare_in1.gtf cuffcompare_in2.gtf --> - <test> + <test expect_num_outputs="7"> <param name="inputs" value="cuffcompare_in1.gtf,cuffcompare_in2.gtf" ftype="gtf"/> <param name="use_ref_annotation" value="Yes"/> <param name="reference_annotation" value="cuffcompare_in3.gtf" ftype="gtf"/> @@ -136,10 +135,10 @@ <param name="discard_single_exon" value="" /> <param name="max_dist_exon" value="100" /> <param name="max_dist_group" value="100" /> - <param name="discard_intron_redundant_transfrags" value="No" /> + <param name="discard_intron_redundant_transfrags" value="false" /> <!-- Line diffs are the result of different locations for input files; this cannot be fixed as cuffcompare outputs full input path for each input. --> - <output name="transcripts_accuracy" file="cuffcompare_out7.txt" lines_diff="2"/> + <output name="transcripts_accuracy" file="cuffcompare_out7.txt" lines_diff="6"/> <output name="input1_tmap" file="cuffcompare_out1.tmap"/> <output name="input1_refmap" file="cuffcompare_out2.refmap"/> <output name="input2_tmap" file="cuffcompare_out3.tmap"/> @@ -248,7 +247,5 @@ -r An optional "reference" annotation GTF. Each sample is matched against this file, and sample isoforms are tagged as overlapping, matching, or novel where appropriate. See the refmap and tmap output file descriptions below. -R If -r was specified, this option causes cuffcompare to ignore reference transcripts that are not overlapped by any transcript in one of cuff1.gtf,...,cuffN.gtf. Useful for ignoring annotated transcripts that are not present in your RNA-Seq samples and thus adjusting the "sensitivity" calculation in the accuracy report written in the transcripts_accuracy file </help> - <citations> - <citation type="doi">10.1038/nbt.1621</citation> - </citations> + <expand macro="citations"/> </tool>
--- a/test-data/cuffcompare_out7.txt Sun Feb 19 12:11:14 2017 -0500 +++ b/test-data/cuffcompare_out7.txt Tue Jun 16 16:56:55 2020 +0000 @@ -1,8 +1,8 @@ # Cuffcompare v2.2.1 | Command line was: -#cuffcompare -o cc_output -r /tmp/tmpZ9KXPVfiles/000/dataset_3.dat -R -e 100 -d 100 ./input1 ./input2 +#cuffcompare -o output -r /tmp/tmpt36dzh0b/files/7/5/1/dataset_75119a58-1488-40cb-bd1d-4e92052062f5.dat -R -e 100 -d 100 /tmp/tmpt36dzh0b/files/1/c/9/dataset_1c9fee4e-77d4-48eb-95ec-43841c5b63b1.dat /tmp/tmpt36dzh0b/files/0/d/d/dataset_0ddbc5aa-4094-447c-ae9e-f30a1f7af6a9.dat # -#= Summary for dataset: ./input1 : +#= Summary for dataset: /tmp/tmpt36dzh0b/files/1/c/9/dataset_1c9fee4e-77d4-48eb-95ec-43841c5b63b1.dat : # Query mRNAs : 50 in 50 loci (0 multi-exon transcripts) # (0 multi-transcript loci, ~1.0 transcripts per locus) # Reference mRNAs : 1 in 1 loci (1 multi-exon) @@ -24,7 +24,7 @@ Missed loci: 0/1 ( 0.0%) Novel loci: 49/50 ( 98.0%) -#= Summary for dataset: ./input2 : +#= Summary for dataset: /tmp/tmpt36dzh0b/files/0/d/d/dataset_0ddbc5aa-4094-447c-ae9e-f30a1f7af6a9.dat : # Query mRNAs : 50 in 50 loci (0 multi-exon transcripts) # (0 multi-transcript loci, ~1.0 transcripts per locus) # Reference mRNAs : 0 in 0 loci (0 multi-exon)
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.combined.gtf Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,90 @@ +chr1 Cufflinks exon 3204755 3204833 . - . gene_id "XLOC_000001"; transcript_id "TCONS_00000001"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.13.1"; nearest_ref "Xkr4"; class_code "c"; tss_id "TSS1"; +chr1 Cufflinks exon 3111450 3111490 . . . gene_id "XLOC_000002"; transcript_id "TCONS_00000002"; exon_number "1"; oId "CUFF.1.1"; class_code "u"; tss_id "TSS2"; +chr1 Cufflinks exon 3111546 3111576 . . . gene_id "XLOC_000003"; transcript_id "TCONS_00000003"; exon_number "1"; oId "CUFF.3.1"; class_code "u"; tss_id "TSS3"; +chr1 Cufflinks exon 3174766 3174792 . . . gene_id "XLOC_000004"; transcript_id "TCONS_00000051"; exon_number "1"; oId "CUFF.1.1"; class_code "u"; tss_id "TSS4"; +chr1 Cufflinks exon 3187402 3187428 . . . gene_id "XLOC_000005"; transcript_id "TCONS_00000052"; exon_number "1"; oId "CUFF.3.1"; class_code "u"; tss_id "TSS5"; +chr1 Cufflinks exon 3188522 3188548 . . . gene_id "XLOC_000006"; transcript_id "TCONS_00000053"; exon_number "1"; oId "CUFF.5.1"; class_code "u"; tss_id "TSS6"; +chr1 Cufflinks exon 3189811 3190789 . . . gene_id "XLOC_000007"; transcript_id "TCONS_00000054"; exon_number "1"; oId "CUFF.11.1"; class_code "u"; tss_id "TSS7"; +chr1 Cufflinks exon 3190859 3191434 . . . gene_id "XLOC_000008"; transcript_id "TCONS_00000055"; exon_number "1"; oId "CUFF.7.1"; class_code "u"; tss_id "TSS10"; +chr1 Cufflinks exon 3191513 3192077 . . . gene_id "XLOC_000009"; transcript_id "TCONS_00000056"; exon_number "1"; oId "CUFF.9.1"; class_code "u"; tss_id "TSS11"; +chr1 Cufflinks exon 3192251 3192336 . . . gene_id "XLOC_000010"; transcript_id "TCONS_00000057"; exon_number "1"; oId "CUFF.13.1"; class_code "u"; tss_id "TSS13"; +chr1 Cufflinks exon 3192442 3192494 . . . gene_id "XLOC_000011"; transcript_id "TCONS_00000009"; exon_number "1"; oId "CUFF.39.1"; class_code "u"; tss_id "TSS14"; +chr1 Cufflinks exon 3192551 3192629 . . . gene_id "XLOC_000012"; transcript_id "TCONS_00000010"; exon_number "1"; oId "CUFF.41.1"; class_code "u"; tss_id "TSS15"; +chr1 Cufflinks exon 3192650 3192676 . . . gene_id "XLOC_000013"; transcript_id "TCONS_00000058"; exon_number "1"; oId "CUFF.15.1"; class_code "u"; tss_id "TSS16"; +chr1 Cufflinks exon 3192732 3192811 . . . gene_id "XLOC_000014"; transcript_id "TCONS_00000011"; exon_number "1"; oId "CUFF.43.1"; class_code "u"; tss_id "TSS17"; +chr1 Cufflinks exon 3192941 3193042 . . . gene_id "XLOC_000015"; transcript_id "TCONS_00000012"; exon_number "1"; oId "CUFF.45.1"; class_code "u"; tss_id "TSS18"; +chr1 Cufflinks exon 3194186 3194226 . . . gene_id "XLOC_000016"; transcript_id "TCONS_00000013"; exon_number "1"; oId "CUFF.47.1"; class_code "u"; tss_id "TSS19"; +chr1 Cufflinks exon 3194303 3194329 . . . gene_id "XLOC_000017"; transcript_id "TCONS_00000014"; exon_number "1"; oId "CUFF.49.1"; class_code "u"; tss_id "TSS20"; +chr1 Cufflinks exon 3194707 3194733 . . . gene_id "XLOC_000018"; transcript_id "TCONS_00000059"; exon_number "1"; oId "CUFF.17.1"; class_code "u"; tss_id "TSS21"; +chr1 Cufflinks exon 3195084 3195110 . . . gene_id "XLOC_000019"; transcript_id "TCONS_00000015"; exon_number "1"; oId "CUFF.51.1"; class_code "u"; tss_id "TSS22"; +chr1 Cufflinks exon 3195451 3195477 . . . gene_id "XLOC_000020"; transcript_id "TCONS_00000016"; exon_number "1"; oId "CUFF.53.1"; class_code "u"; tss_id "TSS23"; +chr1 Cufflinks exon 3197090 3197116 . . . gene_id "XLOC_000021"; transcript_id "TCONS_00000017"; exon_number "1"; oId "CUFF.55.1"; class_code "u"; tss_id "TSS24"; +chr1 Cufflinks exon 3197247 3197273 . . . gene_id "XLOC_000022"; transcript_id "TCONS_00000018"; exon_number "1"; oId "CUFF.57.1"; class_code "u"; tss_id "TSS25"; +chr1 Cufflinks exon 3197347 3197373 . . . gene_id "XLOC_000023"; transcript_id "TCONS_00000019"; exon_number "1"; oId "CUFF.59.1"; class_code "u"; tss_id "TSS26"; +chr1 Cufflinks exon 3197426 3197452 . . . gene_id "XLOC_000024"; transcript_id "TCONS_00000060"; exon_number "1"; oId "CUFF.19.1"; class_code "u"; tss_id "TSS27"; +chr1 Cufflinks exon 3200023 3200191 . . . gene_id "XLOC_000025"; transcript_id "TCONS_00000020"; exon_number "1"; oId "CUFF.7.1"; class_code "."; tss_id "TSS28"; +chr1 Cufflinks exon 3200326 3200352 . . . gene_id "XLOC_000026"; transcript_id "TCONS_00000021"; exon_number "1"; oId "CUFF.5.1"; class_code "u"; tss_id "TSS29"; +chr1 Cufflinks exon 3200431 3200457 . . . gene_id "XLOC_000027"; transcript_id "TCONS_00000061"; exon_number "1"; oId "CUFF.21.1"; class_code "u"; tss_id "TSS30"; +chr1 Cufflinks exon 3201008 3201039 . . . gene_id "XLOC_000028"; transcript_id "TCONS_00000062"; exon_number "1"; oId "CUFF.26.1"; class_code "u"; tss_id "TSS31"; +chr1 Cufflinks exon 3201078 3201481 . . . gene_id "XLOC_000029"; transcript_id "TCONS_00000022"; exon_number "1"; oId "CUFF.9.1"; class_code "."; tss_id "TSS32"; +chr1 Cufflinks exon 3201597 3201666 . . . gene_id "XLOC_000030"; transcript_id "TCONS_00000063"; exon_number "1"; oId "CUFF.29.1"; class_code "u"; tss_id "TSS33"; +chr1 Cufflinks exon 3201673 3201699 . . . gene_id "XLOC_000031"; transcript_id "TCONS_00000023"; exon_number "1"; oId "CUFF.11.1"; class_code "u"; tss_id "TSS34"; +chr1 Cufflinks exon 3201726 3201809 . . . gene_id "XLOC_000032"; transcript_id "TCONS_00000064"; exon_number "1"; oId "CUFF.31.1"; class_code "u"; tss_id "TSS35"; +chr1 Cufflinks exon 3211522 3211561 . . . gene_id "XLOC_000033"; transcript_id "TCONS_00000065"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.33.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS36"; +chr1 Cufflinks exon 3212214 3212292 . . . gene_id "XLOC_000034"; transcript_id "TCONS_00000024"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.15.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS37"; +chr1 Cufflinks exon 3212368 3212439 . . . gene_id "XLOC_000035"; transcript_id "TCONS_00000025"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.19.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS38"; +chr1 Cufflinks exon 3212718 3212801 . . . gene_id "XLOC_000036"; transcript_id "TCONS_00000066"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.35.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS39"; +chr1 Cufflinks exon 3213096 3213192 . . . gene_id "XLOC_000037"; transcript_id "TCONS_00000026"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.17.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS40"; +chr1 Cufflinks exon 3213119 3213242 . . . gene_id "XLOC_000037"; transcript_id "TCONS_00000067"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.37.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS40"; +chr1 Cufflinks exon 3240607 3240633 . . . gene_id "XLOC_000038"; transcript_id "TCONS_00000068"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.39.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS41"; +chr1 Cufflinks exon 3242480 3242512 . . . gene_id "XLOC_000039"; transcript_id "TCONS_00000069"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.41.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS42"; +chr1 Cufflinks exon 3242634 3242923 . . . gene_id "XLOC_000040"; transcript_id "TCONS_00000027"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.25.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS43"; +chr1 Cufflinks exon 3242925 3243005 . . . gene_id "XLOC_000041"; transcript_id "TCONS_00000070"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.43.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS44"; +chr1 Cufflinks exon 3243019 3243079 . . . gene_id "XLOC_000042"; transcript_id "TCONS_00000028"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.21.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS45"; +chr1 Cufflinks exon 3243109 3243154 . . . gene_id "XLOC_000043"; transcript_id "TCONS_00000071"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.45.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS46"; +chr1 Cufflinks exon 3243348 3243401 . . . gene_id "XLOC_000044"; transcript_id "TCONS_00000029"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.23.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS47"; +chr1 Cufflinks exon 3254080 3254106 . . . gene_id "XLOC_000045"; transcript_id "TCONS_00000072"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.47.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS48"; +chr1 Cufflinks exon 3256975 3257011 . . . gene_id "XLOC_000046"; transcript_id "TCONS_00000030"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.27.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS49"; +chr1 Cufflinks exon 3277156 3277182 . . . gene_id "XLOC_000047"; transcript_id "TCONS_00000073"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.49.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS50"; +chr1 Cufflinks exon 3277191 3277218 . . . gene_id "XLOC_000048"; transcript_id "TCONS_00000031"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.61.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS51"; +chr1 Cufflinks exon 3277914 3278390 . . . gene_id "XLOC_000049"; transcript_id "TCONS_00000074"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.51.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS52"; +chr1 Cufflinks exon 3280118 3280144 . . . gene_id "XLOC_000050"; transcript_id "TCONS_00000075"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.53.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS54"; +chr1 Cufflinks exon 3280499 3280525 . . . gene_id "XLOC_000051"; transcript_id "TCONS_00000076"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.55.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS55"; +chr1 Cufflinks exon 3280687 3280741 . . . gene_id "XLOC_000052"; transcript_id "TCONS_00000033"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.65.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS56"; +chr1 Cufflinks exon 3282505 3282531 . . . gene_id "XLOC_000053"; transcript_id "TCONS_00000077"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.57.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS57"; +chr1 Cufflinks exon 3282651 3282677 . . . gene_id "XLOC_000054"; transcript_id "TCONS_00000078"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.59.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS58"; +chr1 Cufflinks exon 3282761 3282832 . . . gene_id "XLOC_000055"; transcript_id "TCONS_00000079"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.61.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS59"; +chr1 Cufflinks exon 3284967 3284993 . . . gene_id "XLOC_000056"; transcript_id "TCONS_00000080"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.63.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS60"; +chr1 Cufflinks exon 3290489 3290553 . . . gene_id "XLOC_000057"; transcript_id "TCONS_00000034"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.67.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS61"; +chr1 Cufflinks exon 3290799 3290859 . . . gene_id "XLOC_000058"; transcript_id "TCONS_00000081"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.65.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS62"; +chr1 Cufflinks exon 3290920 3291273 . . . gene_id "XLOC_000059"; transcript_id "TCONS_00000082"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.69.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS63"; +chr1 Cufflinks exon 3299444 3299640 . . . gene_id "XLOC_000060"; transcript_id "TCONS_00000083"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.67.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS65"; +chr1 Cufflinks exon 3299610 3299664 . . . gene_id "XLOC_000060"; transcript_id "TCONS_00000037"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.73.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS66"; +chr1 Cufflinks exon 3299692 3299733 . . . gene_id "XLOC_000061"; transcript_id "TCONS_00000084"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.71.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS67"; +chr1 Cufflinks exon 3300052 3300078 . . . gene_id "XLOC_000062"; transcript_id "TCONS_00000038"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.75.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS68"; +chr1 Cufflinks exon 3307749 3307775 . . . gene_id "XLOC_000063"; transcript_id "TCONS_00000085"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.73.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS69"; +chr1 Cufflinks exon 3318621 3318647 . . . gene_id "XLOC_000064"; transcript_id "TCONS_00000086"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.75.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS70"; +chr1 Cufflinks exon 3319000 3319051 . . . gene_id "XLOC_000065"; transcript_id "TCONS_00000039"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.77.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS71"; +chr1 Cufflinks exon 3330528 3330554 . . . gene_id "XLOC_000066"; transcript_id "TCONS_00000087"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.77.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS72"; +chr1 Cufflinks exon 3351241 3351311 . . . gene_id "XLOC_000067"; transcript_id "TCONS_00000088"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.79.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS73"; +chr1 Cufflinks exon 3355888 3355914 . . . gene_id "XLOC_000068"; transcript_id "TCONS_00000040"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.79.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS74"; +chr1 Cufflinks exon 3355908 3356119 . . . gene_id "XLOC_000068"; transcript_id "TCONS_00000089"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.81.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS74"; +chr1 Cufflinks exon 3356181 3356225 . . . gene_id "XLOC_000069"; transcript_id "TCONS_00000090"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.83.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS75"; +chr1 Cufflinks exon 3363077 3363176 . . . gene_id "XLOC_000070"; transcript_id "TCONS_00000091"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.85.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS76"; +chr1 Cufflinks exon 3363215 3363278 . . . gene_id "XLOC_000071"; transcript_id "TCONS_00000041"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.81.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS77"; +chr1 Cufflinks exon 3363388 3363446 . . . gene_id "XLOC_000072"; transcript_id "TCONS_00000092"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.87.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS78"; +chr1 Cufflinks exon 3363754 3363849 . . . gene_id "XLOC_000073"; transcript_id "TCONS_00000042"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.83.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS79"; +chr1 Cufflinks exon 3364872 3364919 . . . gene_id "XLOC_000074"; transcript_id "TCONS_00000093"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.89.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS80"; +chr1 Cufflinks exon 3367136 3367162 . . . gene_id "XLOC_000075"; transcript_id "TCONS_00000043"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.85.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS81"; +chr1 Cufflinks exon 3367211 3367237 . . . gene_id "XLOC_000076"; transcript_id "TCONS_00000094"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.91.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS82"; +chr1 Cufflinks exon 3367334 3367382 . . . gene_id "XLOC_000077"; transcript_id "TCONS_00000044"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.87.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS83"; +chr1 Cufflinks exon 3369581 3369607 . . . gene_id "XLOC_000078"; transcript_id "TCONS_00000095"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.93.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS84"; +chr1 Cufflinks exon 3375002 3375028 . . . gene_id "XLOC_000079"; transcript_id "TCONS_00000096"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.95.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS85"; +chr1 Cufflinks exon 3377212 3377262 . . . gene_id "XLOC_000080"; transcript_id "TCONS_00000045"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.89.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS86"; +chr1 Cufflinks exon 3379889 3379915 . . . gene_id "XLOC_000081"; transcript_id "TCONS_00000097"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.97.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS87"; +chr1 Cufflinks exon 3386740 3386836 . . . gene_id "XLOC_000082"; transcript_id "TCONS_00000098"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.99.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS88"; +chr1 Cufflinks exon 3391326 3391352 . . . gene_id "XLOC_000083"; transcript_id "TCONS_00000046"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.91.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS89"; +chr1 Cufflinks exon 3435842 3435880 . . . gene_id "XLOC_000084"; transcript_id "TCONS_00000047"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.93.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS90"; +chr1 Cufflinks exon 3447762 3447788 . . . gene_id "XLOC_000085"; transcript_id "TCONS_00000048"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.95.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS91"; +chr1 Cufflinks exon 3450907 3450965 . . . gene_id "XLOC_000086"; transcript_id "TCONS_00000049"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.97.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS92"; +chr1 Cufflinks exon 3451052 3451109 . . . gene_id "XLOC_000087"; transcript_id "TCONS_00000050"; exon_number "1"; gene_name "Xkr4"; oId "CUFF.99.1"; nearest_ref "Xkr4"; class_code "i"; tss_id "TSS93";
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.cuffcompare_in1.gtf.refmap Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,2 @@ +ref_gene_id ref_id class_code cuff_id_list +Xkr4 Xkr4 c CUFF.13|CUFF.13.1
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.cuffcompare_in1.gtf.tmap Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,51 @@ +ref_gene_id ref_id class_code cuff_gene_id cuff_id FMI FPKM FPKM_conf_lo FPKM_conf_hi cov len major_iso_id ref_match_len +- - u CUFF.1 CUFF.1.1 100 20.607936 0.000000 49.751960 1.317073 41 CUFF.1.1 - +- - u CUFF.3 CUFF.3.1 100 27.255658 0.000000 65.800979 1.741935 31 CUFF.3.1 - +- - u CUFF.5 CUFF.5.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.5.1 - +- - u CUFF.7 CUFF.7.1 100 9.999117 0.000000 19.998234 0.639053 169 CUFF.7.1 - +- - u CUFF.9 CUFF.9.1 100 17.776896 9.153835 26.399957 1.136139 404 CUFF.9.1 - +- - u CUFF.11 CUFF.11.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.11.1 - +Xkr4 Xkr4 c CUFF.13 CUFF.13.1 100 10.695258 0.000000 25.820637 0.683544 79 CUFF.13.1 3634 +Xkr4 Xkr4 i CUFF.15 CUFF.15.1 100 10.695258 0.000000 25.820637 0.683544 79 CUFF.15.1 3634 +Xkr4 Xkr4 i CUFF.17 CUFF.17.1 100 8.710571 0.000000 21.029179 0.556701 97 CUFF.17.1 3634 +Xkr4 Xkr4 i CUFF.19 CUFF.19.1 100 29.337687 3.097262 55.578113 1.875000 72 CUFF.19.1 3634 +Xkr4 Xkr4 i CUFF.21 CUFF.21.1 100 13.851236 0.000000 33.439842 0.885246 61 CUFF.21.1 3634 +Xkr4 Xkr4 i CUFF.23 CUFF.23.1 100 23.470150 0.000000 50.571145 1.500000 54 CUFF.23.1 3634 +Xkr4 Xkr4 i CUFF.25 CUFF.25.1 100 14.567679 5.354270 23.781089 0.931034 290 CUFF.25.1 3634 +Xkr4 Xkr4 i CUFF.27 CUFF.27.1 100 34.253732 0.000000 73.806535 2.189189 37 CUFF.27.1 3634 +- - u CUFF.29 CUFF.29.1 100 107.103219 71.402146 142.804292 6.845070 142 CUFF.29.1 - +- - u CUFF.31 CUFF.31.1 100 122.650461 40.883487 204.417435 7.838710 31 CUFF.31.1 - +- - u CUFF.33 CUFF.33.1 100 109.527366 26.732460 192.322273 7.000000 27 CUFF.33.1 - +- - u CUFF.35 CUFF.35.1 100 96.747183 61.420107 132.074259 6.183206 131 CUFF.35.1 - +- - u CUFF.37 CUFF.37.1 100 104.085013 53.596365 154.573660 6.652174 69 CUFF.37.1 - +- - u CUFF.39 CUFF.39.1 100 23.912983 0.000000 51.525317 1.528302 53 CUFF.39.1 - +- - u CUFF.41 CUFF.41.1 100 10.695258 0.000000 25.820637 0.683544 79 CUFF.41.1 - +- - u CUFF.43 CUFF.43.1 100 10.561567 0.000000 25.497879 0.675000 80 CUFF.43.1 - +- - u CUFF.45 CUFF.45.1 100 20.708956 2.186303 39.231609 1.323529 102 CUFF.45.1 - +- - u CUFF.47 CUFF.47.1 100 20.607936 0.000000 49.751960 1.317073 41 CUFF.47.1 - +- - u CUFF.49 CUFF.49.1 100 15.646767 0.000000 46.940300 1.000000 27 CUFF.49.1 - +- - u CUFF.51 CUFF.51.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.51.1 - +- - u CUFF.53 CUFF.53.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.53.1 - +- - u CUFF.55 CUFF.55.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.55.1 - +- - u CUFF.57 CUFF.57.1 100 15.646767 0.000000 46.940300 1.000000 27 CUFF.57.1 - +- - u CUFF.59 CUFF.59.1 100 15.646767 0.000000 46.940300 1.000000 27 CUFF.59.1 - +Xkr4 Xkr4 i CUFF.61 CUFF.61.1 100 45.263860 0.000000 97.530065 2.892857 28 CUFF.61.1 3634 +Xkr4 Xkr4 i CUFF.63 CUFF.63.1 100 15.646767 0.000000 46.940300 1.000000 27 CUFF.63.1 3634 +Xkr4 Xkr4 i CUFF.65 CUFF.65.1 100 15.362280 0.000000 37.087825 0.981818 55 CUFF.65.1 3634 +Xkr4 Xkr4 i CUFF.67 CUFF.67.1 100 12.998852 0.000000 31.382005 0.830769 65 CUFF.67.1 3634 +Xkr4 Xkr4 i CUFF.69 CUFF.69.1 100 10.058636 0.000000 24.283695 0.642857 84 CUFF.69.1 3634 +Xkr4 Xkr4 i CUFF.71 CUFF.71.1 100 8.621688 0.000000 20.814595 0.551020 98 CUFF.71.1 3634 +Xkr4 Xkr4 i CUFF.73 CUFF.73.1 100 15.362280 0.000000 37.087825 0.981818 55 CUFF.73.1 3634 +Xkr4 Xkr4 i CUFF.75 CUFF.75.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.75.1 3634 +Xkr4 Xkr4 i CUFF.77 CUFF.77.1 100 16.248565 0.000000 39.227507 1.038462 52 CUFF.77.1 3634 +Xkr4 Xkr4 i CUFF.79 CUFF.79.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.79.1 3634 +Xkr4 Xkr4 i CUFF.81 CUFF.81.1 100 13.201959 0.000000 31.872349 0.843750 64 CUFF.81.1 3634 +Xkr4 Xkr4 i CUFF.83 CUFF.83.1 100 13.201959 0.000000 28.446269 0.843750 96 CUFF.83.1 3634 +Xkr4 Xkr4 i CUFF.85 CUFF.85.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.85.1 3634 +Xkr4 Xkr4 i CUFF.87 CUFF.87.1 100 17.243375 0.000000 41.629191 1.102041 49 CUFF.87.1 3634 +Xkr4 Xkr4 i CUFF.89 CUFF.89.1 100 16.567165 0.000000 39.996674 1.058824 51 CUFF.89.1 3634 +Xkr4 Xkr4 i CUFF.91 CUFF.91.1 100 31.293533 0.000000 75.549272 2.000000 27 CUFF.91.1 3634 +Xkr4 Xkr4 i CUFF.93 CUFF.93.1 100 21.664754 0.000000 52.303342 1.384615 39 CUFF.93.1 3634 +Xkr4 Xkr4 i CUFF.95 CUFF.95.1 100 46.940300 0.000000 101.142289 3.000000 27 CUFF.95.1 3634 +Xkr4 Xkr4 i CUFF.97 CUFF.97.1 100 21.481154 0.000000 46.285454 1.372881 59 CUFF.97.1 3634 +Xkr4 Xkr4 i CUFF.99 CUFF.99.1 100 14.567679 0.000000 35.169489 0.931034 58 CUFF.99.1 3634
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.cuffcompare_in2.gtf.refmap Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,1 @@ +ref_gene_id ref_id class_code cuff_id_list
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.cuffcompare_in2.gtf.tmap Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,51 @@ +ref_gene_id ref_id class_code cuff_gene_id cuff_id FMI FPKM FPKM_conf_lo FPKM_conf_hi cov len major_iso_id ref_match_len +- - u CUFF.1 CUFF.1.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.1.1 - +- - u CUFF.3 CUFF.3.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.3.1 - +- - u CUFF.5 CUFF.5.1 100 21.226627 0.000000 59.889707 1.205672 27 CUFF.5.1 - +- - u CUFF.7 CUFF.7.1 100 29.709524 19.806349 39.612698 1.687500 576 CUFF.7.1 - +- - u CUFF.9 CUFF.9.1 100 34.072933 23.364686 44.781179 1.935341 565 CUFF.9.1 - +- - u CUFF.11 CUFF.11.1 100 32.531777 24.582998 40.480555 1.847804 979 CUFF.11.1 - +- - u CUFF.13 CUFF.13.1 100 16.582060 0.000000 35.729373 0.941860 86 CUFF.13.1 - +- - u CUFF.15 CUFF.15.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.15.1 - +- - u CUFF.17 CUFF.17.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.17.1 - +- - u CUFF.19 CUFF.19.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.19.1 - +- - u CUFF.21 CUFF.21.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.21.1 - +- - u CUFF.23 CUFF.23.1 100 16.205195 0.000000 34.917342 0.920455 88 CUFF.23.1 - +- - u CUFF.25 CUFF.25.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.25.1 - +- - u CUFF.26 CUFF.26.1 100 29.709524 0.000000 71.725135 1.687500 32 CUFF.26.1 - +- - u CUFF.29 CUFF.29.1 100 13.581496 0.000000 32.788633 0.771429 70 CUFF.29.1 - +- - u CUFF.31 CUFF.31.1 100 22.635827 0.000000 45.271655 1.285714 84 CUFF.31.1 - +Xkr4 Xkr4 i CUFF.33 CUFF.33.1 100 23.767619 0.000000 57.380108 1.350000 40 CUFF.33.1 3634 +Xkr4 Xkr4 i CUFF.35 CUFF.35.1 100 11.317914 0.000000 27.323861 0.642857 84 CUFF.35.1 3634 +Xkr4 Xkr4 i CUFF.37 CUFF.37.1 100 11.500461 0.000000 24.780049 0.653226 124 CUFF.37.1 3634 +Xkr4 Xkr4 i CUFF.39 CUFF.39.1 100 52.816931 0.000000 113.804669 3.000000 27 CUFF.39.1 3634 +Xkr4 Xkr4 i CUFF.41 CUFF.41.1 100 43.213852 0.000000 93.112911 2.454545 33 CUFF.41.1 3634 +Xkr4 Xkr4 i CUFF.43 CUFF.43.1 100 23.474191 0.000000 46.948383 1.333333 81 CUFF.43.1 3634 +Xkr4 Xkr4 i CUFF.45 CUFF.45.1 100 20.667495 0.000000 49.895746 1.173913 46 CUFF.45.1 3634 +Xkr4 Xkr4 i CUFF.47 CUFF.47.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.47.1 3634 +Xkr4 Xkr4 i CUFF.49 CUFF.49.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.49.1 3634 +Xkr4 Xkr4 i CUFF.51 CUFF.51.1 100 14.948188 7.228977 22.667399 0.849057 477 CUFF.51.1 3634 +Xkr4 Xkr4 i CUFF.53 CUFF.53.1 100 52.816931 0.000000 113.804669 3.000000 27 CUFF.53.1 3634 +Xkr4 Xkr4 i CUFF.55 CUFF.55.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.55.1 3634 +Xkr4 Xkr4 i CUFF.57 CUFF.57.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.57.1 3634 +Xkr4 Xkr4 i CUFF.59 CUFF.59.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.59.1 3634 +Xkr4 Xkr4 i CUFF.61 CUFF.61.1 100 13.204233 0.000000 31.877838 0.750000 72 CUFF.61.1 3634 +Xkr4 Xkr4 i CUFF.63 CUFF.63.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.63.1 3634 +Xkr4 Xkr4 i CUFF.65 CUFF.65.1 100 31.170648 0.000000 62.341295 1.770492 61 CUFF.65.1 3634 +Xkr4 Xkr4 i CUFF.67 CUFF.67.1 100 15.681351 3.378764 27.983938 0.890700 197 CUFF.67.1 3634 +Xkr4 Xkr4 i CUFF.69 CUFF.69.1 100 18.799247 8.750627 28.847866 1.067797 354 CUFF.69.1 3634 +Xkr4 Xkr4 i CUFF.71 CUFF.71.1 100 22.635827 0.000000 54.647722 1.285714 42 CUFF.71.1 3634 +Xkr4 Xkr4 i CUFF.73 CUFF.73.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.73.1 3634 +Xkr4 Xkr4 i CUFF.75 CUFF.75.1 100 52.816931 0.000000 113.804669 3.000000 27 CUFF.75.1 3634 +Xkr4 Xkr4 i CUFF.77 CUFF.77.1 100 17.605644 0.000000 52.816931 1.000000 27 CUFF.77.1 3634 +Xkr4 Xkr4 i CUFF.79 CUFF.79.1 100 13.390208 0.000000 32.326821 0.760563 71 CUFF.79.1 3634 +Xkr4 Xkr4 i CUFF.81 CUFF.81.1 100 11.211141 1.183592 21.238690 0.636792 212 CUFF.81.1 3634 +Xkr4 Xkr4 i CUFF.83 CUFF.83.1 100 21.126772 0.000000 51.004540 1.200000 45 CUFF.83.1 3634 +Xkr4 Xkr4 i CUFF.85 CUFF.85.1 100 19.014095 0.000000 38.028190 1.080000 100 CUFF.85.1 3634 +Xkr4 Xkr4 i CUFF.87 CUFF.87.1 100 24.170460 0.000000 52.080103 1.372881 59 CUFF.87.1 3634 +Xkr4 Xkr4 i CUFF.89 CUFF.89.1 100 29.709524 0.000000 64.015126 1.687500 48 CUFF.89.1 3634 +Xkr4 Xkr4 i CUFF.91 CUFF.91.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.91.1 3634 +Xkr4 Xkr4 i CUFF.93 CUFF.93.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.93.1 3634 +Xkr4 Xkr4 i CUFF.95 CUFF.95.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.95.1 3634 +Xkr4 Xkr4 i CUFF.97 CUFF.97.1 100 35.211287 0.000000 85.007567 2.000000 27 CUFF.97.1 3634 +Xkr4 Xkr4 i CUFF.99 CUFF.99.1 100 19.602160 0.000000 39.204320 1.113402 97 CUFF.99.1 3634
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.loci Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,87 @@ +XLOC_000001 chr1[-]3204563-3661579 Xkr4|Xkr4 CUFF.13.1 - +XLOC_000002 chr1[.]3111450-3111490 - CUFF.1.1 - +XLOC_000003 chr1[.]3111546-3111576 - CUFF.3.1 - +XLOC_000004 chr1[.]3174766-3174792 - - CUFF.1.1 +XLOC_000005 chr1[.]3187402-3187428 - - CUFF.3.1 +XLOC_000006 chr1[.]3188522-3188548 - - CUFF.5.1 +XLOC_000007 chr1[.]3189811-3190789 - CUFF.31.1,CUFF.33.1,CUFF.29.1 CUFF.11.1 +XLOC_000008 chr1[.]3190859-3191434 - - CUFF.7.1 +XLOC_000009 chr1[.]3191513-3192077 - CUFF.35.1,CUFF.37.1 CUFF.9.1 +XLOC_000010 chr1[.]3192251-3192336 - - CUFF.13.1 +XLOC_000011 chr1[.]3192442-3192494 - CUFF.39.1 - +XLOC_000012 chr1[.]3192551-3192629 - CUFF.41.1 - +XLOC_000013 chr1[.]3192650-3192676 - - CUFF.15.1 +XLOC_000014 chr1[.]3192732-3192811 - CUFF.43.1 - +XLOC_000015 chr1[.]3192941-3193042 - CUFF.45.1 - +XLOC_000016 chr1[.]3194186-3194226 - CUFF.47.1 - +XLOC_000017 chr1[.]3194303-3194329 - CUFF.49.1 - +XLOC_000018 chr1[.]3194707-3194733 - - CUFF.17.1 +XLOC_000019 chr1[.]3195084-3195110 - CUFF.51.1 - +XLOC_000020 chr1[.]3195451-3195477 - CUFF.53.1 - +XLOC_000021 chr1[.]3197090-3197116 - CUFF.55.1 - +XLOC_000022 chr1[.]3197247-3197273 - CUFF.57.1 - +XLOC_000023 chr1[.]3197347-3197373 - CUFF.59.1 - +XLOC_000024 chr1[.]3197426-3197452 - - CUFF.19.1 +XLOC_000025 chr1[.]3200023-3200191 - CUFF.7.1 CUFF.23.1 +XLOC_000026 chr1[.]3200326-3200352 - CUFF.5.1 - +XLOC_000027 chr1[.]3200431-3200457 - - CUFF.21.1 +XLOC_000028 chr1[.]3201008-3201039 - - CUFF.26.1 +XLOC_000029 chr1[.]3201078-3201481 - CUFF.9.1 CUFF.25.1 +XLOC_000030 chr1[.]3201597-3201666 - - CUFF.29.1 +XLOC_000031 chr1[.]3201673-3201699 - CUFF.11.1 - +XLOC_000032 chr1[.]3201726-3201809 - - CUFF.31.1 +XLOC_000033 chr1[.]3211522-3211561 - - CUFF.33.1 +XLOC_000034 chr1[.]3212214-3212292 - CUFF.15.1 - +XLOC_000035 chr1[.]3212368-3212439 - CUFF.19.1 - +XLOC_000036 chr1[.]3212718-3212801 - - CUFF.35.1 +XLOC_000037 chr1[.]3213096-3213242 - CUFF.17.1 CUFF.37.1 +XLOC_000038 chr1[.]3240607-3240633 - - CUFF.39.1 +XLOC_000039 chr1[.]3242480-3242512 - - CUFF.41.1 +XLOC_000040 chr1[.]3242634-3242923 - CUFF.25.1 - +XLOC_000041 chr1[.]3242925-3243005 - - CUFF.43.1 +XLOC_000042 chr1[.]3243019-3243079 - CUFF.21.1 - +XLOC_000043 chr1[.]3243109-3243154 - - CUFF.45.1 +XLOC_000044 chr1[.]3243348-3243401 - CUFF.23.1 - +XLOC_000045 chr1[.]3254080-3254106 - - CUFF.47.1 +XLOC_000046 chr1[.]3256975-3257011 - CUFF.27.1 - +XLOC_000047 chr1[.]3277156-3277182 - - CUFF.49.1 +XLOC_000048 chr1[.]3277191-3277218 - CUFF.61.1 - +XLOC_000049 chr1[.]3277914-3278390 - CUFF.63.1 CUFF.51.1 +XLOC_000050 chr1[.]3280118-3280144 - - CUFF.53.1 +XLOC_000051 chr1[.]3280499-3280525 - - CUFF.55.1 +XLOC_000052 chr1[.]3280687-3280741 - CUFF.65.1 - +XLOC_000053 chr1[.]3282505-3282531 - - CUFF.57.1 +XLOC_000054 chr1[.]3282651-3282677 - - CUFF.59.1 +XLOC_000055 chr1[.]3282761-3282832 - - CUFF.61.1 +XLOC_000056 chr1[.]3284967-3284993 - - CUFF.63.1 +XLOC_000057 chr1[.]3290489-3290553 - CUFF.67.1 - +XLOC_000058 chr1[.]3290799-3290859 - - CUFF.65.1 +XLOC_000059 chr1[.]3290920-3291273 - CUFF.69.1,CUFF.71.1 CUFF.69.1 +XLOC_000060 chr1[.]3299444-3299664 - CUFF.73.1 CUFF.67.1 +XLOC_000061 chr1[.]3299692-3299733 - - CUFF.71.1 +XLOC_000062 chr1[.]3300052-3300078 - CUFF.75.1 - +XLOC_000063 chr1[.]3307749-3307775 - - CUFF.73.1 +XLOC_000064 chr1[.]3318621-3318647 - - CUFF.75.1 +XLOC_000065 chr1[.]3319000-3319051 - CUFF.77.1 - +XLOC_000066 chr1[.]3330528-3330554 - - CUFF.77.1 +XLOC_000067 chr1[.]3351241-3351311 - - CUFF.79.1 +XLOC_000068 chr1[.]3355888-3356119 - CUFF.79.1 CUFF.81.1 +XLOC_000069 chr1[.]3356181-3356225 - - CUFF.83.1 +XLOC_000070 chr1[.]3363077-3363176 - - CUFF.85.1 +XLOC_000071 chr1[.]3363215-3363278 - CUFF.81.1 - +XLOC_000072 chr1[.]3363388-3363446 - - CUFF.87.1 +XLOC_000073 chr1[.]3363754-3363849 - CUFF.83.1 - +XLOC_000074 chr1[.]3364872-3364919 - - CUFF.89.1 +XLOC_000075 chr1[.]3367136-3367162 - CUFF.85.1 - +XLOC_000076 chr1[.]3367211-3367237 - - CUFF.91.1 +XLOC_000077 chr1[.]3367334-3367382 - CUFF.87.1 - +XLOC_000078 chr1[.]3369581-3369607 - - CUFF.93.1 +XLOC_000079 chr1[.]3375002-3375028 - - CUFF.95.1 +XLOC_000080 chr1[.]3377212-3377262 - CUFF.89.1 - +XLOC_000081 chr1[.]3379889-3379915 - - CUFF.97.1 +XLOC_000082 chr1[.]3386740-3386836 - - CUFF.99.1 +XLOC_000083 chr1[.]3391326-3391352 - CUFF.91.1 - +XLOC_000084 chr1[.]3435842-3435880 - CUFF.93.1 - +XLOC_000085 chr1[.]3447762-3447788 - CUFF.95.1 - +XLOC_000086 chr1[.]3450907-3450965 - CUFF.97.1 - +XLOC_000087 chr1[.]3451052-3451109 - CUFF.99.1 -
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.stats Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,33 @@ +# Cuffcompare v2.2.1 | Command line was: +#cuffcompare -o output -r cuffcompare_in3.gtf -R -e 100 -d 100 cuffcompare_in1.gtf cuffcompare_in2.gtf +# + +#= Summary for dataset: cuffcompare_in1.gtf : +# Query mRNAs : 50 in 50 loci (0 multi-exon transcripts) +# (0 multi-transcript loci, ~1.0 transcripts per locus) +# Reference mRNAs : 1 in 1 loci (1 multi-exon) +# Super-loci w/ reference transcripts: 1 +#--------------------| Sn | Sp | fSn | fSp + Base level: 2.2 2.3 - - + Exon level: 0.0 0.0 0.0 0.0 + Intron level: 0.0 -nan 0.0 -nan +Intron chain level: 0.0 -nan 0.0 -nan + Transcript level: 0.0 0.0 0.0 0.0 + Locus level: 0.0 0.0 0.0 0.0 + + Matching intron chains: 0 + Matching loci: 0 + + Missed exons: 2/3 ( 66.7%) + Novel exons: 49/50 ( 98.0%) + Missed introns: 2/2 (100.0%) + Missed loci: 0/1 ( 0.0%) + Novel loci: 49/50 ( 98.0%) + +#= Summary for dataset: cuffcompare_in2.gtf : +# Query mRNAs : 50 in 50 loci (0 multi-exon transcripts) +# (0 multi-transcript loci, ~1.0 transcripts per locus) +# Reference mRNAs : 0 in 0 loci (0 multi-exon) + + Total union super-loci across all input datasets: 87 + (0 multi-transcript, ~1.1 transcripts per locus)
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.tracking Tue Jun 16 16:56:55 2020 +0000 @@ -0,0 +1,98 @@ +TCONS_00000001 XLOC_000001 Xkr4|Xkr4 c q1:CUFF.13|CUFF.13.1|100|10.695258|0.000000|25.820637|0.683544|- - +TCONS_00000002 XLOC_000002 - u q1:CUFF.1|CUFF.1.1|100|20.607936|0.000000|49.751960|1.317073|- - +TCONS_00000003 XLOC_000003 - u q1:CUFF.3|CUFF.3.1|100|27.255658|0.000000|65.800979|1.741935|- - +TCONS_00000004 XLOC_000007 - u q1:CUFF.29|CUFF.29.1|100|107.103219|71.402146|142.804292|6.845070|- - +TCONS_00000005 XLOC_000007 - u q1:CUFF.31|CUFF.31.1|100|122.650461|40.883487|204.417435|7.838710|- - +TCONS_00000006 XLOC_000007 - u q1:CUFF.33|CUFF.33.1|100|109.527366|26.732460|192.322273|7.000000|- - +TCONS_00000007 XLOC_000009 - u q1:CUFF.35|CUFF.35.1|100|96.747183|61.420107|132.074259|6.183206|- - +TCONS_00000008 XLOC_000009 - u q1:CUFF.37|CUFF.37.1|100|104.085013|53.596365|154.573660|6.652174|- - +TCONS_00000009 XLOC_000011 - u q1:CUFF.39|CUFF.39.1|100|23.912983|0.000000|51.525317|1.528302|- - +TCONS_00000010 XLOC_000012 - u q1:CUFF.41|CUFF.41.1|100|10.695258|0.000000|25.820637|0.683544|- - +TCONS_00000011 XLOC_000014 - u q1:CUFF.43|CUFF.43.1|100|10.561567|0.000000|25.497879|0.675000|- - +TCONS_00000012 XLOC_000015 - u q1:CUFF.45|CUFF.45.1|100|20.708956|2.186303|39.231609|1.323529|- - +TCONS_00000013 XLOC_000016 - u q1:CUFF.47|CUFF.47.1|100|20.607936|0.000000|49.751960|1.317073|- - +TCONS_00000014 XLOC_000017 - u q1:CUFF.49|CUFF.49.1|100|15.646767|0.000000|46.940300|1.000000|- - +TCONS_00000015 XLOC_000019 - u q1:CUFF.51|CUFF.51.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000016 XLOC_000020 - u q1:CUFF.53|CUFF.53.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000017 XLOC_000021 - u q1:CUFF.55|CUFF.55.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000018 XLOC_000022 - u q1:CUFF.57|CUFF.57.1|100|15.646767|0.000000|46.940300|1.000000|- - +TCONS_00000019 XLOC_000023 - u q1:CUFF.59|CUFF.59.1|100|15.646767|0.000000|46.940300|1.000000|- - +TCONS_00000020 XLOC_000025 - . q1:CUFF.7|CUFF.7.1|100|9.999117|0.000000|19.998234|0.639053|169 q2:CUFF.23|CUFF.23.1|100|16.205195|0.000000|34.917342|0.920455|88 +TCONS_00000021 XLOC_000026 - u q1:CUFF.5|CUFF.5.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000022 XLOC_000029 - . q1:CUFF.9|CUFF.9.1|100|17.776896|9.153835|26.399957|1.136139|404 q2:CUFF.25|CUFF.25.1|100|35.211287|0.000000|85.007567|2.000000|27 +TCONS_00000023 XLOC_000031 - u q1:CUFF.11|CUFF.11.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000024 XLOC_000034 Xkr4|Xkr4 i q1:CUFF.15|CUFF.15.1|100|10.695258|0.000000|25.820637|0.683544|- - +TCONS_00000025 XLOC_000035 Xkr4|Xkr4 i q1:CUFF.19|CUFF.19.1|100|29.337687|3.097262|55.578113|1.875000|- - +TCONS_00000026 XLOC_000037 Xkr4|Xkr4 i q1:CUFF.17|CUFF.17.1|100|8.710571|0.000000|21.029179|0.556701|- - +TCONS_00000027 XLOC_000040 Xkr4|Xkr4 i q1:CUFF.25|CUFF.25.1|100|14.567679|5.354270|23.781089|0.931034|- - +TCONS_00000028 XLOC_000042 Xkr4|Xkr4 i q1:CUFF.21|CUFF.21.1|100|13.851236|0.000000|33.439842|0.885246|- - +TCONS_00000029 XLOC_000044 Xkr4|Xkr4 i q1:CUFF.23|CUFF.23.1|100|23.470150|0.000000|50.571145|1.500000|- - +TCONS_00000030 XLOC_000046 Xkr4|Xkr4 i q1:CUFF.27|CUFF.27.1|100|34.253732|0.000000|73.806535|2.189189|- - +TCONS_00000031 XLOC_000048 Xkr4|Xkr4 i q1:CUFF.61|CUFF.61.1|100|45.263860|0.000000|97.530065|2.892857|- - +TCONS_00000032 XLOC_000049 Xkr4|Xkr4 i q1:CUFF.63|CUFF.63.1|100|15.646767|0.000000|46.940300|1.000000|- - +TCONS_00000033 XLOC_000052 Xkr4|Xkr4 i q1:CUFF.65|CUFF.65.1|100|15.362280|0.000000|37.087825|0.981818|- - +TCONS_00000034 XLOC_000057 Xkr4|Xkr4 i q1:CUFF.67|CUFF.67.1|100|12.998852|0.000000|31.382005|0.830769|- - +TCONS_00000035 XLOC_000059 Xkr4|Xkr4 i q1:CUFF.69|CUFF.69.1|100|10.058636|0.000000|24.283695|0.642857|- - +TCONS_00000036 XLOC_000059 Xkr4|Xkr4 i q1:CUFF.71|CUFF.71.1|100|8.621688|0.000000|20.814595|0.551020|- - +TCONS_00000037 XLOC_000060 Xkr4|Xkr4 i q1:CUFF.73|CUFF.73.1|100|15.362280|0.000000|37.087825|0.981818|- - +TCONS_00000038 XLOC_000062 Xkr4|Xkr4 i q1:CUFF.75|CUFF.75.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000039 XLOC_000065 Xkr4|Xkr4 i q1:CUFF.77|CUFF.77.1|100|16.248565|0.000000|39.227507|1.038462|- - +TCONS_00000040 XLOC_000068 Xkr4|Xkr4 i q1:CUFF.79|CUFF.79.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000041 XLOC_000071 Xkr4|Xkr4 i q1:CUFF.81|CUFF.81.1|100|13.201959|0.000000|31.872349|0.843750|- - +TCONS_00000042 XLOC_000073 Xkr4|Xkr4 i q1:CUFF.83|CUFF.83.1|100|13.201959|0.000000|28.446269|0.843750|- - +TCONS_00000043 XLOC_000075 Xkr4|Xkr4 i q1:CUFF.85|CUFF.85.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000044 XLOC_000077 Xkr4|Xkr4 i q1:CUFF.87|CUFF.87.1|100|17.243375|0.000000|41.629191|1.102041|- - +TCONS_00000045 XLOC_000080 Xkr4|Xkr4 i q1:CUFF.89|CUFF.89.1|100|16.567165|0.000000|39.996674|1.058824|- - +TCONS_00000046 XLOC_000083 Xkr4|Xkr4 i q1:CUFF.91|CUFF.91.1|100|31.293533|0.000000|75.549272|2.000000|- - +TCONS_00000047 XLOC_000084 Xkr4|Xkr4 i q1:CUFF.93|CUFF.93.1|100|21.664754|0.000000|52.303342|1.384615|- - +TCONS_00000048 XLOC_000085 Xkr4|Xkr4 i q1:CUFF.95|CUFF.95.1|100|46.940300|0.000000|101.142289|3.000000|- - +TCONS_00000049 XLOC_000086 Xkr4|Xkr4 i q1:CUFF.97|CUFF.97.1|100|21.481154|0.000000|46.285454|1.372881|- - +TCONS_00000050 XLOC_000087 Xkr4|Xkr4 i q1:CUFF.99|CUFF.99.1|100|14.567679|0.000000|35.169489|0.931034|- - +TCONS_00000051 XLOC_000004 - u - q2:CUFF.1|CUFF.1.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000052 XLOC_000005 - u - q2:CUFF.3|CUFF.3.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000053 XLOC_000006 - u - q2:CUFF.5|CUFF.5.1|100|21.226627|0.000000|59.889707|1.205672|- +TCONS_00000054 XLOC_000007 - u - q2:CUFF.11|CUFF.11.1|100|32.531777|24.582998|40.480555|1.847804|- +TCONS_00000055 XLOC_000008 - u - q2:CUFF.7|CUFF.7.1|100|29.709524|19.806349|39.612698|1.687500|- +TCONS_00000056 XLOC_000009 - u - q2:CUFF.9|CUFF.9.1|100|34.072933|23.364686|44.781179|1.935341|- +TCONS_00000057 XLOC_000010 - u - q2:CUFF.13|CUFF.13.1|100|16.582060|0.000000|35.729373|0.941860|- +TCONS_00000058 XLOC_000013 - u - q2:CUFF.15|CUFF.15.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000059 XLOC_000018 - u - q2:CUFF.17|CUFF.17.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000060 XLOC_000024 - u - q2:CUFF.19|CUFF.19.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000061 XLOC_000027 - u - q2:CUFF.21|CUFF.21.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000062 XLOC_000028 - u - q2:CUFF.26|CUFF.26.1|100|29.709524|0.000000|71.725135|1.687500|- +TCONS_00000063 XLOC_000030 - u - q2:CUFF.29|CUFF.29.1|100|13.581496|0.000000|32.788633|0.771429|- +TCONS_00000064 XLOC_000032 - u - q2:CUFF.31|CUFF.31.1|100|22.635827|0.000000|45.271655|1.285714|- +TCONS_00000065 XLOC_000033 Xkr4|Xkr4 i - q2:CUFF.33|CUFF.33.1|100|23.767619|0.000000|57.380108|1.350000|- +TCONS_00000066 XLOC_000036 Xkr4|Xkr4 i - q2:CUFF.35|CUFF.35.1|100|11.317914|0.000000|27.323861|0.642857|- +TCONS_00000067 XLOC_000037 Xkr4|Xkr4 i - q2:CUFF.37|CUFF.37.1|100|11.500461|0.000000|24.780049|0.653226|- +TCONS_00000068 XLOC_000038 Xkr4|Xkr4 i - q2:CUFF.39|CUFF.39.1|100|52.816931|0.000000|113.804669|3.000000|- +TCONS_00000069 XLOC_000039 Xkr4|Xkr4 i - q2:CUFF.41|CUFF.41.1|100|43.213852|0.000000|93.112911|2.454545|- +TCONS_00000070 XLOC_000041 Xkr4|Xkr4 i - q2:CUFF.43|CUFF.43.1|100|23.474191|0.000000|46.948383|1.333333|- +TCONS_00000071 XLOC_000043 Xkr4|Xkr4 i - q2:CUFF.45|CUFF.45.1|100|20.667495|0.000000|49.895746|1.173913|- +TCONS_00000072 XLOC_000045 Xkr4|Xkr4 i - q2:CUFF.47|CUFF.47.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000073 XLOC_000047 Xkr4|Xkr4 i - q2:CUFF.49|CUFF.49.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000074 XLOC_000049 Xkr4|Xkr4 i - q2:CUFF.51|CUFF.51.1|100|14.948188|7.228977|22.667399|0.849057|- +TCONS_00000075 XLOC_000050 Xkr4|Xkr4 i - q2:CUFF.53|CUFF.53.1|100|52.816931|0.000000|113.804669|3.000000|- +TCONS_00000076 XLOC_000051 Xkr4|Xkr4 i - q2:CUFF.55|CUFF.55.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000077 XLOC_000053 Xkr4|Xkr4 i - q2:CUFF.57|CUFF.57.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000078 XLOC_000054 Xkr4|Xkr4 i - q2:CUFF.59|CUFF.59.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000079 XLOC_000055 Xkr4|Xkr4 i - q2:CUFF.61|CUFF.61.1|100|13.204233|0.000000|31.877838|0.750000|- +TCONS_00000080 XLOC_000056 Xkr4|Xkr4 i - q2:CUFF.63|CUFF.63.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000081 XLOC_000058 Xkr4|Xkr4 i - q2:CUFF.65|CUFF.65.1|100|31.170648|0.000000|62.341295|1.770492|- +TCONS_00000082 XLOC_000059 Xkr4|Xkr4 i - q2:CUFF.69|CUFF.69.1|100|18.799247|8.750627|28.847866|1.067797|- +TCONS_00000083 XLOC_000060 Xkr4|Xkr4 i - q2:CUFF.67|CUFF.67.1|100|15.681351|3.378764|27.983938|0.890700|- +TCONS_00000084 XLOC_000061 Xkr4|Xkr4 i - q2:CUFF.71|CUFF.71.1|100|22.635827|0.000000|54.647722|1.285714|- +TCONS_00000085 XLOC_000063 Xkr4|Xkr4 i - q2:CUFF.73|CUFF.73.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000086 XLOC_000064 Xkr4|Xkr4 i - q2:CUFF.75|CUFF.75.1|100|52.816931|0.000000|113.804669|3.000000|- +TCONS_00000087 XLOC_000066 Xkr4|Xkr4 i - q2:CUFF.77|CUFF.77.1|100|17.605644|0.000000|52.816931|1.000000|- +TCONS_00000088 XLOC_000067 Xkr4|Xkr4 i - q2:CUFF.79|CUFF.79.1|100|13.390208|0.000000|32.326821|0.760563|- +TCONS_00000089 XLOC_000068 Xkr4|Xkr4 i - q2:CUFF.81|CUFF.81.1|100|11.211141|1.183592|21.238690|0.636792|- +TCONS_00000090 XLOC_000069 Xkr4|Xkr4 i - q2:CUFF.83|CUFF.83.1|100|21.126772|0.000000|51.004540|1.200000|- +TCONS_00000091 XLOC_000070 Xkr4|Xkr4 i - q2:CUFF.85|CUFF.85.1|100|19.014095|0.000000|38.028190|1.080000|- +TCONS_00000092 XLOC_000072 Xkr4|Xkr4 i - q2:CUFF.87|CUFF.87.1|100|24.170460|0.000000|52.080103|1.372881|- +TCONS_00000093 XLOC_000074 Xkr4|Xkr4 i - q2:CUFF.89|CUFF.89.1|100|29.709524|0.000000|64.015126|1.687500|- +TCONS_00000094 XLOC_000076 Xkr4|Xkr4 i - q2:CUFF.91|CUFF.91.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000095 XLOC_000078 Xkr4|Xkr4 i - q2:CUFF.93|CUFF.93.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000096 XLOC_000079 Xkr4|Xkr4 i - q2:CUFF.95|CUFF.95.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000097 XLOC_000081 Xkr4|Xkr4 i - q2:CUFF.97|CUFF.97.1|100|35.211287|0.000000|85.007567|2.000000|- +TCONS_00000098 XLOC_000082 Xkr4|Xkr4 i - q2:CUFF.99|CUFF.99.1|100|19.602160|0.000000|39.204320|1.113402|-