2761 days ago
2750 days ago
2750 days ago
Perl script to remove the duplicate sequences from multifasta file
...while(my $seqs = $seqio->next_seq) { my $id = $seqs->display_id; my $seq = $seqs->seq; unless(exists($unique{$seq})) { $outseq->write_seq($seqs); $unique{$se...2747 days ago
Perl script to insert the DNA string in genome
...put fasta file my $chr="test"; #insertion chromosome my $pos=10; # position of the insertion my $seqI...lay_id => $seq->primary_id,-alphabet => "dna" ); $seq_out->write_seq($seq_obj); } e...2716 days ago
Calculate ATGC percentage in parallel with perl
...# Setup a callback for when a child finishes up so we can # get it's exit code $pm->run_on_finish (...$pm->finish($child); # pass an exit code to finish } print "Waiting for Children...\n"; $pm...2713 days ago
BASH script for SelfBLAST a genome
...Expecting you have blast and samtools installed in your system #Author: Jitendra Narayan #USAGE: ./selfB...t" SEQ=$FASTAFILE else echo "Something went wrong $USER - Contact jitendra" fi echo "Doing alig...2709 days ago
Extract a range from genome file with perl.
...-format => 'fasta'); while (my $seq = $in->next_seq() ) { $seq->display_id( $seq->display_id() . "_$start_pos-$end_pos" ); $out->write_seq( $seq->trunc($start_pos,...2684 days ago
Calculate some statistics for a DNA alignment with Perl
...-method => 'Jukes-Cantor'); print $jcmatrix->print_matrix; ## and for measurements of synonymous /nonsynonymous substitutions ## my $in = Bio::A...2706 days ago
2693 days ago