1992 days ago
Split the multifasta in separate files !
cat Avaga_allPalindrome.fa | awk '{ if (substr($0, 1, 1)==">") {filename=(substr($0,2) ".fa")} print $0 > filename }'1968 days ago
Backward Elimination with Python
...n(x, sl): numVars = len(x[0]) for i in range(0, numVars): regr...) if maxVar > sl: for j in range(0, numVars - i):...emp = np.zeros((50,6)).astype(int) for i in range(0, numVars): regr...1965 days ago
Perl script to count occurrence of a character !
#!/usr/bin/env perl # -*- coding: utf-8 -*- #!/usr/bin/perl use strict; use warnings; my %count_of; while ( ) { my @val = split "\t", $_; #my...1911 days ago
1811 days ago
Palindrome Simulation commands !
...nt153CH101.fasta > allPalindromeSimulated.fa 3319 ~/Tools/art_bin_MountRainier/art_illumina -ss MSv3 -s...-p -l 251 -f 100 -m 300 -s 10 -o paired_dat 3320 ~/Tools/art_bin_MountRainier/art_illumina -ss MSv3 -s...1746 days ago
Resume the MIRA assembler run !
mira -r manifest_file Usage: mira [options] manifest_file [manifest_file ...] Options: -c / --cwd= directory Change working directory -r / --resume...1745 days ago
1733 days ago
Perl script to run in parellel !
...sr/bin/perl use strict; use warnings; use Parallel::ForkManager; use Bio::S...e_data_ref}; my $n_processes = 4; my $pm = Parallel::ForkManager->new( $n_pro...= @_; print "starting XFOIL instance with parameters $chr_set!\n"; touc...1697 days ago
Samtools commands for bioinformatician !
...ds ### calling snp and indels with no frequency threshold # omit -v parameter bcftools call -c sal_va...ults into txt file samtools depth /path/to/sorted_bam.bam > /path/to/coverage_results.txt ### note: retu...1613 days ago