Results for "genome letters"
Bio-Scripts
#!/bin/bash #set the path export PATH=$PATH:/home/urbe/Tools/augustus-3.2.3/bin ex...anuCorrect_sspaceLong/GRAAL_Haploid/adineta_p_lvl6/genome_CORRECTED2.fasta -c 30 -l metazoa_odb9 -o GRAALHaploid_BUSCO -m genome...2138 days ago
Perl subroutine to read genome/reads fasta file !
sub readSeqFromFasta{ my $file = $_[0]; my (%map,$id,$seq); open(IN,$file) or die("Cannot open file for reading $file:$!\n"); while(){ chomp; if($_ =~ /^>...2134 days ago
- ➜ bin git:(master) ✗ ./fastq-dump --help Usage: ./fastq-dump [options] [...]...l-filter Filter used in early 1000 Genomes data: no...ter-1 Filter used in current 1000 Genomes...
2127 days ago
- jitendra@jitendra-UNLOCK-INSTALL[SNP] wget https://github.com/broadinstitute/gatk/relea...awWgsMetrics (Picard) Collect whole genome sequencing-related metrics....ct metrics about coverage and performance of whole genome...
2076 days ago
Generates a genome coverage plot with R
library(CoverageView) ##draw a coverage plot for a test case BAM file #get a BAM test file treatBAMfile2056 days ago
Generate simulated polyploid genome !
#Generate 3% divergence msbar -point 4 -count 16558 toy.fasta > toyheterozygous3percent.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheteroz...2054 days ago
Setting up falconUnzip conda environments for genome assembly !
➜ Analysis_Results conda create -n denovo_asm Solving environment: done ## Package Plan ## environment location: /home/urbe/anaconda3/envs/denovo_asm Pro...1991 days ago
- Map them agaist reference avaga genome using following codes git clone https://github.com/lh3/bwa.git cd bwa; make bwa index ref.fa bwa mem -x pacbio ref.fa pacbio.fq > aln.sam bwa mem -x ont2d ref.fa ont-2D.fq > aln.sam
1726 days ago
Perl script to run in parellel !
#!/usr/bin/perl use strict; use warnings; use Parallel::ForkManag...my ($sequence_data_ref) = parse_genome_files($ARGV[0]); my %genome=%{$sequence_data_ref}; my...; foreach my $chr_set (keys %genome) { $count++...e FILE or die $!; } sub parse_genome_...1687 days ago
Find and replace in multifasta or fasta header with perl onliner
You have a fasta file and you want to replace: "|" You are told to replace that by "_" perl -i -p -e "s/\|/_/g" genome.fasta -i = inplace editing -p = loop over lines and print each line (after processing) -e = command line script1616 days ago
