2143 days ago
Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences # https://www.gisaid.org/ has...eqs = seqs['genbank-sequences'] print("got %d sequences" % len(seqs)) from Bio impo...1135 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::Seq...2915 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";2916 days ago
Blast result parser with Perl and Bioperl
#!/usr/local/bin/perl # # Dr. Xiaodong Bai # It may be freely distributed under...her analysis. # use strict; use warnings; use Bio::SearchIO; # Usage informa...V; print "Parsing the BLAST result ..."; my $in = Bio::Sea...2913 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"....DIN as the input\n"; our($opt_h, $opt_m); use Bio::SeqIO; use Getopt::Std;...V; while (my $file = shift) { my $seqio_obj = Bio::Seq...2913 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::Seq...2626 days ago
2703 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::Seq...2689 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 f...; #sequence of the insertion my $seq_in = Bio::SeqIO->new( -format => 'fast...$pos,$length); my $seq_obj = Bio::Seq...2658 days ago