Structure of Binary files used for storing sequencing data-bam and sff
Many times bioinformatician needs to parse binary files like bam and sff. Advantage of binary files is that they occupy less space in memory with maximum information content. Link for those who looking for structure of Bam and sff file: Bam: http://samtools.sourceforge.net/SAMv1.pdf (from...Tags: bam, sff, sam, iontorrent data, 454 data, pyrosequencing data, hxd, hexadecimal, binary files, sequencing
3917 days ago
Tags: sam, bam, conversion
2288 days ago
How to extract all the reads mapped at specific region in a genome?
I want to extract all the reads mapped at chr5: 200-6000.Tags: extract, reads, mapped, bam, sam, ngs
2285 days ago
Qualimap2: Evaluating next generation sequencing alignment data
Qualimap 2 is a platform-independent application written in Java and R that provides both a Graphical User Inteface (GUI) and a command-line interface to facilitate the quality control of alignment sequencing data and its derivatives like feature counts. Supported types of experiments i...Tags: Qualimap2, Evaluating, next, generation, sequencing, alignment, coverage, plot, ngs, bam
2060 days ago
samtool error: open: fail to open file out.bam
samtools sort -o out.bam in.bamopen: No such file or directory[bam_sort_core] fail to open file out.bamTags: samtool, error, open, out, bam, ngs
1431 days ago
How to resolve samtools index error "Chromosome blocks not continuous" ?
samtools index mappedPairs.sorted.bam return the following error ... [bam_sort_core] merging from 0 files and 2 in-memory blocks...[E::hts_idx_push] Chromosome blocks not continuoussamtools index: failed to create index for "mappedPairs.sorted.bam"Tags: samtools, index, sort, bam, reads
994 days ago
BBTools for bioinformatician !
BBMap.sh Mapping Nanopore reads BBMap.sh has a length cap of 6kbp. Reads longer than this will be broken into 6kbp pieces and mapped independently. Code: $ mapPacBio.sh -Xmx20g k=7 in=reads.fastq ref=reference.fa maxlen=1000 minlen=200 idtag ow int=f qin=33 out=mapped1.sam minratio=0.15...Tags: BBMap, Map, NGS, Reads, Sam, Bam, Contamination, Filtration
2268 days ago