Results for "RA"

Bio-Scripts

  • Sync github repo

    ...-local-repo-with-remote-one This makes your local repo exactly like your remote repo. Remember to replace origin and master with the remote and branch that you want to synchroni...

    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

  • Split array in Perl !

    my @tests = ( ["FOO", "BAR"], ["CRATE", "TRACE"], ["CRATE", "CRATE"], ["TRACE", "CRATE"], ["CREATE", "TRACT"], ["DWAYNE", "DUANE"], ); for my $word_pair (@tests) { #Nice way to split the...

    1810 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

  • Run sspace !

    #!/bin/bash cd `pwd` perl ~/apps/SSPACE-1.2_linux-x86_64/SSPACE_v1-2.pl \ -l libraries.txt \ -s Contigs_over200_nocp.fasta \ -k 5 \ -a 0.7 \ -x 1 \ -m 30 \ -o 20 \ -b Rayk31_scaffolds_extension

    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