This is a test version of Biostars. For the public version, visit https://www.biostars.org.
First thing to do when you receive your sequenced genome of interest.

As it is clear in the title I am asking for advice

What is the first thing that I shall do when I receive the sequenced data for the genome am interested in,

like for example is there is a tool to check that the whole genome where sequenced probably? "there was no missing parts that was not sequenced". is the sequence good and I can proceed to the other steps? or I need to repeat something "like if some regions are not covered for example", and what is the other thing that I need to but in my consideration from your expertise and point of view.

Thanks in advance,

assembly genome sequence next-gen
• 4,842 views • link
updated by Ram Ram • written by Medhat Medhat

6 answers

map it (or a subset of it) against the genome that will tell you right away what the data looks like

• 0 views • link
written by Istvan Albert Istvan Albert

What if I do not have reference genome?

• 0 views • link
written by Medhat Medhat

well that's trickier, you can always try a closely related species,

also look for and remove contaminants, we just had two situations where not dealing with contamination right away led to setbacks.

• 0 views • link
written by Istvan Albert Istvan Albert

about removing contamination did you meant using blast for example to decide if there is other sequences than the expected specious and then remove it? or I miss understood?

If I rightly understood "How to remove it if there is such thing?"

• 0 views • link
updated by Ram Ram • written by Medhat Medhat

yes align it with a short read aligner (but not blast because you want to remove only reads that match very closely) to the contaminant then export from the alignment the unaligned reads into a new fastq file.

• 0 views • link
written by Istvan Albert Istvan Albert

I will do a quality check first with FastQC or simple stats with R, sometimes data is bad from the beginning. After that, as Istvan said, you can map your sequences to the genome and compute the average coverage (BWA, STAR or your prefer mapper) and check uncovered or suspicious regions such as high coverage in repetitive regions.

• 0 views • link
written by JC JC

Check the quality of the sequencing reads before mapping(?), for instance, using fastx and trim (adapters/primers, if any) accordingly. Then map!

• 0 views • link
updated by Ram Ram • written by xb xb
FastQC first, then if it is genome re-sequencing map and check coverage metrics, for de-novo genome assemble contigs and check metrics like N50 contig size
• 0 views • link
written by mikhail.shugay mikhail.shugay

If you do not have a reference genome, it is hard to find regions not covered. A few tips:

  • run SGA's preqc, this will tell you quite a bit about your genome and dataset
  • run assemblies and use tools like blobology to assess species content
  • if your species is a vertebrate, run CEGMA to check whether the gene space of your assembly seems complete (helps also to choose between assemblies)
  • if it is a bacteria, on the other hand, run iMetAMOS, it does much of the above in an automated fashion.
• 0 views • link
updated by Ram Ram • written by lexnederbragt lexnederbragt

+1 very informative but, do you have any idea if I'm dealing with plant genome?

• 0 views • link
written by Medhat Medhat

Adapters can make an assembly a real mess. Sometimes, there might be partial adapters also present in the reads. So, trimming them atleast within our scope is better. If adapters are not known, requesting them from the sequencing center and running experiments is fruitful.

• 0 views • link
written by Prakki Rama Prakki Rama

Log in to answer this question.

More posts like this