Hello all
I have a list containing left mate sample reads and a list for right mate sample reads like: list1 = sample11.fastq , ... samplen1.fastq list2 = sample12.fastq , ... samplen2.fastq
I know I should prepare a .sh file for running tophat over my list at the same time, but do not know what is the efficient way. I mean should I write a script like: for i in list1 for j in list2 tophat bwtie_index i j done
and then call my script, is it the best way of doing it? many thanks in advance
2 answers
You can try our scripts at:
Yes this is the easiest way I think ;) don't forget to specify a different output directory for each of your job. And also be carefull of your memory (tophat use several Gb per run depending on the size of the genome and the size of your fastq files)
what is your hardwar configuration ?
Log in to answer this question.
-
Why extracting reads from BAM to fastq produces small amount of reads?
written by artemd •
-
Creating text file as a temporary file in python
written by mdsiddra •
-
Count right and left mate pairs from Sam files
written by PAn •
-
Pipeline for >100 RNA-Seq PE samples
written by umn_bist •
-
Venn Diagram Help
written by kevluv93
-
Expected number of read duplicates
written by dariober
-
Automating Paired End mating
written by Furor •
-
Does The New Tophat2 - Bowtie2 Pipeline Really Map 100% Of The Rna-Seq Reads ? Is This Real ?
written by biorepine
-
Make Paired File Task-Scripture
written by GPR
-
Databases For Gene Family?
written by Dejian
