Gist 6344d0fe9ded24d86cd9e580c72e7101

100%

Gist 6344d0fe9ded24d86cd9e580c72e7101

✓ Published0🌍 Public

Last edited Feb 8, 2019

Created on Feb 8, 2019

This visualization compares read-length distributions before and after error correction of MinION long reads from Silene latifolia. The code uses the LoRDEC tool to correct and trim reads, then assembly-stats to generate summary statistics, likely producing histograms or cumulative frequency plots. The rendering relies on command-line output from assembly-stats, displaying changes in sequence length and count across three datasets: raw, corrected, and trimmed reads.

AI-generated description

Similar vizzes

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

73% match

Loading thumbnail…

Gist fb39aee08b636aec6b8a68d4011241cf

This visualization counts the number of contigs in a genome assembly file. It parses the FASTA-formatted file `gapCloserMinia.break.broken_assembly.fasta` using a Perl script, incrementing a counter each time a line starts with the `>` header character. The final count is printed as a plain integer to the terminal. The script relies solely on Perl's built-in file-handling and regular-expression capabilities, with no external visualization libraries, and does not generate any graphical output.

NNpfn

63% match

Loading thumbnail…

Gist 6cde24482216f46654efb9e0cef0f166

This example shows a Unix shell script for assembling and evaluating MinION nanopore sequencing data. It configures environment variables for Java, Perl, and the Canu assembler, then runs Canu with parameters for a 10kb genome, specifying an error rate and corrected nanopore reads as input. The script uses the `assembly-stats` tool to generate statistics for the resulting contigs and unassembled sequences, printing summary metrics to the console. The code relies on the Canu command-line interface, the `assembly-stats` binary, and shell environment exports for dependency management.

NNpfn

63% 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

62% 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

60% match

Loading thumbnail…

Gist 9a3defa4596db38a11fc0ec9549b4834

This example demonstrates a preprocessing workflow for raw paired-end Illumina sequencing data. It shows the trimming of low-quality bases and removal of adapter sequences using the `fastq-mcf` tool from the ea-utils suite, followed by quality control reporting with `FastQC`. The script runs `fastq-mcf` with parameters for a minimum quality score of 20, a mean quality threshold, and a 150-base pair read length, producing trimmed FASTQ files. These outputs are then processed by `FastQC` to generate quality assessment reports for both the forward and reverse reads, illustrating a standard bioinformatics pipeline for high-throughput sequence data preparation.

NNpfn

59% 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

58% 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

56% match