For RNA seq data analysis
Hi all, I need to identify a software to be used for aligning my reads to the reference genome and getting out the reads that do not align to the reference genome. I have no idea for the softwares available. I hope I can get help. Regards, Bandana
• 3,391 views
•
link
updated
by
Ming Tommy Tang
•
written
by
bandanaschapagain •
•
written
by
bandanaschapagain •
1 answer
my github repo may help you https://github.com/crazyhottommy/RNA-seq-analysis
• 0 views
•
link
written
by
Ming Tommy Tang
Log in to answer this question.
More posts like this
-
Aligning Contigs of various Mtb strains to a reference genome and get the variants
written by Paul
-
Understanding Lowfreq Output
written by bandanaschapagain •
-
For Viral Recombination Identification
written by bandanaschapagain •
-
Bowtie Sam output analysis
written by bandanaschapagain •
-
To find the correct paramter
written by bandanaschapagain •
-
How to get percentage of reads mapped to each Reference Genome
written by bandanaschapagain •
-
To download the sequences by using NCBI entrenz direct
written by bandanaschapagain •
-
Using Cufflinks with SRA Reads
written by bandanaschapagain •
-
For Basespace data download
written by bandanaschapagain •
-
For cufflinks error of using gff files
written by bandanaschapagain •
Have a look here -
Current state-of-the-art for RNA-Seq gene / transcript expression quantification OR just the tools of preference
and specifically for getting reads that do not align, use the
--outReadsUnmappedparameter if you use STAR.Thank you for providing such useful link.
I am using the below command for STAR. Two surprising things for me are: 1.Time taken by the STAR for read aligment 2. I am not getting the unmapped reads even by using
--outReadsUnmapped. I am confused if there something I need to modify for my command here.Please help me with this. Regards, Bandana
Can you be more specific about the time taken by STAR? Was it surprisingly short or long? How long? How many reads is the input?
I am using paired reads to be aligned to the reference genome.My reads are of 2.9 GB. It is very fast for aligning the reads and also when I check the log file, it says that there are no reads that did not match to the reference genome. for building the reference index I used
./STAR --runMode genomeGenerate --genomeFastaFiles /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/Brome_mosaic_virus/Brome_mosaic_virus_uid15052/NC_002026.fna --sjdbGTFfeatureExon /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/Brome_mosaic_virus/Brome_mosaic_virus_uid15052/NC_002026.gtf --genomeDir /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/BMV_002026_incommand. Please Help me with this. Regards, BandanaSince you are just generating the STAR index with this command there were no reads to align. Now you are ready to do the actual alignment using STAR (step 3 in STAR manual).
I am doing the read alignment using step 3 but its giving me error as: segmentation fault core dump/ I am really getting confused since so long. Please anyone can help me with this.
How much memory do you have available? Since you were able to generate genome indexes your STAR install appears to be working correctly. Can you provide the command you are using for the alignment?
I am using this command for the alignment.Please help with this.I have no idea how to correct it.
./STAR --genomeDir /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/Brome_mosaic_virus/Brome_mosaic_virus_INDEX --readFilesIn /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/RNA_Virus_Josef/reads/BMV/JB-BMV_S3_L001_R1_001.fastq /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/RNA_Virus_Josef/reads/BMV/JB-BMV_S3_L001_R2_001.fastq --outReadsUnmapped Fastx
That command looks ok, assuming your indexes have been created properly. Can you list output of
ls -sh1 /home/yyin/RNA_Virus_Josef/RNA_Virus_Josef/Brome_mosaic_virus/Brome_mosaic_virus_INDEXThats the way
--outReadsUnmappedparam is used. You have to sayFastx, and not an output file name. Refer to STAR manual plzbbmap (splice aware aligner) from BBMap suite with the
outu=option will capture reads that do not align.Thank you for your time and response.
Please try to use a more informative title for your question, e.g.
RNA-seq aligner and extracting unaligned reads.Thank you for your suggestion.I will improve my post,for sure!