Results for "5"

Bio-Scripts

  • Coverage / Depth of reads !

    ...$4>=X' | wc -l 32876 # get length of reference genome bowtie2-inspect -s refgenome | awk '{ FS = "\t" } ; BEGIN{L=0}; {L=L+$3}; END{print L}' 45678

    2266 days ago

  • Perl script to check fastq reads qualities !

    #!/usr/bin/env perl use strict; use warnings; sub readfq { my ($fh, $aux) = @_; @$aux = [undef, 0] if (!defined(@$aux)); return if ($aux->[1]); if (!d...

    2210 days ago

  • BloomFilter

    ...867 (delta 7), reused 12 (delta 4), pack-reused 852 Receiving objects: 100% (86...KiB | 0 bytes/s, done. Resolving deltas: 100% (512/512), done. Checking conne...➜ bloomfilter git:(master) swig -Wall -c++ -perl5 BloomFilter.i Unable to find...

    2207 days ago

  • Install Perl Locally !

    ...bfolder named src, and untar: curl --create-dirs -L -o src/perl-5.24.1.tar.gz http://www.cpan.org/src/5.0/perl-5.24.1.tar.gz cd src...L::Simple Test::XPath IO::String Bio::Perl version # Set $PERL5LIB to find these libraries, a...

    2200 days ago

  • Perl script to count number of Ns in a multifasta file !

    ...countN.pl out.fasta scaffold1_size15575755 15575755 71824 0.0046326...0.001605866729631 scaffold10_size2523421 2523421 142217 0.0597248...19159 scaffold27_size971184 971184 55725 0.0608711040035654 scaff...2100447444 scaffold71_size97059 97059 9770 0.11192704693604 scaff...

    2200 days ago

  • Perl script to find the distance beetween all the contigs and scaffolds

    #!/usr/bin/perl use strict; use warnings; use Bio::SeqIO; $| = 1; #Script to see the distance beetween all the contigs and scaffolds #Usage: perl clustalRe...

    2196 days ago

  • Perl script to run SATSUMA in loop !

    ...enamedAdinetaV2.fa"; #Ur target file my $satsumaLoc="/home/urbe/Tools/SATSUMA/satsuma-code-0"; # Location of ur SATSUMA my $maxSize = 5000; my $resolution = 5000; my $dotsize = 1; my $cp...

    2196 days ago

  • Long reads mapper bash script !

    ...lr-0.2.3/ngmlr lamsaLoc=/home/urbe/Tools/LAMSA/lamsa #Parameters accepted toolName=$1 refFasta=$2 longReads=$3 thread=$4 readsType=$5 if [ $# -lt 5 ]; then echo "No or less...

    2179 days ago

  • Running Trinity on RNA-seq !

    ...-k5,5 -T . -S 125G right.fa.K25.stats > right.fa.K25.stats.so....6/trinity_out_dir/inchworm.K25.L25.DS.fa 100 10 > /home/urbe...dir/read_partitions/Fb_0/CBin_53 Monday, July 9, 2018: 11:17...r/read_partitions/Fb_0/CBin_215 Monday, July 9, 2018: 11:18:...

    2183 days ago

  • Fastq-dump for SRA download

    ...-L|--log-level Logging level as number or enum string One of (fatal|sys|int|err|warn|info) or (0-5)...

    2182 days ago