quant coverage
Thanks flxlex,
good suggestions. I indeed intended tolok at the assembly and coverage of that specific contigs. I just find the numbers for repeats so high. In the latter 800 fold coverage was a transposase IS4.....
Didn't know newbler could do that. I will install the software shortly on our linux machine and see how it goes.
Alex
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Your extra-high coverage contig should indeed represent a repeated piece of genome, as the 454 assembly program (newbler) collapses repeats. To find out where it belongs, you can have newbler generate the contig graph file (using the -g option on the command line). This file contains information on the nodes, i.e. the contigs (length, coverage) and on the edges, i.e. the (number of) reads that cross the boundary between contigs (start in one contig and end in another). For your high-coverage contig, you will find that it has several neighbors on both ends. Your next task is to figure out which neighbors belong together (with the high-coverage contig in between), by for instance PCR. There might be sequence variants among the different copies of the repeated region, and sequencing the PCR products will tell youwhich variant belongs where. Of course, this works best for contigs that are within the length that can be sequenced by Sanger reads...
Or, you might want to invest in paired-end sequencing to place the repeats... Apparently, one run from an 8kb jumping paired end library will get you one scaffold (ordered, oriented contigs with estimates of gap sizes) for bacterial genomes (we are about to test this for one of our strains):
Good luck!
Leave a comment:
-
coverage quantitative?
Hi kmcarr,
yes indeed its a bacterial (streptococcus) genome. I guess that since for transcriptomics the technique is quantitative so my guess it should be for genomes as well. Just wondered whether such "artefacts" might be the result of library generation and or the used chemistry for sequence analysis (I am not so much in those details).
Regarding the 4-5 fold more coverage I indeed expect the example to be a multiple duplicated gene (have to look for transposon elements and so) or indeed as you say extra chromosomal like plasmid. The example (I looked it up was a transposase of 1163bp and flanking 270 post and 100 pre). Another contains IS630-spnl. Most indeed transposable elements.... But how to find where they belong????? Looking at the assembly and flanking sequences/SNPs???
In addition I just finished some pre-processing of another strain this afternoon that also has an average coverage of roughly 20-35 times but has several larger regions that have coverage over multiple contigs in synteny with about 100 fold. More extreme even a region with even 800coverage (really smells like plasmid...?). Don't know yet what they contain but I hope to know soon.
So is it generally true that the coverage ratio approx resembles copy number per genome???
Thanks
ALex
Leave a comment:
-
Alex,
What is the ORF? This is obviously a bacterium (2 Mbp genome) and there are duplicated genes observed in bacteria, most notably the rRNA genes. It could also be extrachromosomal (i.e. plasmid) DNA.
Leave a comment:
-
How quantitative is contig coverage?
Hi All,
I am working on draft (resequencing) genomes of about 2Mb. In one run we multiplexed three strains and that works quite ok (as far as i can see now).
I noticed one thing which I find odd and which might be an artefact or meaningful observation.
When I plotted the sequence coverage (contrib to consensus) I noticed that in my draft which has on average about 20-25 fold coveragethat some contigs have extreme high coverage (over that complete contig). this can be as high as 125 fold
coverage. .... Aaah repeats..... no.
I looked at it and one example for instance is of size 1kb containing one complete ORF. But there are more regions similar to this.
Does this mean this ORF is about 4-5 fold present in the genome or is this a too simple conclusion. Could there be amplification or assembly artefacts???
Thanks for sharing any thoughts on this!
AlexTags: None
Latest Articles
Collapse
-
by seqadmin
The first FDA-approved CRISPR-based therapy marked the transition of therapeutic gene editing from a dream to reality1. CRISPR technologies have streamlined gene editing, and CRISPR screens have become an important approach for identifying genes involved in disease processes2. This technique introduces targeted mutations across numerous genes, enabling large-scale identification of gene functions, interactions, and pathways3. Identifying the full range...-
Channel: Articles
08-27-2024, 04:44 AM -
-
by seqadmin
Sequencing mRNA provides a snapshot of cellular activity, allowing researchers to study the dynamics of cellular processes, compare gene expression across different tissue types, and gain insights into the mechanisms of complex diseases. “mRNA’s central role in the dogma of molecular biology makes it a logical and relevant focus for transcriptomic studies,” stated Sebastian Aguilar Pierlé, Ph.D., Application Development Lead at Inorevia. “One of the major hurdles for...-
Channel: Articles
08-07-2024, 12:11 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 08-27-2024, 04:40 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
08-27-2024, 04:40 AM
|
||
New Single-Molecule Sequencing Platform Introduces Advanced Features for High-Throughput Genomics
by seqadmin
Started by seqadmin, 08-22-2024, 05:00 AM
|
0 responses
293 views
0 likes
|
Last Post
by seqadmin
08-22-2024, 05:00 AM
|
||
Started by seqadmin, 08-21-2024, 10:49 AM
|
0 responses
135 views
0 likes
|
Last Post
by seqadmin
08-21-2024, 10:49 AM
|
||
Started by seqadmin, 08-19-2024, 05:12 AM
|
0 responses
124 views
0 likes
|
Last Post
by seqadmin
08-19-2024, 05:12 AM
|
Leave a comment: