Can someone explain the differences between sequence, reads, and contigs of genetic material such as DNA, if possible with an example?
I am new to bioinformatics, and I have not found any conclusive answers for all these concepts on the web.
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It only takes a minute to sign up.
Sign up to join this communityCan someone explain the differences between sequence, reads, and contigs of genetic material such as DNA, if possible with an example?
I am new to bioinformatics, and I have not found any conclusive answers for all these concepts on the web.
My understanding of those three words as follows:
sequence is a generic name describing order of biological letters (DNA/RNA or amino acids). Both contigs and reads are DNA/RNA or aa sequences
reads are just a short hand for sequenced reads. Usually sequenced reads refer to somewhat digital information obtained from the sequencing machine (for example Illumina MySeq) and stored in the fastq
file with quality scores per base. Reads are usually short. However "short" changes rapidly. Right now MySeq produces reads anywhere between 50-150 base pairs long (bp). From a single run (it will really depends on the run) you can get millions of reads, where each read will be set bp size e.g 100bp long. All reads are stored in a single fastq
file per replicate, where all reads in that file are usually of uniform size e.g all 5 million reads are 100bp long.
As a bioinformatician your first job is to identify where about those reads come from. Depending on the experimental goal and on what sort of sequencing you were doing e.g DNA-seq or RNA-seq you may or may not encounter contigs.
contigs are simply reads that have been assembled together. For example if you are doing de novo transcriptomics. Then you would:
I'm going to say the same thing as @Serine but in a slightly different context. Let's take an example where you want to compare smoking persons against non-smokers.
In this context, you'd want to take a DNA sequence of smoking persons. However, due to technology limitation you won't get a single DNA sequence from the sequencing machine. You'll get millions of short overlapping DNA sequences known as reads.
We need an assembler to "map" the reads and compare them with a reference genome. In this example, the reference genome could have been the human HG38.
The assembler would need to merge the overlapping-reads into a set of non-overalapping regions, known as contigs.