2 answers
I want to use bwa mem to align many RNA-seq reads into a reference genome.
I don't think that's the recommended aligner for this purpose. Please have a look at STAR.
I am a beginner and your support is appreciated.
Please show us what you tried. Many tools are well documented, and we rather help with a specific question. Let us know if something doesn't work, and show the command you used and error you get.
Hello, I am actually tempted to close this one since there seems to be no effort shown into reading existing threads and literature. Instead I will just link some good literature that you can go through, then hopefully being able to solve this yourself or ask a specific question:
RNA-seq for starters: https://www.annualreviews.org/doi/abs/10.1146/annurev-biodatasci-072018-021255
Edit: As OP seems to aim for assembly, see for example trinity: https://github.com/trinityrnaseq/trinityrnaseq/wiki
RNA-seq workflow at Bioconductor: https://bioconductor.org/packages/release/workflows/vignettes/rnaseqGene/inst/doc/rnaseqGene.html
Some buzzwords to google: Traditional alignment vs pseudo- and selective alignment. Involved tools: hisat2, STAR, featureCounts, HTSeq, salmon, kallisto, tximport. hisat2 and STAR are aligners, salmon and kallisto are selective- or pseudoaligners. All have extensive documentation including recommended syntax to run from command line.
Some tools for differential expression: DESeq2, edgeR, limma, sleuth (the latter downstream of kallisto), all have extensive manuals, the first three are part of the Bioconductor project, all are in R.
Some other buzzwords: Splice-awareness (this is respecting splice junctions during alignment, bwa-mem does not do that). Gene-level versus transcript-level differential analysis.
What you don't need unless you want to do transcript-level differential analysis (and even then there are probably more recent workflows): https://www.nature.com/articles/nprot.2016.095
As your post does not contain details, please first go through all that content, and if then something is unclear please comment.
Log in to answer this question.
-
BWA mem align reads with MQ = 0 to only one region (reference genome has two similar sequences).
written by captainlabman •
-
How to find reads that "jump" in the reference sequence
written by Gradient •
-
mapping long-reads to a reference library
written by mthm •
-
How do I make bwa-mem and bowtie2 less stringent?
written by gt •
-
Bacterial Seq Alignment
written by ravihansa82
-
How get read group information from fastq files in SRA
written by zamani_2012 •
-
SNP calling using samtools
written by ankit hinsu •
-
problem for setting BWA arguments for alignment
written by modarzi
-
Will the ordering of fasta references influence bwa mem alignment result?
written by hanschen •
-
Tophat error on using NONCODE database gtf file.
written by cfarmeri
Thank you so much.. I will check more information. I am trying to make a draft genome by bwa mem using RNA.seq data.. I was trying this command; "bwa mem ref.fa reads.fq > aln.sam" I was confused with "reads.fq", just only 1 read or can I put many reads here..? However, I hope I got it... I will try to make a 1 read file combining all reads..
You cannot make a genome out of RNA-seq, only a transcriptome. bwa-mem is in any case the wrong tool, it is a non-splice aware aligner. Aligners so to say "put back reads to a reference genome", but you need de novo assembly. You should google for de novo transcriptome assembly pipelines (not genome). Not my field at all but I hear people using something like trinity: https://github.com/trinityrnaseq/trinityrnaseq/wiki Be sure to read previous threads as well, this has for sure been asked before.
thank you... your comment is really helpful