my_first_tool

This program consists of 3 tools:

• 'dna_rna_tools' - List of procedures: transcribe (return the transcribed sequence), reverse (return the reversed sequence), complement (return the complementary sequence), reverse_complement (return the reverse complementary sequence);
• 'protein_tools'- List of procedures: length (return the number of amino acids in protein sequence(s)), percentage (return percentage of each amino acid in sequence), pattern (return all non-overlaping instances of a given pattern in sequences), 3Letter_name (return three-letter amino acids into a three-letter amino acids), DNA_code ( return transformed protein sequence(s) to DNA sequence(s)), fastq_tools (return percentage of each amino acid in sequence);
• 'fastq_tools' - Procedure: filtering of the dictionary of sequences by the length, quality of sequencing of each nucleotides and GC% content.
• 'bio_files_processor' - Procedure: converting the multilined fasta to one line fasta

run_dna_rna_tools.py

description of how the run_dna_rna_tools.py program works

This program contains the function run_dna_rna_tools. The run_dna_rna_tools function takes as input an arbitrary number of arguments containing DNA or RNA sequences in the form (str), as well as the name of the procedure to be performed, specified as the last argument. After this, the command performs the specified action on all transmitted sequences. If one sequence is supplied, a string with the result is returned. If several are submitted, a list of strings is returned.

Use example

run_dna_rna_tools('ATG', 'transcribe') # 'AUG'
run_dna_rna_tools('ATG', 'reverse') # 'GTA'
run_dna_rna_tools('AtG', 'complement') # 'TaC'
run_dna_rna_tools('ATg', 'reverse_complement') # 'cAT'
run_dna_rna_tools('ATG', 'aT', 'reverse') # ['GTA', 'Ta']

protein_tools.py

Discription how the protein_tools.py works:
This program contains the function protein_tool. The protein_tool function takes as input an arbitrary number of arguments in the form of amino acid (aa)/protein sequences of type str, as well as the name for the procedure to be performed. After this, the function performs the specified action on all provided sequences. Carefully read the rules of usage for each option, because they specify correct ways of entering arguments, as well as the output and the type of data in the output.

⚠️ Attention: 1) The programm is register-dependent; 2) Before using some of the options read 'Procedures description' carefully. 3) If you input sequenses or 'options' incorrectly, the program will provide you with helpful error messages.

compare Introduction The compare procedure compares the first amino acid sequence provided with the following ones. Inputs To start using the length procedure, enter sevreal arguments:

• an arbitrary number of sequences, where the first sequence is a reference to which the following sequences are compared; each argument should be of type 'str'.
• second-to-last argument is the number of decimals to round the number to; type 'int'
• last argument determines whether percentages are returned instead of fractions; type 'bool' Outputs It returns a 'dict' object where:
• keys are compared-to sequences (type str)
• values are either fractions or percentages (type float). Usage example

protein_tool('LAlLAlwWGPdPA', 'LAlLAl', 3, False, options = 'compare') # {'LAlLAl': 1.0}
protein_tool('LAlLAlwWGPdPA', 'LAlLAl', 'GPdPA', 3, True, options = 'compare')) # {'LAlLAl': 100.0, 'GPdPA': 20.0}

length Introduction The length procedure calculates the length of protein sequence(s) (equal to the number of amino acids). Inputs To start using the length procedure, enter one or more protein sequences for which you want to get a summary, and at the end add options = ‘length’. Outputs The result of the procedure is a list with the numbers of amino acids in each sequence. The list contains only numbers of amico acids in the sequence. Usage example*

protein_tool('LAlLAlwWGPdPA', options = 'length') # [13]
protein_tool('RRRrrrR', 'WGPdPA', 'LAlLAlw', options = 'length') # [7, 6, 7]

percentage Introduction The percentage procedure calculates the percentage of all 20 proteinogenic amino acid residues, case-sensitive in the protein sequences Input To start using the count_percentage procedure, enter one or more protein sequences for which you want to get a summary, and at the end add options = ‘percentage’. Outputs The result of the procedure is a list of dictionaries with the percentages of the corresponding amino acids in each sequence. The dictionary contains only amino acid residues whose percentage in the sequence is not equal to 0 (which are contained in the sequence at all). Also, the dictionary is ordered from the largest percentage of content to the smallest. Cases of amino acid residues are taken into account.

⚠️ Attention: We use rounding to 2 decimal places. In some cases, the sum of percentages of all amino acid residues for sequence may not be exactly 100% due to rounding. Usage example

protein_tool('LAlLAlwWGPdPA', options = 'percentage') # [{'A': 23.08, 'L': 15.38, 'l': 15.38, 'P': 15.38, 'w': 7.69, 'W': 7.69, 'G': 7.69, 'd': 7.69}]
protein_tool('RRRrrrR', 'WGPdPA', 'LAlLAlw', options = 'percentage') # [{'R': 57.14, 'r': 42.86}, {'P': 33.33, 'W': 16.67, 'G': 16.67, 'd': 16.67, 'A': 16.67}, {'L': 28.57, 'A': 28.57, 'l': 28.57, 'w': 14.29}]

