Multiple Sequence Alignment of Smaller Groups Within a Larger File
Hello,
I'm trying to align small groups of RNA-Seq reads within a larger FASTQ file. I need a separate alignment for each group. I can separate the read groups with an ID number or something similar, but it's not feasible to separate the groups into their own files, since there are around 10^6. I only have one reference sequence I'm trying to align to, and there isn't any alternative splicing.
Is there something that I can use for this? I've looked at the manuals for EMBOSS and Clustal O, but they don't seem to have anything appropriate. BWA has the option "-R" for setting the read group id, but I think that's just for output.
Thanks in advance for any help.
• 1,380 views
•
link
written
by
newsome6 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Finding the expression levels of non-canonical RNA splicing
written by Afnan Nuruzzaman •
-
FastQ ColorSpace to BaseSpace format in 2018
written by ZheFrench
-
Kissplice: Problem getting SNPs in test data
written by WeepingMeadow •
-
Altanalyze for Arabidopsis thaliana
written by mail2steff
-
BQSR on one sample vs. all samples in a run
written by ari.nazarian •
-
How do I find input arguments (RGPU, RGLB, etc) for AddOrReplaceReadGroups
written by Max
-
Splicing programs to find alternative polyadenylation
written by kiminsigne •
-
Downloading multiple sequence alignment as clustal format file from clustal omega.
written by siddharth.avadhanam •
-
Publicly Available Rna-Seq Data For Study Splicing Change In Human Brain
written by michealsmith
-
Mostly Used (Fast) Global Pairwise Alignment Tool
written by Pdubois •
You could use
filterbyname.shfrom BBMap, search with that name on this page to see usage OR faSomeRecords utility from Jim Kent. Both will allow you to pull out subsets of records from your large file on demand.I will try that, thank you so much!