I am new to bioinformatics. Currently, I am working with software that generates artificial FASTQ files from a given reference genome. These FASTQ files are supposed to resemble the reads that would come out of a modern next gen sequencer. Moreover, the user can customize the length of reads (the default is 76bp). What I have understood so far, is that each of the reads in a FASTQ file is a sequenced fragment of the genome. Why is that the genome needs to be fragmented into short pieces in order to be sequenced by a modern sequencer?
1 answer
Sequencing can only be performed for fairly short strands (100 to 5000 basepairs) and longer sequences must be subdivided into smaller fragments to sequence them.
The main reason being that the quality of the base (confidence with which a photo or chemical signal can be interpreted into a nucleotide identity) decreases with length and after a point it becomes hard to identify the actual base or nucleotide call.
See these links:
Log in to answer this question.
You may also be interested in this tool: https://github.com/lh3/wgsim