1109 days ago
1109 days ago
1109 days ago
1109 days ago
Commandline to Extract a list of specific read IDs from a bam file
#Save all the IDs in IDs.txt file #Run it on BAM file samtools view file.bam | fgrep -w -f IDs.txt1059 days ago
Perl onliner to print fasta headers !
#Save all your fasta in seq.fa and run the following ... perl -ne 'print if /^>/' seq.fa #Print header with line number perl -ne 'print "$. $_" if /^>/ ' seq.fa1059 days ago
1053 days ago
Commandline for paired end reads simulation with BBMap !
..., unique=true insRate=0.0, max=0, len=(0-0) delRate=0.0, max=0, len=(0-0) subRate=0.0, max=0, len=(0-0) nRate =0.0, max=0, len=(0-0) genome=1 PERFECT_READ_RATIO=0.0 ADD_ERRO...1053 days ago
1044 days ago
Print in terminal with python !
#!/usr/bin/env python import time import curses def pbar(window): height, width = window.getmaxyx() for i in range(10): window.addstr(height -1, 0, "[" + ("=" * i) + ">" + (" " * (10 - i )) + "]") window.refresh() time.sleep(0.5) curses.wrapper(pbar)1041 days ago