# HG changeset patch # User tyty # Date 1419015530 18000 # Node ID 361dc1047dc0040eb1fe7aaccf9313cf5fdb75c1 # Parent 40b8071a44f7318d41c8859f0fe7da5334c14152 Uploaded diff -r 40b8071a44f7 -r 361dc1047dc0 get_reads/get_read.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/get_reads/get_read.xml Fri Dec 19 13:58:50 2014 -0500 @@ -0,0 +1,46 @@ + + derives the reverse transcriptase (RT) stop count on each nucleotide from a mapped file provided by the Iterative Mapping module + get_read.py $lib_file $map_file $output + + biopython + numpy + samtools + + + + + + + + + + + + + + + + + + +**Function** + +Get RT Stop Counts derives the RT stop counts on each nucleotide of each transcript in the reference transcriptome based on a mapped file (.bam), typically the output from the Iterative Mapping module. + +----- + +**Input**: + +* 1. A mapped (.bam) file from Bowtie (or any other mapping program) +* 2. Reference library sequences (fasta) used to map the reads to + +----- + +**Output**: + +A text file with reverse transcription stop counts mapped to each nucleotide (RTSC file) + + + + +