You never know whether DNA sequences are forward or reverse until after you map them; the mapper will figure that out on a per-read basis. You don't need to do anything special.
The usage is pretty straightforward:
Index:
bwa index ref.fa
Map single-ended reads:
bwa mem ref.fa reads.fq > aln-se.sam
...or, map paired reads:
bwa mem ref.fa read1.fq read2.fq > aln-pe.sam
-
bwa-mem Roche454
Goodmorning,
(sorry for my bad english), I would use bwa-mem to map sequences against human genome. I have the file fast.q of Roche454 Junior, in this there is not specified if the sequence is in forward or in reverse. In this case, what is the option which allows to use BWA-MEM?
I really hope someone from seqanswer can help me.
Thank you in advance.
-
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...07-08-2024, 03:19 PM
Leave a comment: