Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Layla
    replied
    Hi Tim,

    It was a presentation I came across on the web by Stephan Trong on microbial genomes..viewed by downloading as a pdf. I must say this error was found using next gen sequencing on microbial genomes.

    Layla

    Leave a comment:


  • timread
    replied
    Originally posted by Layla View Post
    Dear 454 data analysers...
    ... but I then read in a paper that stated 64% of errors were assigned a Newbler QS of 64!

    Layla
    Layla - can you give the reference for this statement? Thanks in advance,

    Tim

    Leave a comment:


  • Layla
    replied
    Dear 454 data analysers...
    We have data generated from 454 using nimblegens capture array experiment. I have read that 454 has issues with AT rich genomes and homopolymers. Having looked at the HCdiffs file I was going to focus on variances with Quality scores > 40, but I then read in a paper that stated 64% of errors were assigned a Newbler QS of 64! Does anyone know how to calculate the error rate from 454 data using the files created from gsMapper or if this error rate can be seen in any of the generated files? Ideally I would like to see how many of the variances seen are "real".

    New to sequence data analysis!
    Thankyou
    Layla

    Leave a comment:


  • timread
    replied
    The thresholds are calculated for each run by the analysis software based on the signal distribution. If you look at the 454BaseCallerThresholds.txt file in the Analysis (D_) directory you can see some stats for the run. e.g.

    distributionPeaks
    {
    signalPeak = 0, 0.12, Found;
    signalPeak = 1, 0.98, Found;
    signalPeak = 2, 2.00, Found;
    signalPeak = 3, 2.98, Found;
    signalPeak = 4, 3.92, Found;
    signalPeak = 5, 4.84, Found;
    signalPeak = 6, 5.68, Found;
    }

    thresholdsUsed
    {
    threshold = 0, 1, 0.68;
    threshold = 1, 2, 1.58;
    threshold = 2, 3, 2.52;
    threshold = 3, 4, 3.48;
    threshold = 4, 5, 4.44;
    threshold = 5, 6, 5.32;

    interpolationAmount = 0.94
    }


    For a good description of 454 errors see this manuscript ..

    Leave a comment:


  • Alex Clop
    replied
    454 sequences by Pyrosequencing. As defined by Nature Reviews Genetics Glossary, Pyrosequencing is A DNA sequencing technique that relies on detection of pyrophosphate release on nucleotide incorporation rather than chain termination with dideoxynucleotides.

    This chemistry i) incorporates a known nucleotide to the sequencing reaction at a time, ii) then eliminates the remaining nucleotides and iii) then incorporates another known nucleotide again. It does this for the 4 nucleotides sequentially and then the cycle is repeated again. Each time a nucleotide is incorporated at the elongating chain, a luminiscent reaction (luciferase) will occur and detected by a CCD camera. Thus, and will permit the basecalling and the sequence "reading".

    Regarding your question, the intensity of the luciferase ("I") activity is proportional to the number of - identical - nucleotides incorporated ("N") in the synthesized chain. However, if "N" is too high (I don't think there is a well defined threshold), "I" will reach saturation which will lead to wrong "quantification" of the number of contiguous A, T, C, or G that have been incorporated.

    It is not simple to explain but it may help.

    Leave a comment:


  • bioinfosm
    replied
    As I understand..
    454 does not have reverse end-blockers to nucleotides and attach multiple bases at once, using the total intensity to determine how many bases may have been incorporated.

    Others use blockers to make sure only one nucleotide goes at a time, and make it a simple yes/no question with less chance for error

    Leave a comment:


  • vasvale
    replied
    below is a reference for the PCR error, can somebody eplain why homopolymers error are seen with Roche and not with other platforms (as I heard)

    Clarke LA, Rebelo CS, Gonçalves J, Boavida MG, Jordan P.
    PCR amplification introduces errors into mononucleotide and dinucleotide repeat
    sequences.
    Mol Pathol. 2001 Oct;54(5):351-3.
    PMID: 11577179 [PubMed - indexed for MEDLINE]

    Leave a comment:


  • glacerda
    replied
    I have compared a segment of 100K that we have sequenced by sanger (at > 20X) and by 454 ( > 50X)

    Even at 50X, 454 made indel errors (all of them in homopolymers greater than 5 bp)

    in you case, you have a 4bp homopolymer in which an indel error would be very very rare. I also thin this is a pcr issue or a real frameshift.

    Leave a comment:


  • timread
    replied
    Ian -

    IMHO, when you have a ratio of 200 read to 1 it is very unlikely that you are seeing homoploymer miscalls. Typically this occurs when the base calling software has difficulty accurately setting the threshold between say 3C and 4C for a particular extension and you will see both reads well represented. I'd go for a PCR problem, or more interestingly, a real frameshift. Verify the result with a Sanger read.

    tim

    Leave a comment:


  • ian Adams
    started a topic 454 homopolymer errors or????

    454 homopolymer errors or????

    I have some 454 data produced from cDNA randomly amplified using proofreading taq, sheared and then sequenced. Once the amplification tags were removed, we'ved assembled it using CLC genomic work bench. Looking at the protein translations of two particular genes we have very good homology to genes from related organisms but in two places we appear to have frame shifts in the DNA sequences which introduce stop codons leading to incomplete proteins. Correcting the frameshifts would be CCC to CCCC in one case and CCCC to CCC in the other and completely restore the protein homology. Both regions have over 200 fold coverage reading the frameshift with only one read giving the "corrected" unshifted sequence. I have read that the 454 does have problems with homopolymers. What I am wondering is if this is biological (a mutation leading to an unusual terminated protein--I think unlikley these proteins are essential) a sequencing error (introduced by the 454 / emPCR) or sample prep (PCR error-but I was using proofreading taq). Does anyone have any experience of this sort of thing and any comments on the "introduced by the sequencing / emPCR option"?
    We are going to Sanger sequence across these regions and see what we get but that will take time.

    thx

    Ian

Latest Articles

Collapse

  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM
  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:12 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-23-2024, 04:11 PM
0 responses
67 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
73 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
62 views
0 likes
Last Post seqadmin  
Working...
X