Results for "pyrosequencing data"
Bio-Scripts
Estimate Genome Size with Jellyfish and R
jellyfish count -t 8 -C -m 19 -s 5G -o 19mer_out --min-qual-char=? /common/Tutorial/Genome_e...ore, where "?" = 30 jellyfish histo -o 19mer_out.histo 19mer_out #Plot dataframe19...2245 days ago
Perl script to find palindromic regions in DNA sequences
use strict; use warnings; my $pp = qr/(?: (\w) (?1) \g{-1} | \w? )/ix; my $filename = $ARGV[0]; open(my $fh, '2154 days ago
Perl script to check fastq reads qualities !
#!/usr/bin/env perl use strict; use warnings; sub readfq { my ($fh, $aux) = @_; @...qlen), "\n"; __END__ my @aux = undef; # this is for keeping intermediate data while (my ($name, $seq, $q...2154 days ago
Installing Busco version 4.0.6
(base) jit@jit-HP-Pro-3335-MT: conda install -c bi...busco=4.0.6 Collecting package metadata (repodata.json): done Solving environm...pan-meta-yaml-0.018-pl526_0 perl-data-dumper bioconda/linux-64::p...bioconda/linux-64::perl-module-metadata-1.000036-pl526_0 perl-para...1452 days ago
- ➜ Tools git:(master) ✗ git clone git://github.com/bcgsc/bloomfil...xec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving direct...xec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving direct...
2151 days ago
Biological Sequence handling with Perl !
package Sequence::Generic; # File: Sequence/Generic.pm use st...} else { croak "Doesn't look like sequence data" unless $sequence=~/^[gactnu\s]+$/i; $self->{'data'} = $self->_canonicalize($seq...2181 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 her...2140 days ago
- ➜ bin git:(master) ✗ ./fastq-dump --help Usage:...plitting Sequence data may be used in raw form or...Filter used in early 1000 Genomes data: no...Output to stdout, all split data become...
2126 days ago
Install Parrot Virtual Machine !
#Parrot is a virtual machine designed to efficiently compile and execute bytecode for dynami...other build files...........done. gen::config_pm - Record configuration data for later retrieval......done...1537 days ago
- # list sample names & download URLs. sample_links = {"ERR458493"...ES=sample_links.keys() # download yeast rna-seq data from Schurch et al, 2016 stud...nload_all: input: expand("rnaseq/raw_data/{sample}.fq.gz", sample=SAMPL...
1516 days ago
