Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.or...en(seqs)) from Bio import Entrez allseq = {} for x in seqs: if 'gene-region' in x and x['gene-region...1134 days ago
Read lines from input file – print lines that match a regular expression
#!/usr/bin/perl use strict; use warnings; my $line; #read the line-by-line for each line ask if it matches the regex print it if it matches while($line = ){...2914 days ago
Parse a genbank file using regular expressions
#! /usr/local/bin/perl -w $genbank = "genbank_file.txt"; open (GENBANK, $genbank) || die "canno...ANISM\s*)(.*)/) { print "Organism: $2\n"; } elsif(/(gene)(\s*)(\d*)(\.\.)(\d*)/) { print "Gene length...2914 days ago
Retrieve NCBI GenBank records with a range of accession numbers
#!/usr/bin/perl #FILE: ncbi_search.pl #AUTH: Paul Stothard (paul.stothard@gmail.com) use warnings; use strict; use Getopt::Long; use LWP::Simple; use URI::Esca...2913 days ago
2510 days ago
Plot the density of genes in R
#column1 = chromosome name and column2 = start position of the gene # check if ggplot2 is installed, if so, lo...rary(ggplot2) } # import a text file with gene positions # columns should b...chr, position (no end or gene name required) genes...2282 days ago
Plot custom gene density with R
library(karyoploteR) pp2247 days ago
Perl script to find coding regions in DNA sequences
#!/usr/bin/perl -w use strict; # if the number of input arguments is lower than 2 # return a message showing the error if (scalar(@ARGV) < 2) { print "dnal...2153 days ago
Perl script to convert GFF 2 FASTA !
#!/usr/bin/perl use strict; use warnings; use Bio::Seq; use Bio::Seq...>$ARGV[2].cdna.fasta" ); my $outfile_gene = Bio::SeqIO->new( -format => 'fasta', -file => ">$ARGV[2].gene.fasta" ); my $outfile_upstre...the stop codon) # gene - the entire gene sequen...2136 days ago
2073 days ago