anvi-script-gen-short-reads [program]

Generate short reads from contigs. Useful to reconstruct mock data sets from already assembled contigs.

🔙 To the main page of anvi’o programs and artifacts.

Authors

A. Murat Eren (Meren)

Can consume

configuration-ini

Can provide

short-reads-fasta

Usage

This program uses already assembled contigs to create a mock list of short reads. You can then use these short reads to reassemble your data in order to test alternative assembly programs or analysis methods as a positive control.

Basically, this attempts to undo the assembly and produce a data set that could have been directly received from laboratory sequencing. While the computer’s mock short reads won’t be perfect, they can be used to make sure your analysis pipeline is working from step 1.

Example Usage

This program takes an INI file - a form of text file containing various information. For this program, the example provided in the anvi’o test suite looks like this:

[general]
short_read_length = 10
error_rate = 0.05
coverage = 100
contig = CTGTGGTTACGCCACCTTGAGAGATATTAGTCGCGTATTGCATCCGTGCCGACAAATTGCCCAACGCATCGTTCCTTCTCCTAAGTAATTTAACATGCGT

Note that this file contains both the contig that you want to break down, and various information about the short reads that you want to create. To run this program, just call

anvi-script-gen-short-reads configuration-ini \ --output-file-path short-reads-fasta

The resulting FASTA file with short reads will cover the contig with short reads that are 10 nts long at 100X coverage. There will also be an error-rate of 0.05, to mimic the sequencing errors you would get from sequencing in the wet lab.

Edit this file to update this information.

Additional Resources

Are you aware of resources that may help users better understand the utility of this program? Please feel free to edit this file on GitHub. If you are not sure how to do that, find the __resources__ tag in this file to see an example.