pattern Introduction The pattern procedure finds all non-overlaping cases of a given pattern in amino acid sequence(s) provided. Inputs To start using the pattern procedure, enter one or more protein sequences for which you want to get a summary, where the first sequence is a pattern, which is searched for in the following sequences; each argument should be of type 'str' and at the end add options = ‘pattern’. The find_pattern() function goes through a sequence in the following way: it takes a subsequence of amino acids in front of an index equal in length to the pattern and compares it to the pattern. If there is no match, index is moved one amino acid to the end of the sequence. If there is a match, the index is saved, and the function jumps to an aminoacid next to the end of the subsequence, then the algorithm repeats. Comparison is performed by compare_pattern subfunction.
The image explanation of that function.

Outputs The result of this procedure is a 'dict' object where:

• keys are amino acid sequences (type 'str')
• values are lists where the first element is a number of pattern instances in a given sequence, and the following elements are indexes of these occurances Usage example

protein_tool('LAlLAlwWGPdPA', 'LAlLAl', 'GPdPA', options = 'pattern') # {'LAlLAl': [2, 0, 3], 'GPdPA': [0]}
protein_tool('LAlLAlwWGPdPA', 'AlLAl', options = 'pattern') # {'AlLAl': [1, 2]}

3Letter_name Introduction The 3Letter_name procedure transforms one-letter amino acid entry sequences to three-letter amino acid sequences, separated by a specified separator. It is a case-sensitive procedure. Inputs To start using the rename_three_letter_name procedure, enter one or more protein sequences for which you want to get three-letter sequences. After the protein sequences put a symbol (type 'str') that will be a separator. And specify the options = ‘3Letter_name’. Outputs The result of the procedure is a list of three-letter sequences. Each amino acid is separated by the specified separator. The case of the three-letter amino acid coincides with the case of the one-letter designation at the input. Usage example

protein_tool('wWGPdPA', '', options = '3Letter_name') # ['trpTRPGLYPROaspPROALA']
protein_tool('LAlLAlwWGPdPA', '-', options = '3Letter_name') # ['LEU-ALA-leu-LEU-ALA-leu-trp-TRP-GLY-PRO-asp-PRO-ALA']
protein_tool('RRRrrrR', 'WGPdPA', 'LAlLAlw', options = 'percentage') # [{'R': 57.14, 'r': 42.86}, {'P': 33.33, 'W': 16.67, 'G': 16.67, 'd': 16.67, 'A': 16.67}, {'L': 28.57, 'A': 28.57, 'l': 28.57, 'w': 14.29}]
protein_tool('qwerty', 'G', options = '3Letter_name') # ['glnGtrpGgluGargGthrGtyr']

DNA_code Introduction The DNA_code procedure transforms a protein into a DNA sequence that may encode it (this can be used in genetic ingeneering). P.S. codons chosen at the discretion of the tool authors. Inputs To start using the DNA_code procedure, enter one or more protein sequences for which you want to get a summary, and at the end add options = ‘DNA_code’. Outputs The result of the procedure is a list with type 'str' elements - nucleotide sequence that corresponds to the amino acid sequence. Usage example

protein_tool('LAlLAlwWGPdPA', options = 'DNA_code') # ['TTAGCAttaTTAGCAttatggTGGGGGCCCgcaCCCGCA']
protein_tool('RRRrrrR', 'WGPdPA', 'LAlLAlw', options = 'DNA_code') # ['CGACGACGAcgacgacgaCGA', 'TGGGGGCCCgcaCCCGCA', 'TTAGCAttaTTAGCAttatgg']

fastq_tools.py

description of how the run_dna_rna_tools.py program works

This program contains the function fastq_tool. The fastq_tools function takes as input a dictionary of the fastaq data and sort it by the length, quality of sequencing of each nucleotides and GC% content. The result of the running of the tool is sorted dictionary of the dictionary from the input. 'gc_bounds' - GC composition interval (in percentage) for filtering (by default results (0, 100), i.e. all reads in the direction). If you pass one number as an argument, it is considered to be the upper limit. Examples: gc_bounds = (20, 80) - save only reads with GC content from 20 to 80%, gc_bounds = 44.4 - save reads with GC content less than 44.4%. 'length_bounds' - length of the interval for filtering, still gc_bounds, but by default it is (0, 2**32). 'quality_threshold' - threshold value of average read quality for the filter, default is 0 (phred33 scale). Reads contribute to quality for all nucleotides below the threshold are discarded.

⚠️ Attention: please be careful. currently that program for the input argument takes the namu of analising file only in fastq format. So your agrument shold be like 'example_fastq.fastq'. The and fot the output name write only the name of the file; the prefix .fatsq will be add automaticly.

Use example

fasta_filtering('example_fastq.fastq') # will be created the folder with the file inside example_fastq.fastq with filtered sequenses. 

bio_files_processor.py

description of how the run_dna_rna_tools.py program works That programm create from the file with multilined sequenses one line form of sequse.

Use example

convert_multiline_fasta_to_oneline('example_fastq.fastq') # will be created the folder with the file inside example_fastq.fastq with sequenses written in one line istead of multiline in input fils. 

⚠️ Attention: please be careful. currently that program for the input argument takes the namu of analising file only in fastq format. So your agrument shold be like 'example_fastq.fastq'. The and fot the output name write only the name of the file; the prefix .fatsq will be add automaticly.