BOL: Related items

Hercules: a profile HMM-based hybrid error correction algorithm for long reads

Jit — Mon, 20 Aug 2018 14:14:11 -0500

Choosing whether to use second or third generation sequencing platforms can lead to trade-offs between accuracy and read length. Several studies require long and accurate reads including de novo assembly, fusion and structural variation detection. In such cases researchers often combine both technologies and the more erroneous long reads are corrected using the short reads. Current approaches rely on various graph based alignment techniques and do not take the error profile of the underlying technology into account. Memory- and time- efficient machine learning algorithms that address these shortcomings have the potential to achieve better and more accurate integration of these two technologies. Results: We designed and developed Hercules, the first machine learning-based long read error correction algorithm. The algorithm models every long read as a profile Hidden Markov Model with respect to the underlying platformtextquoterights error profile. The algorithm learns a posterior transition/emission probability distribution for each long read and uses this to correct errors in these reads. Using datasets from two DNA-seq BAC clones (CH17-157L1 and CH17-227A2), and human brain cerebellum polyA RNA-seq, we show that Hercules-corrected reads have the highest mapping rate among all competing algorithms and highest accuracy when most of the basepairs of a long read are covered with short reads. Availability:

Hercules source code is available at https://github.com/BilkentCompGen/Hercules

Address of the bookmark: https://github.com/BilkentCompGen/Hercules

How to sequence the human genome - Mark J. Kiel

Fri, 30 May 2014 13:24:11 -0500

View full lesson: http://ed.ted.com/lessons/how-to-sequence-the-human-genome-mark-j-kiel Your genome, every human's genome, consists of a unique DNA sequence of A's, T's, C's and G's that tell your cells how to operate. Thanks to technological advances, scientists are now able to know the sequence of letters that makes up an individual genome relatively quickly and inexpensively. Mark J. Kiel takes an in-depth look at the science behind the sequence. Lesson by Mark J. Kiel, animation by Marc Christoforidis.

Pacasus: Correction of palindromes in long reads from PacBio and Nanopore

BioStar — Mon, 12 Nov 2018 05:26:48 -0600

Tool for detecting and cleaning PacBio / Nanopore long reads after whole genome amplification. Check the poster from the Revolutionizing Next-Generation Sequencing (2nd edition) conference in the source folder: https://github.com/swarris/Pacasus/blob/master/vib2017.pdf.

The prepint version is found on http://www.biorxiv.org/content/early/2017/08/09/173872

It uses the pyPaSWAS framework for sequence alignment (https://github.com/swarris/pyPaSWAS)

Address of the bookmark: https://github.com/swarris/Pacasus

Genomics and Personalized Medicine

Sun, 01 Jun 2014 23:38:42 -0500

(October 20, 2009) Michael Snyder, Professor of Genetics and Chair of the Department of Genetics at Stanford, discusses advances in gene sequencing, the impact of genomics on medicine, the potential for personalized medicine. and efforts at Stanford to further study these issues. Stanford Mini Med School is a series arranged and directed by Stanford's School of Medicine, and presented by the Stanford Continuing Studies program. Featuring more than thirty distinguished, faculty, scientists and physicians from Stanford's medical school, the series offers students a dynamic introduction to the world of human biology, health and disease, and the groundbreaking changes taking place in medical research and health care. Stanford University http://www.stanford.edu Stanford University School of Medicine http://med.stanford.edu Stanford Continuing Studies http://continuingstudies.stanford.edu Stanford University Channel on YouTube: http://www.youtube.com/stanford

AlignGraph2: similar genome-assisted reassembly pipeline for PacBio long reads

Rahul Nayak — Sun, 14 Mar 2021 09:42:47 -0500

AlignGraph2 is the second version of AlignGraph for PacBio long reads. It extends and refines contigs assembled from the long reads with a published genome similar to the sequencing genome.

More at https://academic.oup.com/bib/advance-article-abstract/doi/10.1093/bib/bbab022/6146772

Address of the bookmark: https://github.com/huangs001/AlignGraph2

Adhoc Bioinformatics Faculty Position @ NIT

Tue, 03 Jun 2014 16:19:52 -0500

NATIONAL INSTITUTE OF TECHNOLOGY, DEPARTMENT OF BIOTECHNOLOGY, WARANGAL – 506 021, Andhra Pradesh

No.NITW/BT/2014/adhoc

APPLICATIONS ARE INVITED FOR THE APPOINTMENT OF ADHOC FACULTY ON CONTRACT BASIS IN THE DEAPARTMENT OF BIOTECHNOLOGY

Period of Contract: Initially the appointment is for one semester i.e., from July 2014 up to December 2014 only.

Essential Qualifications:

i) B. Tech or equivalent in Biotechnology/ Industrial Biotechnology/ Biochemical Engineering / Chemical Engg. Or M. Sc in Microbiology/ Botany/ Zoology/ Biochemistry/Biotechnology and ii) M. Tech or equivalent in Biotechnology/Industrial Biotechnology/Bioinformatics

