changing quality score of the read
Hi,
I want to change the quality score of the read, any suggestions would be appreciated.
Thanks
sequence
alignment
• 2,720 views
•
link
written
by
madhu.9124 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Coverage information of the read using samtool
written by madhu.9124 •Hi, I want to obtain coverage information of the read using samtool. I gave the following command- ./samtools stats -c mappedreversecomp.bam . I got an …
-
Meaning of the Output from BBMap
written by madhu.9124 •Hi, I ran command , ./bbmap.sh in=ecoliReads1.fastq outm=mappedadd_ecoli.bam outu=unmapped2.fa ref=ecoliRef.fasta. The output that was printed was: BBMap version 37.28 Retaining first best site only for …
-
Extension of reads
written by madhu.9124 •Hi, I want to extend the read to the 3' or 5' end of the read. Kindly suggest some command.
-
Edit bam files using Pysam
written by zaidnab •Hi there! I am trying to edit values in a bam file, using Pysam, or any other interface. I want to change the base quality, …
-
Extension of reads
written by madhu.9124 •Hi, I just want to extend the read by 20bp to the 3' or 5' end of the read, with high quality score, based on …
-
Correction of errors or mismatches in the read file
written by madhu.9124 •Hi, I am testing BBMap tool, so I want to correct any errors or mismatch in the read file. How do I go about this …
-
Obtaining Quality score of the mapped sequences - BBMap
written by madhu.9124 •Hi, I am new to BBMap, so I kindly request to help me to find the **quality score of the mapped sequences. I also want …
-
MEGAN - changing query coverage? HELP
written by tjf516 •Hello, Does anyone know if it is possible to change the parameters in MEGAN to only consider fragments with 100% query coverage? When attempting the …
-
Generating Read Length?
written by nikelle.petrilloHello! I am new to RNA seq. I quality trimmed my fastq sequences via fastX toolkit using a phred score of 30. I would like …
-
Fastq Quality Read And Score Length Check
written by thecuriousbiologist<p>I have a huge Fastq file and I want to check if the sequence length and the quality score length is the same for each …
Perhaps you can clarify this? Explaining what kind of data you have, as well as why and how you want to change the quality would be a good start. Also, what you mean by "quality score" (sometimes people mean base quality, and sometimes they mean mapq).
Previous questions by OP have been about MAPQ (Obtaining Mapping Quality ( MAPQ) value from samtool ). It is not clear what OP exactly wants to do.
I am basically testing BBmap application, so I am checking whether increase in quality score of the mapped sequence influence the output mapping. I have read and mapped with the reference ecoli read. I want to know increase the quality score of the mapped sequence.
Can I just ask - how is the quality score increased?
Increasing the quality score of reads will lead to an increased mapping rate, because BBMap ignores sequences with really low quality, and does not even try to map reads in which most of the kmers are almost certainly incorrect, according to the quality scores. Thus, you can inflate the mapping rate by assigning Q41 to all bases. The difference will usually be slight, but interesting if you have data with incorrect quality scores. Unfortunately, modern Illumina software does not report correct quality scores. The quality scores from modern Illumina software are garbage and should never be used without recalibration. In fact, I have tested their NovaSeq platform and the quality scores they assign to reads are actually worse than picking numbers at random.
thanks, which tool would you recommend to retrieve quality score and I want to manipulate the score as well.
"The quality scores from modern Illumina software are garbage" - do you have any review paper that looks into this?
found it - NovaSeq quality analysis