Bowtie alignment for PE SOLiD data
Sorry my mistake: solution of simple:
bowtie -t -p 4 -C --sam Genome/mouce_database_10/mm10_c -f -1 /SOLiD_data/F3/reads/5500xl_F3.csfasta -2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.csfasta --Q1 /SOLiD_data/F3/reads/5500xl_F3.QV.qual --Q2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.QV.qual output1.sam 2> stdout.err
I think the question to use Bowtie to align PE SOLiD data should be easy, however, I cannot understand the error.
To perform PE alignment I am using the following code:
bowtie -t -p 4 -C --sam Genome/mouce_database_10/mm10_c -f -1 /SOLiD_data/F3/reads/5500xl_F3.csfasta -2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.csfasta -Q1 /SOLiD_data/F3/reads/5500xl_F3.QV.qual -Q2 /SOLiD_data/F5-RNA/reads/5500xl_F5-RNA.QV.qual output1.sam 2> stdout.err
• 2,332 views
•
link
written
by
nikkihathi •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Viral positive and negative strand with paired sequencing and bowtie
written by Lluís R.
-
Beagle: Index -l out of bounds for length 2
written by Tom •
-
error with BWA
written by storm1907 •
-
Trimming SOLiD adapters
written by bioinformaticsfilesdrive •
-
Tophat2 alignment for csfasta file
written by nikkihathi •
-
strandedness - RNASeq library preparation
written by Sashi Kiran Challa
-
Solid / Colorspace with bowtie : ERRO -> Extra parameter(s) specified
written by essepf •
-
Count number of mapped raw reads for each transcript in SAM file
written by manekineko
-
How can I extract and count unmapped reads using htseq count?
written by farah.dahalan •
-
How do I obtain uniquely mappable reads from bowtie on galaxy?
written by Jason2 •
