Python script to download covid genome !
...Entrez.efetch(db='nucleotide',id=nm, rettype = 'fasta', retmode= 'text').read().split("\n")[1:] allseq[nm] = ''.join(dna) import json with open("data/allseq.json", "w") as f...1119 days ago
2899 days ago
Parse a genbank file using regular expressions
#! /usr/local/bin/perl -w $genbank = "genbank_file.txt"; open (GENBANK, $genbank) || die "cannot open $gb_report for reading: $!"; # Flag for multiline translation; 1 means translation "in prog...2899 days ago
Perl script to Mutate a DNA Sequence
...use File::Path; print "ENTER THE FILENAME OF THE DNA SEQUENCE:= "; $dnafilename = ; chomp $dnafilename; unless ( open(DNAFILE, $dnafilename) ) { print "Cannot open file \"$dnafilename\"\n\n";...2898 days ago
Perl script to count the number of Adenine, Thymine, Guanine and Cytosine in your DNA Sequence
...THE DNA SEQUENCE:= "; $dna_filename = ; chomp $dna_filename; unless ( open(DNAFILE, $dna_filename) ) { print "Sorry the file does not exist!!! \n"; print "Cannot open file \"$dna_filename\"\n";...2898 days ago
Retrieve NCBI GenBank records with a range of accession numbers
...if ( $retmax > $count ) { $retmax = $count; } open( my $OUTFILE, ">" . $param{outputFile} ) or die("Error: Cannot open $param{outputFile} : $!");...2898 days ago
Extract a random sequence from a file
#!/usr/local/bin/perl -w use strict; use warnings; use autodie; use List::Util qw/ shuffle /; my $outputfile = 'randomoutput.txt'; open my $in_fh, '', $outputfile;...2897 days ago
Blast result parser with Perl and Bioperl
...s,$outfile) = @ARGV; print "Parsing the BLAST result ..."; my $in = Bio::SearchIO->new(-format => 'blast', -file => $infile); open (OUT,">$outfile") or die "Cannot open $outfile: $!"; # print the...2897 days ago
Perl program to implement sliding window !
#!/usr/bin/perl -w my $filename = 'data.txt'; open(my TR, '2897 days ago
2897 days ago