BOL: Related items

Funny software engineer

Ram Yash Pal — Fri, 09 May 2014 06:57:04 -0500

Alien Genome !!!

Jit — Sat, 13 Dec 2014 00:24:32 -0600

Genome sequencing, analysis and expression of Alien genome.

Note: This image/cartoon is create only for fun. It has nothing to do with any scientific findings.

A Bioinformatician’s Lament

LEGE — Thu, 29 May 2025 01:33:31 -0500

"I have a presentation tomorrow," they say,

With hopeful eyes, like it’s all child's play.
As if results bloom overnight, full-grown—
Not wrangled from chaos, and error-prone.

Oh brave soul, sit, let’s walk through the tale,
Of pipelines broken and servers that fail.
The journey starts: “The data? It’s there—
Just fetch it from S3, easy, I swear.”

Now I summon awscli with dread,
Reset my keys, credentials fed.
Configure regions, IAM roles too—
All this, and still no peek at the view.

Next up, the tool: “It’s open source!”
On GitHub, rotting, no sign of remorse.
Python 2.7, some GCC trick—
The install alone might make you sick.

Finally, progress! The pipeline runs…
Till RAM collapses and error stuns.
Oh, and the metadata? A crime,
Merged cells, font soup, out of time.

Sample IDs—what a cryptic game:
Sample_1, S1, sample-1... the same?
Controls mislabeled, cases flipped,
No wonder my sanity's starting to slip.

Then QC plots, PCA joy—
Wait, that’s a tumor labeled as a boy?
Clusters cross, and axes lie,
And I still don’t know which sample’s "guy."

But the clock ticks on, and it’s half-past doom,
They want the final UMAP soon.
With pastel colors, labeled clear—
"Can we move that legend to right here?"

Tweak by tweak, I adjust each frame,
Resize Panel B, annotate a name.
Export the plot—it starts to gleam…
Then my laptop crashes. I scream.

This is the grind, the long-haul game,
Where science hides behind code and flame.
No “Export to Nature” button to press,
Just toil and logic and hope for success.

So next time you whisper that fated line—
“I have a talk, can you make it shine?”
Know: bioinformatics is craft, not a click,
It’s science with scars, not just a quick fix.

To all who debug at 3AM light,
Who ghostwrite figures through sleepless night—
You are the backbone, silent and true,
First-author-worthy, if only they knew.

"कल मेरी प्रेज़ेंटेशन है," वो कहते हैं,

आशा भरी आँखों से, जैसे सब सहज है।
जैसे परिणाम रातोंरात प्रकट हो जाएं—
ना कि डेटा की भूलभुलैया से उखाड़े जाएं।

आओ बैठो, एक किस्सा सुनाता हूँ,
जहाँ पाइपलाइन टूटती है, और सर्वर भी थक जाते हैं।
कहानी शुरू होती है: “डेटा तो है—
बस S3 बकेट में, एकदम पास में कहीं।”

अब awscli बुलाता हूँ डरते हुए,
कुंजी सेट करूँ, क्रेडेंशियल जोड़ूं, रीजन भरूँ।
इतनी मशक्कत, फिर भी डेटा नहीं मिला,
बस सेटअप में ही पूरा दिन चला।

फिर आता है टूल: “ओपन-सोर्स है!”
GitHub पर है, 2019 से सूखा पड़ा है।
Python 2.7 चाहिए, एक पुराना कम्पाइलर,
और साथ में थोड़ी सी दुआ की ताकत।

आख़िरकार टूल चला, खुशी सी हुई,
लेकिन रन करते ही, मेमोरी ने हार मानी।
और मेटाडेटा? एक एक्सेल की आफ़त,
मर्ज़ किए हुए सेल, बस और क्या चाहिए काफ़ियत?

सैंपल आईडी? बस भगवान ही जाने—
Sample_1, sample-1, S1, और control1—
ये सब एक ही सैंपल हैं क्या?
पता तब चलता है जब पूछो दो-तीन बार।

काउंट मैट्रिक्स तैयार, अब R या Python की बारी,
QC करो, PCA प्लॉट—पर कुछ गड़बड़ भारी।
ट्यूमर और नॉर्मल का अदला-बदली खेल,
बार-बार, वही पुरानी झमेल।

आख़िर में आया मॉडलिंग का समय,
स्टैट्स, प्लॉट्स, डिफरेंशियल एक्सप्रेशन का श्रम।
लेकिन घड़ी में 5 बज चुके हैं जनाब,
और 8 बजे तक UMAP चाहिए, साफ़-सुथरा जबाब।

तो मैं कोड लिखता हूँ रात भर बैठ कर,
कलर पैलेट, जीन लेबल, लीजेंड बाहर रख कर।
फ़ॉन्ट, पैनल, एक्सिस सब सुधार,
एक्सपोर्ट करता हूँ... और लैपटॉप कहता है—"अब नहीं यार!"

इसीलिए बायोइन्फॉर्मेटिक्स में लगता है समय,
ये “बस सीरत चलाओ” या “वोल्कैनो प्लॉट बनाओ” नहीं है।
ये है सिस्टम एडमिन का काम, डेटा की सफ़ाई,
QC, डिबगिंग, और सांइस की सच्ची लड़ाई।

तो कुछ सीखें इस व्यथा से आप भी आज:
24 घंटे पहले चमत्कार मत माँगिए।
अच्छे फ़िगर साफ़ डेटा से बनते हैं।
बायोइन्फॉर्मेटिक्स जादू नहीं, विज्ञान है।
समय से बात कीजिए, प्रक्रिया का सम्मान कीजिए।

