2072 days ago
Perl script to reverse complement a DNA sequence !
#!/usr/bin/perl -w $DNA = 'ACGGGAGGACGGGAAAATTACTACGGCATTAGC'; print "Here is the starting DNA:\n\n";...A; # See the text for a discussion of tr/// $revcom =~ tr/ACGTacgt/TGCAtgca/; print "Here is the...2040 days ago
Installing pb-assembly on Linux !
...M4/1.4.17-GCCcore-5.4.0 bzip2/1.0.6-foss-...NLopt/2.4.2-foss-2016b fftw2/gcc/64/double/2.1.5...LAPACK-3.6.0 help2man/1.47.4-GCCcore-6.4.0...hfc679d8_5 conda-forge libgcc 7.2.0...2031 days ago
2016 days ago
Installing Platypus on Ubuntu !
...HAVE_DRAND48 1' >> config.h gcc -g -Wall -O2 -I. -c -o kfunc.o kfunc.c gcc -g -Wall -O2 -I. -c -o knet...al.pico cram/cram_external.c gcc -g -Wall -O2 -I. -fpic -c -...S=64 -g -Wno-unused-function gcc -pthread -shared -B /home/ur...ng 'variantcaller' extension gcc -pthread -B /home/urbe/anaco...2016 days ago
CollectGcBiasMetrics.jar will generate a GC bias plot for each contig
...ted.bam $i > aln-pe.mapped.sorted.$i.bam; java -Xmx2g -jar $(which CollectGcBiasMetrics.jar) R=data/Cdiff078.fa I=aln-pe.mapped.sorted.$i.bam O=aln-pe.mapped.sorted.${i}_GCBias.txt CHART=aln-pe.mapped.s...2001 days ago
Setting up falconUnzip conda environments for genome assembly !
...libffi: 3.2.1-hfc679d8_5 conda-forge libgcc: 7.2.0-h69d50b8_2 conda-forge libgcc-ng: 7.2.0-hdf6...1987 days ago
Finding Kmers from fasta sequence file
Save it in sample.fa >test TAATGCCATGGGATGTT jellyfish count -m 3 -s 100000 sample.fa -o sample.jf jellyfish dump -c sample.jf It return TGT 1 GAT 1 GGG 1 GGA 1 CAT 1 TGC 1 TAA 1 GCC 1 CCA 1 GTT 1 TGG 1 ATG 3 AAT 11969 days ago
Perl script to run in parellel !
...y $sequence = $gene_info->seq(); my $accession_number = $gene_info->display_id; my $len = $gene_info->length; my $GCcount = $sequence =~ tr/GC|gc//; my $GCcontent =...1683 days ago
1556 days ago