Results for "low-coverage"
Bio-Scripts
Create genome scaffolding with Perl
#!/usr/bin/perl use warnings; use strict; use English; use Pod::Usage; ## uses pod documentation in usage code use Getopt::Long qw(:config auto_version auto_help...2312 days ago
- # 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...
2214 days ago
Perl script to find coding regions in DNA sequences
#!/usr/bin/perl -w use strict; # if the number of input arguments is lower than 2 # return a message showing the error if (scalar(@ARGV) < 2) { print "dnalo...2159 days ago
- ➜ trinityrnaseq-Trinity-v2.6.6 git:(master) ✗ ./Trinity --seqType fq --max_memory 250G --left /media/urbe/MyPassport/rnaseq_boris/AIG_AAAOSW_8_1_C1CBGACXX.IND11_noribo_cl...
2131 days ago
Install Parrot Virtual Machine !
#Parrot is a virtual machine designed to efficiently compile and execute bytecode for dynamic languages jit@jit-HP-Pro-3335-MT:~/Downloads/parrot-8.1.0$ wget ftp://ftp....1541 days ago
- 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/b...
2079 days ago
Generates a genome coverage plot with R
library(CoverageView) ##draw a coverage plot for a test case BAM file #get a BAM test file treatBAMfile2059 days ago
- #!/bin/bash Plot the coverage script chr=$1 start=$2 end=$3 samtools depth deduped_MA605.bam > deduped_MA605.coverage awk '$1 == $chr {print $0}' deduped_MA605...
2008 days ago
Samtools commands for bioinformatician !
## count mapped reads samtools view -c -F 260 mapping_file.bam ### converting sam file into fasta samtools fasta reads_mapped.sam > reads.fasta ### converting sa...1607 days ago
Installing docker for Bioinformatics on Ubuntu !
jit@jit-HP-Pro-3335-MT:~/Downloads$ sudo apt-get remove docker docker-engine docker.io containerd runc Reading package lists... Done Building dependency tree Rea...1517 days ago
