BOL: Wire posts tagged with 'samtools'

Wire post by Aaryan Lokwani

Aaryan Lokwani — Fri, 22 Jun 2018 04:41:13 -0500

Count the contigs length in index fasta file with samtools. $ awk '{s+=$2}END{print s}' scaffolds.fasta.fai #Fasta #Index #Length #Count #Size #Samtools

Wire post by Shruti Paniwala

Shruti Paniwala — Wed, 26 Oct 2016 09:24:08 -0500

Extract a fasta sequence from multifasta file: samtools faidx multitest.fa and then samtools faidx multitest.fa fasta_id > out.fa #Samtools #Fasta #Extract #Multifasta

Wire post by Jit

Jit — Thu, 30 Jun 2016 04:48:56 -0500

Extract reference nucleotide from BAM file: samtools view your.bam | awk '{print substr( $10, 100, 1)}' #SAM #BAM #Samtools #Tricks #NGS

Wire post by Poonam Mahapatra

Poonam Mahapatra — Fri, 24 Jun 2016 09:16:32 -0500

#Extract all #Reads from #BAM file for a #region Chr10:18000-45500 using #samtools: samtools view input.bam "Chr10:18000-45500" > output.bam