Results for "in-depth"
Bio-Scripts
# 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...2221 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...
2087 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 treatBAMfile2067 days ago
Generate simulated polyploid genome !
#Generate 3% divergence msbar -point 4 -count 16558 toy.fasta > toyheterozygous3percent.fasta #Cat both files cat toy.fasta toymutated3percent.fasta > toyheterozy...2065 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_MA6...
2016 days ago
- # $1 = vcf file # $2 = input bam file # $3 = output QV file module load samtools NUM_BP=`samtools depth $2 | perl -e '$c = 0; while(){chomp; @s = split(/\t/); if...
1616 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...1615 days ago
Tadpole is 250x faster than SPADes assembler !
lege@jit-Lenovo-ideapad-320-15ISK:~/Downloads/MyTools/Vir$ tadpole.sh Written by Brian Bushnell Last modified July 16, 2018 Description: Uses kmer counts to ass...988 days ago
