# HG changeset patch # User tyty # Date 1419013326 18000 # Node ID 6b7429975fdefb7a41138456a3da2f541ee44951 # Parent c3092769835edfad550de18c76890ab2dc84656f Uploaded diff -r c3092769835e -r 6b7429975fde 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:22:06 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** + +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) + + + + +