और उन सभी बायोइन्फॉर्मेटिशियनों को सलाम,
जो दूसरों की प्रेज़ेंटेशन के लिए रातों में जागते हैं—
तुम हो फ़िगर्स के भूत लेखक,
तुम हो बिना नाम के सह-लेखक।
तुम पहले लेखक बनने के हक़दार हो—
और एक लंबी नींद के भी।

Note: Written with the help of AI/LLM Tools !

Bioinformatician on Valentine's Day

BioQueen — Sun, 14 Feb 2021 11:36:32 -0600

Once asked to a bioinformatician "How is ur Valentine's Day?"

Bioinformatician replied:
if ($date == "Valentine's Day" && $me =! Bioinformatician) {
rose_day(); promise_day(); kiss_day();
}
else {
hack_genome(); ko-fi(); youtube(); do_scripting(); sleep();
)

BEAP: Blast Extension and Assembly Program

Shruti Paniwala — Mon, 11 Jun 2018 04:52:56 -0500

The Blast Extension and Assembly Program (BEAP) is a computer program that uses a short starting DNA fragment, often a EST or partial gene segment, as "primer", to recursively blast nucleotide databases in an attempt to obtain all sequences that overlaps, directly or indirectly, with the "primer" therefore help to "extend" the length of the original sequence for constructing a "full length" sequence for functional analysis, or at least to obtain neighboring regions of the segment for SNP discovery and linkage disequilibrium analysis. The confidence of assembling the resulting sequences is achieved by using a known genome, such as human genome, as a reference. https://www.animalgenome.org/tools/beap/

Address of the bookmark: https://www.animalgenome.org/tools/beap/

Charpak exchange program

Fri, 10 Nov 2023 02:33:38 -0600

The scholarship is designed for Indian students from all fields and streams of study, enrolled in an Indian institution at the Bachelors or Master’s degree level, who wish to undertake a study exchange semester programme in France (for a period of one to six months).

BENEFITS
The Charpak exchange program offers the following benefits to the awardees based on merit: monthly stipend of 860 euros social security student visa and Campus France fee waiver assistance in finding an affordable student accommodation (subject to availability)

https://www.inde.campusfrance.org/charpak-exchange-scholarship-spring-session-jan-june

HISAT2: a fast and sensitive alignment program for mapping next-generation sequencing reads

Rahul Nayak — Tue, 08 May 2018 04:27:22 -0500

HISAT2 is a fast and sensitive alignment program for mapping next-generation sequencing reads (both DNA and RNA) to a population of human genomes (as well as to a single reference genome). Based on an extension of BWT for graphs [Sirén et al. 2014], we designed and implemented a graph FM index (GFM), an original approach and its first implementation to the best of our knowledge. In addition to using one global GFM index that represents a population of human genomes, HISAT2 uses a large set of small GFM indexes that collectively cover the whole genome (each index representing a genomic region of 56 Kbp, with 55,000 indexes needed to cover the human population). These small indexes (called local indexes), combined with several alignment strategies, enable rapid and accurate alignment of sequencing reads. This new indexing scheme is called a Hierarchical Graph FM index (HGFM).

more at https://ccb.jhu.edu/software/hisat2/index.shtml

Address of the bookmark: https://github.com/infphilo/hisat2

PANDASEQ is a program to align Illumina reads, optionally with PCR primers embedded in the sequence, and reconstruct an overlapping sequence.

BioStar — Fri, 21 Sep 2018 10:19:52 -0500

Development packages for zlib and libbz2 are needed, as well as a standard compiler environment. On Ubuntu, this can be installed via:

sudo apt-get install build-essential libtool automake zlib1g-dev libbz2-dev pkg-config

On MacOS, the Apple Developer tools and Fink (or MacPorts or Brew) must be installed, then:

sudo fink install bzip2-dev pkgconfig

Address of the bookmark: https://github.com/neufeld/pandaseq

LTR_Finder: an efficient program for finding full-length LTR retrotranspsons in genome sequences.

Neel — Sun, 13 Jan 2019 07:05:53 -0600

LTR_Finder is an efficient program for finding full-length LTR retrotranspsons in genome sequences.

The Program first constructs all exact match pairs by a suffix-array based algorithm and extends them to long highly similar pairs. Then Smith-Waterman algorithm is used to adjust the ends of LTR pair candidates to get alignment boundaries. These boundaries are subject to re-adjustment using supporting information of TG..CA box and TSRs and reliable LTRs are selected. Next, LTR_FINDER tries to identify PBS, PPT and RT inside LTR pairs by build-in aligning and counting modules. RT identification includes a dynamic programming to process frame shift. For other protein domains, LTR_FINDER calls ps_scan (from PROSITE, http://www.expasy.org/prosite/) to locate cores of important enzymes if they occur.

Address of the bookmark: https://github.com/xzhub/LTR_Finder

Gooey: Turn (almost) any Python command line program into a full GUI application with one line

BioStar — Fri, 01 Nov 2019 00:29:27 -0500

Turn (almost) any Python command line program into a full GUI application with one line

The easiest way to install Gooey is via pip

pip install Gooey

Alternatively, you can install Gooey by cloning the project to your local directory

git clone https://github.com/chriskiehl/Gooey.git

run setup.py

python setup.py install

Address of the bookmark: https://github.com/chriskiehl/Gooey