Results for "genome%20overlapping%20features"

Bio-Scripts

  • Estimate Genome Size with Jellyfish and R

    jellyfish count -t 8 -C -m 19 -s 5G -o 19mer_out --min-qual-char=? /common/Tutorial/Genome_estimation/sample_read_1.fastq /common/Tutorial/Genome_estimation/sample_read_2.fastq...

    2259 days ago

  • Genome Covered !

    zero=$(bedtools genomecov -ibam BAM -g hg38.fasta -bga | awk '$4==0 {bpCountZero+=($3-$2)} {print bpCountZero}' | tail -1) nonzero=$(bedtools genomecov -ibam BAM -g hg38.fasta -bga | awk '$4>0 {bpCountNonZero+=($3-$2)} {print bpCountNonZero}' | tail -1) percent=$(bc

    2223 days ago

  • Coverage / Depth of reads !

    # get total number of bases covered at MIN_COVERAGE_DEPTH or higher samtools mpileup mapping_result_sorted.bam | awk -v X="${MIN_COVERAGE_DEPTH}" '$4>=X' | wc -l 32876 # get l...

    2223 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 clustalReads.pl genome.f...

    2153 days ago

  • Perl script to run SATSUMA in loop !

    #!/usr/bin/perl -w use strict; use File::Temp qw(tempfile); # Usage perl 1by1.pl for SATSUMA analysis # User need to set the reference multifasta file name here my $seqfil...

    2153 days ago

  • Perl script to convert GFF 2 FASTA !

    #!/usr/bin/perl use strict; use warnings; use Bio::Seq; use Bio::SeqIO; use Bio::DB::Fasta; $| = 1; # Flush output my $outfile_cds = Bio::SeqIO->new( -format => 'fasta'...

    2152 days ago

  • Bash script to run Busco2 !

    #!/bin/bash #set the path export PATH=$PATH:/home/urbe/Tools/augustus-3.2.3/bin export PATH=$PATH:/home/urbe/Tools/augustus-3.2.3/scripts export AUGUSTUS_CONFIG_PATH=/home/...

    2150 days ago

  • Perl subroutine to read genome/reads fasta file !

    sub readSeqFromFasta{ my $file = $_[0]; my (%map,$id,$seq); open(IN,$file) or die("Cannot open file for reading $file:$!\n"); while(){ chomp; if($_ =~ /^>/){ if(d...

    2146 days ago

  • Fastq-dump for SRA download

    ➜ bin git:(master) ✗ ./fastq-dump --help Usage: ./fastq-dump [options] [...] ./fastq-dump [options] INPUT -A|--accession Replaces accession derived from i...

    2139 days ago

  • Downloading GATK !

    jitendra@jitendra-UNLOCK-INSTALL[SNP] wget https://github.com/broadinstitute/gatk/releases/download/4.0.2.1/gatk-4.0.2.1.zip --2018-08-30 11:47:20-- https://github.com/broadinsti...

    2088 days ago