2872 days ago
picard tools command to get some insert statistics
...lace #module load picard/2.0.1 java -Xmx16g -XX:PermSize=8g -jar $PICARD_HOME/picard.jar CollectInsertSizeMetrics MINIMUM_PCT=0 HISTOGRAM_FILE=.pdf INPUT= OUTPUT=.sorted.colle...1337 days ago
Download the genome from NCBI using bash script/command
...ent file name wget --input $f #cat $f cd .. done #Reference #ftp://ftp.ncbi.nlm.nih.gov/pub/factsheets/HowTo_Downloading_Genomic_Data.pdf2532 days ago
Perl script to read multi fasta sequence one by one
...E #perl rohanRun.pl seq.fa my $outfile='tmp.fa'; my $fastaSeq_ref = readfasta ("$ARGV[0]"); my %fastaS...file $outfile $!"; print OUT "$key\n$fastaSeq{$key}\n"; } sub readfasta { (my $file)=@_; m...2221 days ago
Perl script to check fastq reads qualities !
...env perl use strict; use warnings; sub readfq { my ($fh, $aux) = @_; @...) = (0, 0, 0); while (($name, $seq, $qual) = readfq(\*STDIN, \@aux)) { ++$n;...iate data while (my ($name, $seq, $qual) = readfq(\*STDIN, \@aux)) { if...2155 days ago
Perl subroutine for reading multifasta file !
sub readfasta { (my $file)=@_; my %sequence; my $header; my $temp_seq; #suppose fasta files contains multiple sequences; open (IN, "2145 days ago
Perl script to run SATSUMA in loop !
...r SATSUMA my $maxSize = 5000; my $resolution = 5000; my $dotsize = 1; my $cpu=40; my @ids; #Store the ids if ($ARGV[1] eq "ids") { my $idFile="palindrome_ids.txt"; ope...2141 days ago
Installing Busco version 4.0.6
...fortran_impl_linux-64-7.3.0| hdf63c60_5 9.4 MB conda-...gxx_impl_linux-64-7.3.0 | hdf63c60_1 18.7 MB gx...x-64::gfortran_impl_linux-64-7.3.0-hdf63c60_5 gfortran_linux-64...da-forge::certifi-2020.4.5.1-py37hc8dfbb8_0 Proceed ([y]/n)? y1453 days ago
2076 days ago
2075 days ago