Python script to download covid genome !
...nm = x['accession'] print("downloading", nm) dna = Entrez.efetch(db='nucleotide',id=nm, rettype = 'fasta', retmode= 'text').read().split("\n")[1:] allseq[nm] = '...1154 days ago
Perl script to count the number of Adenine, Thymine, Guanine and Cytosine in your DNA Sequence
...me\"\n"; die; } @DNA = ; close DNAFILE; $DNA = join( '', @DNA); print " \n The original DNA file is:\n $DNA \n"; $DNA =~ s/\s//g; @DNA = split( '', $DNA ); $count_of_A = 0...2934 days ago
Perl program to implement sliding window !
#!/usr/bin/perl -w my $filename = 'data.txt'; open(my TR, '2932 days ago
Perl to print indivisual nucleotide from a sequence!
#!/usr/bin/perl use strict; use warnings; my $string = "ATGCTTGCGT?AAATG??CT?GCGTA"; my @chars = split("", $string); print "First character: $chars[0]\n";2932 days ago
Needleman-Wunsch Algorithm in Perl
#!/usr/bin/perl # USAGE: perl nw.pl HEAGAWGHEE PAWHEAE BLOSUM50.txt -8 # See: "Biological sequence anaysis" Durbin et al. ed. CUP 1998, Pg. 19 # Needleman...2924 days ago
Count GC Content in nucleotide sequence with Perl
...#finish up last line. &process_it; close(IN); close(OUT); sub usage { $0 get_gc_content.pl } sub process_it { @letters = split(//, $seq); $gccount = 0;...2924 days ago
Perl script to extract lines with matching ids !!
#!/usr/bin/perl use strict; use warnings; my %patterns; #USAGE: perl extactByIds.pl Idsfile1 file2 > Result # Open file and get patterns to search for open(my $fh2,"2893 days ago
Perl script to extract fasta sequence by matching name/ids !!
...while () { chomp; next if /^\s*$/; s/>//g; my @ids=split (/\t/, $_); $select{$ids...$fasta" or die; while () { trim($_); s/>//g; my ($id) = split (/\n/, $_); #my @i=split...2893 days ago
2890 days ago
2711 days ago