Results for "Utilities #Tools #Scripts #Pierre"
Bio-Scripts
(JitMetaENV) ➜ JitMeta conda install grabseqs -c louiejtaylor -c bioconda -c conda-forge Collecting package metadata (current_repodata.json): done Solving environment: done ==> WARNING: A newer version of conda exists.1211 days ago
- # To follow these instructions, we'll assume your system admins have already installed these essentials: # Debian/Ubuntu system admins should run: sudo apt-get instal...
2138 days ago
Palindrome Simulation commands !
(base) ➜ palindromeAssemblySim more allCommands 3315 mutate.sh in= mutant15CH101.fasta out=mutant153CH101.fasta id=97 3316 mutate.sh in=mutant15CH101.fasta out=...1729 days ago
Read a tab delimited file and search with perl
use strict; use warnings; use Data::Dumper; use Text::CSV; use IO::Handle; my $file = "/home/urbe/Tools/Alienomics_v0.1/Alienomics/output/intermediate_files/rRNA/refGene.megablast"; open my $fh, "[0]\n"; warn Dumper $row; # To see the structure }2519 days ago
Loop over with all files in a directory in bash
#!/bin/bash FILES=/media/ComparativeGenomics/ncbi-genomes-2017-11-13/* ref=/media/ComparativeGenomics/ncbi-genomes-2017-11-13/GCA_000196735.1_ASM19673v1_genomic.fna...2325 days ago
- # Count k-mer occurrence using Jellyfish 2.2.6 jellyfish count -t 8 -C -m 19 -s 5G -o 19mer_out --min-qual-char=? sread_1.fastq sread_2.fastq # points for a histogram jellyfish histo -o 19mer_out.histo 19mer_out #Plot results using R ##load the data into dataframe19 dataframe19
2258 days ago
Perl script to convert fastq to fasta file
#!/usr/bin/env perl use strict; use warnings; use Bio::Factory::EMBOSS; my $usage = "perl $0 in.fq out.fa"; my $infile = shift or die $usage; my $outfile =...2245 days ago
Mapping with BWA-mem or BWA-sampe in one go with python script !
BAM files and mapping BESST requires sorted and indexed BAM files as input. Any read aligner + samtools can be used to obtain such files. Read pairs needs to be aligne...2118 days ago
- ➜ Tools git:(master) ✗ git clone git://github.com/bcgsc/bloomfilter.git Cloning into 'bloomfilter'... remote: Counting objects: 867, done. remote: Compressing objec...
2145 days ago
Perl script to run SATSUMA in loop !
#!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here...2134 days ago
