To convert just one specific read group to fastq
...ta for each cat selected.txt | parallel fastq-dump -O reads -X 1000 --split-files {} # Create a direct...eader:" samtools view -H all.bam echo "" echo "Number of alignments with read group: GROUP-SRR1972919...1603 days ago
1602 days ago
Install python-numpy python-scipy python-matplotlib on Ubuntu !
...ython-tz tk8.6-blt2.5 ttf-bitstream-vera Suggested package...-sip python-tornado python-traits python-wxgtk3.0 texlive-e...tu bionic/universe amd64 ttf-bitstream-vera all 1.10-8 [352 kB...ously unselected package ttf-bitstream-vera. Preparing to unp...1602 days ago
1600 days ago
Bash oneliner to extract all ids from a multifasta file
#List of ids - one per line in allIds.txt $ awk 'BEGIN{while((getline0)l[">"$1]=1}/^>/{f=!l[$1]}f' seq.fa # You can play with this f=!l[$1 ] if wanted to extract or not extract the ids1599 days ago
Install Parrot Virtual Machine !
...ode for dynamic languages jit@jit-HP-Pro-3335-MT:~/Downloads/pa...e pulling some defaults from its configuration. do "lib/Parr...e byteorder for wordsize.....little-endian. auto::va_ptr -...`make test' to run the test suite. Happy Hacking,...1598 days ago
1586 days ago
1586 days ago
1585 days ago
Exit Bash Script When Any Command Fails
This can actually be done with a single line using the set builtin command with the -e option. #1--------...------------------------ # exit when any command fails #Putt...commands return a non-zero exit code. set -e #all other...cho an error message before exiting trap 'echo "\"${last_comm...1581 days ago