2128 days ago
Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/ has...seqs = seqs['genbank-sequences'] print("got %d sequences" % len(seqs)) from Bio i...1120 days ago
BioPerl to convert between sequence formats from Fasta to Genbank
#!/usr/local/bin/perl -w # Sequence formats to choose: Fasta, EMBL. GenBank, Swissprot, PIR and GCG use Bio::SeqIO; $inFile = "BRCA2.fa"; $in = Bio::...2899 days ago
Check if your coputer ready to use BioPerl
#!/usr/bin/perl use strict; use warnings; #bioperl example code use strict; use warnings; #make the bioperl module (class) accessible to your program use Bio::Seq; print"ok - ready to use Bio::Seq";2900 days ago
Retrieve NCBI GenBank records with a range of accession numbers
#!/usr/bin/perl #FILE: ncbi_search.pl #AUTH: Paul Stothard (p...HTTP::Request::Common; my %param = ( query => undef,...( 'q|query=s' => \$param{query}, 'o|output_file=s..., 'h|help' => \$param{help} ); if ( defined( $p...t) ) { message( $param...2899 days ago
Blast result parser with Perl and Bioperl
#!/usr/local/bin/perl # # Dr. Xiaodong Bai # It may be freely distributed under G...ther analysis. # use strict; use warnings; use Bio::SearchIO; # Usage informa...GV; print "Parsing the BLAST result ..."; my $in = Bio::...2898 days ago
Find and replace ambiguous characters in fasta file with Perl and Bioperl
#!/usr/bin/perl -w my $usage="\nUsage: $0 [-h] [-m char] [fastaFileName1 ...]\n"....TDIN as the input\n"; our($opt_h, $opt_m); use Bio::SeqIO; use Getopt::Std;...GV; while (my $file = shift) { my $seqio_obj = Bio::...2898 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_po...= $ARGV[1]; my $end_pos = $ARGV[2]; my $in = Bio::SeqIO->new ( -file => $in_file, -format => 'fasta'); my $out = Bio::...2610 days ago
2688 days ago
Perl script to remove the duplicate sequences from multifasta file
use strict; use Bio::SeqIO; my %unique; my $file = "myseqs.fa"; my $seqio = Bio::SeqIO->new(-file => $file, -format => "fasta"); my $outseq = Bio::...2673 days ago