Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 4 additions & 4 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -47,8 +47,8 @@ Depending on the size of sequence set and on desired number of seeds to be analy
Input files:
- Necessary:
- `rna_fastafile`: a fasta file with RNA sequences of interest
- `exp_values_file`: expression values in a csv format. It should have 2 or more columns, one column with names of of sequences, another column with values. The names of these two columns have to be specified with the `--anno_name_column` and `--measur_column` arguments. All the names of the sequences listed in this file must also have a corresponding sequence record in the fasta file provided with `rna_fastafile`. The measurement values might be integer or float numbers.
- `seeds_file`: a binary file containing the seeds to search through. Such file can be created with `seed_generator.py` script
- `exp_values_file`: expression values in a tsv format. It should have 2 or more columns, one column with names of of sequences, another column with values. The names of these two columns have to be specified with the `--anno_name_column` and `--measur_column` arguments. All the names of the sequences listed in this file must also have a corresponding sequence record in the fasta file provided with `rna_fastafile`. The measurement values might be integer or float numbers.
- `seeds_file`: a binary file containing the seeds to search through. Such file can be created with `seeds_generator.py` script
- Optional:
- user can include a file with RNA structure probing data (SHAPE or DMS-seq) to guide the possible match selection. There is no commonly used standard format for SHAPE RNA reactivity data; therefore, we are using the two-column SHAPE file format used by RNAstructure package ([link](https://rna.urmc.rochester.edu/Text/File_Formats.html#SHAPE)). SHAPE file provided by user should contain SHAPE profiles for multiple sequences, separated with `>`, like in fasta file. SHAPE file can be provided to the `filter_profiles_by_folding.py` script with the `--shape_profile` argument

Expand All @@ -69,7 +69,7 @@ pyteiser_pipeline --rna_fastafile <path to test_seqs.fa> --exp_values_file <path
Arguments for the automatic pipeline (parameters for all the individual steps included):
- input / output files:
- `rna_fastafile`: fasta file with RNA sequences, see above
- `exp_values_file`: expression values in a csv format, see above
- `exp_values_file`: expression values in a tsv format, see above
- `anno_name_column`: column name in exp_values file that contains annotations, see above
- `measur_column`: column name in exp_values file that contains expression measurements, see above
- `seeds_file`: file with seeds in binary format, see above
Expand Down Expand Up @@ -97,7 +97,7 @@ Arguments for the automatic pipeline (parameters for all the individual steps in
- `jackknife_n_samples`: how many permutations to do in jackknife test
- `jackknife_fraction_retain`: what fraction of the sample to retain for each test
- `jackknife_min_fraction_passed`: what fraction of all iterations should
- arguments for `seed_generator.py` script
- arguments for `seeds_generator.py` script
- `outfolder`: output folder
- `prefix`: prefix for naming the seed file
- `num_motifs_per_file`: maximal number of seeds to write into a single file
Expand Down
6 changes: 4 additions & 2 deletions pyteiser/seeds_generator.py
Original file line number Diff line number Diff line change
Expand Up @@ -3,8 +3,10 @@
import os
import sys

from . import glob_var
from . import structures
import glob_var
import structures



# this function generates all possible seeds with specified length of stem and loop
# to reduce the search space, it only keeps the seeds that have a pre-specified number of informative bases (non-Ns)
Expand Down
4 changes: 2 additions & 2 deletions pyteiser/structures.py
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@
import struct
import sys

from . import glob_var
import glob_var


class w_motif:
Expand Down Expand Up @@ -531,4 +531,4 @@ def copy_n_motif(motif):

motif_copy = n_motif(motif.stem_length, motif.loop_length, copy_sequence, copy_structure)

return motif_copy
return motif_copy