Results for "No"

Bio-Scripts

  • Installing nupack3.0.6

    ...~~~ init.c:795:5: warning: ignoring return value of ‘fgets’,...~~ init.c:1077:5: warning: ignoring return value of ‘fgets’,...variable ‘pair_flag’ set but not used [-Wunused-but-set-varia...’: subopt.c:74:5: warning: ignoring return value of ‘scanf’,...

    2081 days ago

  • Install HTSLIB on Ubuntu !

    ...b-1.9/test/c1#bounds.sam htslib-1.9/test/c1#clip.sam htslib-1.9/test/c1#noseq.sam htslib-1.9/test/c1#pa...-1.9/test/c1#pad2.sam htslib-1.9/test/c1#pad3.sam htslib-1.9/test/c1#unknown.sam htslib-1.9/test/c1.fa...

    2074 days ago

  • Installing Platypus on Ubuntu !

    ...... Processing triggers for gnome-menus (3.13.3-11ubuntu1.1)...ocal/lib/pkgconfig mkdir: cannot create directory ‘/usr/local...E -D_FILE_OFFSET_BITS=64 -g -Wno-unused-function creating bui...tabix/index.c.pysam.c:37:1: note: in expansion of macro ‘KHA...

    2074 days ago

  • Collecting arguments with R

    #! /usr/bin/Rscript ## Collect arguments args

    2059 days ago

  • Setting up falconUnzip conda environments for genome assembly !

    ...25 KB bioconda genomicconsensus-2.3.2 |...conda-forge genomicconsensus: 2.3.2-py2...################## | 100% genomicconsensus-2.3 | 97 KB |...# # We make no warranty that any such issue...

    2045 days ago

  • One liner for reads mapping !

    bwa aln -t 8 targetGenome.fa reads.fastq | bwa samse targetGenome.fa - reads.fastq\ | samtools view -bt targetGenome.fa - | samtools sort - reads.bwa.targetGenome samtools index reads.bwa.targetGenome.bam

    1889 days ago

  • Map the long reads

    Map them agaist reference avaga genome using following codes git clone https://github.com/lh3/bwa.git cd bwa; make bwa index ref.fa bwa mem -x pacbio ref.fa pacbio.fq > aln.sam bwa mem -x ont2d ref.fa ont-2D.fq > aln.sam

    1780 days ago

  • Run sspace !

    #!/bin/bash cd `pwd` perl ~/apps/SSPACE-1.2_linux-x86_64/SSPACE_v1-2.pl \ -l libraries.txt \ -s Contigs_over200_nocp.fasta \ -k 5 \ -a 0.7 \ -x 1 \ -m 30 \ -o 20 \ -b Rayk31_scaffolds_extension

    1778 days ago

  • Perl script to run in parellel !

    ...SeqIO; my ($sequence_data_ref) = parse_genome_files($ARGV[0]); my %genome=%{$sequence_data_ref}; m...ount = 0; foreach my $chr_set (keys %genome) { $count...die $!; system ("augustus --species=caenorhabditis --outfile=$fn $seq -...

    1741 days ago

  • Find and replace in multifasta or fasta header with perl onliner

    You have a fasta file and you want to replace: "|" You are told to replace that by "_" perl -i -p -e "s/\|/_/g" genome.fasta -i = inplace editing -p = loop over lines and print each line (after processing) -e = command line script

    1670 days ago