| Previous changeset 3:e320ef2d105a (2019-09-05) Next changeset 5:378565f5a875 (2019-11-22) |
|
Commit message:
Uploaded |
|
modified:
fasta_interlacer.py |
|
added:
Galaxy_integration.org README.html __init__.py clean_repository_tarball.sh example1.sh example2.sh name_affixer.py test_data/test_run1 tmp/.dummy |
|
removed:
test_data/VTS_contigs.info.minRD5 test_data/prefix_suffix.fasta |
| b |
| diff -r e320ef2d105a -r d397f5a85464 Galaxy_integration.org --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/Galaxy_integration.org Wed Sep 18 06:30:04 2019 -0400 |
| [ |
| @@ -0,0 +1,58 @@ + +#+BEGIN_SRC sh +/home/petr/anaconda3/bin/planemo shed_init --name=repeatexplorer_utilities \ + --owner=repeatexplorer \ + --description="some utilities for data preprocessing" \ + --long_description="some utilities for data preprocessing" \ + --category="Fasta Manipulation" +#+END_SRC +# this create file .shed.yml + + +# create repository manualy: +#+BEGIN_SRC sh +planemo shed_build . +#+END_SRC + +#+RESULTS: +: Created: /mnt/raid/users/petr/workspace/re_utilities.tar.gz + +this creates archive: ../re_utilities.tar.gz +unecessary file from archive are deleted; + +#+BEGIN_SRC bash :tangle clean_repository_tarball.sh + #!/bin/bash + containsElement () { + local e match="$1" + shift + for e; do [[ "$e" == "$match" ]] && return 0; done + return 1 + } + ARCHIVE_GZ=../re_utilities.tar.gz + TMP_TAR=`mktemp` + echo $TMP_TAR + ARCHIVE_GZ_CLEAN=../re_utilities_clean.tar.gz + zcat $ARCHIVE_GZ > $TMP_TAR + ARCHIVE_FILE_LIST=`tar -tz -f ../re_utilities.tar.gz` + ls -l $TMP_TAR + GIT_LIST=`git ls-files` + + for FILE in $ARCHIVE_FILE_LIST + do + containsElement $FILE ${GIT_LIST[@]} + if [ $? != 0 ] + then + echo "Deleting ${FILE}" + tar --delete -f $TMP_TAR $FILE + fi + done + ls -l $TMP_TAR + echo "compressing.." + gzip -c $TMP_TAR > $ARCHIVE_GZ_CLEAN + echo "output in ${ARCHIVE_GZ_CLEAN}" + +#+END_SRC + + + + |
| b |
| diff -r e320ef2d105a -r d397f5a85464 README.html --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/README.html Wed Sep 18 06:30:04 2019 -0400 |
| b |
| @@ -0,0 +1,70 @@ +<?xml version="1.0" encoding="UTF-8" ?> +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" + "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> + +<html xmlns="http://www.w3.org/1999/xhtml"> + +<head> +<title>README.html</title> + +</head> + +<body> + +<h1>RepeatExplorer utilities</h1> +<p>This repository include utilities for preprocessing of NGS data to suitable format for RepeatExplorer and TAREAN +analysis. Each tool include also XML file which define tool interface for Galaxy environment</p> +<h2>Available tools</h2> +<h3>Paired fastq reads filtering and interlacing</h3> +<p>tool definition file: <code>paired_fastq_filtering.xml</code></p> +<p>This tool is designed to make memory efficient preprocessing of two fastq files. Output of this file can be used as input of RepeatExplorer clustering. Input files can be in GNU zipped archive (.gz extension). Reads are filtered based on the quality, presence of N bases and adapters. Two input fastq files are procesed in parallel. Only complete pair are kept. As the input files are process in chunks, it is required that pair reads are complete and in the same order in both input files. All reads which pass the quality filter fill be writen into output files. If sampling is specified, only sample of sequences will be returned. Cutadapt us run with this options:</p> +<p><code>--anywhere='AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT' +--anywhere='AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTAGATCTCGGTGGTCGCCGTATCATT' +--anywhere='GATCGGAAGAGCACACGTCTGAACTCCAGTCAC' +--anywhere='ATCTCGTATGCCGTCTTCTGCTTG' +--anywhere='CAAGCAGAAGACGGCATACGAGAT' +--anywhere='GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC' +--error-rate=0.05 +--times=1 --overlap=15 --discard</code></p> +<p>Order of fastq files processing</p> +<ol> +<li>Trimming (optional)</li> +<li>Filter by quality</li> +<li>Discard single reads, keep complete pairs</li> +<li>Cutadapt filtering</li> +<li>Discard single reads, keep complete pairs</li> +<li>Sampling (optional)</li> +<li>Interlacing two fasta files</li> +</ol> +<h3>single fastq reads filtering</h3> +<p>tool definition file: <code>single_fastq_filtering.xml</code></p> +<p>This tool is designed to perform preprocessing +of fastq file. Input files can be in GNU zipped archive (.gz extension). Reads +are filtered based on the quality, presence of N bases and adapters. All reads +which pass the quality filter fill be writen into output files. If sampling is +specified, only sample of sequences will be returned. </p> +<h3>fasta afixer</h3> +<p>tool definition file: <code>fasta_affixer.xml</code></p> +<p>Tool for appending prefix and suffix to sequences names in fasta formated sequences. This tool is useful +if you want to do comparative analysis with RepeatExplorer and need to +append sample codes to sequence identifiers</p> +<h2>Dependencies</h2> +<p>R programming environment with installed packages <em>optparse</em> and <em>ShortRead</em> (Bioconductor) +python3 +cutadapt</p> +<h2>License</h2> +<p>Copyright (c) 2012 Petr Novak (petr@umbr.cas.cz), Jiri Macas and Pavel Neumann, +Laboratory of Molecular Cytogenetics(http://w3lamc.umbr.cas.cz/lamc/) +Institute of Plant Molecular Biology, Biology Centre AS CR, Ceske Budejovice, Czech Republic</p> +<p>This program is free software: you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation, either version 3 of the License, or +(at your option) any later version.</p> +<p>This program is distributed in the hope that it will be useful, +but WITHOUT ANY WARRANTY; without even the implied warranty of +MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +GNU General Public License for more details. +You should have received a copy of the GNU General Public License +along with this program. If not, see <a href="http://www.gnu.org/licenses/">http://www.gnu.org/licenses/</a>.</p> +</body> +</html> |
| b |
| diff -r e320ef2d105a -r d397f5a85464 clean_repository_tarball.sh --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/clean_repository_tarball.sh Wed Sep 18 06:30:04 2019 -0400 |
| [ |
| @@ -0,0 +1,29 @@ +#!/bin/bash +containsElement () { + local e match="$1" + shift + for e; do [[ "$e" == "$match" ]] && return 0; done + return 1 +} +ARCHIVE_GZ=../re_utilities.tar.gz +TMP_TAR=`mktemp` +echo $TMP_TAR +ARCHIVE_GZ_CLEAN=../re_utilities_clean.tar.gz +zcat $ARCHIVE_GZ > $TMP_TAR +ARCHIVE_FILE_LIST=`tar -tz -f ../re_utilities.tar.gz` +ls -l $TMP_TAR +GIT_LIST=`git ls-files` + +for FILE in $ARCHIVE_FILE_LIST +do + containsElement $FILE ${GIT_LIST[@]} + if [ $? != 0 ] + then + echo "Deleting ${FILE}" + tar --delete -f $TMP_TAR $FILE + fi +done +ls -l $TMP_TAR +echo "compressing.." +gzip -c $TMP_TAR > $ARCHIVE_GZ_CLEAN +echo "output in ${ARCHIVE_GZ_CLEAN}" |
| b |
| diff -r e320ef2d105a -r d397f5a85464 example1.sh --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/example1.sh Wed Sep 18 06:30:04 2019 -0400 |
| b |
| @@ -0,0 +1,4 @@ +#!/bin/sh +# minimal example : +./ChipSeqRatioAnalysis.py -c test_data/seq_C_10k -i test_data/seq_I_10k -k test_data/VTS_contigs.info.minRD5 -o tmp/chi_seq_ratio.csv -ht tmp/html_report.html +# output goes to current directory |
| b |
| diff -r e320ef2d105a -r d397f5a85464 example2.sh --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/example2.sh Wed Sep 18 06:30:04 2019 -0400 |
| b |
| @@ -0,0 +1,4 @@ +#!/bin/sh +# minimal example2 : +./ChipSeqRatioAnalysis.py -c test_data/seq_C_20k -i test_data/seq_I_10k -k test_data/VTS_contigs.info.minRD5 -o tmp/chi_seq_ratio2.csv -ht tmp/html_report2.html --max_cl 300 +# output goes to current directory |
| b |
| diff -r e320ef2d105a -r d397f5a85464 fasta_interlacer.py --- a/fasta_interlacer.py Thu Sep 05 09:04:56 2019 -0400 +++ b/fasta_interlacer.py Wed Sep 18 06:30:04 2019 -0400 |
| [ |
| @@ -79,7 +79,7 @@ "last character of sequence id must be used for distinguishing pairs!") exit(1) # check first thousand! - for i in range(1000): + for i in range(3): seqA = readSingleSeq(fA) seqB = readSingleSeq(fB) if (not seqA) or (not seqB): @@ -105,10 +105,8 @@ buffB_names = [] while True: - seqA = readSingleSeq(fA) seqB = readSingleSeq(fB) - if not seqA and not seqB: break # end of file @@ -138,8 +136,8 @@ for i in buffA_names: seqtmp = {"name": i + charA, "sequence": buffA[i]} writeSingleSeq(single, seqtmp) - buffA = {} - buffA_names = [] + buffA = {} + buffA_names = [] j = 0 for i in buffB_names: @@ -166,8 +164,8 @@ for i in buffB_names: seqtmp = {"name": i + charB, "sequence": buffB[i]} writeSingleSeq(single, seqtmp) - buffB = {} - buffB_names = [] + buffB = {} + buffB_names = [] j = 0 for i in buffA_names: @@ -183,17 +181,19 @@ else: buffB[seqB["name"][:-1]] = seqB['sequence'] buffB_names.append(seqB["name"][:-1]) - fA.close() - fB.close() - fPairs.close() - # write rest of singles: + + fA.close() + fB.close() + fPairs.close() + + # write rest of singles: for i in buffA: seqtmp = {"name": i + charA, "sequence": buffA[i]} writeSingleSeq(single, seqtmp) for i in buffB: seqtmp = {"name": i + charB, "sequence": buffB[i]} writeSingleSeq(single, seqtmp) - single.close() + single.close() if __name__ == "__main__": |
| b |
| diff -r e320ef2d105a -r d397f5a85464 name_affixer.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/name_affixer.py Wed Sep 18 06:30:04 2019 -0400 |
| [ |
| @@ -0,0 +1,31 @@ +#!/usr/bin/env python +import sys + +from optparse import OptionParser + +parser = OptionParser() +parser.add_option("-f", "--fastq", dest="fastq", help="fastq file") +parser.add_option("-p", "--prefix", dest="prefix", help="prefix to be added to names") +parser.add_option("-s", "--suffix", dest="suffix", help="suffix to be added",default='') +parser.add_option("-n", "--nspace", dest="nspace", help="number of spaces to ignore",default='0') +options, args = parser.parse_args() +nspace=int(options.nspace) + +f=open(options.fastq,"r") +j=0 +for oneline in f: + if oneline=="": + continue + j+=1 + if j==5: + j=1 + if not oneline: + break + + if (oneline[0]=="@" and j==1) or (oneline[0]=="+" and len(oneline)>2 and j==3): + header=" ".join(oneline.split()[:1+nspace]) + header_out=header[0]+options.prefix+header[1:]+options.suffix+"\n" + sys.stdout.write(header_out) + else: + sys.stdout.write(oneline) + |
| b |
| diff -r e320ef2d105a -r d397f5a85464 test_data/VTS_contigs.info.minRD5 --- a/test_data/VTS_contigs.info.minRD5 Thu Sep 05 09:04:56 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 |
| b |
| b'@@ -1,428103 +0,0 @@\n->CL41Contig243 (655-2038.9-1335472)\n-ATTGCCAGCTTACACCTCAGAGAGCCGCAAAACTTATTCCCCCAGTGGAGTTGCCGGCTC\n-ATTACCCACAAAATTCAGGTCATGTTCACAACATGTTGCACTCACATGTTCATCATATAT\n-AAAGCACTTTCATTGTGGTATCTTAGGGGCTACAAATTGTGCATTCTGTTGTATTTAAGT\n-CTACCTTCATCCCTCCGGATCGAAGAAGCCTTAAATAGGGGCATCTGTCATACCCAATTT\n-TTTGACCCCCCTAAGATCCCACATGCATCATATCAAGCATTGCATTCACATCAACCATAA\n-GATCATACTTATGTTGATAGTTGATTCTCTCTCCTAAGGATTGGCCTCCTCTGAGGACAC\n-TTAGGCATTTGAGACTCCATGATCATCACTAACCACTAATCATATCCATCATACAAAGAT\n-ATGCTTGTGCTTTGTTCATCTCACCCATTGCAGGTAAAGTTGGCTCAAAGGATCATCAAA\n-GCATCCCCAACTAGGGTTTAGTCCCTCAGATGAATTTCTTCCTCAAAATGGTGAAGCATC\n-AATGCTAATCAACCACCATGGACCACAAGTCACTCTATGCATCAAGGACAACCTTCACAA\n-CATAGCTTTTGACCTACAATGGCTCAGAAGTTTTGAAGCTTGCTCCCCAAGGCAT\n->CL64Contig342 (450-1716.7-772504)\n-TTATGATATGGACCTCATTCTTTCACTCAAGTAATCTTCAATGAGTCCTCTACAAATGTT\n-ATGTTGAGAGCAAAGATCAATTATGGATGGAAGCATTTCAATTTCATCAAGACATTACAG\n-GTCATTTTGTCTGTCTATAGGAGCCTAAGTTCACAAGATCATATCTCTCTCAATTTTCAT\n-GATATGGAGGTGATTCCTTTTCTCAACATCTTGTATTAACCTCCTCTTCAATTCATCTTC\n-ACAAAGTTTGGGCCAAAAATGCATAGAAAAGGTAGCTATCCATGGGAGCTAGGCTTATTT\n-CTTAACTTAGTGAAATTTTCAAGGCATAAAAACCCTAAAATCTGCCAACTTCATGGCCAA\n-ATTCCACCAAGATGCAAGGTGAATCAAGTCAAGTCATCAACTACAAAGTTGAAGAGCATC\n-CCCTCACCTTTCTAAAAAGTCCAAGATCAC\n->CL56Contig206 (406-1569.1-637072)\n-TCAACCACCATGGACCACAAGTCACTCTATGCATTAAGGACAATCTTCACAACATAGCTT\n-GTTGACCTACAATGGCTCAAAAGCTTGAAGCTTGCTCCTCAAGGCATCATCAGCTAGGGT\n-TTGATCCTCTGATGAATTGGATCATCAGTCTAAGGACTTCCAACACAAGACAATCATAAA\n-TGGTCTTAGGACTCACTTGTCTTAAGCATTAACACCATTTCATGACCATTCAAGTCATTC\n-TCAGCTCAAGGGTTAGGAAAACCCAGTCATCAAGAGATCTCCTAAGTCCTCAGCATCACC\n-TGAGACTTGGTCCACAAGCAATCTCAATGACTACATCAAATTATCCAAGGATTCCTTCAT\n-GGATCATTCCATTGACTTCATGTCATCTCATTCAATTGAGAAAAGT\n->CL81Contig246 (411-1562.1-642013)\n-TCCGTAGTACTCTACGGGAGCCCGTTAGTTTTTATTCTGTCTTTTACTGCTTTAACTGTT\n-TGTGTGTTTGCCGTTCATGGGCAAAGTTTCTCCTTTGTGTTTTTTATATGGGTGAGGACC\n-TAAAGATGTGTGGATCTTATCATGCTCGCAGAGATCTCTGGATCTCATGCCTACATATCC\n-CTAGAGGAAGTCAGAGCATCTGTAGTTCGGGGACTACGGGTATCTAGGGTTTGTTTTGAT\n-GGTGATTATGCTCGCAGAGGATAGACCTCGTGCCTACGTATTCTCAAAAGGAAGTTGAGA\n-AATCAGAGCAAACGTAGTTCCCACTTATATGCCAGTGGAGGCAGAGTTGTGTTTTTGTTT\n-TTGAAGAGACGAAGACAAGTCGTCATTTTAATAAAATATAAATCACCTACT\n->CL85Contig254 (680-1497.5-1018285)\n-TATCGATTGGATCGCCTCTATACACGGCAAGGATCCCACTTAAGGATCCAGGGAAGGCAT\n-CTTCTCGATAGGCGACAAGGATCCGACATAAGGATTCAATGAATGGATTTCCTCTTTAGA\n-CGGCAAGGATCCAGCGAATGGCTCGCATTTATAAGCAGCAAGGATCCCATTAAAGGATCT\n-ATCGAATGGATCGCCTCTATAGGCGGCAAGGATCCCACTTAAGGATCCAGGGAATGGATC\n-TTCTCGATAGGCGACAAGGATCCCACATAAGGATTCAATGAATGGATTTCCTCTTTAGAC\n-GGCAAGGATCCCACTTAAGGATCCAGGGAATGGATCTTCTCGATAGGCGACAAGGATCCG\n-ACATAAGGATTCAATGAATGGATTTCCTCTTTAGACGGCAAGGATCCAGCGAATGGCTCG\n-CATTTATAAGCAGCAAGGATCCCATTAAAGGATCTATCGAATGGATCGCCTCTATAGGCG\n-GCAAGGATCCCACTTAAGGATCCAGGGAATGGATCTTCTCGATACGCGACAAGGATCCGA\n-CATAAGGATTCAATGAATGGATTTCCTCTTTAGTCGGCAAGGATCTAGCGAATCGATTAC\n-ATTTGTAAGCCACAAGGATCCCACTAAAGAATACAGGGAATGGATCGCCTCAATAGGCGA\n-CAAGGATCCCACTTAAGGAT\n->CL136Contig14 (845-1299.1-1097761)\n-TACATACACAATTAAAGCGTAGTACTTAACAGCTGCGATCATACCAGCACTAATGTACTG\n-GATCCCATCAGAACTCCGCAGTTAAGCGTGCTTGGGCGAGAGTAGTACTAGGATGGGTGA\n-CCTCCTGGGAAGTCCTTGTGTTGTTATTTTCGTTAATCGGCGTGAAATGTTATCGCTAGT\n-AAGGACTGATTTTTTTTATATTTTTGTTAATCGGCGTGAAATATTATCGGTAGTGAGGTT\n-AATCGGCGTGAAATATTATCGGTAGTGAGGATAAAACAAAAATAGTTGCGTTAATTTTAT\n-TGAGTGGTGTGATAAATGTCGTAAGTAAGGAGAGAAGATATATACGGTTCGTTTTGTTTA\n-CCGGTAAAAGAAAATGGAGTTAGTAAGGAGAAAAGGAATATACATACACAATTAAAGCGT\n-AGTACTTATCAGGTGCGATCATACCAGCACTAATGCACCGGATCCCATCAGAACTCCGCA\n-GTTAAGCGTGCTTGGGCGAGAGTAGTACTAGGATGGGTGACCTCCTGGGAAGTCCTTGTG\n-TTGCACCTCTTTTTTTGTTTTTGCAAGGAATTCTTTTTATTTTCGTTAATCGGCGTGAAA\n-TGTTATCGCTAGTAAGGACTGATTTTTTTATATTTTTGTTAATCGGCGTGAAATATTATC\n-GGTAGTGAGGTTAATCGGCGTGAAATATTATCGGTAGTGAGGTTAATCGGCGTGAAATAT\n-TATCGGTAGTGAGGATAAAACAAAAATAGTTGCGTTAATTTTATTGAGTGGTAGATAAAT\n-GTCGTAAGTAAGGACAGAAGATATATACGTTTCGTTTTGTTTACCGGTAAAAGAAAATGG\n-AGTTA\n->CL33Contig1110 (580-1253.5-727047)\n-TGCTGGTTTGCTTACTCAAAGGTATGTACAAAGGTTAAAGGATAAAAAAGGGTAAGTATT\n-TTATTTACAATTTTGATTCGATAGAGTTTAAAGACTCTATGCCTACGTACTCCCTCGTGC\n-AATGGGAAAGTCAGAATCACGTAGTTCTGC'..b'TT\n-AACTGGGAGAAAGTTTGAATTGATTTGATTTGAAAAGATATTTTTGAATGTGGTGTGATA\n-GAAGAAAATATTTTATAGAGATGAATCATGTCTAAATCTCTAAGTGTTTTTGATGATTTT\n-TTGAATGGGTGTTGAGGTGCTTTTAAGAAGCCTGAGGGTCCTTGTATGAA\n->CL11Contig153 (258-5.0-1283)\n-TCAATAACTATAAATTGAGTAAACTCTTTAGTTTCATTGAGACAAAACCTTTTGGAAAAG\n-GTTCTGACTTTTGGATGCTTCTGATTTATGAGAACTTCACTGTTCAAAAGTCTCACACTT\n-CAGAAACGCTTAAGCCTCAGAATCTCATTTTATCGTTCTGAAAGATAAAAGTTATGAAAT\n-TCTCATTTCTTTCATTATGGACATAATTCTATGTTCAAATGTTTGAGAATGAAAACAAAT\n-CTATCTTCAGCTATGGAT\n->CL207Contig19 (159-5.0-795)\n-CAAGTTATCGTGATTTCTTAAATAATTTTTGATCAATTCATCCTTAATCCATATTTCCTT\n-GAAGAATAGGTCATGACAAATGCCTTTGGCCCAAAAATTCCAATTTTGTCCCTTTTGAAT\n-CAGATGAAATTTCATCAATCAAAACCCTAATTTGACCGG\n->CL33Contig936 (112-5.0-558)\n-GATAGAAAAAAGGGTTAGGGATCACTTGACGGAGTCAAGGTTTATTATGAAAGCATTTGA\n-AGGGTTAAGGTAAAGGTTAAGGATATATGTACAAGGATAGGGTTGAGATGTT\n->CL25Contig1626 (159-5.0-797)\n-CATTCCACCCGCCAAGGGGATTTCGCCTCCTTCAATGAGGTCTTAATCCACTATAACCAA\n-ACTGATTACAAACACAAAGACCTACTGTCAATGTCTTCTTGAGAATCTGACTATACCCTA\n-GTCTCTCAAGGAATATCAACTCACAAGTTAAGATACAAA\n->CL175Contig200 (149-5.0-741)\n-TGATTGTGTTTGTGAAGTTTAAAATGGGAGTTTTGGGTAAAAAGGGTTGAAAACTCGTAT\n-TTGGAAAAATGGTGATTTTTGGGTTCTGAACTGGTGGGAACCGGTTACCACCCTGGAGGG\n-AACCGGTTAGCACAAGAGACAGTAGCGTG\n->CL86Contig1868 (186-5.0-924)\n-TGGATTATTGACACGGGTAACTGATCCTTGTATCATGCTTCAACTTTCTTATGGATTTCA\n-GATTCTCCATGGGACTCGTGCCCAATTATAAGAATTATCGCTCCCTATGGGACTCAAACC\n-CATTTCGAGCCTTCTTCCTCGTATGGGACTCTAACCCACTTTGACGTCCACCGCTCTATG\n-GGACTC\n->CL62Contig2990 (143-5.0-719)\n-TGAAGAGAAGCACAAAGGCTCAGATTTTGAGTTGTCAAAATCTAAAGTTCCAGACTCTGA\n-AGGTTCATCAGACTCTAAAGTCTATGCAACTTCTGAAGGAGGCTCAGACTCTGAAGAGTG\n-AGACTTTATCTCAACCTCTGAAG\n->CL6320Contig1 (222-5.0-1102)\n-TCCTTTCTGAACAAATAATATAACTGAACATGGATTAAAAAAAATCTGCAACTCATAGCT\n-AGTACATTCCGAGTTAACATATCCGACTATCATGTTTAAATTACTAACTTTTACTGTTAA\n-AACTAGTCTAACTTAAAATTGGCACAAGAGGATTGACATTAACTTAGCTTAATTTGCTAC\n-AACCAATTACTAATACTCTAACTGTTTTGTAACTAACACTAA\n->CL9Contig1287 (136-5.0-675)\n-TAACTCTTTGGTACGTTACCTCGGTGTTCTTTAAACCAAAAGGCATAACCTTGTAGCAAA\n-AAGTTCCCCATGGGGTAATAAATGTGGTCTTTTCCATATCTTCTGGCGCCATCCTGATTT\n-GATTGTATCCCAAAAA\n->CL71Contig2784 (156-5.0-780)\n-GTTTTTTCTAGTTCGGGTATGTTTTCGGATATACATATCTGAATTAACCTCACGTATTTA\n-AAAAAAATCAATCAGGGATAGTTTCGGATATTAATCTCCGAAGCTACCCTTTTTTGTAAA\n-AAAAGTTTTGGTTTATTCGGACATGAATATCCGAAA\n->CL76Contig1338 (133-5.0-663)\n-GTTTTAACCTCTTAGATGTATGGCTAGGTTCTCACTCTTAGATGTATGCTTAAGATATTT\n-TATTTCCATGTCCATCTGATGATTATTTCTCTTGATTATGATTATCATGTGAATGATGCT\n-TGTATGTATGAAT\n->CL4207Contig1 (348-5.0-1724)\n-ATTGGGCATGGGTAATATAAGATTGTAAGTAAATGGTAGGGCACACATCATCACTAGAGT\n-TTTCTTTGTGCCTGAATTGAAGAATAATTTACTAAGCATTGGCCAACTGCAAGCCAAAAA\n-CCTAACTATCATGTTTCAAGGTGGAAAATGCAAAGTATTTCATCCTGAAAAAGGTCTCCT\n-CATGGAAACAGATATGGCTGCCAATCGAATGTTCATTCTGAATGTTGTGGCTCAACCTAA\n-GTCTTCCACGTGTTTTAACACTGTTGCTGAGGATGAATGGTATCTTTGGCATTGTCGGTA\n-CGGACATTTAGGTTTCAAGGGACTCCAAACTCTTGTTCAGAAAAAGAT\n->CL1Contig321 (162-5.0-806)\n-CCATGTACTTGTATACTACAGAAACTTGTTTTTCTGCACATAGTGCGATCCTGCACTCCA\n-ATTTTCCATCCCACGCATCCATAATCAATCCAAAACCTTACATACATGTTATATATACCT\n-CTATAACATCTAATTGCATCAAAACAACACCTTAATTCATGG\n->CL44Contig2458 (116-5.0-583)\n-TCGAAATTGTTGCCATTTTCTTTAAACGATAAAAAACAAAACGAAAATGTAACAAAACGC\n-TTTTATTAATGATGATATTCATTTAAAACAAAAGGGGCCCTACATATGATCGTTTA\n->CL76Contig880 (167-5.0-838)\n-CTAAACACTTGTCATCTTTCTAAACCCGATCCTTTTTTCAAGATCTTCGCCCCATTTAAA\n-CATTTTCAAGGTATCAGTCGGTCTTTTTATCCAATTTTAGTTAGAGTCTTTGATTTTCTT\n-AAACAAATGAGTCTTAATCAATCTTAAATAAGTGAATTGATCGTCAC\n->CL67Contig523 (262-5.0-1308)\n-ATCATTGCTAGTGTTATGTTTAAGTCTTTTCCCATTTACAATAAAAGGACCCGACGATTG\n-GTTTTGAATCTCAACAGCACCAGATTGCATAACTTTTGTAACTCTGAATGGCCCAGTCCA\n-TCTAGAGCGTAACTTACCAGGGAAAAGGCGCAATCTGGAGTTGAAAAGGAGCACTAAGTC\n-GCCTGCAGTAAAGTCTTTCTTCGTAATTTGCCTATCATGCCATTTCTCAGTGCGTTCCTT\n-ATAGATTACAGTGTTTTCGTAG\n->CL60Contig1946 (146-5.0-733)\n-AAGTCCAGAAGATTCCCAAAGAAGTCTTCTGGGATCAGAATTCAGGAACCATCTTCTAAA\n-GTTCCTTCAGGTATGCCACCTATACCTTATGTTACTTTTGTTTCTATTTCAGTTTATATG\n-CCAACATCTGAACCAACTCCCCCAAT\n->CL32Contig2720 (132-5.0-664)\n-ATTATTTGTTGGTCGTTTGTTTTCAAAAAAACATTTCCGTGGTTGATTCCGGATCTGTTA\n-TTCCTCAGTGTTTCCGTTGATGATTCCAGATTCACTTGTTTCCAGCTTATAGGATATTTC\n-TGTTGATGATTC\n' |
| b |
| diff -r e320ef2d105a -r d397f5a85464 test_data/prefix_suffix.fasta --- a/test_data/prefix_suffix.fasta Thu Sep 05 09:04:56 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 |
| b |
| b'@@ -1,267832 +0,0 @@\n->TESTERR215189.1OK\n-GTGAGGAGTAGTCATCGGCCGCCAACTGATCCTGGAGTGAAGTGTGTCTCGACCCTGTCCTCAATTCCCATGGTAGTACCGCCTCGTGACCATAC\n->TESTERR215189.1OK\n-CCGATCAAGGAACGATGTTTACGTCAGGAGAATTCGAAGAGTTCACGGCCGATATGGGAATCAAATTATTAAACTCCTCTCCGTACTATGCCCAA\n->TESTERR215189.2OK\n-ATATGCTTAAACTCAGCGGGTAGTCCCGTCTGACCTGGGGTCGCGGTCCGAGCGCCTAAACGCTTCGGTCGTGTATGGGTCCTTGAGGCCAATAA\n->TESTERR215189.2OK\n-AAGTTTGGGCTGCCGGCATAACTTGTCGGGCACCGCACGTGGTGGGCGACATCTAGTTGTTCTCGGTGCAGTGTCCTGGCATGCAGCTGGCTTAT\n->TESTERR215189.3OK\n-ATCTGGGGCCAACCTCTTGGTTTGGGGCCGGTATGTGCACTTGGAGGGCTCGGGGCGGCCGGGGTTGGAGAGGGGAGAGCACGATGGAGGGCAGT\n->TESTERR215189.3OK\n-CTGCCGCCGCCCCTCTTTCCTCCCCTGCGCCAAGCCTGCCGCAGCCCCCGGTGACGAGGACCTCAGCGATAGTCCAGTTCCCCTCTTGCGACGAT\n->TESTERR215189.4OK\n-AATACCCGAAATGCATGCCGCCTCACATTTCCGCATACGCAGTCTTATTACTAGTACTGGTAGTAGCAGTAGTAGGACTGCAGGAGACCTGCAGC\n->TESTERR215189.4OK\n-CGGACAGGCTAGCCAGGCAAAGTCATTAAAGGCCAACGACGACGACGACGTGGGCCGTGGCCGAGTCGCACGTTGGTAGCTGGGCCGTACCGAGG\n->TESTERR215189.6OK\n-ATTTTGTTAAGGCCCGACACTGGTAGTGTTCGTGATCAAGTGTTTGAAAGTACTAATCTCATACCTAGTATGGGATGGGGAAGCCTAGTACCTGA\n->TESTERR215189.6OK\n-CTAATCTCACGACATCCGCGCACCACAAGGGTTGCTTCCTGTGTCGGCCGTCCCCATCAATTCCCAGGCACGTGTCAGGTCCAAACTTCCCTTGG\n->TESTERR215189.7OK\n-GCCCACCGACGGGTGGTCCTCTGGCCAATGGGCCCCGCGACGCAGCAGCTCTGTACGTGGCGGCCCCGGGTGCGCGGAGAGGGGCTCATAAATTT\n->TESTERR215189.7OK\n-GTCAGAAGACGCACGGACGCAAAGGTGCGAGCGCGAGTGGCACAAACACGTTTCCCTCCGTCCTCCCGCGCGCCCCACCGGAAAAGTCCAGACCC\n->TESTERR215189.8OK\n-AGTGAATTGAAGGCCCAGATGAAATCTGGGTTATTTCTTAATGAACATGAGGTTTGTGCTACTCAAGTTCTAATTCATTTGCTTAAAGCGTTATT\n->TESTERR215189.8OK\n-TGCTTTTTGTCATGTGTGGATACACTAAATCCAGATATGTACAGCTAAGACATGTTAAAATAAAAAATTTACATTAGTTTATGATTTATTAAATT\n->TESTERR215189.9OK\n-CTAACTTCCAGCACTTCACCACGTGTATCATATGTGTGCCGAGTTCCGCATCACGAGTATGCGAGGGTTGATAGGGAGTTCTATTCAAAGGACCC\n->TESTERR215189.9OK\n-CCCTTCGTGCCAGCACTCACGAAAGCGTCCCCAAGCATTACCGTTCACTATAGAAACGACACTCATGCATTTAACGCTGCATGGACAGCACAACA\n->TESTERR215189.10OK\n-GATCGTGTTACTAACATGTAGTGCCTAAAAAATTACAAAATGTAGGAAAGGCATAATCATTGGTGGATATTGCGCCCGGGTATATTCGGAATCAT\n->TESTERR215189.10OK\n-CGAAACATTGTAAGAAAGAGTTAATAGCGACAAATATCGAACCCGAGGGAGTACGTTACTACTGAAATGCAACGACGAATACACTTGAGCCAAAA\n->TESTERR215189.11OK\n-AGCACCGAGCCGGTCGCCACAATGCGTGGCTGATTGCCGCTAAGCCAACATGCAGTTGCAAACGCACTGCCGTACTCTAACCAGTCAACACTCCA\n->TESTERR215189.11OK\n-TCTCCGGGTGAAAACCATGCTTGAGTCCTCAGGCGGGCGATGGCGGCGCTTTAGCATCGTTTCTTTGGAGATCATGTTCCGGCCTTTTTCTCGGC\n->TESTERR215189.12OK\n-GGCTCCTCCGAGCCCATCGCGACGCGAAGGTAGAAGTCTCCGACCACGCGTTGCGTGCCGTAGCGGATCCACGCGTCGATCTGGTCGGGACGGGG\n->TESTERR215189.12OK\n-ATTCATCGGCGGCTTCCTGGACTGGGTGCTCTCCCAGCGCGGCGACGTCGACATGGACTCTCTCCTCATCTCCATGAAGACAGGGGACTGCCCCC\n->TESTERR215189.13OK\n-CATCGCCTACGCCGCCCAATAAAGGAGGATGCGCCGCCCGGTCTCCCGCTCGGCCTAGCGGTAGGGGCGTGTAGCCGCGCCTCCTGGTCAGCACG\n->TESTERR215189.13OK\n-GAAGGACCGAGCCCCTGGGGTCGCCGCCTTCCCTCCCGGTCTTTGCCTCCTGTCAGCATGACCGCTGAAGGGGGATGGTTGTCGTTGGGCCCCGT\n->TESTERR215189.14OK\n-GACTTGATTGATGGTAATAAGCTATCAGTGCTTAGATTCATTACAGCTAGAAGTAGAAATTGCTTAATTTGCTGCTTACGTTATCTGATTTGTGT\n->TESTERR215189.14OK\n-GAATAATTCCACAAGAATAGTACATGCAATTAGAGAAAGGGTTTTAATTGAAGCCTAATAAACACAGAAACAAACCAAAGCATGCAATTGAACAT\n->TESTERR215189.15OK\n-AATTTGCCTTGGAAATAACGTGAAGAAGAAAATGCACTTTCAAAGCTAGAAAACGCAAAACAGCAAGATACGCCTTCGTCAGACGTATACGTCTC\n->TESTERR215189.15OK\n-GTCCTCGCTTGCAGGCGTGCACTTGTGGCCACCGATGCAACCCAATATCCTTCTGCTCGGCGTCCGTCAGGTATCTCTGAGAGGATCCATGAACC\n->TESTERR215189.16OK\n-TCAGTCCTTGATGCAGGTCAACATCTCCACCGTCTCCGGCAACAACCGCTAGCGCCACTCCTCGATCACCCTGCAGTAAGACTAAAGCAAGTTTC\n->TESTERR215189.16OK\n-CTTTGGTGCTCCATCTAGTTCAGGTGATTCTCCTATTCTTGCATCTACACCTCCTATATCTCCTGCTTCTGAATTGTCATCCTACTTGGACAATG\n->TESTERR215189.17OK\n-CGTGTCCCCGTTGTACAGAGCAGCCGCCGCCGCGTGTGTTTACTGTTTAGTGCGATGGGTGTGTAGTGTTAGGATCATCTCTTCTGTTTGTACTG\n->TESTERR215189.17OK\n-CTGGCACATGAACCTAGCTGTCAGATCCACGTCGTAGCACAACCAGCGTGCCACACAAGCGGGTTAACATAAGTTGTGAGGAACTAGAACTTGCA\n->TESTERR215189.18OK\n-CTTCTGGTACTCGTCCGCTGACGTCATTGAAGAGCTGCCGAATCGGTGAAGCGTCTTCTACTAATCGGCCGATGTCCCCTTGAAGGAGGGATCGG\n->TESTERR215189.18OK\n-CCAGCATCCGATCCCTCCTTCAAGGGGACATCGGCCGATTAGTAGAAGACGCTTCACCGATTCGGCAGCTCTTCAATG'..b'TAGGA\n->TESTERR215189.74979OK\n-TTTTGGAGGCCCTTATCCTCATTACGCCTAGCATTGAGTGGGCTGGATATTTACCTTATCAACTAGCAAATCAATAGGGGTTCTATTTGATTAAG\n->TESTERR215189.74981OK\n-CATTTGCAGACCCAACTGGACCATCTGGACAAACAAATGTGAATTGTATGGTTCTTGTTCTTGTAGTAAAATAGTTTCCTTTTGTTGTACACAAC\n->TESTERR215189.74981OK\n-CTTACACCACTTCCTCCAAAGCATCTGGTAAAAACAAATGAGCAGACAAGCAACCAAACACATTGCCTACCCAGCTACATGATCGTGCCAGTGCT\n->TESTERR215189.74982OK\n-CGGCAACGTCGTCGAGGCCGTGCGCCACATCCGCTCCGTCATGGGGGACGTCCGCGCGCTGCGCAACATGGACGACGACGAGGTCTTCGCCTACG\n->TESTERR215189.74982OK\n-CACGATGGCTCGCGCCCGCCGCGCCGGGTCGCCGCTCTTGAAGATGCCGGAGCCGACGAAGACGCCGTCGCATCCGAGCTGCATCATGAGCGCGG\n->TESTERR215189.74984OK\n-ATTGGTTAACTTATCTTTCAAGGTGCTGAACGCAACCTCTTGTGAATCACTCCAAGCAAAGGGAACATCCTTCTTTGTAAGCTCATGTAGAGGCG\n->TESTERR215189.74984OK\n-AAATTCGAAGCTTTCTTGGCCTTGCAGGTTTCTACCGTCGGTTTGTTCGTGATTTTAGCTCCATTGCAGCGCCTCTACATGAGCTTACAAAGAAG\n->TESTERR215189.74986OK\n-CAACGTTGGTAGCAACAGGCCGCTACCAGGAGGACCGGCTGAGATGTGGTGTCCAAGGGCACAGGTTTTTGCTCCACTGTCGATCTTTGCATGAA\n->TESTERR215189.74986OK\n-CCCTTTTCCTAGAACTTGTGTGCTGCCGAATCTTGACAACACTATCTGTAGAAAGTTTGTTTTCATGCAAAGATCGACAGTGGAGCAAAAACCTG\n->TESTERR215189.74987OK\n-GCCAGGCTTTGTTGATCCATATGATCTGTACGATATGGAAGAAGCGGACGAGGATGAAATAGAGGAGCTGAAGAAATATGAGTTCGCTTTGCTGC\n->TESTERR215189.74987OK\n-TTCAGCATACTTGCGACAGCACGGAATGCAGCAAAACGAGGGTCGTCGCAGAGATAGTCCTCCAATGCACGATGCCTACCCTCGAGTCCAGGGTG\n->TESTERR215189.74989OK\n-ATGACTTGGCCTCATCCTCTCCTTCCTCCGGCTTAACACCGGCGGTCTGTTCAGGGTTCCAAACTCATAGTGGCAACTAAACACGAGGGTTGCGC\n->TESTERR215189.74989OK\n-GAGGGGTGCCCTCGGGAACGCGGACACAGGTGGTGCATGGCTGTCGTCAGCTCGTGCCGTAAGGTGTTGGGTTAAGTCTCGCAACGAGCGCAACC\n->TESTERR215189.74990OK\n-AAAGATCGTCGTACGATGTCCGATAATTTTATACGCGTAGCGGCGTAGGCCGTGTTCACGGAACGTTCCGTCACATCCCCCGTATCAACCTGAGC\n->TESTERR215189.74990OK\n-AGCGGCCGCGGGCACGGACAATCCGGGACAGCAGGGCCTCGCCGGCATTATTACCCCCGCCGTGCCGCCGCTTAATTGCCACCGTGGAGCCGCGG\n->TESTERR215189.74991OK\n-GGAGAAAGAGAATAATTGTTGGCACAAGTCGATAAGGTCAGATAGTCCTTACGTTGTCTAACGTATGGTGCAGGATGTCATGACTAGTCCAACTT\n->TESTERR215189.74991OK\n-GTTGTTCCGACAGTGGCGGTCGTGAACATCCATTACTAGTCTTTAACGGCCTATTCTAAAATGCCTTTTCAGTCTTAATTCCCAACCTTTGGAGG\n->TESTERR215189.74992OK\n-GCTGGTTAAGCAGAAGCTTTGCCGTTTTGCTTAATATTGCAAGGAGGCCATTAGCATATAAATAAATAAGCTTTTAGCAGCAGGTTTCATCCGTG\n->TESTERR215189.74992OK\n-TTTTTTACAAGGACTGGATTTGCTAACCAGTCGGGGTGCTTACGTTCACGGATGAAACCTGCTGCTAAAAGCTTATTTATTTATATGCTAATGGC\n->TESTERR215189.74993OK\n-GGGAGAGCAGAGCTGGGCTGCTGCCACACAGCTGTTTGCTGCTGATGGCAAAAGCGGTGTCACCGTGCCACGTGCCGGCGCAGGTGTGCGACGTG\n->TESTERR215189.74993OK\n-CCTTTCCCGGATATGATATTTAGTAGAGTGAGTTGACAGAGCCAGGAAGAGCTAAGAGAAAAAGCACAGTTAGGGCGATCTAGTTTCTTTTTCTT\n->TESTERR215189.74994OK\n-CTTTGGTTTTCAGCCCATCTCTCCCCATCACTCTACAATGCTTCTGCTTTAGCGTGTCTTTCATTGAACCACGCAACAAACCTGTAGAAGGTGAA\n->TESTERR215189.74994OK\n-GGTGTGAAGGATAATTTTGAACTGTTTCAGGTTCCTGCTCGTCAGCTTGAGGACTATTACTTTCAGCACACGTTGGAGCAGGTCAGGTCCACAAC\n->TESTERR215189.74996OK\n-TTGTTCAGCCACACCATAAGCCACCATGATAAGCTTGTTTTTCACGGGAGTAACAACCTCACAAGGAACCCTCCCAGTAATGTTGTCCACGTGGT\n->TESTERR215189.74996OK\n-TGCAAGAACACGTAGCTTTGACAGAGGCAAGAATGGAGGAAGCAATCGACCGGCGTGTGGCTCTAGCTCTTAGCAAACAAGCCATCAGGGGCTAA\n->TESTERR215189.74997OK\n-TGGCGAGGAGGCGCTCGCGCGTCTGGCCGACATGGGCCGGCTAGGCGTCCCCGTCGACAGCTACAGCTACGCGCACGGGCTTAAGGCCTGCATTG\n->TESTERR215189.74997OK\n-CAGGGCGCTCATTCTTGGCGTAGCATCCAATCATGGCACTCCATGACACAACATTCCTCTCGGGCATCCAAGCAAACACACGCTCTGCATAGGTG\n->TESTERR215189.74998OK\n-GTGTTCCAGGCGGCGTAAGTTGTGGAACAATTCCGCAACTCTTTAGATGATTGCTAAACTCGTGGCTCAAATACTCGCCTCCACGATCAGATCGT\n->TESTERR215189.74998OK\n-GTCTGAAACCTTTGAAAAGTTCAAGGAATTTCAGAATGAAGTTGAAAATCAGCGTGGCAAGAAAATTAAGGCCTTACGATCTGATCGTGGAGGCG\n->TESTERR215189.74999OK\n-CTTATCCGACAAAGCAAACATGCGCCTAGCGAACTGGGCTGGCAAGCTCCTCACGCCCGCAGGAAGGAGAACCCTAGTCCGGGCAGTCCTATCCT\n->TESTERR215189.74999OK\n-AGAAGCGAATAGGTTCTTGTCTTTTTCAGTGCACGGAGTTGGCATTTCGACCCAAGGTTTGCGTGGCTCCTTCCATTCATTCCAAAGCCAGCGTA\n->TESTERR215189.75000OK\n-CCACTCCTCCAAGCTATATGTGTGGGGAGAGAAGGAAAGCTCCACACACGCGCACTCGCTCACCACGCCACGCCTTGCCTCGCCGGGCCGGGCGG\n->TESTERR215189.75000OK\n-GGCAAAAACCCTAGCTTGATTAGATCATCTTTTCGGTAGACGCCAATAGGTCGTATATCTAGAAGAACCCATGATCTCACGCCACGATCGTGATC\n' |
| b |
| diff -r e320ef2d105a -r d397f5a85464 test_data/test_run1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test_data/test_run1 Wed Sep 18 06:30:04 2019 -0400 |
| b |
| @@ -0,0 +1,202 @@ +"Cluster" "Chip_Hits" "Input_Hits" "Ratio" +0 7945058 1071746 7.41319118522486 +1 278875 215538 1.2938553758502 +2 43504 94919 0.458327626713303 +3 37000 82538 0.448278368751363 +4 37253 104022 0.358126165618811 +5 32188 94960 0.338963774220724 +6 41039 81447 0.503873684727491 +7 39166 81623 0.479840241108511 +8 19339 66583 0.290449514140246 +9 29697 61497 0.482901604956339 +10 34273 82270 0.416591710222438 +11 27681 96899 0.285668582751112 +12 40545 96101 0.421899876171944 +13 31935 83895 0.380654389415341 +14 33932 90760 0.373865138827677 +15 29663 63240 0.469054395951929 +16 29084 91679 0.317237317160964 +17 17457 55788 0.312916756291676 +18 27249 55273 0.49298934380258 +19 21490 77952 0.275682471264368 +20 25851 67348 0.383842133396686 +21 27818 76020 0.365930018416206 +22 23859 72563 0.328803935890192 +23 85592 65619 1.30437830506408 +24 20936 59046 0.354571012430986 +25 21352 79958 0.267040196102954 +26 28586 62558 0.456951948591707 +27 24346 56413 0.4315671919593 +28 31336 64221 0.487940081904673 +29 16956 47383 0.357849861764768 +30 25372 59508 0.426362841970827 +31 19306 74537 0.259012302614809 +32 12498 50658 0.24671325358285 +33 22528 57473 0.391975362344057 +34 11578 38760 0.298710010319917 +35 17775 53446 0.332578677543689 +36 20979 71877 0.291873617429776 +37 24357 43553 0.559249649851905 +38 19401 65495 0.296221085579052 +39 23919 53615 0.446125151543411 +40 17944 37992 0.472309959991577 +41 21754 49360 0.440721231766613 +42 19911 59247 0.336067648994886 +43 17630 44057 0.400163424654425 +44 14639 47108 0.3107540120574 +45 10675 40313 0.26480291717312 +46 40498 55429 0.730628371430118 +47 15591 32845 0.474684122393058 +48 13229 45448 0.291079915507833 +49 15226 41012 0.371257193016678 +50 23851 49910 0.477880184331797 +51 11010 42552 0.258742244782854 +52 14233 41167 0.345738091189545 +53 17944 35795 0.5012990641151 +54 20288 38659 0.524793709097493 +55 13152 38394 0.342553523988123 +56 16631 38914 0.427378321426736 +57 249976 205389 1.21708562776001 +58 13347 28879 0.462169742719623 +59 16808 65546 0.256430598358405 +60 13355 42455 0.314568366505712 +61 16294 29082 0.560277835086995 +62 9790 52857 0.185216716801937 +63 17801 47726 0.372983279554121 +64 17157 32211 0.532644127782435 +65 20927 47539 0.440206987946738 +66 14999 49875 0.300731829573935 +67 16260 36851 0.441236330086022 +68 15925 42157 0.377754584054843 +69 11127 28004 0.397336094843594 +70 7739 25074 0.308646406636356 +71 6621 29947 0.221090593381641 +72 10923 40365 0.270605722779636 +73 9914675 84502 117.330654895742 +74 8758 38412 0.228001666145996 +75 9310 24364 0.382121162370711 +76 10614 47205 0.224849062599301 +77 13063 32597 0.400742399607326 +78 13404 29084 0.460871957089809 +79 9299 25667 0.362293996181868 +80 10109 26945 0.375171645945444 +81 15658 31730 0.493476205483769 +82 7483 32182 0.232521285190479 +83 16345 30634 0.533557485147222 +84 13293 27484 0.483663222238393 +85 5322 38749 0.137345479883352 +86 21424 33627 0.637107086567342 +87 11548 29314 0.393941461417753 +88 7762 29456 0.263511678435633 +89 12153 46433 0.261731957874787 +90 13081 22936 0.570326124869201 +91 6183 33349 0.185402860655492 +92 6357 35098 0.181121431420594 +93 48995 31410 1.5598535498249 +94 8371 31091 0.269241902801454 +95 13228 20080 0.658764940239044 +96 8453 27056 0.312426079243051 +97 9568 19505 0.490540886952064 +98 7938 32634 0.243243243243243 +99 7824 16977 0.460858808976851 +100 76514 28367 2.6972891035358 +101 10176 25158 0.404483663248271 +102 11009 30853 0.356821054678637 +103 6696 20379 0.328573531576623 +104 6767 26109 0.25918265732123 +105 12572 66266 0.189720218513265 +106 9522 20128 0.473072337042925 +107 10249 29304 0.349747474747475 +108 10362 34525 0.300130340333092 +109 9815 35803 0.274139038628048 +110 4387 16229 0.27031856553084 +111 10448 17881 0.584307365359879 +112 8117443 57680 140.732368238558 +113 9645 19854 0.485796313085524 +114 6879 27269 0.252264476145073 +115 8197 18651 0.439493860918986 +116 7111 21775 0.326567164179104 +117 16273 33720 0.482591933570581 +118 10242 29605 0.345955075156224 +119 12252 37892 0.323340019001372 +120 10751 27590 0.389670170351577 +121 6448 18451 0.349466153596011 +122 9191 25908 0.354755287941948 +123 5235 24699 0.211951900886676 +124 6465 20952 0.308562428407789 +125 7721 17814 0.4334231503312 +126 6593 26577 0.248071640892501 +127 12142 18977 0.639827159192707 +128 6505 18595 0.349825221833826 +129 3710 21375 0.173567251461988 +130 6348 16298 0.389495643637256 +131 4208 19672 0.213908092720618 +132 4136 9478 0.436378982907786 +133 7285 17145 0.424905220180811 +134 7642 14034 0.544534701439362 +135 3953 13049 0.302935090811556 +136 7778 14820 0.524831309041835 +137 7588 20607 0.368224389770466 +138 1208 3015 0.400663349917081 +139 7826 20711 0.377866834049539 +140 7551 13678 0.552054393917239 +141 4938 15140 0.326155878467635 +142 3788 16800 0.22547619047619 +143 4537 15722 0.288576516982572 +144 20198 10665 1.8938584153774 +145 28177 18344 1.53603358046228 +146 5743 19920 0.288303212851406 +147 3763 14154 0.265861240638689 +148 56754 17007 3.33709648968072 +149 5988 11740 0.510051107325383 +150 784 1172 0.668941979522184 +151 2718 13509 0.201199200532978 +152 3807 14812 0.25702133405347 +153 8382 45786 0.183069060411479 +154 3647 15198 0.239965784971707 +155 8500 32566 0.261008413682982 +156 6815 14389 0.473625686288137 +157 2307 13795 0.167234505255527 +158 2912 11803 0.246716936372109 +159 4316 11848 0.364280891289669 +160 992 2170 0.457142857142857 +161 5305 13186 0.402320643106325 +162 4569 17654 0.258808202107171 +163 3935 16460 0.23906439854192 +164 2725 13244 0.205753548776805 +165 3772 11890 0.317241379310345 +166 3081 15437 0.199585411673253 +167 837 2591 0.323041296796604 +168 2438412 16297 149.62336626373 +169 6131 16203 0.378386718508918 +170 3283 13179 0.249108430078155 +171 3106 12624 0.246039290240811 +172 9370 45868 0.204281852271736 +173 3537 15712 0.225114562118126 +174 747 2267 0.329510366122629 +175 4326 13004 0.332666871731775 +176 6665 9325 0.714745308310992 +177 4657 12573 0.370396882207906 +178 2538 10016 0.253394568690096 +179 2662 9086 0.292978208232446 +180 2441 10297 0.237059337671166 +181 544 1658 0.32810615199035 +182 664 1018 0.652259332023576 +183 31045 22110 1.4041157847128 +184 1365 4387 0.311146569409619 +185 2159 11375 0.189802197802198 +186 1933 7834 0.246744957875925 +187 2062 6753 0.305345772249371 +188 3962 11049 0.358584487283917 +189 2234 9201 0.242799695685252 +190 834 6610 0.126172465960666 +191 520 1629 0.319214241866176 +192 2677 10891 0.245799283812322 +193 509 805 0.632298136645963 +194 2155 4528 0.475927561837456 +195 2477 7482 0.331061213579257 +196 14445 15443 0.935375250922748 +197 3703 8763 0.42257217847769 +198 2873 28533 0.100690428626503 +199 1727 9618 0.179559159908505 +200 3448 6234 0.553095925569458 |