I have 100bp paired end RNA-Seq data from a mixed human/bacterial sample, I was wondering how RNA-Star would deal with this? Can it identify human reads with a high degree of accuracy from a mixed sample such as mine?I have been using tophat and bowtie2 to map to human and bacterial database respectively, but there are millions of sequences that are identified as being human if tophat is run first, and bacterial if bowtie2 is run first and given how tophat works, running tophat with one large database of human and bacterial would return the same result as doing bowtie first and then tophat.
--
You received this message because you are subscribed to a topic in the Google Groups "rna-star" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/rna-star/TqOdXiEFYrI/unsubscribe.
To unsubscribe from this group and all its topics, send an email to rna-star+u...@googlegroups.com.
Visit this group at http://groups.google.com/group/rna-star.
I should mention that I intend to analyze these resulting separated BAM files in HTSeq then DeSeq2.
On Friday, June 22, 2018 at 4:49:15 PM UTC-4, George Tollefson wrote:
Hi Alex,I see this is an old thread but don't see a newer one anywhere. I have a follow up question. I am following your instructions to generate a combined two species genome index file. My question is: once I have mapped my mixed sample rna reads to this combined genome, what is the best way to separate reads by species into new BAM files for future analysis. Does a tool exist? Will I have to write a script to split and save a new BAM file where the mapped reads for the second species begin?
Thanks,George
On Tuesday, April 14, 2015 at 6:18:27 PM UTC-4, Alexander Dobin wrote:Hi Jing,you need to make sure that chromosome names are different for the two species in the FASTA files, and the transcript names are different in the GTF files. The chromosome naming should be consistent between FASTA and GTF files for each species.After that, the FASTA files can be concatenated or simply listed one after another (space-separated) in the --genomeFastaFilesThe GTF files have to be concatenated into one GTF file.CheersAlex
On Monday, April 13, 2015 at 2:19:38 PM UTC-4, Jing Yang wrote:Hi Alex,If the reference genome is built by combing human and mouse genome. How should I create the index file with the annotation files? Should I combine two fasta files directly and two gtf files directly? Or I need two build a total new gtf file?
On Thursday, February 5, 2015 at 11:57:44 PM UTC+1, Alexander Dobin wrote:Hi Martin,I would use the same parameters you are using for the normal mouse mapping, since you cannot separately control the parameters for different "chromosomes". Any splice junctions you will see in the bacterial genome will be either mapping or library prep artifacts.CheersAlex
On Tuesday, February 3, 2015 at 9:39:18 AM UTC-5, martin....@gmail.com wrote:Hi Alex,Any pointers on how one should parameterize STAR for mapping mixed mouse+bacteria (the bacteria is known) samples to a combined genome given the difference in splicing?Thanks in advance!
--
You received this message because you are subscribed to the Google Groups "rna-star" group.
To unsubscribe from this group and stop receiving emails from it, send an email to rna-star+unsubscribe@googlegroups.com.
Visit this group at https://groups.google.com/group/rna-star.