BOL: Related items

Gblocks: eliminates poorly aligned positions and divergent regions of a DNA or protein alignment

Poonam Mahapatra — Sat, 02 Jun 2018 07:36:05 -0500

Gblocks eliminates poorly aligned positions and divergent regions of a DNA or protein alignment so that it becomes more suitable for phylogenetic analysis. This server implements the most important features of the Gblocks program to make its use as simple as possible without loosing the functionality that it is necessary in most of the cases. Other options can be changed in the stand-alone program. You can see here an example output file showing the blocks selected from a protein alignment. Further information can be found in the online documentation.

Address of the bookmark: http://molevol.cmima.csic.es/castresana/Gblocks_server.html

FGENESH - Program for predicting multiple genes in genomic DNA sequences

BioStar — Thu, 20 Dec 2018 11:55:08 -0600

FGENESH is the fastest (50-100 times faster than GenScan) and most accurate gene finder available - see the figure and the table below. In recent rice genome sequencing projects, it was cited "the most successful (gene finding) program (Yu et al. (2002) Science 296:79) and was used to produce 87% of all high-evidence predicted genes (Goff et al. (2002) Science 296:79).

Address of the bookmark: http://www.softberry.com/berry.phtml?topic=fgenesh&group=help&subgroup=gfind

3D de novo assembly (3D DNA) pipeline

Jit — Sun, 02 Feb 2020 13:41:55 -0600

For a detailed description of the pipeline and how it integrates with other tools designed by the Aiden Lab see Genome Assembly Cookbook on http://aidenlab.org/assembly.

For the original version of the pipeline and to reproduce the Hs2-HiC and the AaegL4 genomes reported in (Dudchenko et al., Science, 2017) see the original commit.

For the detailed description of the merge section see https://github.com/theaidenlab/AGWG-merge.

Address of the bookmark: https://github.com/theaidenlab/3d-dna

Metabuli 분리 improves metagenomic read classification

Abhi — Sat, 03 Jun 2023 20:15:04 -0500

Metabuli 분리 improves metagenomic read classification through metamers, DNA-AA k-mers, to be sensitive and specific, recovering 99% and 98% of DNA or AA classifiers.

Metabuli is metagenomic classifier that jointly analyze both DNA and amino acid (AA) sequences. DNA-based classifiers can make specific classifications, exploiting point mutations to distinguish close taxa. AA-based classifiers have higher sensitivity in detecting homology between query and reference sequences, leverageing higher conservation of AA sequences. Metabuli combines the information of both sequence types using a novel k-mer structure, metamer, to enable both specific and sensitive characterization of metagenomic samples. In addition, it can classify reads against a database of any size as long as it fits in the hard disk.

Address of the bookmark: https://github.com/steineggerlab/Metabuli

SVbyEye: R Package to visualize alignments between two or multiple DNA sequences

LEGE — Tue, 17 Sep 2024 02:34:57 -0500

R Package to visualize alignments between two or multiple DNA sequences including
a number of functionalities to facilitate processing of alignments in PAF format.

SVbyEye, an open-source R package to visualize and annotate sequence-to-sequence alignments along with various functionalities to process alignments in PAF format. The tool facilitates the characterization of complex SVs in the context of sequence homology helping resolve the mechanisms underlying their formation. Availability and implementation SVbyEye is available at https://github.com/daewoooo/SVbyEye.

Author: David Porubsky

Address of the bookmark: https://github.com/daewoooo/SVbyEye

BioGeek Fun

Jit — Sun, 16 Mar 2014 06:33:31 -0500

1. A futuristic computational biology student was told to write "It is in my gene!!!" on the board 100 times as a punishment. here's his response -

use warnings;
for ($count=1; $count <=100; $count++) { print "It is in my gene!!!";}

I guess, he is gonna to be a real biogeek. Nice try though. Smart kid.

2. In some perl script I found this
. . . . . .
. . . . . .
# It works for me, only God understood how it is working
while (/(<\/[^>]+>)|(<[^>]+>)|(<[^>]+>)$|([^><]+)/go) {
            $startGene=$1;
            $beginChromosome=$2;

. . . . . .
.. . . . . .
}

3. One more interesting message in Perl found …. It will must tickle you bone :)
open(my $fh, "<", "gene.txt") or kill " Me if you think this is a mistake :$!";

4. From the Perl

while () { # "The Mothership Connection is here!"
print “$_\n”; # Printing the offspring :)

5. Perl message
if ($1) { print “Just found a the error in chromosome !!!, yahoo…”; else { “That is not error, but mutation you moron!”;

6. One genome database curator walk in wine bar asked the bartender:
CREATE TABLE gene IF NOT EXISTS SexOnTheBeach;

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();
)

GOLD:Genomes Online Database

Jit — Wed, 26 Jul 2017 07:49:29 -0500

GOLD:Genomes Online Database, is a World Wide Web resource for comprehensive access to information regarding genome and metagenome sequencing projects, and their associated metadata, around the world.

https://gold.jgi.doe.gov/

Address of the bookmark: https://gold.jgi.doe.gov/

coursera genome assembly tutorial

Jit — Sat, 25 Nov 2017 08:57:25 -0600

Solutions to Coursera Genome Sequencing (Bioinformatics II)

Address of the bookmark: https://github.com/iansealy/coursera-assembly