Mercurial > repos > tyty > structurefold

changeset 4:a292aaf51735 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author tyty
date Mon, 15 Sep 2014 14:52:43 -0400
parents f4cc06e92530
children 7a8ddf1819b1
files Iterative_mapping/iterative_map.py Iterative_mapping/iterative_map.xml Iterative_mapping/map_ex.py Iterative_mapping/read_file.py Iterative_mapping/read_s_file.py Iterative_mapping/remove_map.py Iterative_mapping/seq_track.py Iterative_mapping/tool_dependencies.xml Iterative_mapping/truncate.py Iterative_mapping/unmap.py predict/.DS_Store predict/parse_dis_pac.py predict/predict_RNAs.py predict/predict_RNAs.xml predict/rRNA.txt predict/rRNA_react.txt predict/read_file.py predict/rtts_plot.py predict/tool_dependencies.xml
diffstat 19 files changed, 301 insertions(+), 402 deletions(-) [+]
line wrap: on
line diff
--- a/Iterative_mapping/iterative_map.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,103 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-import os
-from read_file import *
-from read_s_file import *
-
-type_input = sys.argv[1]
-seq_file = sys.argv[2]
-ref_file = sys.argv[3]
-shift = sys.argv[4]
-length = sys.argv[5]
-map_type = sys.argv[6]
-output_file = sys.argv[7]
-
-
-if map_type!="default":
-    s = ""
-    s = s+"-v "+sys.argv[8]
-    s = s+" -5 "+sys.argv[9]
-    s = s+" -3 "+sys.argv[10]
-    s = s+" -k "+sys.argv[11]
-    if sys.argv[12]:
-        s = s+" -a"
-    if int(sys.argv[13])>=1:
-        s = s+" -m "+sys.argv[13]
-    if sys.argv[14]:
-        s = s+" --best --strata "
-    
-else:
-    s = "-v 3 -a --best --strata "
-
-ospath = os.path.realpath(sys.argv[0])
-ost = ospath.split('/')
-syspath = ""
-for i in range(len(ost)-1):
-    syspath = syspath+ost[i].strip()
-    syspath = syspath+'/'
-
-os.system("bowtie-build -f "+ref_file+" "+syspath+"ref > "+syspath+"log.txt")
-
-os.system("cp "+seq_file+" "+syspath+"seq0.fa")
-
-if type_input == "fasta":
-    tp = 'fasta'
-if type_input == "fastq":
-    tp = 'fastq'
-
-k = 0
-print(type_input)
-while(True):
-    if type_input == "fasta":
-        os.system("bowtie "+s+"-f "+syspath+"ref"+" "+syspath+"seq"+str(k)+".fa --quiet -S > "+syspath+"map"+str(k)+".sam")
-    if type_input == "fastq":
-        os.system("bowtie "+s+"-q "+syspath+"ref"+" "+syspath+"seq"+str(k)+".fa --quiet -S > "+syspath+"map"+str(k)+".sam")
-    os.system("samtools view -Sb -F 0xfff "+syspath+"map"+str(k)+".sam > "+syspath+"mapped"+str(k)+".bam 2>"+syspath+"log.txt") #get mapped reads
-    os.system("samtools view -Sb -f 0x4 "+syspath+"map"+str(k)+".sam > "+syspath+"umapped"+str(k)+".bam 2>"+syspath+"log.txt") #get unmapped reads
-    os.system("samtools view -Sb -f 0x10 "+syspath+"map"+str(k)+".sam > "+syspath+"rmapped"+str(k)+".bam 2>"+syspath+"log.txt") #get reversed mapped reads
-    os.system("samtools merge -f "+syspath+"unmapped"+str(k)+".bam "+syspath+"umapped"+str(k)+".bam "+syspath+"rmapped"+str(k)+".bam") #get reversed mapped reads
-    os.system("samtools view -h -o "+syspath+"unmapped"+str(k)+".sam "+syspath+"unmapped"+str(k)+".bam") #get reversed mapped reads
-    if k>0:
-        os.system("samtools view -h -o "+syspath+"mapped"+str(k)+".sam "+syspath+"mapped"+str(k)+".bam") #get reversed mapped reads
-        os.system("cut -f 1 "+syspath+"unmapped"+str(k)+".sam > "+syspath+"unmapped"+str(k)+".txt")
-        os.system("cut -f 1 "+syspath+"mapped"+str(k)+".sam > "+syspath+"mapped"+str(k)+".txt")
-        os.system("python "+syspath+"remove_map.py "+syspath+"unmapped"+str(k)+".txt "+syspath+"mapped"+str(k)+".txt "+syspath+"runmapped"+str(k)+".txt")
-        os.system("rm "+syspath+"mapped"+str(k)+".sam")
-        os.system("rm "+syspath+"mapped"+str(k)+".txt")
-        os.system("rm "+syspath+"unmapped"+str(k)+".txt")
-    else:
-        os.system("cut -f 1 "+syspath+"unmapped"+str(k)+".sam > "+syspath+"runmapped"+str(k)+".txt")
-    
-    os.system("rm "+syspath+"unmapped"+str(k)+".bam")
-    os.system("rm "+syspath+"umapped"+str(k)+".bam")
-    os.system("rm "+syspath+"rmapped"+str(k)+".bam")
-    os.system("python "+syspath+"seq_track.py "+syspath+"runmapped"+str(k)+".txt "+syspath+"seq"+str(k)+".fa "+syspath+"unmap_seq"+str(k)+".fa "+tp) #get unmapped sequence
-    os.system("python "+syspath+"truncate.py "+syspath+"unmap_seq"+str(k)+".fa "+shift+" "+syspath+"seq"+str(k+1)+".fa "+length) #truncate unmapped sequence
-    os.system("rm "+syspath+"seq"+str(k)+".fa") #Remove sequences being mapped
-    os.system("rm "+syspath+"map"+str(k)+".sam") #Remove mapping file
-    os.system("rm "+syspath+"unmap_seq"+str(k)+".fa") #Remove unmapped sequnce
-    os.system("rm "+syspath+"runmapped"+str(k)+".txt")
-    os.system("rm "+syspath+"unmapped"+str(k)+".sam")
-    
-    os.system("wc -l "+syspath+"seq"+str(k+1)+".fa > "+syspath+"count"+str(k+1)+".txt")
-    c = read_sp_file(syspath+"count"+str(k+1)+".txt")
-    if c[0][0] == '0': #If no reads is in the sequence file, stop
-        os.system("rm "+syspath+"count"+str(k+1)+".txt")
-        os.system("rm "+syspath+"seq"+str(k+1)+".fa")
-        break
-    os.system("rm "+syspath+"count"+str(k+1)+".txt")
-    k = k+1
-
-ss = ""
-for i in range(0,k+1):
-    ss = ss+" "+syspath+"mapped"+str(i)+".bam"
-
-
-os.system("samtools merge -f "+output_file+" "+ss)
-#print("samtools merge mapped_all.bam"+ss)
-os.system("rm "+syspath+"mapped*.bam")
-os.system("rm "+syspath+"ref*")
-
-
--- a/Iterative_mapping/iterative_map.xml	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,80 +0,0 @@
-<tool id="iterative_map_pipeline" name="Iterative mapping" version="1.0">
-	<description></description>
-	<command interpreter="python">
-        #if $mapping_file.type == "user"
-            iterative_map.py $file_format.type $file_format.seq_file $reference_file $shift $length $mapping_file.type $output $mapping_file.param_v $mapping_file.param_five $mapping_file.param_three $mapping_file.param_k $mapping_file.param_a $mapping_file.param_m $mapping_file.param_best  
-        #else
-            iterative_map.py $file_format.type $file_format.seq_file $reference_file $shift $length $mapping_file.type $output
-        #end if
-    </command>
-        <requirements>
-                <requirement type="package" version="1.61">biopython</requirement>
-                <requirement type="package" version="1.7">numpy</requirement>
-                <requirement type="package" version="0.1.18">samtools</requirement>
-                <requirement type="package" version="0.12.7">bowtie</requirement>
-        </requirements>
-	<inputs>
-                <conditional name="file_format">
-                  <param name="type" type="select" label="Format of the file of the reads (Default FASTQ)">
-                    <option value="fastq">FASTQ</option>
-                    <option value="fasta">FASTA</option>
-                  </param>
-                  <when value="fastq">
-                    <param name="seq_file" type="data" format="fastq" label="Fastq file"/>
-                  </when>
-                  <when value="fasta">
-                    <param name="seq_file" type="data" format="fasta" label="Fasta file"/>
-                  </when>
-                </conditional>
-		<param name="reference_file" type="data" format="fasta" label="Reference genome/transcriptome"/>
-                <param name="shift" type="integer" value="1" label="Number of nucleotide trimmed each round"/>
-                <param name="length" type="integer" value="21" label="Minimum requirement of read length for mapping"/>
-                <conditional name="mapping_file">
-                  <param name="type" type="select" label="Bowtie mapping flags (Default -v 0 -a --best --strata)">
-                    <option value="default">Default</option>
-                    <option value="user">User specified</option>
-                  </param>
-                  <when value="default"/>
-                  <when value="user"> 
-                    <param name="param_v" type="integer" value="0" label="Number of mismatches for SOAP-like alignment policy (-v)"/>
-                    <param name="param_five" type="integer" value="0" label="Trim n bases from high-quality (left) end of each read before alignment (-5)"/>
-                    <param name="param_three" type="integer" value="0" label="Trim n bases from high-quality (right) end of each read before alignment (-3)"/>
-                    <param name="param_k" type="integer" value="1" label="Report up to n valid alignments per read (-k)"/>
-                    <param name="param_a" type="boolean" checked="False" truevalue = "1" falsevalue = "0" label="Whether or not to report all valid alignments per read (-a)"/>
-                    <param name="param_m" type="integer" value="-1" label="Suppress all alignments for a read if more than n reportable alignments exist (-m), -1 for unlimited"/>
-                    <param name="param_best" type="boolean" checked="False" truevalue = "1" falsevalue = "0" label="Whether or not to make Bowtie guarantee that reported singleton alignments are 'best' in terms of stratum and in terms of the quality values at the mismatched positions (--best --strata)"/>
-                  </when>
-                </conditional>
-
-	</inputs>
-	<outputs>
-		<data name="output" type="data" format="bam"/>
-	</outputs>
-
-	<help>
-
-
-**TIPS**:
-
------
-
-**Input**:
-
-* 1. Sequence file type (FASTA/FASTQ)
-* 2. Sequence file (fasta/fastq format) {Default: fastq file}
-* 3. Reference file (e.g. cDNA library [fasta])
-* 4. “Shift” (The length of the sequence that will be trimmed at the 3’end of the reads before each round of mapping)
-* 5. “Length” (The minimum length of the reads for mapping after trimming)
-* [Optional]
-* 1. Bowtie mapping flags (options) [Default: -v 0 -a --best --strata] (-v flag indicates the number of allowed mismatches. use -5/-3 flag to trim nucleotides from 5'/3' end of the reads)
-
------
-
-**Output**:
-
-A bam file with all of the reads that are mapped	
-
-
-
-	</help>
-</tool>
--- a/Iterative_mapping/map_ex.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,31 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-from read_file import *
-from Bio import SeqIO
-
-map_file = sys.argv[1]
-result_file = sys.argv[2]
-
-
-#reads = read_t_file(read_file);
-
-f = open(map_file);
-h = file(result_file, 'w')
-
-for aline in f.readlines():
-    tline = aline.strip();
-    tl = tline.split('\t');
-    if len(tl)>4:
-        if int(tl[1].strip())== 0:
-            h.write(tline)
-            h.write('\n')
-
-
-f.close();
-h.close()
-
-
-
-
--- a/Iterative_mapping/read_file.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,21 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-
-
-
-def read_t_file(in_file):
-    f = open(in_file);
-    result = [];
-    for aline in f.readlines():
-        temp = [];
-        tline = aline.strip();
-        tl = tline.split('\t');
-        for i in range(0, len(tl)):
-            temp.append(tl[i].strip());
-        result.append(temp);
-    f.close();
-    return result;
-
-
--- a/Iterative_mapping/read_s_file.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,22 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-
-
-
-def read_sp_file(in_file):
-    f = open(in_file);
-    result = [];
-    for aline in f.readlines():
-        temp = [];
-        tline = aline.strip();
-        tl = tline.split(' ');
-        for i in range(0, len(tl)):
-            if len(tl[i].strip())>0:
-                temp.append(tl[i].strip());
-        result.append(temp);
-    f.close();
-    return result;
-
-
--- a/Iterative_mapping/remove_map.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,29 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-from read_file import *
-
-
-unmap_file = sys.argv[1]
-map_file = sys.argv[2]
-result_file = sys.argv[3]
-
-
-unmap = read_t_file(unmap_file)
-mapped = read_t_file(map_file)
-h = file(result_file, 'w')
-
-maps = set()
-for i in range(len(mapped)):
-    maps.add(mapped[i][0])
-
-
-for i in range(len(unmap)):
-    name = unmap[i][0]
-    if name not in maps:
-        h.write(name)
-        h.write('\n')
-
-
-h.close()
--- a/Iterative_mapping/seq_track.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,38 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-from read_file import *
-from Bio import SeqIO
-
-unmap_file = sys.argv[1]
-reads_file = sys.argv[2]
-result_file = sys.argv[3]
-tp = sys.argv[4]
-
-
-unmap = read_t_file(unmap_file);
-
-h = file(result_file, 'w')
-
-reads = SeqIO.parse(reads_file,tp)
-um = set()
-for i in range(0, len(unmap)):
-    id_r = unmap[i][0]
-    um.add(id_r)
-
-for read in reads:
-    if read.id in um:
-        h.write('>')
-        h.write(read.id)
-        h.write('\n')
-        h.write(read.seq.tostring())
-        h.write('\n')
-    
-
-
-h.close()
-
-
-
-
--- a/Iterative_mapping/tool_dependencies.xml	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,15 +0,0 @@
-<?xml version="1.0"?>
-<tool_dependency>
-    <package name="biopython" version="1.61">
-        <repository changeset_revision="ae9dda584395" name="package_biopython_1_61" owner="biopython" prior_installation_required="True" toolshed="http://toolshed.g2.bx.psu.edu" />
-    </package>
-    <package name="numpy" version="1.7">
-        <repository changeset_revision="ef12a3a11d5b" name="package_numpy_1_7" owner="iuc" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
-    </package>
-    <package name="samtools" version="0.1.18">
-        <repository changeset_revision="171cd8bc208d" name="package_samtools_0_1_18" owner="devteam" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
-    </package>
-    <package name="bowtie" version="0.12.7">
-        <repository changeset_revision="9f9f38617a98" name="package_bowtie_0_12_7" owner="devteam" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
-    </package>
-</tool_dependency>
--- a/Iterative_mapping/truncate.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,32 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-from Bio import SeqIO
-
-fasta_file = sys.argv[1]
-shift_in = sys.argv[2]
-result_file = sys.argv[3]
-length = sys.argv[4]
-
-shift = int(shift_in)
-    
-fasta_sequences = SeqIO.parse(open(fasta_file),'fasta');
-h = file(result_file,'w')
-for seq in fasta_sequences:
-        nuc = seq.id;
-        sequence = seq.seq.tostring();
-        if (len(sequence)-shift)>=int(length):
-                h.write('>'+nuc)
-                h.write('\n')
-                h.write(sequence[0:(len(sequence)-shift)])
-                h.write('\n')
-
-
-
-
-h.close()
-
-
-
-
--- a/Iterative_mapping/unmap.py	Mon Sep 15 14:52:20 2014 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,31 +0,0 @@
-#!/usr/bin/env python
-# -*- coding: utf-8 -*-
-
-import sys
-from read_file import *
-from Bio import SeqIO
-
-map_file = sys.argv[1]
-result_file = sys.argv[2]
-
-
-#reads = read_t_file(read_file);
-
-f = open(map_file);
-h = file(result_file, 'w')
-
-for aline in f.readlines():
-    tline = aline.strip();
-    tl = tline.split('\t');
-    if len(tl)>4:
-        if int(tl[1].strip()) != 0:
-            h.write(tl[0].strip());
-            h.write('\n');
-
-
-f.close();
-h.close()
-
-
-
-
Binary file predict/.DS_Store has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/parse_dis_pac.py	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,43 @@
+#parse reactivity file into a dictionary
+
+import sys
+
+def parse_dist(in_file):
+    result = []
+    distribution = {}
+    name = []
+    f = open(in_file)
+    for aline in f.readlines():
+        line = aline.strip()
+        dis = line.strip()
+        dist = dis.split('\t') #split the line and the reactivites or reads are in a list
+        if len(dist) > 0:
+            if len(dist) == 1:
+                if dist[0].strip().find('coverage')==-1:
+                    name.append(line) #add the name in the name list
+                    flag = 1
+                    t_name = line
+            else:
+                distri = []
+                for i in range(0, len(dist)):
+                    distri.append(dist[i].strip())
+                distribution[t_name] = distri #add the list of reactivities into a dictionary
+    result.append(name)
+    result.append(distribution) #Output the dictionary
+    f.close()
+    return result
+                
+                
+
+
+
+
+
+
+
+        
+
+
+
+
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/predict_RNAs.py	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,93 @@
+#RNA structure prediction & Output and illustrate reactivities
+
+import sys
+from parse_dis_pac import *
+from read_file import *
+from Bio import SeqIO
+import os
+from rtts_plot import *
+
+
+id_file = sys.argv[1]
+seq_file = sys.argv[2]
+output_file = sys.argv[4]
+
+
+flag = 0
+if sys.argv[3]!='None': #input reactivity file if provided
+    react_file = sys.argv[3]
+    react = parse_dist(react_file)
+    react = react[1]
+    flag = 1
+
+ospath = os.path.realpath(sys.argv[0])
+ost = ospath.split('/')
+syspath = ""
+for i in range(len(ost)-1):
+    syspath = syspath+ost[i].strip()
+    syspath = syspath+'/' 
+
+ids = read_t_file(id_file)
+sequences = SeqIO.parse(seq_file, 'fasta')
+
+
+seqs = {}
+for seq in sequences:
+    seqs[seq.id] = seq.seq.tostring()
+
+if len(ids)>10: #setup a limit of the number of sequence to be predicted
+    print("Number of sequences exceeds limitation!")
+    sys.exit(0)
+    
+
+#predict RNA structures
+os.system("mkdir "+syspath+"output_f")
+for i in range(len(ids)):
+    id_s = ids[i][0]
+    print(id_s)
+    #Put RNA sequence and reactivities into files
+    if id_s in seqs:
+        f = file(syspath+"temp.txt", 'w')        
+        f.write('>'+id_s)
+        f.write('\n')
+        f.write(seqs[id_s])
+        f.close()
+        if flag == 0:
+            os.system("Fold "+syspath+"temp.txt"+" "+syspath+"output_f/"+id_s+".ct")
+        if flag == 1:
+            if id_s in react:
+                f = file(syspath+"constraint.txt",'w')
+                make_plot(react[id_s],id_s,(syspath+"output_f/")) #make a plot of the distribution of the reactivites of the input RNA
+                #h = file(syspath+"output_f/transcript_reactivities.txt", 'w')
+                #h.write(id_s)
+                #h.write('\n')
+                for j in range(0, (len(react[id_s]))):
+                    if react[id_s][j]!='NA':
+                        f.write(str(j+1))
+                        f.write('\t')
+                        f.write(str(react[id_s][j]))
+                        f.write('\n')
+                    #h.write(str(react[id_s][j])) #Output the reactivities
+                    #h.write('\t')
+                f.close()
+                #h.write('\n')
+                #h.write('\n')
+                os.system("Fold "+syspath+"temp.txt"+" -sh"+" "+syspath+"constraint.txt"+" "+syspath+"output_f/"+id_s+".ct")
+            else:
+                print(id_s+" not in the data of react!")
+        os.system("draw "+syspath+"output_f/"+id_s+".ct "+syspath+"output_f/"+id_s+".ps")
+    else:
+        print(id_s+" not in the data of sequences!")
+
+#Remove the unnecessary files
+os.system("tar -zcvPf "+output_file+" "+syspath+"output_f/"+"*.*")
+os.system("rm -f "+syspath+"temp.txt")
+os.system("rm -r "+syspath+"output_f")
+if flag == 1:
+    os.system("rm -f "+syspath+"constraint.txt")
+ #   h.close()
+    
+    
+
+
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/predict_RNAs.xml	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,63 @@
+<tool id="predict_pipeline" name="RNA structure prediction" version="1.0">
+	<description></description>
+	<command interpreter="python">predict_RNAs.py $rna_list $reference_file $reactivity_file $output </command>
+        <requirements>
+                <requirement type="package" version="1.61">biopython</requirement>
+                <requirement type="package" version="1.7">numpy</requirement>
+        </requirements>
+	<inputs>
+        <param name="rna_list" type="data" format="txt" label="List of RNA ids to predict"/>
+        <param name="reference_file" type="data" format="fasta" label="Reference genome/transcriptome"/>
+        <param name="reactivity_file" type="data" optional = "true" label="Reactivity"/>
+<!---
+        <conditional name="reactivity">
+            <param name="choose" type="select" label="Reactivity files exists?">
+                <option value="N">No</option>
+                <option value="Y">Yes</option>
+            </param>
+            <when value="Y">
+                <param name="reactivity_file" type="data" optional = "true" label="Reactivity"/>
+            </when>
+            <when value="N"/>
+        </conditional>  
+
+-->		
+	</inputs>
+	<outputs>
+		<data name="output" format=".tgz"/>
+	</outputs>
+
+	<help>
+
+
+**TIPS**:
+
+-----
+
+**Input**:
+
+* 1. A file with transcript Ids (Max num. 20), (each ID one line)
+* 2. Reference file (fasta) used to map the reads
+* [Optional]:
+* 1. A reactivity file with structural reactivity for each nucleotide on the sequence provided
+
+-----
+
+**Output**:
+
+* 1. .ct files with predicted RNA structures [transciptID.ct]
+* 2. .ps files which depict the predicted RNA structures [[transciptID.ps]
+* [Optional]
+* 3. .png files that shows the distribution of the reactivity of each nucleotide on the transcripts of interest. [transciptID.png]
+* 4. A .txt file that includes the reactivities of all the nucleotides on the transcripts of interest. [transciptID.txt]
+
+-----
+
+**Attention**
+
+Make sure any of the transcript Ids does not contain "|" or space!	
+
+
+
+	</help>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/rRNA.txt	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,8 @@
+>25s rRNA 3375nts
+GCGACCCCAGGTCAGGCGGGATTACCCGCTGAGTTTAAGCATATCAATAAGCGGAGGAAAAGAAACTAACAAGGATTCCCTTAGTAACGGCGAGCGAACCGGGAAGAGCCCAGCTTGAAAATCGGACGTCTTCGGCGTTCGAATTGTAGTCTGGAGAAGCGTCCTCAGCGACGGACCGGGCCTAAGTTCCCTGGAAAGGGGCGCCAGAGAGGGTGAGAGCCCGTCGTGCCCGGACCCTGTCGCACCACGAGGCGCTGTCTACGAGTCGGGTTGTTTGGGAATGCAGCCCCAATCGGGCGGTAAATTCCGTCCAAGGCTAAATACGGGCGAGAGACCGATAGCGAACAAGTACCGCGAGGTAAAGATGAAAAGGACTTTGAAAAGAGAGTCAAAGAGTGCTTGAAATTGTCGGGAGGGAAGCGGATGGGGGCCGGCGATGCGTCCTGGTCGGATGCGGAACGGAGCAATCCGGTCCGCCGATCGATTCGGGGCGTGGACCGACGCGGATTACGGTGGCGGCCTAAGCCCGGGCTTTTGATACGCTTGTGGAGACGTCGCTGCCGTGATCGTGGTCTGCAGCACGCGCCTAACGGCGTGCCTCGGCATCAGCGTGCTCCGGGCGTCGGCCTGTGGGCTCCCCATTCGACCCGTCTTGAAACACGGACCAAGGAGTCTGACATGTGTGCGAGTCAACGGGTGAGTAAACCCGTAAGGCGCAAGGAAGCTGATTGGCGGGATCCTCGCGGGTGCACCGCCGACCGACCTTGATCTTCTGAGAAGGGTTCGAGTGTGAGCATGCCTGTCGGGACCCGAAAGATGGTGAACTATGCCTGAGCGGGGTAAAGCCAGAGGAAACTCTGGTGGAAGCCCGCAGCGATACTGACGTGCAAATCGTTCGTCTGACTTGGGTATAGGGGCGAAAGACTAATCGAACCATCTAGTAGCTGGTTCCCTCCGAAGTTTCCCTCAGGATAGCTGGAGCTCGGACGCGAGTTCTATCGGGTAAAGCCAATGATTAGAGGCATTGGGGGCGCAACGCCTCGACCTATTCTCAAACTTTAAATAGGTAGGACGTGTCGGCTGCTTTGTTGAGCCGTCACACGGAATCGAGAGCTCCAAGTGGGCCATTTTTGGTAAGCAGAACTGGCGATGCGGGATGAACCGGAAGCCGGGTTACGGTGCCCAACTGCGCGCTAACCTAGAACCCACAAAGGGTGTTGGTCGATTAAGACAGCAGGACGGTGGTCATGGAAGTCGAAATCCGCTAAGGAGTGTGTAACAACTCACCTGCCGAATCAACTAGCCCCGAAAATGGATGGCGCTTAAGCGCGACCTATACCCGGCCGTCGGGGCAAGAGCCAGGCCTCGATGAGTAGGAGGGCGCGGCGGTCGCTGCAAAACCTAGGGCGCGAGGCGCGGAGCGGCCGTCGGTGCAGATCTTGGTGGTAGTAGCAAATATTCAAATGAGAACTTTGAAGGCCGAAGAGGGGAAAGGTTCCATGTGAACGGCACTTGCACATGGGTTAGTCGATCCTAAGAGTCGGGGGAAACCCGTCTGATAGCGCTTAAGCGAACTTCGAAAGGGGATCCGGTTAAAATTCCGGAACCGGGACGTGGCGGTTGACGGCAACGTTAGGGAGTCCGGAGACGTCGGCGGGGGCCTCGGGAAGAGTTATCTTTTCTGTTTAACAGCCTGCCCACCCTGGAAACGGCTCAGCCGGAGGTAGGGTCCAGCGGCTGGAAGAGCACCGCACGTCGCGTGGTGTCCGGTGCGCCCCCGGGCGCCCTTGAAAATCCGGAGGACCGAGTGCCGCTCACGCCCGGTCGTACTCATAACCGCATCAGGTCTCCAAGGTGAACAGCCTCTGGTCGATGGAACAATGTAGGCAAGGGAAGTCGGCAAAATGGATCCGTAACTTCGGGAAAAGGATTGGCTCTGAGGGCTGGGCTCGGGGGTCCCAGTTCCGAACCCGTCGGCTGTCAGCGGACTGCTCGAGCTGCTTCCGCGGCGAGAGCGGGTCGCCGGCTGCCGGCCGGGGGACGACTGGGAACGGCTCTCTCGGGAGCTTTCCCCGGGCGTCGAACAGTCAGCTCAGAACTGGTACGGACAAGGGGAATCCGACTGTTTAATTAAAACAAAGCATTGCGATGGTCCCTGCGGATGCTAACGCAATGTGATTTCTGCCCAGTGCTCTGAATGTCAAAGTGAAGAAATTCAACCAAGCGCGGGTAAACGGCGGGAGTAACTATGACTCTCTTAAGGTAGCCAAATGCCTCGTCATCTAATTAGTGACGCGCATGAATGGATTAACGAGATTCCCACTGTCCCTGTCTACTATCCAGCGAAACCACAGCCAAGGGAACGGGCTTGGCAGAATCAGCGGGGAAAGAAGACCCTGTTGAGCTTGACTCTAGTCCGACTTTGTGAAATGACTTGAGAGGTGTAGGATAAGTGGGAGCTTCGGCGCAAGTGAAATACCACTACTTTTAACGTTATTTTACTTACTCCGTGAATCGGAGGCCGGGGTACAACCCCTGTTTTTGGTCCCAAGGCTCGCTTCGGCGGGTCGATCCGGGCGGAGGACATTGTCAGGTGGGGAGTTTGGCTGGGGCGGCACATCTGTTAAAAGATAACGCAGGTGTCCTAAGATGAGCTCAACGAGAACAGAAATCTCGTGTGGAACAAAAGGGTAAAAGCTCGTTTGATTCTGATTTTCAGTACGAATACGAACCGTGAAAGCGTGGCCTATCGATCCTTTAGACTTCGGAATTTGAAGCTAGAGGTGTCAGAAAAGTTACCACAGGGATAACTGGCTTGTGGCAGCCAAGCGTTCATAGCGACGTTGCTTTTTGATCCTTCGATGTCGGCTCTTCCTATCATTGTGAAGCAGAATTCACCAAGTGTTGGATTGTTCACCCACCAATAGGGAACGTGAGCTGGGTTTAGACCGTCGTGAGACAGGTTAGTTTTACCCTACTGATGCCCGCGTCGCGATAGTAATTCAACCTAGTACGAGAGGAACCGTTGATTCGCACAATTGGTCATCGCGCTTGGTTGAAAAGCCAGTGGCGCGAAGCTACCGTGCGCTGGATTATGACTGAACGCCTCTAAGTCAGAATCCGGGCTAGAAGCGACGCATGCGCCCGCCGCCCGATTGCCGACCCTCAGTAGGAGCTTAGGCTCCAAAGGCACGTGTCGTTGGCTAAGTCCGTTCGGCGGAACGGTCGTTCGGACCGCCTTGAATTATAATTACCACCGAGCGGCGGGTAGAATCCTTTGCAGACGACTTAAATACGCGACGGGGTATTGTAAGTGGCAGAGTGGCCTTGCTGCCACGATCCACTGAGATTCAGCCCTTTGTCGCTAAGATTCGA
+>gi|20197903:2706-4513 Arabidopsis thaliana chromosome 2 BAC F23H14 genomic sequence, complete sequence
+TACCTGGTTGATCCTGCCAGTAGTCATATGCTTGTCTCAAAGATTAAGCCATGCATGTGTAAGTATGAACGAATTCAGACTGTGAAACTGCGAATGGCTCATTAAATCAGTTATAGTTTGTTTGATGGTAACTACTACTCGGATAACCGTAGTAATTCTAGAGCTAATACGTGCAACAAACCCCGACTTATGGAAGGGACGCATTTATTAGATAAAAGGTCGACGCGGGCTCTGCCCGTTGCTCTGATGATTCATGATAACTCGACGGATCGCATGGCCTCTGTGCTGGCGACGCATCATTCAAATTTCTGCCCTATCAACTTTCGATGGTAGGATAGTGGCCTACCATGGTGGTAACGGGTGACGGAGAATTAGGGTTCGATTCCGGAGAGGGAGCCTGAGAAACGGCTACCACATCCAAGGAAGGCAGCAGGCGCGCAAATTACCCAATCCTGACACGGGGAGGTAGTGACAATAAATAACAATACTGGGCTCTTTCGAGTCTGGTAATTGGAATGAGTACAATCTAAATCCCTTAACGAGGATCCATTGGAGGGCAAGTCTGGTGCCAGCAGCCGCGGTAATTCCAGCTCCAATAGCGTATATTTAAGTTGTTGCAGTTAAAAAGCTCGTAGTTGAACCTTGGGATGGGTCGGCCGGTCCGCCTTTGGTGTGCATTGGTCGGCTTGTCCCTTCGGTCGGCGATACGCTCCTGGTCTTAATTGGCCGGGTCGTGCCTCCGGCGCTGTTACTTTGAAGAAATTAGAGTGCTCAAAGCAAGCCTACGCTCTGGATACATTAGCATGGGATAACATCATAGGATTTCGATCCTATTGTGTTGGCCTTCGGGATCGGAGTAATGATTAACAGGGACAGTCGGGGGCATTCGTATTTCATAGTCAGAGGTGAAATTCTTGGATTTATGAAAGACGAACAACTGCGAAAGCATTTGCCAAGGATGTTTTCATTAATCAAGAACGAAAGTTGGGGGCTCGAAGACGATCAGATACCGTCCTAGTCTCAACCATAAACGATGCCGACCAGGGATCAGCGGATGTTGCTTATAGGACTCCGCTGGCACCTTATGAGAAATCAAAGTTTTTGGGTTCCGGGGGGAGTATGGTCGCAAGGCTGAAACTTAAAGGAATTGACGGAAGGGCACCACCAGGAGTGGAGCCTGCGGCTTAATTTGACTCAACACGGGGAAACTTACCAGGTCCAGACATAGTAAGGATTGACAGACTGAGAGCTCTTTCTTGATTCTATGGGTGGTGGTGCATGGCCGTTCTTAGTTGGTGGAGCGATTTGTCTGGTTAATTCCGTTAATGAACGAGACCTCAGCCTGCTAACTAGCTACGTGGAGGCATCCCTTCACGGCCGGCTTCTTAGAGGGACTATGGCCGTTTAGGCCAAGGAAGTTTGAGGCAATAACAGGTCTGTGATGCCCTTAGATGTTCTGGGCCGCACGCGCGCTACACTGATGTATTCAACGAGTTCACACCTTGGCCGACAGGCCCGGGTAATCTTTGAAATTTCATCGTGATGGGGATAGATCATTGCAATTGTTGGTCTTCAACGAGGAATTCCTAGTAAGCGCGAGTCATCAGCTCGCGTTGACTACGTCCCTGCCCTTTGTACACACCGCCCGTCGCTCCTACCGATTGAATGATCCGGTGAAGTGTTCGGATCGCGGCGACGTGGGTGGTTCGCCGCCCGCGACGTCGCGAGAAGTCCACTAAACCTTATCATTTAGAGGAAGGAGAAGTCGTAACAAGGTTTCCGTAGGTGAACCTGCGGAAGGATCATTG
+>Arabidopsis thaliana 1
+GGATGCGATCATACCAGCACTAATGCACCGGATCCCATCAGAACTCCGCAGTTAAGCGTGCTTGGGCGAGAGTAGTACTAGGATGGGTGACCTCCTGGGAAGTCCTCGTGTTGCATCCCTC
+>gi|186498419|ref|NR_022453.1| Arabidopsis thaliana (AT2G01020) rRNA
+AAAACGACTCTCGGCAACGGATATCTCGGCTCTCGCATCGATGAAGAACGTAGCGAAATGCGATACTTGGTGTGAATTGCAGAATCCCGTGAACCATCGAGTCTTTGAACGCAAGTTGCGCCCCAAGCCTTCTGGCCGAGGGCACGTCTGCCTGGGTGTCACAA
\ No newline at end of file
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/rRNA_react.txt	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,8 @@
+gi|186498419|ref|NR_022453.1|
+0.0349942738544	0.219596062195	0.0958706408233	0.194395004716	0.165304592733	NA	0.649363521545	0.004703371147	NA	0.0	NA	0.0	NA	NA	0.226052814765	0.145104959196	0.235114323931	0.590973505942	NA	NA	0.671823965086	NA	0.600180467352	NA	0.0	NA	0.67204697444	NA	NA	0.0	NA	0.0	NA	0.0	NA	0.0	0.0905633447413	NA	0.0844580185128	NA	0.0	NA	NA	0.387056722073	0.063190678881	NA	0.302648296659	0.123623615737	0.0	NA	NA	0.288529445998	NA	0.0	NA	0.272682267617	1.18142319227	0.222698183429	NA	NA	0.0	NA	0.308507027145	NA	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	NA	1.09905990052	0.331081493258	NA	NA	NA	0.0	0.316741455	NA	0.152869347497	0.358830282006	NA	0.172524825658	0.0232297542077	0.397612802938	NA	NA	NA	0.591410529537	0.531312624449	0.0105725487651	0.0	0.544371401206	NA	0.764412925075	NA	1.31480559071	NA	NA	0.338476420086	NA	NA	NA	NA	0.315926121887	0.0	0.0	NA	0.0	0.708180867703	1.4961467391	NA	NA	NA	NA	0.664465128059	NA	0.0	0.0	0.562535200045	0.476678390556	2.19691885985	0.524311802546	NA	1.61857603596	0.0	NA	NA	0.0	NA	NA	NA	0.0	0.0	NA	1.46994373288	NA	NA	NA	0.0	0.0	0.0	NA	NA	0.0	NA	NA	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	0.0	0.0	0.0	0.0
+25s
+NA	0.445057897744	NA	0.0	0.0	0.0	0.616883939924	0.0	0.0544377002455	NA	NA	NA	0.0	0.274299815387	NA	NA	0.0161734735511	NA	NA	NA	0.113806481313	NA	NA	0.0	0.0	0.0	0.0	NA	0.0902360616756	NA	NA	0.0259497293124	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.128445973988	NA	0.0	NA	0.0	0.0	0.0561933770321	NA	0.0	0.292152642133	NA	0.0	NA	NA	0.0	NA	NA	0.0	0.0	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	0.0	0.0	0.0	1.05193677381	0.112019960608	NA	NA	0.0	NA	NA	0.0	0.0	0.0	NA	NA	0.454696508183	NA	NA	0.0	0.0	0.0541585126648	NA	NA	0.0	NA	0.0	NA	0.0	NA	0.0	0.0	0.0	0.00805821063273	NA	NA	NA	0.0	0.160579429597	NA	0.0212575846522	NA	0.0	0.0	0.0	0.0	NA	0.544027973389	NA	NA	NA	0.425230790667	0.733207830254	0.0	0.0	NA	0.0	NA	NA	0.0	0.0	NA	NA	0.0	NA	NA	0.236885996269	NA	NA	0.0	NA	NA	NA	0.0	NA	0.00917986657696	0.903567358643	NA	NA	NA	NA	0.358438939384	NA	NA	0.0	NA	NA	NA	0.378435644084	NA	0.565879841924	0.255131654267	NA	0.570022482976	NA	NA	0.38092251132	0.0	NA	0.0	0.456921019156	NA	0.34916390672	NA	0.401409258122	0.0396326598324	NA	NA	0.807425004098	0.00846128317379	0.0394311220708	NA	NA	NA	0.21059141834	0.239063190705	NA	0.0979608844846	0.0	NA	NA	NA	0.0474692522788	0.0	0.0	NA	NA	NA	0.656707278731	0.0	0.0761444662556	NA	NA	NA	NA	0.391726941792	NA	0.0	0.0	1.30489991374	NA	0.945708714773	NA	0.72781357993	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	0.0	0.132493170615	0.34803155101	NA	NA	0.129549012388	NA	NA	NA	0.0	0.0	0.29614477645	NA	NA	0.704763040047	0.56956988712	0.0	0.0600574362663	NA	NA	NA	0.362742350085	NA	0.56405635138	0.645775512364	0.896823517909	0.0	0.0540423619053	0.52126275377	NA	0.315047770556	NA	NA	0.0	NA	0.0	NA	NA	NA	0.0	NA	0.156974005797	0.0	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.531410545675	1.18326261838	NA	NA	0.0152722862549	0.127338939497	NA	0.0353599494945	0.0	0.0142090682171	0.679734480705	0.0	0.263031834225	NA	0.985375666218	NA	NA	NA	0.0972886961604	NA	NA	NA	0.847064861501	0.968544261356	0.639603497573	NA	NA	0.366911053247	0.0	NA	NA	0.0	0.0925945255189	0.842712849319	0.788198195629	NA	NA	0.0	NA	0.207548618929	0.343212298205	0.0	NA	0.0	0.832367460529	NA	NA	NA	0.0	NA	0.181366734314	NA	0.0	NA	0.496484060853	0.0	0.0	NA	0.0	NA	0.0	NA	0.0	NA	0.887394863232	0.0	0.0	0.294131087652	0.0	NA	NA	0.0	0.0	0.0	NA	0.304310857899	NA	0.260357268467	NA	NA	NA	0.955488223116	1.26688198445	0.875229842333	NA	0.0	NA	NA	1.63278092007	0.514088552204	0.408365845387	0.652619391573	NA	NA	0.667954630332	0.0	NA	NA	NA	NA	1.01245281715	0.908098734184	0.834658046014	0.374243366292	NA	0.296597361408	NA	0.521754407574	NA	NA	0.412732837452	0.826636773533	0.0	0.325315528841	NA	0.0	NA	NA	NA	0.0	NA	NA	NA	0.297045398287	0.119693277826	0.0	NA	NA	NA	NA	0.0	NA	NA	NA	0.413444210292	NA	NA	NA	0.0	0.614760267946	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	NA	0.0203387523684	0.0	NA	NA	0.662888317611	NA	0.810411698898	NA	NA	0.0	NA	NA	0.0	0.0	NA	NA	NA	NA	0.0388212343243	NA	NA	0.648759772383	NA	NA	0.0	NA	NA	0.303717184162	0.0	0.0	NA	NA	0.535165912151	NA	0.597271148946	0.362095296283	0.0858411830887	NA	0.143598296218	0.0381308456186	NA	NA	NA	0.0215817669597	1.40687529297	NA	0.0	0.0	NA	0.0303834231878	NA	0.216921517271	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	NA	NA	NA	0.0	1.06699879852	0.0	NA	0.0	0.252863258544	NA	0.0	NA	NA	0.292213631555	NA	NA	0.249414499458	0.0	NA	NA	NA	NA	NA	0.473128266176	NA	NA	0.591834136228	0.327887954759	NA	0.73713903679	0.688145314485	NA	0.753552501186	0.0	0.0	NA	NA	NA	0.0445763726679	NA	NA	NA	NA	NA	0.484668453557	NA	0.972325323458	0.204039342759	NA	0.026875539766	NA	NA	NA	NA	NA	NA	1.27487785668	NA	1.01701187793	0.282312329051	NA	NA	0.137105964441	NA	0.0	NA	NA	0.0	0.0	NA	NA	NA	0.0	NA	1.28138528782	NA	NA	NA	NA	NA	0.380108239113	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	NA	0.0	NA	0.0	0.110639273201	NA	0.156758986775	0.232110678863	0.256643329142	NA	NA	1.25139241899	NA	NA	NA	0.844770381547	0.168404091838	NA	0.227649354233	NA	NA	0.101789989504	0.365846871061	NA	0.0	0.0	NA	0.0	NA	NA	NA	0.369098007488	NA	1.20446170804	0.0	NA	NA	NA	0.0	NA	NA	0.749456424764	NA	NA	0.0	0.0	NA	NA	NA	NA	NA	NA	0.614882956175	NA	0.0	0.0	0.0	0.880504978446	0.0	NA	NA	0.461424354396	NA	0.610311556777	0.0	0.0	0.0	NA	NA	0.0	NA	NA	NA	0.0	0.0	0.0	0.0	0.0	0.0	NA	NA	0.0	0.503920101526	0.0	0.249077057033	0.0	NA	NA	0.0	NA	NA	0.0	NA	NA	0.0	0.0	0.0	NA	NA	NA	NA	NA	NA	0.0	NA	0.152786269974	NA	NA	0.0	1.47121659957	1.06909650565	0.700144661674	NA	NA	NA	NA	NA	0.350911905592	NA	NA	0.312682673566	1.6894131616	1.01938882717	0.404565875038	1.51974693619	1.12623464472	NA	NA	1.37154283401	1.2688958914	NA	NA	0.26040033479	NA	0.600549521307	0.539166211903	0.333588872611	NA	NA	0.506909443298	1.41856307216	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	0.132579540803	NA	NA	NA	0.0	NA	0.0	0.221359920289	NA	0.0345263309576	NA	0.432527467916	NA	NA	NA	NA	NA	0.0	0.900492034843	0.119809739847	0.725470179848	NA	0.30434074559	0.889141767065	NA	0.101279262275	0.567410860836	0.0	NA	0.261404253614	0.0	0.0	NA	NA	NA	1.24681193361	NA	0.684881579094	NA	NA	1.17264517378	NA	NA	1.08565570148	NA	1.13223364991	0.553855213729	NA	NA	NA	NA	NA	0.132071626699	NA	0.485748763644	NA	NA	NA	NA	NA	0.0	NA	0.0	0.0	NA	NA	0.359498260565	0.240057267542	NA	NA	NA	0.0504800340658	NA	NA	NA	0.0	0.0	0.0	0.197420570382	NA	0.0137503922192	0.552197932985	0.00654157195938	NA	0.554028257933	NA	NA	NA	NA	NA	0.0115936188099	0.0	0.0	NA	0.0	NA	NA	0.331852524525	0.0	NA	NA	0.224970334328	NA	0.348557316092	NA	NA	NA	NA	NA	1.09270921566	1.22259027297	0.602107735616	NA	0.155315571382	0.0	0.0	NA	0.317736158732	NA	NA	0.401170148918	0.0	0.609158426745	0.0	NA	0.239006446129	NA	NA	NA	NA	NA	NA	0.404763847023	1.46883204568	NA	0.0	0.0	0.0	NA	0.0	0.3617221521	NA	0.278893870402	NA	0.850989559775	NA	0.410969919324	0.393388305487	NA	NA	1.0090276252	0.802590372932	NA	NA	NA	0.937695050443	1.6304479889	1.43736635079	0.792504827655	NA	0.212438225392	NA	NA	NA	0.825240261064	NA	NA	0.519076933164	NA	NA	0.182490592056	0.0	NA	NA	NA	NA	NA	NA	0.525494638628	NA	0.0424499960303	NA	NA	NA	NA	0.0	NA	0.382662671779	1.03662578195	0.478051534735	NA	0.900318367021	0.0	NA	0.0283154099686	0.302896986365	NA	0.0457690881864	NA	0.359568906648	0.101968582452	0.0	0.0	0.419460405434	NA	0.0	NA	0.0	NA	NA	0.259228552941	NA	0.0	NA	NA	NA	NA	NA	0.0	0.257818736114	0.0	NA	0.0	0.0766320737715	NA	0.131365657118	0.0329649802112	NA	NA	NA	NA	0.0	0.157019453832	0.0	NA	0.0452289137179	0.355878493915	NA	NA	0.0	NA	0.0	NA	0.0	NA	NA	NA	0.0	NA	0.0	NA	0.0	NA	NA	0.721547552016	0.0	NA	0.0	NA	0.0506067078662	NA	NA	NA	0.0	NA	0.0	NA	0.0	NA	NA	NA	NA	0.240751379674	0.296570626543	0.228497209418	NA	0.0	0.0	0.202045599166	0.403533950241	NA	NA	0.0277353113723	NA	NA	0.0661029025382	NA	0.20681965929	NA	NA	0.199947408468	0.215123796802	NA	NA	NA	NA	NA	NA	NA	0.822633165225	NA	0.0	0.0	0.0	0.0	NA	0.0	0.0	NA	0.0	NA	0.0635676990266	0.0	0.0	NA	0.243936204533	NA	NA	0.0	NA	0.0	0.0767961818183	0.0	0.0	0.0	NA	NA	NA	0.0	0.0	0.489262555503	NA	0.0	NA	NA	NA	0.0	NA	NA	0.263720493712	0.186243556243	NA	NA	NA	NA	0.0	NA	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	NA	NA	0.0212100286931	NA	0.0	0.0	NA	NA	0.0	0.0	0.0	0.191781514732	0.323808925525	NA	NA	0.0	0.0	NA	0.0	NA	0.0942877016485	NA	0.0	NA	0.0	NA	0.0	0.0	0.0	0.0	NA	NA	NA	NA	NA	0.0	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	NA	0.0	0.0515659980884	NA	0.0	0.0	NA	0.0	0.0728981237199	0.0	NA	NA	NA	0.0	NA	0.0497386479542	NA	NA	0.344131947399	NA	NA	NA	0.0	NA	NA	0.0	0.0	0.0	0.0	NA	NA	0.832992398226	0.0	NA	0.297655103355	0.0	NA	NA	NA	NA	NA	0.0	0.0	NA	NA	NA	NA	0.0	0.0	0.0	0.0	0.0	0.0	NA	NA	0.0	NA	0.0	NA	0.0	NA	0.0	0.563992297514	0.0746194468951	0.239855394121	NA	0.0	NA	0.0413315152918	0.0789091146893	0.0	0.450063900032	0.0	0.0187354091156	0.0	0.805138725071	0.377808848639	0.0	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.422221172566	NA	0.083850305445	NA	NA	1.69187275574	0.797712628988	NA	0.254165673967	0.389917879097	1.29764224015	NA	0.146198168789	1.08772896875	NA	NA	0.25177095078	0.132249177723	NA	NA	NA	NA	NA	NA	0.0	0.513332504368	NA	NA	NA	1.35016047735	0.918964696424	NA	NA	0.638955628195	NA	0.796127761392	0.674153030358	0.279244384849	NA	0.374262073004	0.141562990812	NA	0.0	NA	0.695289932334	0.549909312314	NA	NA	0.242268748968	NA	NA	NA	NA	NA	NA	1.42861304531	0.96200637613	0.723945126863	1.14184709367	0.974881083766	0.160783218804	NA	0.0	0.785116068271	0.129550326685	0.0742103337741	NA	NA	0.0887103295255	0.293326828932	NA	0.0609302788116	0.299174331913	NA	0.148074442128	0.776482856082	0.282593720521	0.292328881825	NA	0.182891160252	NA	0.00902344477709	0.0	0.0	0.0	NA	0.714717858445	0.543998368327	0.284979851972	0.424633242858	NA	NA	NA	0.237626583335	NA	NA	NA	1.49891284952	NA	0.796471346303	NA	NA	0.265860524032	0.0	NA	0.0	NA	0.198705264957	NA	0.194885087665	0.0424079813341	0.0	NA	0.0749223934241	NA	0.0	0.0	0.0	0.0	NA	NA	0.0	0.0	NA	NA	0.00802841612499	NA	NA	NA	NA	0.0	0.31999992992	0.0	NA	0.0694765030851	NA	0.0	0.0308050109141	0.189156337939	NA	NA	0.0	0.0	NA	0.157957315137	NA	0.327825282739	NA	NA	0.181627968758	NA	NA	0.278239963013	NA	NA	0.802814747133	NA	NA	NA	0.0	NA	0.168130304182	NA	NA	0.0093429342157	NA	NA	NA	0.0	NA	0.0214410702375	NA	NA	0.904803223216	0.639136999102	0.935987858518	0.18690133213	0.0	0.0	0.0	NA	0.0	NA	NA	NA	0.0	NA	0.749456424764	NA	0.0	NA	NA	0.0	NA	0.0	NA	NA	0.0145916681828	NA	0.0	NA	NA	0.0	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	1.42644076068	NA	0.223706487904	NA	0.0	NA	NA	NA	NA	NA	NA	NA	NA	0.0	NA	NA	0.0130617898734	NA	0.0	0.0898999792737	0.0546930962057	0.177149998102	NA	0.0741706624637	NA	NA	0.0	0.277016298373	0.0607612564539	0.613418445477	NA	NA	0.582872285854	NA	0.408040736596	0.492172366594	0.0	NA	NA	NA	NA	0.299449789138	0.77269892527	NA	NA	0.0	0.0	NA	0.382989370304	0.0	NA	0.452094370139	NA	NA	NA	NA	0.792394759449	0.744680150853	0.259704725548	NA	NA	NA	NA	0.0	0.0	0.33788703532	NA	NA	NA	NA	1.62621954143	1.04027934063	0.431206826489	NA	NA	0.714243138538	1.34097011843	1.01730053896	NA	NA	NA	0.0	0.172844565107	0.0	0.0328826909876	NA	NA	NA	NA	NA	NA	0.283882554849	NA	NA	0.0	NA	0.0	NA	0.0	0.0	NA	0.287635419384	0.110647612752	NA	0.202972997278	NA	NA	0.0833531998795	NA	NA	NA	NA	NA	0.214849954963	0.464265623903	0.974767502827	0.0	0.588378927844	0.0	NA	NA	1.49891284952	NA	NA	1.32747000763	NA	2.04364253198	NA	1.09013162636	NA	0.0	NA	NA	0.0	0.0	NA	0.0	NA	0.592807348042	0.314970673548	0.0	NA	NA	0.0	NA	0.280396915869	0.379614403384	0.364879338494	NA	NA	NA	NA	0.056664582617	NA	0.0	0.0	NA	NA	NA	NA	0.622332337653	0.255556133808	1.40117902037	0.944664923194	NA	NA	0.450287253117	0.0	NA	NA	0.307473484428	0.27287125866	0.164471335469	0.0	NA	NA	NA	0.157270120359	0.299740996864	NA	NA	NA	NA	0.0413075139923	NA	NA	NA	NA	NA	0.0	0.0	NA	NA	0.0	0.0	0.0	0.0	NA	NA	NA	0.303568964489	NA	NA	NA	0.531856039163	NA	NA	0.0201986551663	0.435382169479	NA	NA	0.0987027824716	NA	0.31342651955	0.00343474041446	NA	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0387128132167	NA	NA	NA	0.57720593517	0.0	NA	0.0988224987239	NA	NA	NA	0.26035391713	NA	0.0	NA	NA	NA	NA	0.143064536743	NA	NA	NA	NA	NA	0.029330256003	0.0	0.0	0.0	NA	0.0	0.256081822821	NA	NA	0.0	0.0	0.0	0.0	0.0182523670474	0.0	0.0	NA	NA	NA	0.0	0.774600225958	0.26663821268	0.0	NA	NA	0.0924618711091	NA	0.0	0.0	NA	0.0	0.0	NA	NA	0.307646651131	NA	NA	NA	0.313241564771	NA	NA	NA	NA	0.0	0.0	0.403326031766	NA	0.0	NA	NA	0.0	NA	NA	NA	0.0549226886851	0.0	NA	0.0	NA	0.0	0.183196997031	0.0630570247126	0.185268839967	NA	0.0	0.174964572408	0.0	NA	NA	0.0	NA	0.0	NA	NA	NA	NA	NA	NA	NA	0.0241371990068	0.0	NA	NA	NA	NA	0.504888168933	NA	0.0	0.213183553079	0.0	0.0	0.0	NA	NA	NA	0.506931246785	NA	0.0	0.0	0.0	NA	NA	NA	0.0	0.322641160202	0.141652321934	0.0	NA	0.252264521562	0.0	NA	NA	0.47525284449	NA	NA	0.942520573299	0.304319464713	0.219700385375	NA	0.0	NA	NA	NA	0.0	0.0	NA	0.0	NA	0.0	0.673987482878	0.0	NA	0.0741770095132	0.0372303546807	0.0	NA	NA	NA	0.0247189995351	NA	NA	0.0	0.0	NA	0.0	0.133958409668	NA	0.606476608455	0.941802285886	0.48275655054	0.454217968731	NA	1.44362500116	1.73619693469	NA	0.0	1.88486859437	NA	NA	NA	0.0	NA	0.0354162387935	0.0	0.241105533729	0.57779902804	NA	NA	NA	NA	0.12919921605	0.0	0.0	0.562347584973	NA	0.0	0.719744679504	NA	0.0	NA	NA	NA	NA	0.709876516692	NA	1.55959736769	NA	NA	NA	0.297851670351	0.0	0.186655077141	0.0289179295796	1.38889882368	NA	NA	NA	0.762646831607	NA	NA	0.0	0.652857745214	0.687614402238	NA	NA	NA	0.712520865755	0.0	NA	NA	0.109301474262	NA	NA	0.147989151834	0.0	0.0804848478165	0.392549716256	0.508391954576	NA	NA	NA	0.0	NA	0.165743013216	0.0	NA	NA	0.298884242816	0.0	0.0	NA	NA	0.683826721934	NA	NA	NA	1.047682093	0.513582124643	1.27357043773	1.31598581539	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	0.0470783180706	NA	NA	0.277242232928	NA	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	0.15071339473	0.0	0.324671043927	NA	NA	NA	0.0402310633871	0.269111626505	NA	1.00950433989	0.616185671929	0.0	0.626989519149	0.637145772493	NA	NA	0.0	NA	NA	0.0	NA	NA	NA	0.0	0.459563800145	NA	0.445174910398	NA	NA	0.573674350403	0.0	NA	NA	0.142500112434	NA	0.0	NA	0.763077300905	NA	0.0	NA	NA	0.364308231584	NA	NA	0.175737376483	0.460129803976	NA	0.265308812468	NA	NA	0.498238947092	NA	0.951525704053	NA	0.320744812576	NA	0.263022765696	NA	NA	NA	NA	0.0	NA	0.0	0.0	NA	NA	0.0	NA	NA	0.0	0.0	NA	NA	1.02113585779	0.0	NA	NA	NA	NA	NA	0.353334515549	0.218282500099	NA	0.309348847372	0.160238205731	NA	NA	NA	NA	0.53337093343	0.545847288001	0.724783656107	NA	NA	0.0	NA	0.383269965231	NA	0.547383577675	NA	0.717142438509	NA	NA	NA	0.373790691961	NA	0.0	NA	NA	NA	0.0	0.0	0.0	1.36961325023	NA	NA	NA	0.0	NA	NA	0.704623310128	NA	0.802736920461	1.02041633588	0.0	0.751436521461	NA	NA	0.430990670005	0.858662368122	NA	0.0	NA	0.0	0.742978589189	NA	0.139090719898	0.253993630943	0.221483843841	NA	NA	NA	NA	0.846560527242	0.813471732436	NA	NA	0.0	0.268751171191	0.562103557971	0.127880516129	NA	NA	NA	NA	0.770878720713	0.4526261378	NA	0.133319038193	0.0	NA	0.391126042068	0.0	NA	NA	NA	NA	NA	0.427740357065	0.142048733194	NA	NA	0.750848291344	0.151052218888	0.0	0.0	0.0	0.142865219696	0.0779391056953	0.362880826221	NA	0.0	0.173655539399	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	NA	0.0	0.0	0.28438912012	NA	NA	0.767157756193	NA	NA	0.0302528408745	NA	NA	0.873517219644	NA	0.0857635382405	0.44958306086	0.280078763569	NA	0.534544257192	0.0	0.123436051388	NA	NA	NA	NA	0.10854493184	NA	NA	NA	0.0	NA	NA	0.0	0.0	0.688340229423	1.01457894475	NA	NA	NA	0.0	NA	0.0	NA	NA	0.393415255785	0.410023899091	NA	NA	NA	0.0	0.0	0.59401830049	0.928373357296	NA	NA	NA	0.582069932474	0.588659773799	NA	1.20348415374	0.0850980939041	0.263774553372	NA	NA	0.0	0.0932578899745	0.230311934834	0.0	0.0	0.208589996817	0.0	NA	0.0	NA	0.0	NA	NA	NA	NA	0.0771848784039	0.0	0.0	0.0	NA	NA	0.127282431712	NA	NA	NA	0.306664035399	NA	NA	0.0	1.26462679711	1.02287296995	NA	0.289475252316	NA	NA	0.0	0.247260442475	NA	0.173659501846	NA	0.138763812029	NA	NA	0.566430199234	0.299412518113	NA	NA	NA	0.498217188314	NA	0.0	0.226901041464	1.21176762632	0.761625388306	0.60166959396	NA	NA	0.164936025044	0.0	NA	0.508914790911	NA	NA	0.0	0.0	NA	0.0	NA	0.0	0.0	NA	NA	0.120782117778	NA	NA	NA	0.403573895969	0.234762284786	NA	0.172461584345	NA	1.17084784873	0.59469814674	NA	NA	0.0	0.421132491162	NA	NA	NA	0.1306794702	NA	NA	0.0	0.0	0.0	NA	0.127513437691	NA	0.0	NA	NA	0.0	0.0	0.0	0.231737505521	0.0	NA	NA	NA	0.0	0.0	0.0	NA	NA	NA	0.0	NA	0.0	0.0	NA	0.204813197391	NA	0.167954420125	0.0	0.0	NA	0.0	NA	1.25963308137	1.65490763074	0.349665039522	0.0	0.0	0.0	0.0	0.088135396391	NA	0.0	0.0	0.052063263328	1.10030351945	NA	NA	NA	0.0	0.0	0.0	NA	NA	NA	0.0	NA	NA	NA	NA	0.0	1.19282107102	NA	0.351500033721	0.0	NA	0.0	0.0	NA	0.0	NA	NA	NA	NA	0.639486692424	1.72246658214	1.27928669427	NA	1.79434484187	1.43487503629	NA	1.11459423067	1.28050635355	0.0	0.033808248598	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	0.037351845069	0.0	NA	0.0	NA	0.0	NA	NA	0.0	0.0	NA	0.0680509646718	0.0	NA	NA	NA	NA	NA	NA	0.216734779888	0.238453108898	0.0047777392186	NA	NA	0.0	0.0	NA	NA	NA	0.24091195962	NA	0.367136308324	NA	NA	NA	NA	NA	0.0364529921341	NA	NA	0.290430964122	NA	0.0292200666431	0.0	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	NA	0.376846428488	NA	NA	0.0538172630509	NA	0.478261010535	0.721287586992	0.582408590163	NA	NA	NA	0.41727516366	0.401127778986	0.723354419621	NA	0.095607199307	0.258806359595	0.0835986123205	0.0	0.0	NA	0.0	0.0	NA	NA	NA	NA	0.440304834094	0.265761431623	0.519660229692	NA	NA	NA	0.168302233778	NA	NA	NA	NA	0.157327271761	0.0	NA	NA	0.0	0.0	NA	0.0	0.0	NA	NA	NA	0.0924344002851	0.0	NA	0.0	NA	NA	0.375769028351	NA	NA	0.0446187487194	0.0	NA	NA	NA	NA	NA	0.516864204729	0.172416681072	0.0567192634098	0.0378721224169	0.00299445087817	0.0932354718163	0.0743292528517	0.0912629153496	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.992289703074	0.0	0.316404114911	0.107107875064	0.520226987768	NA	NA	0.0	NA	0.361658745396	NA	0.0	NA	NA	0.651767453937	NA	NA	0.373975972699	NA	NA	NA	NA	0.0	NA	1.38394570998	NA	0.0	1.00495180372	NA	NA	NA	0.284672456867	NA	NA	0.0	NA	NA	0.431826210178	0.0	0.0	NA	NA	NA	NA	0.0	0.0901844111313	NA	NA	NA	NA	NA	NA	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	NA	NA	NA	NA	NA	1.33336901534	NA	NA	0.191191895748	1.4841175405	0.392189953008	1.45736007574	NA	0.311778885904	NA	NA	NA	NA	0.176831354715	0.673338414422	0.0443852672859	0.0929438936537	NA	0.0	NA	0.0	0.0	0.0	NA	0.0	0.0694646091221	NA	NA	NA	NA	NA	0.0	0.0	NA	0.192693932018	0.288541449812	NA	0.475655710856	NA	NA	0.388660038013	NA	0.00984462464751	NA	0.0	1.81119542895	1.58177307546	0.571881702456	NA	0.841799847816	NA	0.0	1.8803751819	1.55005205033	0.817218685186	NA	1.88597635494	1.3912539199	1.62016777444	NA	0.218813430261	NA	0.209125600132	NA	NA	NA	NA	NA	NA	0.00152129255625	0.0692149342689	0.482781374021	0.286646187881	0.116245486156	0.258681032169	0.0426374923982	NA	NA	NA	NA	1.19938393364	0.632277359129	0.551425712127	0.824991309934	NA	0.00546462277811	NA	0.0	NA	NA	NA	NA	NA	0.190663657047	NA	NA	0.0552614588428	NA	NA	0.340937408688	NA	NA	NA	NA	0.13847541988	0.20559094886	NA	NA	0.781555754406	0.0	NA	0.525171284317	0.436527075283	NA	0.0421028137972	0.0307349474561	NA	0.281933713892	0.0	0.629440274879	0.179786036677	NA	NA	NA	0.11835642342	0.198474100584	0.261444284574	NA	0.0	NA	NA	NA	NA	0.0	0.0	NA	0.26387724218	NA	0.0	NA	0.472854357389	NA	1.74976904497	0.87854869935	NA	NA	NA	0.0	NA	0.151450492281	0.0	NA	NA	0.0	NA	NA	0.351662153841	0.0	NA	NA	NA	NA	0.0701877141513	0.429452307682	NA	0.0	NA	0.0	NA	0.0	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	NA	NA	0.0533163327192	0.0	0.0	0.0	0.0	0.0	NA	NA	NA	0.508254834906	NA	0.0	0.0	0.0	NA	NA	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	0.853405247499	NA	0.0	NA	NA	NA	0.0	0.0	NA	0.800049449505	NA	0.0	NA	0.0381971897849	0.0	NA	NA	NA	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	NA	0.0	0.0	NA	NA	0.0	NA	0.653339817705	NA	NA	NA	0.0	NA	NA	0.0	NA	0.0	NA	NA	0.0	0.0	NA	0.0	NA	0.0	0.0	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0	NA	0.678467855295	0.0	NA	NA	0.0	0.0	0.0	0.0	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	0.0	NA	NA	NA	NA	NA	0.0	0.0	0.0	0.0	0.0	0.018664182086	0.0	0.0	0.0	0.0	NA	0.0	NA	NA	NA	0.0	0.121678706886	0.355299241427	NA	NA	NA	1.73293458871	NA	0.0493727138012	NA	NA	NA	NA	NA	NA	NA	0.0568785504168	NA	0.134312284659	0.0	0.0	NA	NA	0.0	NA	NA	NA	0.433579288967	NA	0.445397859363	0.0	0.975058935559	NA	NA	NA	NA	0.334019900058	NA	NA	NA	NA	NA	0.819899325822	0.0	0.0	0.0	NA	0.0	0.0	NA	NA	0.0	NA	NA	0.14790191642	0.355048973175	0.632457142469	NA	0.0885424973481	NA	NA	0.0	NA	0.245308306361	NA	1.45476259015	NA	0.461488242458	NA	NA	0.0403145742559	0.0669040008976	NA	NA	0.189059523364	0.223383021291	0.443345417968	0.302699986613	0.261935974538	NA	0.453463826854	NA	NA	0.281780796046	0.0	NA	0.222728714314	NA	0.496511941732	NA	NA	0.648092143251	0.867025205093	0.698984683105	0.342390556037	NA	NA	NA	NA	0.0	NA	NA	0.0	NA	0.0	0.0	0.0	0.8816043907	0.0	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0	NA	0.143927736109	NA	0.0	NA	NA	NA	NA	NA	NA	NA	0.0	0.643731475471	1.04357444649	0.662993815823	NA	0.0	0.0461426678103	0.0	NA	NA	NA	NA	0.0	NA	0.0	NA	0.0	0.033505630595	NA	0.0	NA	0.133469644664	0.0	0.234879130863	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	0.688130547851	NA	NA	0.732542656381	NA	NA	0.0	0.0	NA	NA	0.438119294374	0.558376720682	0.0	NA	0.918102306041	0.56240697358	NA	0.0	NA	0.237088351477	0.0	NA	NA	0.0	0.20233663462	NA	0.600412306412	0.183536339486	NA	0.0	0.0	NA	NA	NA	0.0	NA	0.0	NA	0.732387193297	0.63658249172	NA	0.010121757995	NA	0.0	0.0	NA	0.0	0.641739077143	NA	NA	0.0731629955416	NA	0.0	0.0	0.0	NA	0.0	0.0	NA	0.0	0.0	0.0	NA	0.0	NA	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	0.0	0.0	NA	NA	0.389549029322	NA	NA	0.0	NA	0.472854357389	NA	NA	0.0	NA	NA	0.0	NA	0.0	0.0	0.0	0.0	0.37028839264	NA	NA	0.0	0.0481524671794	0.387465196563	NA	NA	NA	NA	0.0	NA	NA	NA	NA	NA	0.0	NA	0.640486051737	0.354013107665	NA	NA	0.337971646441	0.0	NA	NA	NA	0.0	NA	NA	0.0	NA	NA	0.447020729153	0.309541139682	0.0	NA	NA	NA	0.686037910468	NA	NA	NA	0.3406752016	NA	NA	0.146691758883	0.0	0.0	NA	0.0	0.0	NA	NA	NA	0.913095760914	1.03984524992	NA	NA	0.868123808953	NA	0.273507880758	0.0	NA	NA	0.0	0.0562709433796	0.0	0.117755450683	1.09301118286	0.3406752016	NA	0.279937783278	NA	0.472854357389	NA	NA	0.0	NA	NA	NA	NA	1.1279640247	NA	2.021448478	0.0	NA	0.940786080429	1.03814311506	NA	NA	NA	NA	0.844770381547	0.661304456557	NA	0.741654233914	0.0	NA	0.0975430676565	0.0437438686722	NA	NA	0.594144441129	0.350319089768	1.93721008621	NA	1.22231078215	0.0	NA	0.0	NA	0.0	0.538436231703	NA	NA	NA	NA	NA	0.0	NA	NA	NA	NA	0.0508327001275	0.279131420348	NA	NA	NA	NA	0.0	0.109423953351	NA	0.0	NA	NA	NA	NA	0.271679433027	0.0	NA	NA	NA	0.0	NA	NA	0.0	0.0	0.0	0.46793172304	NA	0.0	NA	0.0	0.0	0.0	0.0	NA	NA	0.749456424764	NA	0.0	NA	NA	0.472854357389	1.22231078215	NA	0.0	0.0	0.0	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	0.0	0.0	NA	0.0	NA	NA	0.0	NA
+Arabidopsis
+NA	NA	0.0	NA	NA	3.53700365182	NA	0.0	NA	0.0	0.0	NA	0.0	0.0	0.0	5.10017382764	NA	0.0	0.0	0.0	NA	0.0	0.0	NA	NA	0.0	0.0	0.0	0.0	NA	NA	0.0	NA	0.0	0.0	0.0	1.04740785071	NA	0.751507319045	0.679065423588	NA	0.0	1.80093151001	0.0	NA	0.0	0.0	NA	0.0	0.0	NA	NA	NA	0.0	0.0	NA	2.35404575646	NA	NA	NA	4.93385220878	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	4.45571252703	NA	NA	0.0	NA	NA	0.0	0.0	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	NA	4.93385220878	0.0	0.0	NA	0.0	0.0	NA	NA	NA	NA	0.0	0.0	NA	NA	0.0	0.0	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	NA
+gi|20197903:2706-4513
+NA	0.012685795205	0.0	0.0	NA	NA	NA	NA	NA	NA	0.369865767367	NA	0.0	0.0	NA	NA	0.0	0.0	0.0	NA	NA	0.291074839291	NA	NA	0.0	0.0779489717102	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	0.388739681544	0.607358801578	0.86505141743	0.419532419715	NA	0.775525397956	NA	NA	0.418634405395	0.310232647615	NA	0.0	0.0	0.51629506558	NA	NA	0.0	0.435316902631	NA	NA	NA	NA	NA	0.406071541412	0.102260948087	NA	NA	1.02827139604	NA	NA	0.555559168062	0.301220165622	0.0	NA	0.428351005927	0.0	NA	NA	0.0942381418988	0.340249958839	NA	0.402129377649	0.639234686516	NA	NA	NA	NA	0.935138133481	0.0592665799408	0.143289116981	0.0473507766862	NA	NA	0.0	NA	0.0266155577151	0.0982764742903	NA	NA	NA	0.0	NA	0.0	0.11309561083	NA	NA	0.628690973329	0.479906074642	0.102631619096	NA	0.0	0.434378189774	NA	NA	NA	0.0807132307206	NA	0.0	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.359209156119	NA	NA	NA	NA	0.648363725715	1.00111919062	1.04262692243	NA	0.755897037378	0.819932067068	NA	0.455162779021	1.20773264872	NA	0.0	NA	NA	0.296441851982	NA	0.37547834061	0.357493014857	0.0258046929258	0.0	NA	NA	0.505600332852	NA	NA	0.215493760216	0.137308939841	NA	NA	0.0	NA	0.327508086131	NA	0.0	NA	0.0	NA	0.110060629173	0.144925722316	NA	0.181195205113	0.198697622929	NA	NA	NA	0.430508890438	0.615380473507	0.413085955445	0.213399108191	0.970869271469	0.940875618616	0.816055880816	0.400666932367	0.0	0.0368366876912	0.0	NA	0.982940202181	0.0152093127598	NA	NA	0.772061687224	NA	NA	NA	0.524486960169	1.14258931334	NA	NA	NA	0.772909687503	0.350161317332	NA	0.0	0.339428903455	NA	NA	NA	0.289638535352	NA	NA	0.301380368516	NA	0.186709719344	NA	0.723493077315	0.549043071946	0.368933272512	0.346567668516	NA	NA	NA	0.0736517765796	NA	0.225216490874	0.0	NA	0.0	NA	NA	NA	0.282065581517	NA	0.0622206275738	NA	NA	0.0	0.0	0.0	NA	NA	NA	NA	0.0979798698289	NA	0.418716246672	NA	NA	0.127229218645	NA	NA	0.122048875919	NA	NA	0.0	0.459466284108	NA	NA	0.22611651355	NA	0.257988938793	0.242043547687	0.0445554892629	NA	0.0	NA	0.211377998276	0.0539458902463	NA	NA	0.705263055871	NA	0.122896017727	NA	0.0	0.364816537092	NA	NA	NA	0.0	0.0	NA	0.0	NA	NA	NA	NA	0.0	NA	NA	NA	0.134627940359	NA	0.384818842789	0.0	NA	0.0	0.298065186972	NA	0.0	0.148933065236	NA	NA	0.108571364914	0.245360133204	0.412047373188	0.0776738288945	NA	NA	NA	0.0	NA	NA	0.385371255467	0.0134780068585	0.186653352588	NA	0.670459919659	NA	0.0	0.163864042171	0.13720252139	0.0	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	NA	0.0	NA	NA	0.0478979106433	NA	0.0	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	0.215656364652	NA	NA	NA	NA	NA	NA	NA	0.676399501898	0.444238412823	0.0	NA	NA	NA	NA	NA	0.402502781349	0.324546337434	NA	NA	0.72428064676	NA	0.365349818072	0.282000801415	NA	NA	0.825136824379	NA	NA	NA	NA	NA	0.0	NA	0.285295564149	NA	NA	0.215476057191	0.182212745127	NA	NA	0.146330023967	NA	0.421564261825	NA	NA	NA	0.353084055047	NA	0.0	0.0	NA	NA	0.0	NA	0.236357075379	0.0	0.244357718994	0.128121995322	NA	NA	0.175198188965	NA	0.218494857155	0.0	0.0	0.00451933044375	0.0	0.927212739519	NA	0.132103213724	0.0	0.740927800969	0.957180040116	NA	NA	0.656796842311	0.580572417551	NA	NA	0.0	0.0	NA	0.0	0.0563230505075	NA	NA	0.120005485202	NA	0.226509753369	NA	0.0	0.471160925233	0.0877285693903	0.768176413935	NA	NA	0.235484190643	0.0	0.0	0.118062685793	0.419161295576	0.135347828419	NA	0.219059490441	0.426910429649	NA	NA	0.971233728998	0.0	0.0	0.48732449835	NA	NA	NA	NA	0.24138759394	NA	NA	NA	0.253899493465	NA	NA	NA	0.825841393383	0.82831906089	0.0	0.0717935208264	NA	0.279778367449	0.282811723237	0.229633772945	NA	0.0	0.0	0.0	0.509924844919	0.429039400771	NA	0.283392930386	0.0	NA	NA	NA	NA	0.0	NA	0.0	NA	NA	NA	0.117058667345	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0965635309754	0.092863243134	NA	NA	NA	NA	1.04941610567	1.19161444221	NA	NA	1.7376246406	NA	NA	0.391195828095	0.375473658421	0.433647125904	0.279359521091	NA	0.501405696729	NA	0.687660733308	0.864175250339	0.79158652207	NA	1.20320006049	0.56311253719	0.253888387685	NA	NA	1.83397160852	1.53706044242	0.961584043251	NA	0.784802938308	NA	NA	1.28330857645	NA	0.268894203299	0.314451464739	0.423300010695	NA	NA	NA	NA	0.259483681982	NA	NA	NA	1.26295217444	0.661502700034	0.652566369353	NA	NA	0.500429839944	NA	NA	NA	NA	NA	0.103898205626	0.255468262008	0.163276078797	NA	0.48540533193	0.395593446683	NA	0.2174452409	0.101616378857	NA	0.205131663835	NA	NA	NA	1.31997685379	0.753084647947	NA	NA	0.0519500996846	0.387428790866	0.34011021097	NA	0.0748065258161	NA	0.0580149462222	0.0768278927631	0.158589807897	0.25491474666	NA	0.347714135251	NA	0.11241524085	NA	NA	0.336941565822	NA	0.0887724160932	NA	NA	NA	0.916896362256	0.278780105213	NA	NA	NA	NA	NA	NA	NA	0.115173560016	0.31061772735	NA	NA	NA	0.486301566768	0.347179567168	0.617356828937	0.444803810174	0.408511448093	NA	0.0782866766856	NA	0.0	NA	NA	0.515740275652	NA	NA	NA	NA	0.61541889439	0.97298662914	0.462906210807	0.542827568418	NA	NA	NA	NA	NA	0.694038563065	NA	NA	NA	NA	NA	0.109092812624	NA	NA	0.0	0.117730829926	NA	NA	NA	0.134948233631	0.152419220299	NA	0.0511684766977	0.0	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.445995841454	0.47369948038	NA	NA	NA	NA	NA	0.040184020931	NA	NA	0.0	NA	NA	NA	NA	0.183024542293	0.464959665882	0.0240688858908	NA	NA	0.261038892925	NA	NA	NA	0.0905114061626	NA	NA	0.330123652207	NA	0.547654296191	NA	0.148523577717	0.0353180022104	NA	0.710181117619	NA	0.0	0.0498588966609	NA	NA	NA	NA	0.037765972332	NA	NA	0.293764451656	0.0	NA	NA	NA	NA	0.0	0.0	NA	NA	NA	NA	0.2319218601	NA	NA	NA	0.0	0.0007880491039	NA	0.0	0.0	NA	NA	0.484697941103	NA	0.48957925096	NA	NA	NA	NA	0.178775276637	0.0	NA	NA	NA	NA	0.204280518962	0.138427143763	NA	0.226834291993	0.582955027497	0.221502626079	NA	NA	0.172379156747	NA	0.00774474492492	NA	NA	NA	0.0	NA	0.0	0.0812527197534	0.139070064757	0.826428013686	NA	0.0	0.566373614269	0.297167065861	NA	0.0	0.43208771255	NA	0.585964397663	0.264915308651	NA	0.173247086337	NA	0.764112125968	NA	NA	NA	0.965842425791	NA	0.346877570077	0.0	0.400547148974	NA	NA	0.0400887452555	NA	0.0	0.0	NA	NA	NA	NA	0.233593147394	NA	0.238066822223	0.35044287177	0.0	1.22637073435	NA	0.588880241101	1.28674125846	NA	0.816194027366	NA	NA	0.21624168024	NA	NA	NA	0.0	NA	0.460702630062	NA	0.0844162750446	0.213999087673	NA	0.543841759443	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.0682752748169	0.0	NA	NA	0.0	NA	NA	NA	0.234279937188	NA	0.174056981406	NA	NA	0.085017264826	NA	NA	0.153333534594	0.0	NA	NA	0.703295059249	NA	NA	0.0	0.0	0.0	0.0216527294522	NA	NA	NA	0.0	0.0	0.0	NA	NA	0.0190758579281	NA	NA	NA	NA	NA	0.0	0.167793515723	NA	NA	0.0	NA	NA	0.0	NA	NA	NA	0.0	0.0	NA	0.0	NA	NA	0.0	0.0562847889823	NA	0.460573683496	NA	NA	NA	NA	0.223670417655	0.233256106881	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	NA	NA	0.0	NA	NA	0.084028304493	0.505561356557	0.0623049297976	NA	0.0	0.0	NA	0.833699137957	0.403279901024	0.0	0.0	0.0322768659755	0.0	NA	NA	0.0	NA	0.0	0.0	0.0	NA	0.0	0.00045836350462	NA	NA	NA	NA	0.0	0.0	0.0	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	NA	0.0	0.0	NA	NA	0.193128657755	0.0	NA	0.0	0.899536353647	0.363201566969	NA	0.0	0.0	0.0	NA	0.0	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	NA	0.0	NA	0.0	NA	0.0	0.110477671057	NA	0.0	0.0	NA	0.0	NA	0.0	0.0	NA	0.182880927986	NA	0.108818565811	0.0	0.0	NA	NA	0.0	0.0	NA	0.0	NA	NA	0.0	NA	0.0	0.0	0.0	0.0	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	0.0210625394541	NA	NA	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	NA	NA	0.0	NA	0.0	0.0	NA	0.0	NA	NA	0.0	NA	NA	NA	NA	NA	0.0385781896847	NA	NA	0.0	NA	0.0	NA	NA	0.0	0.0	NA	0.0	0.0	NA	0.0	NA	NA	NA	0.0	0.0	0.0	0.0	NA	NA	0.0554462223962	NA	NA	0.0	NA	0.0	0.0	0.0	NA	0.0	0.0	0.0	0.0	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.0	0.0	NA	NA	NA	NA	NA	NA	0.0	NA	NA	0.446085727425	NA	NA	NA	NA	0.0	NA	0.0	0.0	0.0	NA	NA	0.0	NA	NA	0.146663316857	0.0	0.0	0.0	NA	NA	0.0	0.285522876287	1.41360776485	NA	NA	0.188205541415	0.0	NA	NA	NA	0.0453327030088	0.0	NA	NA	0.0	0.0492143521014	NA	NA	NA	0.0	0.0	0.0291224196667	0.0	0.0	0.0	0.0	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	NA	0.206914122658	0.0	NA	NA	0.0839474875882	NA	NA	0.0	NA	NA	0.134142893649	0.0	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.820586591704	0.0	0.0	0.0	0.386179006774	NA	NA	NA	NA	0.103032523431	0.0518166812743	0.0	0.120233931127	NA	NA	0.0	0.0	0.0	0.366525621541	NA	NA	NA	0.164870879261	0.0	1.38881890967	NA	0.404286868276	0.0	0.0	NA	0.178472663189	NA	NA	1.31563377475	1.30805369088	NA	NA	0.46978614339	NA	NA	NA	0.390001422112	0.333118860081	0.562108250761	NA	0.946892023632	0.448824683197	NA	NA	0.635745699413	NA	0.46622236075	NA	1.57309356534	NA	0.784800876851	NA	NA	NA	0.516184760729	NA	NA	NA	0.434129860905	NA	NA	0.0250374027677	NA	0.403626419597	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.374841885012	0.340220603837	NA	NA	NA	0.0	0.263065361814	NA	NA	NA	0.269610329411	NA	NA	0.317151321718	NA	NA	NA	NA	NA	NA	NA	NA	0.258599594016	NA	0.110617236497	NA	0.0	NA	NA	NA	NA	NA	0.324479568493	NA	NA	NA	NA	NA	0.325987512869	0.250209550258	NA	NA	0.0774563184444	0.112366803266	NA	NA	NA	0.102919447072	0.123330374332	NA	NA	0.348326050915	0.224773446775	0.0	NA	0.0808889707107	NA	0.214762378779	0.0287355857736	0.0	NA	0.252725257437	0.314236062616	NA	0.0554561654111	0.056815664406	NA	NA	0.345500702375	NA	0.0588579499196	0.0105952720591	0.0456298120613	NA	0.115106069305	NA	0.0199334602017	NA	0.577423401418	0.0462313906594	NA	NA	NA	NA	0.265048014067	NA	NA	0.419418099209	0.113243020283	NA	0.376066546077	0.208296306638	0.0	NA	NA	0.0	0.363219015946	0.198400194261	NA	NA	0.14356423081	0.379033229976	NA	NA	0.0	NA	NA	0.0	NA	NA	0.285411491071	NA	0.343270380421	NA	NA	NA	0.154003422396	0.0913176962236	NA	0.414093438881	NA	NA	NA	0.0	0.216108105084	NA	NA	NA	NA	0.405657537386	NA	NA	0.0	0.0	0.398647870694	0.466721779026	NA	NA	0.362549468307	0.0	NA	NA	NA	NA	NA	0.078297350411	NA	NA	0.0	0.52712134107	0.210547080514	NA	0.084029706675	0.0851834255013	0.0	0.546016042147	NA	NA	NA	0.0	NA	NA	NA	NA	0.210098470701	NA	NA	0.00513323027424	0.0	0.0	NA	NA	0.0872105330697	NA	0.173770776458	NA	NA	NA	NA	0.0	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.0530451107643	0.0	NA	0.0	NA	0.0	NA	0.0	NA	0.0	0.0	0.0314542857359	0.0	NA	NA	0.0	NA	NA	NA	0.0	NA	NA	0.0	0.0	0.184888572061	0.0	NA	0.299087013424	NA	NA	NA	0.253824308291	0.476498867804	0.637367560491	0.557762304345	0.0424647165889	0.0227038975612	NA	NA	NA	NA	0.0	0.164048813678	NA	0.174332609454	0.0887971281539	0.517338202643	NA	NA	0.0	0.0	0.0356275913231	NA	NA	NA	NA	0.219870784412	0.16181267273	NA	0.0254540277981	NA	NA	NA	NA	0.0136899888528	0.0	0.0	NA	NA	NA	0.0	0.0	NA	0.0	NA	NA	NA	0.124832560383	NA	NA	NA	NA	NA	0.291168607862	NA	0.0	NA	0.0	NA	0.0	0.0	NA	NA	NA	0.0	0.0927979100847	0.300974630983	NA	NA	NA	NA	NA	NA	NA	NA	0.0	NA	NA	0.0704858304923	0.0	0.0	0.0	NA	0.0	NA	NA	0.0835696851504	0.0518565943435	NA	NA	0.0	0.0335320265405	NA	0.0	NA	NA	0.12361132429	0.0692121040509	NA	0.0	NA	0.260064432493	NA	0.0959604666378	NA	NA	0.0	0.0904897461535	NA	0.096024873638	0.551806072603	NA	0.0	NA	0.0	NA	0.282845890039	NA	NA	NA	NA	0.0	0.0	NA	0.0	0.206168383343	NA	NA	0.0337910137675	0.237818480755	0.0	NA	NA	0.0	0.0	0.0	NA	NA	NA	NA	NA	0.177374574	0.0	0.0	0.100596272282	0.340124774949	0.336749926715	0.0	NA	0.0	0.0	0.0	NA	NA	0.0	NA	0.0	NA	0.0	0.0	NA	0.0	0.0	0.0	NA	0.0	NA	NA	NA	0.0	0.0	NA	NA	0.0	NA	0.0	0.0	NA	NA	NA	NA	0.226897983481	0.0	NA	NA	NA	NA	NA	0.0	NA	NA	0.0734972121117	NA	0.0	NA	0.572649935308	NA	NA	0.383893601435	NA	0.0628733418693	0.296329616229	NA	NA	NA	NA	NA	NA	NA	NA	NA	NA	0.255561953725	NA	0.0	0.157858376637	NA	0.393736372309	0.0	0.119742668772	NA	0.393736372309	NA	1.01428769186	0.0	NA	NA	0.272044062708	NA	0.733760005611	NA	0.0660086032582	NA	0.900327296663	0.879392059123	NA	NA	0.0898683975317	0.0	0.0	0.0	NA	0.377579916714	1.12452142695	0.0	1.64185114287	0.029392400791	NA	NA	0.0	NA	0.490308287953	0.0	NA	NA	NA	0.0	NA	0.980828030569	NA	NA	0.746093235966	0.0	NA	NA	0.0	NA	0.54743287524	0.821826359884	NA	NA	0.720669950832	NA	NA	0.746435235532	0.0	0.0	0.852194994485	0.496018177117	NA	NA	NA	NA	NA	0.914227544844	0.0	NA	NA	0.0	NA	NA	NA	NA	0.0	0.0	0.0	0.0	NA	NA	0.0	NA	NA	0.0	0.0	NA	NA	0.0	NA	0.0	0.0	NA	NA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/read_file.py	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,21 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+import sys
+
+
+
+def read_t_file(in_file):
+    f = open(in_file);
+    result = [];
+    for aline in f.readlines():
+        temp = [];
+        tline = aline.strip();
+        tl = tline.split('\t');
+        for i in range(0, len(tl)):
+            temp.append(tl[i].strip());
+        result.append(temp);
+    f.close();
+    return result;
+
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/rtts_plot.py	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,56 @@
+#!/usr/bin/env python
+#Make a plot of reactivity distribution
+
+import sys
+import numpy as np
+import matplotlib
+from pylab import *
+import math
+
+#Convert the reactivities (Make NA to 0)
+def convert_react(a):
+    r = []
+    for i in range(len(a)):
+        if a[i]!='NA':
+            r.append(float(a[i]))
+        else:
+            r.append(float(0))
+    return r
+        
+
+#Make a plot of the distribution
+def make_plot(ar,id_s,path):
+    N = len(ar)
+    a = convert_react(ar)
+    w = 1
+    ind = np.arange(N)
+
+    fig = figure()
+    fig, ax = subplots()
+    ax.bar(ind+w, a, width = w, color = 'r',edgecolor = 'r')
+    ax.set_ylabel('DMS Reactivity')
+    ax.set_xlabel('Nucleotide Index')
+
+    
+    mag = int(math.log(N,10))-1
+    tail = 10**mag
+
+    intervel = int(math.ceil(float(N)/tail)/5)
+    print(N)
+    print(intervel)
+    tl = []
+    k = 0
+    ax.set_xticks(np.arange(0,N,intervel*tail))
+    print(np.arange(0,N,intervel*tail))
+    ax.set_xticklabels(np.arange(0,N,intervel*tail))
+
+    ax.set_title(id_s+" reactivity distribution")
+    savefig(path+id_s+'.tif')
+
+
+
+    
+    
+    
+
+
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/predict/tool_dependencies.xml	Mon Sep 15 14:52:43 2014 -0400
@@ -0,0 +1,9 @@
+<?xml version="1.0"?>
+<tool_dependency>
+    <package name="biopython" version="1.61">
+        <repository changeset_revision="ae9dda584395" name="package_biopython_1_61" owner="biopython" prior_installation_required="True" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+    <package name="numpy" version="1.7">
+        <repository changeset_revision="ef12a3a11d5b" name="package_numpy_1_7" owner="iuc" prior_installation_required="False" toolshed="http://toolshed.g2.bx.psu.edu" />
+    </package>
+</tool_dependency>