How to write a program in Perl that read a FASTQ file? after then, there is a method to convert the quality score from ASCII to numeric?
Hi dears, i'm a new student of bioinformatic, this is my first approach on this field, i will be very grateful if someone of you can help me with this question. If you have problem to write the entire exercise i will be grateful if you can write me some tips. King Regards, Tommaso from Italy.
• 2,849 views
•
link
written
by
prototommaso •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How should we replace NAs from a genelight object?
written by jaafari.omid
-
Allele frequency comparison between pairs of populations
written by jaafari.omid
-
GC bias correction for Chip-SEQ
written by Farron16 •
-
How to merge different BAM files?
written by jaafari.omid
-
Blog: Micro Binfie Podcast
written by Lee Katz
-
Fst Visulization from WGS data
written by jaafari.omid
-
vcfR for Whole Genome Data
written by jaafari.omid
-
Differences between lfc=log2(1.5) and lcf=1.5
written by Joe Kherery
-
Forum: Statistics graduate wants to do some bioinformatics analysis
written by shanjida986 •
-
Bowtie - Mapping
written by bioinfo_student •
Could you be more descriptive about your question. a example of Input and Output would suffice. As you're saying that it's your beginner, I would recommend you to use bioperl.
Buongiorno Tommaso, this should have been a Question, not a Tutorial. I have moved it to a Question.
Thank for the answers, The input is a nucleotide sequence, as output i need a rode FASTQ file [example at the link]
https://www.safaribooksonline.com/library/view/bioinformatics-with-r/9781783283132/graphics/3132OS_08_03.jpg
You can't get a FASTQ from an input nucleotide sequence. There is no quality information in that input (unless you have data you haven't told us about). You could create a script which would emulate a FASTQ, and give you some dummy quality scores, but why you'd want to do that I don't know (even as a test assignment :/ )
This looks like an assignment question and will likely be closed. We discourage these because while someone will give you a quick answer it will prevent you from learning the basics.
The actual conversion of the scores can be done based on the information you can find in the WikiPedia page for FASTQ format. You should know how to write the rest of the stuff (reading files, parsing lines, checking to see if the line is has quality scores, doing conversion and writing result out).
If this is not an assignment, consider bioperl libraries in parsing files
If this is not an assignment, consider biopython libraries in parsing files
;-)
Perl might have been an amazing language for bioinformatics, but frankly, I don't think there is a good enough reason to get started with it now.