Integrated M. Tech in Biotechnology/Industrial Biotechnology/ Bioinformatics

Candidates must possess First class (60% aggregate marks or 6.5 CGPA) at B. Tech/ M. Sc and M. Tech.

Desirable: Ph. D Pay Package: All selected candidates shall be eligible for a consolidated pay of Rs.30, 000/- per month. Candidates with Ph. D shall be eligible for an additional amount of Rs.5, 000/- per month.

How to apply : Applications on plain paper with attested photocopies of certificate and bio data along with justification for eligibility should reach to the Head, Department of Biotechnology, National Institute of Technology, Warangal AP 506004 in the form of soft or hard copy on or before 21st June 2014 email : biotech_hod@nitw.ac.in

Intimation: No separate call letters will be sent to the candidates. All the eligible candidates will be notified in the institute web site on 23rd June 2014. All the eligible candidates are requested to report for the interview to the Head, Department of Biotechnology at 9:00 AM on 27th June 2014

Joining: Selected candidates will be informed and they are expected to join immediately.

http://www.nitw.ac.in/nitw/announcements/2014/Bio-Adhoc%20Advt.%20May-2014.pdf

Cerulean: A hybrid assembly using high throughput short and long reads

Rahul Nayak — Tue, 05 Jun 2018 10:10:15 -0500

Cerulean extends contigs assembled using short read datasets like Illumina paired-end reads using long reads like PacBio RS long reads. Cerulean v0.1 has been implemented with bacterial genomes in mind. The method is fully described in Deshpande, V., Fung, E. D., Pham, S., & Bafna, V. (2013). Cerulean: A hybrid assembly using high throughput short and long reads. arXiv preprint arXiv:1307.7933. http://arxiv.org/abs/1307.7933

Address of the bookmark: https://sourceforge.net/projects/ceruleanassembler/

Search Shell Command History

Rahul Nayak — Thu, 12 Jun 2014 17:43:34 -0500

We use couple of hundreads of command in daily basis. Most of them are actually repeated several time. The question remain open how do I search old command history under bash shell and modify or reuse it?

Now a days almost all modern shell allows you to search command history if enabled by user. Use history command to display the history list with line numbers. Lines listed with with a * have been modified by user.

Shell history search command

Type history at a shell prompt:
$ history

It will display the list of all used commandline history with an serial number.

To search particular command, enter:
$ history | grep command-name
$ history | egrep -i 'scp|ssh|ftp'
Emacs Line-Edit Mode Command History Searching

To get previous command containing string, hit [CTRL]+[r] followed by search string:

(reverse-i-search):

To get previous command, hit [CTRL]+[p]. You can also use up arrow key.

CTRL-p

To get next command, hit [CTRL]+[n]. You can also use down arrow key.

CTRL-n

fc command

Apart from hostory command there are fc command to extract the command from history. The fc stands for either "find command" or "fix command.

For example list last 10 command, enter:
$ fc -l 10
To list commands 130 through 150, enter:
$ fc -l 130 150
To list all commands since the last command beginning with ssh, enter:
$ fc -l ssh
You can edit commands 1 through 5 using vi text editor, enter:
$ fc -e vi 1 5

Delete command history

The -c option causes the history list to be cleared by deleting all of the entries:
$ history -c

Understating pacbio reads name !

BioJoker — Fri, 23 Nov 2018 07:36:46 -0600

m140415_143853_42175_c100635972550000001823121909121417_s1_p0/553/3100_11230 0.99 24
└1┘└─────2─────┘└──3─┘└────────────────4────────────────┘└5┘└6┘└7┘└────8────┘└─9─┘└10┘

"m" = movie
Time of Run Start (yymmdd_hhmmss)
Instrument Serial Number
SMRT Cell Barcode
Set Number (a.k.a. "Look Number". Deprecated field, used in earlier version of RS)
Part Number (usually "p0", "X0" when using expired reagents)
ZMW hole number
Subread Region (start_stop using polymerase read coordinates)
readScore
barcodeScore

Bioinformatician’s Pocket Reference !!

RAJESH DETROJA — Sun, 08 Jun 2014 09:56:58 -0500

It is amusing how brain of bioinformaticians work! Learning a new programming language for days feels so much of fun that making 5 minute discussion with neighbours (unless under special circumstances!) in our own mother-tongue. Today every bioinformatician keeps more than few languages and core IT toolkits on their plate. It has become mandatory to be able to mould different code snippets to build our own custom workflows, and thus keeping syntax at our fingertips has become essential.Although Google is best way to get syntax problem solved, it is not a bad idea to keep reference sheets is our smartphones or stick out some printed sheets on the back of your door, in the old fashion way!!

Address of the bookmark: http://infoplatter.wordpress.com/2014/04/06/bioinformaticians-pocket-reference/