tophat 2 rna seq
hi all. i have a fastq file with total sequences 28020920 but i only want to do 10.000 f.e. . which option is there because i cant seem to find that in the manual.
• 1,806 views
•
link
updated
by
Oskar •
•
written
by
dimitrischat
•
written
by
dimitrischat
1 answer
Hi Dimitris - I guess you want a small sample of your reads for testing and debugging purposes. If so, you can create a file containing a small number of reads. For example:
$ zcat myreads.fastq.gz | head -400000 | gzip > Test100k.fastq.gz
which takes the first 100k reads from “myreads” and stores them in “Test100k”
Hope it helps!
Log in to answer this question.
More posts like this
-
Extracting VCF of all variants in TCGA exomes from all cancer types
written by j.lunger18 •
-
remove duplicates with markdup
written by dimitrischat
-
Blast many read sequences
written by samuel
-
Remove hundreds of probes
written by Floydian_slip
-
cuffdiff / htseq-count
written by dimitrischat
-
Reduce/Select reads/lines ffrom files using terminal
written by dimitrischat
-
Retrieve Trembl/Uniprot previous version
written by leon •
-
Exclude sites with a certain fraction of heterozygosity from vcf
written by thal •
-
I Have A Problem With My Proteomics Data Analysis.
written by 9921 •
-
What Script (Or Program) Should I Use To Filter Low Complexity And Repetitive Reads From Rna-Seq?
written by Geparada
Not sure why you want to do that but you could use
reformat.shfrom BBMap suite to sample the 10000 reads into a new file and then use that.thank you. but cant you do it using tophat?
You can check the manual but I don't think tophat has an option to sample a fraction of reads.
Hi dimitrischat,
It's worth noting that TopHat2 has been, essentially, deprecated by the developers, who recommend using HISAT2 instead. Unless you have a very specific reason to adopt TopHat2 in your pipeline, it's probably best to follow their advice.