2033 days ago
Finding Kmers from fasta sequence file
Save it in sample.fa >test TAATGCCATGGGATGTT jellyfish count -m 3 -s 100000 sample.fa -o sample.jf jellyfish dump -c sample.jf It return TGT 1 GAT 1 GGG 1 GGA 1 CAT 1 TGC 1 TAA 1 GCC 1 CCA 1 GTT 1 TGG 1 ATG 3 AAT 12024 days ago
Perl script to split fasta sequence and create overlaps
#!/usr/bin/perl use strict; use warnings; my $len = 5000; my $over = 200; my $seq_id=$ARGV[0]; my $seqFile = $ARGV[1]; my $seq; open(my $fh, "2024 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 }'2009 days ago
Backward Elimination with Python
...ion(x, SL): numVars = len(x[0]) temp = np.zeros((50,6)).astype(int) for i in range(0, numVa...x_rollback = np.delete(x_rollback, j, 1) print (regressor_OLS.summary())...2007 days ago
Perl script to count occurrence of a character !
...in/perl use strict; use warnings; my %count_of; while ( ) { my @v...t "\t", $_; #my ( $word) = m/(\w+)/; $count_of{$val[13]}++; } foreach...{$a} $count_of{$b} } keys %count_of ) { print "$word\t$count_of{$word}\n";...1952 days ago
1813 days ago
Palindrome Simulation commands !
...more allCommands 3315 mutate.sh in= mutant15CH101.fasta out=mutant153CH101.fasta id=97 3316 m...101_read33155_template_pass_FAH31515.faa mutant3CH101.fasta mutant15CH101.fas...10 -o paired_dat 3320 ~/Tools/art_bin_MountRainier/art_illumina -ss MSv3...1787 days ago
Resume the MIRA assembler run !
...nge working directory -r / --resume Resume an interupted assembly -h / --help Print short help and exit -v / -...1786 days ago
1777 days ago