I got a few fastq files from RNA-seq experiment and somebody ask me to check the number of reads that come from retroelements L1-ORF1/2p. Is it a simple way to perform blast of the whole fastq file agains the sequence of the elements (instead of align it to the reference genome using hisat/bowtie/bwa etc) and check the number of reads in each sample? I can try a grep but it can take a long time. Thank you for advise.
2 answers
You are better off with NGS aligners like bwa. Reasons:
1) They would be faster than blast / grep.
2) They are quality aware.
3) grep cannot do a fuzzy search. So reads with sequencing error or polymorphism will not be found with grep. grep is terribly slow for NGS searches unless you optimize the search params and restrict the search to specific LOCALE.
4) It's supereasy to do it! (just index the spliced sequenced and run bwa mem, for example)
Since you're interested in a repeat, you might find TETranscript useful. That will allow more accurate quantification. You might also try TECounts from my pull request on that repository.
Log in to answer this question.
-
Uniquely mapped reads on Bowtie for colorspace RNA-seq
written by IsiI •
-
bowtie (1) 0% alignment on paired-end RNAseq data
written by gulamaltab •
-
bowtie reported reads and aligned reads number doesn't match
written by Amaranta_Remedios •
-
Number of reference sequence contigs in sam/bam header does not match the number in my fasta
written by marcusanaymik •
-
miRNA blast to RNA database?
written by Buffo
-
Run mutltiple Fastq files on Hisat2
written by V
-
Reversing R1.fastq and R2.fastq
written by L. A. Liggett
-
Aligning Paired-End Reads In Single-End Mode
written by GPR
-
How To Generate A Paired End File Suitable For Bwa
written by GPR
-
Aligment Of Reads To A Bisulfite Converted Genome
written by Tonig
Since you are talking about RNA-seq you want a splice-aware aligner so bwa and bowtie will not suffice. HISAT is a good one, alternatives are STAR and bbmap.
Thank you. I'll try to run it after indexing. By the way is it possible to run bam files against indexed sequence?
Please use
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized.If you only need a rough idea of the reads that could be from that transcript you may be able to use kallisto or Salmon to do pseudoalignments.
Thanks. Could you please update the link to instructions?
Again. Please use
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized. I moved your post now.