Tag: cram

Read bam/cram file with IGV from aws s3

Hi all, We store our alignment files on aws s3. I would like to be able to open them with IGV without needing to download them completely, but I can’t find an optimal solution. If I get a pre-signed url it works but it’s not convenient. I try to follow…

Continue Reading Read bam/cram file with IGV from aws s3

Ubuntu Manpage: samtools reheader – replaces the header in the input file

Provided by: samtools_1.13-2_amd64 NAME samtools reheader – replaces the header in the input file SYNOPSIS samtools reheader [-iP] [-c CMD | in.header.sam ] in.bam DESCRIPTION Replace the header in in.bam with the header in in.header.sam. This command is much faster than replacing the header with a BAM→SAM→BAM conversion. By default…

Continue Reading Ubuntu Manpage: samtools reheader – replaces the header in the input file

The Biostar Herald for Tuesday, September 21, 2021

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from Istvan Albert, and was edited by Istvan…

Continue Reading The Biostar Herald for Tuesday, September 21, 2021

Bedtools: Merging Many Bed Files

Bedtools: Merging Many Bed Files 2 I am using the algorithm CookHLA for my research. As part of its preparation, I need to feed it a bed file representing at least 100 of my samples. I have made the bed files for 500 samples using samtools and bedtools in a…

Continue Reading Bedtools: Merging Many Bed Files

Best Omic file compressor?

Best Omic file compressor? 1 Our team has been having storage space issues; we predicted that we will not have enough available memory to store the files generated by our pipelines. Standard file compressors (gzip, bzip2, 7zip) weren’t cutting it and I started experimenting with file-specific compressors. This is where…

Continue Reading Best Omic file compressor?

[main_samview] fail to read the header from “-“.

[main_samview] fail to read the header from “-“. 1 I am attempting to run a file through an algorithm I have been using, HLA*LA. On running the samtools command within the algorithm, I have unfortunately been getting this error. After trying to debug this following other guides, I am seeking…

Continue Reading [main_samview] fail to read the header from “-“.

How to extract all sequences mapped to a transcript from Kallisto output

How to extract all sequences mapped to a transcript from Kallisto output 0 I ran Kallisto with the –pseudobam option. How do I extract all the short reads that are mapped to a single transcript (e.g. ENST00000367969.8)? As a person without any previous SAM/BAM experience, I tried the following things…

Continue Reading How to extract all sequences mapped to a transcript from Kallisto output

install GenomicFeatures fail

install GenomicFeatures fail 1 @5b9023e7 Last seen 19 hours ago China BiocManager::install(‘GenomicFeatures’) results show ‘getOption(“repos”)’ replaces Bioconductor standard repositories, see ‘?repositories’ for details replacement repositories: CRAN: mirrors.tuna.tsinghua.edu.cn/CRAN/ Bioconductor version 3.14 (BiocManager 1.30.16), R 4.1.0 (2021-05-18) Installing package(s) ‘GenomicFeatures’ also installing the dependencies ‘Rhtslib’, ‘Rsamtools’, ‘GenomicAlignments’, ‘rtracklayer’ Packages which are only…

Continue Reading install GenomicFeatures fail

.tar.gz = same size as before?

BAM compression: .tar.gz = same size as before? 2 I tried to compress 5 bam files using: tar -czvf original_bams.tar.gz *.bam The resulting file sizes (“ll –block-size=M”) are: 8067M file1.bam 6962M file2.bam 10662M file3.bam 7794M file4.bam 7346M file5.bam 40828M original_bams.tar.gz There’s a difference of 3MB between the archive and the…

Continue Reading .tar.gz = same size as before?