commandlines and sequences
Guy, I have a fastq file and I need to write a command to shorten each sequence to just 10 nucleotides while keeping the sequence ID unchanged. I have been trying this for days but with no success. I need help from some gurus here. For your info, a fastq file consists of the seq ID starting with @, the sequence, then a separator between the sequence and it its quality scores and sequence quality scores which in my file start with 2. Thanks
• 179 views
•
link
written
by
Baby Bioinformatician
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Is it Acceptable to Have Uniform Quality Scores in a FASTQ File?
written by Κοσμάς •
-
Is it normal to obtain a perfectly dark blue "Per tile sequence quality" plot in fastqc?
written by cg1440 •
-
All against all structure alignments
written by Baby Bioinformatician
-
How to extract the last 1000 nt from a group of sequences in a FASTA file?
written by Andrei •
-
commandlines and sequences
written by Baby Bioinformatician
-
Scripting method to parse fastq diploid consensus sequence bins based on quality
written by memory_donk
-
How to visualize sequence quality control distribution
written by faizlotfy •
-
problem with newbler(error)
written by midox
-
How To Extract Set Of Reads From Fastq (Or Eventually Fasta And Qual) Based On List Of Ids?
written by Biomonika (Noolean)
-
Trimming A Fastq File Representing Rna-Seq Data
written by Varun Gupta
And what have you tried? And what was the problem?
fin swimmer
Could you perhaps write a more specific title for your thread? This one doesn't say anything about the content.
Did you delete your previous thread commandlines and sequences after receiving suggestions there?
Hello Charles Olweny !
Questions similar to yours can already be found at:
We have closed your question to allow us to keep similar content in the same thread.
If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.
Cheers!
You can use
reformat.shfrom BBMap suite to do this.If you have paired-end data then you would want to use something like this