Results for "Lateral gene transfer PLOS"
Bio-Scripts
Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/...len(seqs)) from Bio import Entrez allseq = {} for x in seqs: if 'gene-region' in x and x['gene...1128 days ago
Parse a genbank file using regular expressions
#! /usr/local/bin/perl -w $genbank = "genbank_file.txt"; open (GENBANK, $genbank) ||....*)/) { print "Organism: $2\n"; } elsif(/(gene)(\s*)(\d*)(\.\.)(\d*)/) { print "Gene...2908 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::E...2907 days ago
- #!usr/bin/perl use strict; use warnings; my %id2seq = (); my $id = 0;...; my @val = split /\s+/, $_; my @newval=@val; shift @newval; print ">$id [ gene...
2504 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, load it...library(ggplot2) } # import a text file with gene pos...be: chr, position (no end or gene name required) gene...2276 days ago
Plot custom gene density with R
library(karyoploteR) pp2241 days ago
Perl script to convert GFF 2 FASTA !
#!/usr/bin/perl use strict; use warnings; use Bio::Seq; use Bio::SeqIO;...=> ">$ARGV[2].cdna.fasta" ); my $outfile_gene = Bio::SeqIO->new( -format => 'fasta', -file => ">$ARGV[2].gene.fas...* as the stop codon) # gene - the entire gene...2130 days ago
Pack a perl program with their dependencies on Ubuntu !
#Follow steps to create your own executable ./web jit@jit-HP-...m highly‐expressed A. gambiae genes is appended below. Biote...s that about sixty percent of genes are conserved (Adams et al. 2000...Detecting Horizontal Gene Co-transfer between Bacteria Using Gene-g...MatesRpackageHorizontalGeneCo-tran...1503 days ago
Sequence Ids conversion files !
ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/ Name Size Date Modified ARCHIVE/...30:00 ASN_BINARY/ 03/07/2020, 07:49:00 GENE_INFO/ 03/07/2020, 07:48:00...1, 05:30:00 ftp://ftp.ncbi.nlm.nih.gov/gene/DAT...ene2pubmed.gz ftp://ftp.ncbi.nlm.nih.gov/gene...1394 days ago
Bash script to extract intronic fragments !
#To obtain introns, we simply need the gene and exonic coordinates; #by subtracting the exonic regi...the intronic region. gunzip -c genome_file.gtf.gz | awk 'BEGIN{OFS="\t";} $3=="gene...1357 days ago
