Entering edit mode
7.9 years ago
bandanaschapagain
▴
40
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
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
--outReadsUnmapped
parameter 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
Thats the way
--outReadsUnmapped
param is used. You have to sayFastx
, and not an output file name. Refer to STAR manual plzCan 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_in
command. 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_INDEX
bbmap (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!