Cap3 Assembly

Does CAP3 always allocate a Q20 score to all the bases in the final assembled contig when assembling a forward and a reverse read whose original base quality scores were high( above 20)? or I am missing something?

assembly

• 7,157 views • link

updated 13.8 years ago by Ying W

• written 13.8 years ago by hadasa

1 answer

Are you specifying qualify values correctly for cap3?

For fasta files, the quality values must be in a separate file of the type fileoffilenames.qual.

From the cap3 documentation: http://staden.sourceforge.net/manual/gap4_unix_88.html#SEC161

EDIT: actually that link seems a bit outdated but the documentation link on their man site is down. Luckily its been archived: http://web.archive.org/web/20100423111938/http://deepc2.psi.iastate.edu/aat/cap/capdoc.html

• 1 views • link

written 13.8 years ago by Ying W

Yes the quality files were specified as outlined in the manual.

• 0 views • link

written 13.8 years ago by hadasa

More posts like this

FastQC interpretation - 16S sequencing
written 2.2 years ago by Christopher •

I'm new to the field of bioinformatics and have been reading a lot about 16S and RNA-seq analysis, as well as FastQC interpretation. I understand …
Python program that reads a fastq file and calculate Phred score
written 5.3 years ago by sherrygeek •

Hi! I need help with solving the following question with Python3: Write a Python program that reads a fastq file and calculate how many bases …
Use pre-assembled paired-end reads as single-end for metagenome assembly with metaSPAdes
written 6.5 years ago by somnja •

Hey, i am tasked with assembling a metagenome. From my predecessor i got in fasta format one file that contains: 1) normal paired-end reads (forward …
Standard value of quality score and sequence length in QC trimming of RNA-Seq
written 7.2 years ago by takoyaki

I want to know the standard value of quality score or sequence length in trimming by fastx_toolkit. In my RNA-Seq textbook, author trim fastq sequences …
What is a good Trimmomatic parameter?
written 8.2 years ago by Arindam Ghosh

I am trying to clean RNAseq reads with TRIMMOMATIC using the parameters as: java -jar trimmomatic-0.36.jar PE -trimlog file_trim_log input_1.fastq.gz input_2.fastq.gz output_1P_clean.fq.gz output_1U_clean.fq.gz output_2P_clean.fq.gz output_2U_clean.fq.gz …
Changing samtools mpileup I16 calculations?
written 8.7 years ago by emilywynn6 •

Hello, I am trying to filter some vcf files from mpileup by using the average base quality of reference bases and the average base quality …
Vcf files fields: AD vs DP4
written 9.9 years ago by Rob

In Vcf file, generate by samtools mpileup, there is an AD field (Allelic depth). This field is suppose to replace the previous DP4 field (Number …
RNA seq FASTX quality trimming
written 10.4 years ago by Rahul •

Hello, I have filtered my illumina pair-end reads (Forward lib-24 million reads, Reverse Lib 24 million) using FASTX_Quality_Filter by applying the Q20 score to 90 …
Combining quality scores when aligning/assembling
written 12.0 years ago by earonesty

Suppose I have 2 reads, and I align them as a part of a joining program or an assembly: ``` AATGCGTTTGA GCGTTAGATGCTGA ________^ ``` Then …
Read To Consensus Quality Scores
written 15.8 years ago by Lee Katz

<p>Hi, what is the actual formula for finding the quality score in a consensus sequence base? </p> <p>I'll give a hypothetical example: in one base …