We have generalized the processing script to be able to handle multiple scenarios. Originally it was designed for PCR free data only, we added all the steps necessary to start from fastq and process RNA-seq as well as non-human data. This is our go to script in TechDev. * add optional "starting from fastq" path to the pipeline * add mark duplicates (optionally) to the pipeline * add an option to run with the mouse data (without dbsnp and with single ended fastq) * add option to process RNA-seq data from topHat (add RG and reassign mapping quality if necessary) * add option to filter or include reads with N in the cigar string * add parameter to allow keeping the intermediate files |
||
|---|---|---|
| .. | ||
| R | ||
| c | ||
| chainFiles | ||
| doc | ||
| java | ||
| keys | ||
| packages | ||
| perl | ||
| scala | ||
| testdata | ||