Results for "genome"
Bio-Scripts
Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/ has more (1200?), but they require you to sign up import...1125 days ago
Download genome using NCBI esearch command line !
#Download using command line "NC_005956.1" / "CP003784.1" esearch -db nucleotide -query "NC_005956.1" | efetch -format fasta > seq1.fa esearch -db nucleotide -quer...1091 days ago
R script to visualize fastANI core-genome comparison
####### # Purpose: Visualize fastANI core-genome comparison # Usage: Rscript # Output: .pdf # Uses genoPlotR package: http://genoplotr.r-forge.r-project.org #Parse command...1090 days ago
Extract a random sequence from a file
...utputfile; my $size = 21; my $count = 10; while (my $line = ) { next unless $line =~ /^([ATGCN]+)/; my $genome = $1; my $len_genome = length $genome; my @s...2904 days ago
- ...+ ($size/2); # Figure out URL for the DAS server. Example: # http://genome.ucsc.edu/cgi-bin/das/calJac3/dna?segment=chr14:482280,483280 my $URL_gene ="http://genome.ucsc.edu/cgi-bin/das/papAnu2/...
2899 days ago
Perl script to generate a random psuedo DNA sequence !
#!/usr/bin/perl print "Enter a number of nucleotides: \n"; chomp ($N = ); @b=qw/A T G C/;print ">Genome\n";while($l2872 days ago
Blast script to index and extract sequence !!
...coli O157:H7 str. EC4115 chromosome, complete genome. AGCTTTTCATTCTGACTGCAACGGGCA...coli O157:H7 str. EC4115 chromosome, complete genome AGCTTTTCATTCTGACTGCAACGGGCAA...coli O157:H7 str. EC4115 chromosome, complete genome TTAAAA...2841 days ago
Extract a range from genome file with perl.
#!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; my $in_file = $ARGV[0]; my $start_pos = $ARGV[1]; my $end_pos = $ARGV[2]; my $in = Bio::SeqIO->n...2616 days ago
Perl script to insert the DNA string in genome
#!/usr/bin/perl use warnings; use strict; use Bio::SeqIO; use Bio::Seq; my $file = $ARGV[0]; # input fasta file (genome file) my $out = $ARGV[1]; # output fasta file my $...2648 days ago
BASH script for SelfBLAST a genome
#!/bin/bash #self BLAST a genome -- Expecting you have blast and samtools installed in your system #Aut...[ -f $MYDB.nhr ] then echo "BLAST database for MergedContigs.fasta genome exists" else echo "Thanks...2641 days ago
