Mercurial > repos > petr-novak > re_utils
view fastq_name_affixer.xml @ 20:5a05925340b0 draft
Uploaded
author | petr-novak |
---|---|
date | Mon, 07 Jun 2021 08:46:07 +0000 |
parents | c2c69c6090f0 |
children | 58807b35777a |
line wrap: on
line source
<tool id="names_affixer" name="FASTQ Read name affixer" version="1.0.0"> <description> Tool appending suffix and prefix to sequences names </description> <command interpreter="python"> ${__tool_directory__}/name_affixer.py -f $input -p "$prefix" -s "$suffix" -n $nspace > $output </command> <inputs> <param format="fastq" type="data" name="input" label="Choose your FASTQ file" /> <param name="prefix" type="text" size="10" value="" label="Prefix" help="Enter prefix which will be added to all sequences names" /> <param name="suffix" type="text" size="10" value="" label="Suffix" help="Enter suffix which will be added to all sequences names"/> <param name="nspace" type="integer" size="10" value="0" min="0" max="1000" label="Number of spaces in sequence name to ignore" help="Sequence name is a string before the first space. If you want name to include spaces in name, enter positive integer. All other characters beyond ignored spaces are omitted"/> </inputs> <outputs> <data format="fastq" name="output" label="FASTQ dataset ${input.hid} with modified sequence names" /> </outputs> <help> **What is does** Tool for appending prefix and suffix to sequences names in fastq formated sequences. **Example** The following Solexa-FASTQ file: :: @CSHL_4_FC042GAMMII_2_1_517_596 GGTCAATGATGAGTTGGCACTGTAGGCACCATCAAT +CSHL_4_FC042GAMMII_2_1_517_596 40 40 40 40 40 40 40 40 40 40 38 40 40 40 40 40 14 40 40 40 40 40 36 40 13 14 24 24 9 24 9 40 10 10 15 40 is renamed to: :: @prefixCSHL_4_FC042GAMMII_2_1_517_596suffix GGTCAATGATGAGTTGGCACTGTAGGCACCATCAAT +prefixCSHL_4_FC042GAMMII_2_1_517_596suffix 40 40 40 40 40 40 40 40 40 40 38 40 40 40 40 40 14 40 40 40 40 40 36 40 13 14 24 24 9 24 9 40 10 10 15 40 different format: :: @HISEQ1:92:c0190acxx:8:1101:1252:2230 2:N:0:CGATGT AGAGGAAAAAACATAGTTCTTGTCTAAAAAAATCCCTTGAAAAAGGGCAGATGTATAGAAATAGAAAATTTCAAAGAAAAACTCTCTACAAATGGAAGAGA + CCCFFFFFHHHHHJJJJIJJJJJJJJJJJJJJJIJJJJJIIJJJJJJGIJIJIHHHHHHHHFFFFFFDEEEEEDCDDDDDDDCCDDDEDDDDD>CCCCB@9 is renamed to: :: @prefixHISEQ1:92:c0190acxx:8:1101:1252:2230suffix AGAGGAAAAAACATAGTTCTTGTCTAAAAAAATCCCTTGAAAAAGGGCAGATGTATAGAAATAGAAAATTTCAAAGAAAAACTCTCTACAAATGGAAGAGA + CCCFFFFFHHHHHJJJJIJJJJJJJJJJJJJJJIJJJJJIIJJJJJJGIJIJIHHHHHHHHFFFFFFDEEEEEDCDDDDDDDCCDDDEDDDDD>CCCCB@9 note that string after first space is omitted! Because sequence names sometimes containg spaces which delimit the actual name. By default, anything after spaces is excluded from sequences name. In example sequence: :: @SRR352150.23846180 HWUSI-EAS1786:7:119:15910:19280/1 CTGGATTCTATACCTTTGGCAACTACTTCTTGGTTGATCAGGAAATTAACACTAGTAGTTTAGGCAATTTGGAATGGTGCCAAAGATGTATAGAACTTTC + IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIGIIIHIIIIIFIIIIIIHDHBBIHFIHIIBHHDDHIFHIHIIIHIHGGDFDEI@EGEGFGFEFB@ECG when **Number of spaces in name to ignore** is set to 0 (default) the output will be: :: @prefixSRR352150.23846180suffix CTGGATTCTATACCTTTGGCAACTACTTCTTGGTTGATCAGGAAATTAACACTAGTAGTTTAGGCAATTTGGAATGGTGCCAAAGATGTATAGAACTTTC + IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIGIIIHIIIIIFIIIIIIHDHBBIHFIHIIBHHDDHIFHIHIIIHIHGGDFDEI@EGEGFGFEFB@ECG If you want to keep spaces the setting **Number of spaces in name to ignore** to 1 will yield :: @prefixSRR352150.23846180 HWUSI-EAS1786:7:119:15910:19280/1suffix CTGGATTCTATACCTTTGGCAACTACTTCTTGGTTGATCAGGAAATTAACACTAGTAGTTTAGGCAATTTGGAATGGTGCCAAAGATGTATAGAACTTTC + IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIGIIIHIIIIIFIIIIIIHDHBBIHFIHIIBHHDDHIFHIHIIIHIHGGDFDEI@EGEGFGFEFB@ECG </help> </tool>