BOL: Related items

Bioinformatician Dreams

Jitendra Narayan — Wed, 21 Aug 2013 10:50:45 -0500

Bioinformatician life is interconnected, they always dream for a powerful server, little more space on server as they are generating lots of data per run, dream to publish results in good impact journals, meetings reminders :) and research analysis off course!!!

16sRNA Database Download

LEGE — Wed, 24 Apr 2024 04:33:15 -0500

Downloading 16S rRNA databases can be crucial for various bioinformatics analyses, especially in microbiome research. However, it's important to note that databases can vary based on your specific needs, such as the taxonomic coverage you require or the type of analysis you're performing. Here's a general guideline on how you can obtain 16S rRNA databases:

NCBI (National Center for Biotechnology Information):
- NCBI provides various databases related to genetic information, including 16S rRNA sequences.
- You can access the 16S ribosomal RNA sequences from NCBI's Nucleotide database (https://www.ncbi.nlm.nih.gov/nucleotide/).
- Perform a search using keywords like "16S rRNA" or specific bacterial names to find relevant sequences.
- You can download sequences individually or in batches using the provided tools.
GreenGenes:
- GreenGenes is a widely used 16S rRNA gene sequence database.
- You can access it at http://greengenes.secondgenome.com/.
- GreenGenes provides precompiled databases for various purposes, including classification, alignment, and phylogenetic analysis.
SILVA:
- SILVA (https://www.arb-silva.de/) is another comprehensive database for ribosomal RNA (rRNA) sequences.
- It covers not only 16S rRNA but also other ribosomal RNA sequences.
- SILVA provides precompiled databases for various purposes, including taxonomic classification and alignment.
Ribosomal Database Project (RDP):
- RDP (http://rdp.cme.msu.edu/) is a curated database that offers 16S rRNA sequences.
- It provides tools for sequence analysis and classification.
- You can download sequences and taxonomy information from their website.
QIIME (Quantitative Insights Into Microbial Ecology):
- QIIME (https://qiime2.org/) is a widely used bioinformatics platform for microbiome analysis.
- It provides tools for analyzing microbial communities, including processing 16S rRNA sequences.
- QIIME often includes its own preprocessed 16S rRNA databases that can be used for analysis within the platform.

Before downloading any database, make sure to read the terms of use and citation requirements, as some databases may have specific usage policies. Additionally, consider the compatibility of the database with your analysis pipeline and software tools.

NCBI 16s RNA database location ftp://ftp.ncbi.nih.gov/blast/db/16SMicrobial.tar.gz

BIOINFORMATICS

Wed, 28 Aug 2013 19:16:33 -0500

This is a promo video for the brand new cross-boarder branch of study - BIOINFORMATICS. It´s a co-operation between Johannes Kepler University in Linz (Austria) and University of South Bohemia in České Budějovice (Czech Republic). Written, Edited and Directed by, DOP, VFX: Jan Míka Sound by: Mirek Šmilauer Narrator: Jack Bright Produced by: FILMOFON (http://www.filmofon.cz) Released: Nov 2012

mosdepth: fast BAM/CRAM depth calculation for WGS, exome, or targeted sequencing

Jit — Wed, 13 Nov 2019 22:20:19 -0600

mosdepth can output:

per-base depth about 2x as fast samtools depth--about 25 minutes of CPU time for a 30X genome.
mean per-window depth given a window size--as would be used for CNV calling.
the mean per-region given a BED file of regions.
a distribution of proportion of bases covered at or above a given threshold for each chromosome and genome-wide.
quantized output that merges adjacent bases as long as they fall in the same coverage bins e.g. (10-20)
threshold output to indicate how many bases in each region are covered at the given thresholds.
A summary of mean depths per chromosome and within specified regions per chromosome.

Address of the bookmark: https://github.com/brentp/mosdepth

R and Bioconductor Tutorial

Jitendra Narayan — Fri, 23 Aug 2013 08:23:59 -0500

This tutorial is intended to introduce users quickly to the basics of R, focusing on a few common tasks that biologists need to perform some basic analysis: load a table, plot some graphs, and perform some basic statistics. More extensive tutorials can be found on the project website and via bioconductor (not covered here).

You can add more tutorial links in comments if found new pages.

Address of the bookmark: http://manuals.bioinformatics.ucr.edu/home/R_BioCondManual

Public Databases for Bioinformatics !

Jit — Tue, 23 Mar 2021 05:32:15 -0500

https://www.nature.com/articles/s41467-020-17155-y

Server Infrastructure:

File Server:

dhara: Synology 3614 Storage Appliance
4 Core Xeon
108TB disk storage
10Gb ethernet to SCG3
Access atx: dhara:5000
Has btsync server (try it - its much better than dropbox)

Compute Servers:

nandi: Kundaje and Phi Server
24 intel cores
256GB RAM
500GB of SSD storage 
36TB RAID6 local storage
4 Intel Phi's (space for 4 more GPU's)


durga: Montgomery and sensitive data
24 intel cores
256GB RAM
500GB of SSD RAID0 storage 
60TB RAID6 local storage

mitra: Bassik and Web/DB Server
24 core
256GB RAM 
500GB of SSD RAID0 storage 
36TB RAID6 local storage

vayu: Kundaje GPU server
4 core
64GB RAM 
200GB of SSD storage 
8TB RAID10 local storage
4 Nvidia GTX 970 4GB GPUs

amold: Bickel and SGE server
32 AMD core
128GB RAM 
200GB of SSD storage 
12TB RAID5 local storage

wotan: Bickel and SGE server
64 AMD core
256GB RAM 
200GB of SSD storage 
12TB RAID5 local storage

Filesystem:

/users/$USER
default home directory
full backups nightly 
nfs mount to dhara
should store code, papers, and other highly processed data here

/mnt/data/
globally accessible data
should store common data here
e.g. genomes and indexes, annotations, ENCODE data  
if you dont want this to count towards your quote you must chown

/mnt/lab_data/$LAB/
lab accessible data
should store lab project data here 
e.g. ATAC-seq prediction data, enhancer prediction, motif calls

/srv/scratch/$USER
fast local storage
not backed up, but on raid and data will never be deleted
most analysis should be performed here

/srv/persistent/$USER
fast local storage
synced nightly, but not backed up
       ie if the hard drives fail or you delete something and notice 
       within 24 hours we can recover. Otherwise not. (vs home which is 
       properly backed up )  
intermediate analysis products that would be hard to recover should be stored here 
       e.g. stochastic analysis results that need to be kept so that paper 
       results can be reproduced

/srv/www/$LABNAME/
web accessible from mitra.stanford.edu
*NOT BACKED UP*

Some parallel programming patterns:

# gzip a bunch of files
parallel gzip -- *.FILESTOGZIP

# fork example in python:
(for more detailed examples look at 
 https://github.com/nboley/grit/ grit/lib/multiprocessing_utils.py)

import os
import time
import random

import multiprocessing

class ProcessSafeOPStream( object ):
    def __init__( self, writeable_obj ):
        self.writeable_obj = writeable_obj
        self.lock = multiprocessing.Lock()
        self.name = self.writeable_obj.name
        return
    
    def write( self, data ):
        self.lock.acquire()
        self.writeable_obj.write( data )
        self.writeable_obj.flush()
        self.lock.release()
        return
    
    def close( self ):
        self.writeable_obj.close()

def worker(queue, ofp):
    # Try without this
    random.seed()
    while True:
        i = queue.get()
        if i == 'FINISHED': return
        # simulate an expensive function
        x = random.random()
        time.sleep(x/10)
        print i, x
        ofp.write("%i\t%s\n" % (i, x))

NSIMS = 10000
NPROC = 25

# populate queue
todo = multiprocessing.Queue()
for i in xrange(NSIMS): todo.put(i)
for i in xrange(NPROC): todo.put('FINISHED')

ofp = ProcessSafeOPStream( open("output.txt", "w") )

pids = []
for i in xrange(NPROC):
    pid = os.fork()
    if pid == 0:
       worker(todo, ofp)
       os._exit(0)
    else:
       pids.append(pid)  

for pid in pids:
    os.waitpid(pid, 0)

ofp.close()

print "FINISHED"

For use case 1 we obtained the following ENCODE and ROADMAP datasets https://www.encodeproject.org/files/ENCFF446WOD/@@download/ENCFF446WOD.bed.gz, https://www.encodeproject.org/files/ENCFF546PJU/@@download/ENCFF546PJU.bam, https://www.encodeproject.org/files/ENCFF059BEU/@@download/ENCFF059BEU.bam. Blacklisted regions were obtained from http://mitra.stanford.edu/kundaje/akundaje/release/blacklists/hg38-human/hg38.blacklist.bed.gz. The human genome version hg38 was obtained from http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz.

For use case 2 we used the set of narrowPeak files summarized in https://github.com/wkopp/janggu_usecases/tree/master/extra/urls.txt (archived version v1.0.1). The human genome version hg19 was obtained from http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/hg19.fa.gz

For use case 3 we used the ENCODE datasets https://www.encodeproject.org/files/ENCFF591XCX/@@download/ENCFF591XCX.bam, https://www.encodeproject.org/files/ENCFF736LHE/@@download/ENCFF736LHE.bigWig, https://www.encodeproject.org/files/ENCFF177HHM/@@download/ENCFF177HHM.bam as we as the GENCODE annotation v29 from ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_29/gencode.v29.annotation.gtf.gz.

Address of the bookmark: http://mitra.stanford.edu/

What is Bioinformatics?

Wed, 28 Aug 2013 06:53:05 -0500

Illustration and Animation: Rachel Robinson Script: Tiffany Trent Voice-over: Kris Monger Sound: Glisten Carefully by Guennadi Malyshevski

IBM Research Computational Biology Center

Thu, 29 Aug 2013 08:43:59 -0500

The IBM Computational Biology Center embraces activities at Yorktown Heights, with strong affiliations with activities at Almaden and other IBM Research Centers. Computational Biology (CompBio) including bioinformatics is the study of how computer systems can manage, analyze, and simulate the complex structures and processes inherent in living systems. CompBio Research at IBM spans pattern recognition in sequences, structures and processes, the studying of systems ranging from single protein molecules through to complex molecular interactions, and the data analysis, interpretation and reverse-engineering of complex disease-lifestyle-genomic interactions for fuller biological understanding. "CompBio" has a flavor of its own independant of its parents, biology and computer science. Nonetheless, CompBio is inherently a multi- disciplinary field with important applications in biology, chemical physics, materials science, agriculture, chemistry and ultimately nanotechnology.

jobTree based python wrapper to run the genome simulation tool suite Evolver

Jit — Fri, 08 Dec 2017 16:26:32 -0600

evolverSimControl (eSC) can be used to simulate multi-chromosome genome evolution on an arbitrary phylogeny (Newick format). In addition to simply running evolver, eSC also automatically creates statistical summaries of the simulation as it runs including text and image files. Also included are convenience scripts to: check on a running simulation and see detailed status and logging information; extract fasta sequence files from the leaf nodes of a completed simulation; extract pairwise multiple alignment files (.maf) from leaf and branch nodes from a completed simulation and with the help of mafJoin, join them together into a single maf covering the entire simulation.

Address of the bookmark: https://github.com/dentearl/evolverSimControl

May 17, 2013 - Differential Network Analysis

Wed, 04 Sep 2013 16:22:28 -0500

Steve Horvath presents Differential Network Analysis at the UCLA Human Genetics/Biostatistics Network Course