Skip to main content
100%

deltaG p7308 calculations

✓ Published0🌍 Public
Aakeshavan
Last edited Nov 22, 2013
Created on Nov 22, 2013

A histogram displays the distribution of calculated free-energy values (ΔG) for 500-base-pair subsequences of the p7308 genome, with results extracted from NuPack’s minimum free energy (mfe) output files. The plot uses matplotlib via prettyplotlib to render 30 bins, labeling axes with ΔG and sequence count, and titles the figure with the temperature condition. The underlying pipeline processed 7,308 subsequences, as indicated by the file count and extraction function.

AI-generated description

Similar vizzes

Loading thumbnail…

Using nipype for mfe calculation of p7308 sub-sequences.

This example visualizes the distribution of mean free energy (ΔG) values for overlapping 500-base-pair sub-sequences of the p7308 genome. The notebook computes these values using NuPack’s mfe binary through a Nipype pipeline, then plots a histogram of the results with matplotlib and prettyplotlib. The code extracts ΔG values from text output files and displays the histogram with 30 bins, labeled axes, and a title indicating the sequence length and temperature.

Aakeshavan
80% match
Loading thumbnail…

Gist 143778d0c50388a317ebc30c0ec5de86

This Perl-based command-line tool calculates GC content and nucleotide percentages from a FASTA file. It reads the input file line by line, skipping headers, and tallies counts of G, C, A, T, and N bases. The output prints GC%, individual base percentages, and the N count, using basic string operations and conditional logic rather than external visualization libraries. The script demonstrates a simple text-based data summary workflow for bioinformatics sequence analysis.

NNpfn
55% match
Loading thumbnail…

Gist 6b30d347b0402263522888c0a96653ad

This histogram visualizes the distribution of sequence lengths from Oxford Nanopore sequencing data, comparing 1D+2D and 2D reads of *Silene latifolia* DNA. The chart shows how many reads fall into each length bin, highlighting the range and frequency of fragment sizes for both datasets. The underlying shell script uses `bioawk` with the `fastx` input mode to extract the `length()` of each sequence, writing the resulting measurements to separate text files before plotting. The visualization relies on the bioawk library’s FASTA/FASTQ parsing capabilities to generate the raw length counts.

NNpfn
53% match
Loading thumbnail…

Gist ac0e4a9bdb2f03d365c4402d04ecba09

This example shows a command-line workflow for detecting genetic variants from sequencing reads mapped to a reference genome. It first uses `samtools mpileup` to generate a binary call format (BCF) file, then `bcftools call` to produce a VCF file of variants. A second, independent variant call is made using `freebayes` with a ploidy setting of one, outputting its own VCF. The script processes all sorted BAM files in a mapping directory, using the NC_012967.fas reference, demonstrating two standard variant-calling pipelines for haploid samples.

NNpfn
52% match
Loading thumbnail…

Gist f919aa6a77ff3feb4d8a6d77af5dee85

This example computes the number and length distribution of contiguous sequences from a FASTA file, then writes the results to a tab-delimited text file. It uses Perl’s `Math::Complex` module for mathematical operations, though only basic arithmetic appears. The script reads `gapCloserMinia.break.broken_assembly.fasta`, tracks sequence lengths via regex matching for `[ACTGN]+` lines, and outputs a numbered list of contig sizes, along with the mean and standard deviation. The visualization shows a histogram or table of contig lengths, derived directly from the computed file `tailles_contigues.txt`. There is no graphical rendering; the output is purely textual.

NNpfn
52% match
Loading thumbnail…

Gist 3d7c2ec8630fd8451ee2b2ed513c2c4a

This example documents the assembly of a microbial genome from paired-end Illumina sequencing data. It shows the command-line workflow that converts raw FASTQ files to FASTA format and then performs a de novo assembly using the IDBA-UD software. The script demonstrates the use of the `fq2fa` utility to merge and convert paired reads, followed by an `idba_ud` assembly run with a single thread. The final assembly output is moved to a dedicated directory for further analysis, illustrating a typical bioinformatics pipeline for genomic sequence reconstruction.

NNpfn
51% match
Loading thumbnail…

Gist babdf0de5f2b986506f621199e6a642d

This Perl script calculates the total length of a genome assembly by reading a FASTA file named `gapCloserMinia.break.broken_assembly.fasta`. It iterates line by line, summing the number of nucleotide characters (A, C, T, G) found on each line, ignoring headers and other non-sequence content. The script uses a simple regex match to detect sequence lines and accumulates the total count, printing the final genome size in base pairs to standard output. The computation is performed entirely with Perl's built-in file handling and string operators, requiring no external libraries or modules.

NNpfn
51% match
Loading thumbnail…

Gist aa8c974207019f98e6d45a49488b2c4b

This visualization shows the process of generating pileup files from sequenced reads for a set of 29 samples. The Bash script iterates through a list of sample identifiers, using the `samtools mpileup` command to produce alignment summaries for each sample against a reference genome. The script relies on the command-line APIs of `samtools` for pileup generation, operating on sorted BAM files. The output from this script would be a text-based pileup format, which is not a graphical visualization but a data-preparation step for subsequent analysis. The code references a reference file and a list of SRR accession numbers, indicating the source is genomic sequencing data.

NNpfn
51% match