Python script to download covid genome !
#!/usr/bin/env python3 # these are the publicly available "complete" sequences #...= x['accession'] print("downloading", nm) dna = Entrez.efetch(db='nucleotid...').read().split("\n")[1:] allseq[nm] = ''.join(dna) im...1127 days ago
2908 days ago
Count the frequency of base G in a given DNA sequence
#!/usr/bin/perl use strict; use warnings; my $DNA = "GATTACACAT"; #initialize $countG and $cur...y $currentPos = 0; #calculate the length of $DNA my $DNAlength = length($DNA)...f 'yes' increment $countG while($currentPos < $DNAlength)...2908 days ago
2908 days ago
Perl script to Mutate a DNA Sequence
#!/usr/local/bin/perl -w # This script randomly mutates the DNA sequence and generates 10 success...it asks for the file name of the DNA sequence. # If the DNA seque...print "ENTER THE FILENAME OF THE DNA SEQUENCE:= "; $dnafilename =...nt = mutate($DNA); print "Mutate DNA\n\n";...2907 days ago
Perl script to count the number of Adenine, Thymine, Guanine and Cytosine in your DNA Sequence
#!/usr/local/bin/perl -w # While executing this script it asks for the file name of the DNA sequence. If the sequence file...h the path. In windows: c:\dnafile.txt, In Linux: /home/us...AME OF THE DNA SEQUENCE:= "; $dna_filename = ; chomp $dna_file...rors = 0; foreach $base (@DNA) {...2907 days ago
2901 days ago
Implementation of biological random mutation with Perl
#!/usr/bin/perl -w use strict; use warnings; #sequence for a better recognition my $DNA="AAAAAAAAAAAAAAAAAAAAAAAAAAAA...rand (time | $$). $mutant=mutate ($DNA); print "Mutate \ n". $DNA...NA \ n"; print "Here is the mutant DNA: \ n \...2892 days ago
Perl script to generate a random psuedo DNA sequence !
#!/usr/bin/perl print "Enter a number of nucleotides: \n"; chomp ($N = ); @b=qw/A T G C/;print ">Genome\n";while($l2874 days ago
Perl script introduces control structures, arrays and hashes.
#!/usr/bin/env perl use strict; use warnings; my @first_array = ('DNA', 'ATGCGTGC', 5, 'RNA', 'AUGC'); print $first_array[0], "\n\n";...'Number of seqs' => 2 ); print $sequence{'DNA'}, "\n...2853 days ago