2678 days ago
Perl script to find coding regions in DNA sequences
...# return a message showing the error if (scalar(@ARGV) < 2) { print "dnaloglkh.pl codontable DNAsequence\n"; exit(1); } # create two vars contaning the filenames my...2140 days ago
Biological Sequence handling with Perl !
...: Sequence/Nucleotide.pm use Sequence::Generic; use Sequence::Nucleotide::Subsequence; use Sequence::Alignment; use Carp; use strict; use vars '@ISA'; :Generic'; my %CO...2166 days ago
2111 days ago
Samtools commands for bioinformatician !
...n samtools mpileup -g -f sal_ref_sej.fasta sal_sej_sorted.bam.bam > sal_vars.bcf ### calling snp and...find snp # -v : output only potential variants bcftools view -c -v sal_vars.bcf > sal_vars.vcf ### note:...1588 days ago