Sortdep: Difference between revisions
Jump to navigation
Jump to search
No edit summary |
No edit summary |
||
| (One intermediate revision by the same user not shown) | |||
| Line 5: | Line 5: | ||
The functionality is similar to '''sort -n depths.txt |uniq -c'''. | The functionality is similar to '''sort -n depths.txt |uniq -c'''. | ||
Problems with unix sort is that it floods the | Problems with unix sort is that it floods the ''$TMPDIR''. This utility program keeps the data in memory, | ||
=Installation= | =Installation= | ||
| Line 30: | Line 30: | ||
#or | #or | ||
cat depths.txt |cut -f3 |sed 1d|./sortDepth >depths | cat depths.txt |cut -f3 |sed 1d|./sortDepth >depths | ||
</pre> | |||
;Input as filename can be gz compressed, input as stdin are required to be uncompressed. | |||
Latest revision as of 14:45, 20 March 2014
A ultrafast and efficient tabulation of integers.
Very usefull for finding coverage and depth of coverage of sequencing data.
The functionality is similar to sort -n depths.txt |uniq -c.
Problems with unix sort is that it floods the $TMPDIR. This utility program keeps the data in memory,
Installation
wget http://popgen.dk/software/download/sortDep.cpp g++ sortDep.cpp -O3 -o sortDep -Liz
File below called depths.txt
chr pos depth 1 1 5 1 2 10 1 3 1 ....
./sortDep -skip 1 -col 3 depths.txt >depths #or cat depths.txt |./sortDepth -skip 1 -col3 >depths #or cat depths.txt |cut -f3 |sed 1d|./sortDepth >depths
- Input as filename can be gz compressed, input as stdin are required to be uncompressed.