Comments (13)
Hi,
if the process ends without any error and the .mem file is empty, then there could be two reasons behind this behavior:
- none of the input reads align to the splicing graph of that gene, maybe it is necessary to use different values of the parameters -l and -e or maybe most of the reads align to novel exons, that are currently not fully support by ASGAL
- there is some problem in the input data. For example, the format of some input file (such as the GTF) is different from the one expected by ASGAL
If it is possible, can you please share your data with me? It would be easier for me to troubleshoot the problem.
Edit: are you using the version on the master branch?
from galig.
Hi,
Thanks for the reply. I am sharing the data in fastq file reads (for the NEK1 gene sequences that I have extracted from alignment BAM file of the sample) which I have used in your tool as input for the program SpliceAwareAligner with the following annotation.gtf that I have attached whcih is having the coordinates of NEK1 gene exons, which I have extracted from human hg19. gtf.
command I have used:
/bin/SpliceAwareAligner -g hg19.fa -a annotation_NEK1.gtf -s NEK1_NM1.fastq -o asgal/NEK1_NM1_output.mem
python3 scripts/detectEvents.py -g hg19.fa -a annotation_NEK1.gtf -m asgal/NEK1_NM1_output.mem -o asgal/NEK1_NM1_output.events.csv
I need to check the AS events for this gene NEK1 only for the time being.
It would be a great help if you can test with these files.
Thanks in advance
NEK1_NM1.fastq.gz
from galig.
Ok, I think I figured out what the problem is. If you directly use the SpliceAwareAligner executable, you should not use the full genome as reference (ie the hg19.fa) but you should use the sequence of the chromosome the gene belongs to. In our case, you should use this file. I tried with this file and I got a non-empty .mem file (NEK1_NM1.mem.tar.gz).
Let me know if this helps you
from galig.
Thanks for the reply and initiative.
Sorry ,but I didn't see any output in that NEK1_NM1.mem.tar.gz file that you have attached.
It has a file called aligns.mem which is empty :(
I tried with the same chr4.fa after indexing it with samtools faidx, and used the same annotation_NEK1.gtf.gz file as gtf (attached earlier) but still that *.mem file is empty for me.
Not getting what is wrong ??
from galig.
Oh, I'm sorry... I attached the wrong file. This is the correct one.
In any case, it's very strange that you are still getting an empty .mem file. Can you please paste here the exact command you are running?
from galig.
Oh, thanks for the correct file. Here is the command that I am using
/bin/SpliceAwareAligner -g chr4.fa -a annotation_NEK1.gtf -s NEK1_NM1.fastq -o asgal/NEK1_NM1_output.mem
I have tried the next command detectEvents.py using the .mem file that you have sent, there also I have got a lot of error messages !! (attaching file)
error_msg_detectEvents.txt
python3 scripts/detectEvents.py -g chr4.fa -a annotation_NEK1.gtf -m asgal/NEK1_NM1.mem -o asgal/NEK1_NM1_output.events.csv
I hope I am using all the correct commands.
from galig.
It's strange that you are still getting an empty .mem file. Does SpliceAwareAligner return any error?
Regarding the detectEvents.py script, if a file annotation_NEK1.gtf.db
exists, try to remove it and rerun the script.
from galig.
No, SpliceAwareAligner did not return any error in nohup also. I have oe more sample to run with same command that you have used for SpliceAwareAligner.
For detectEvents now after removing gtf.db that .csv generated with headers line only.
Here is one more sample file I have to run for the entire set.
from galig.
This is the .mem file for the new sample. If the CSV is empty, it means that no events have been found. Maybe using a lower value for the parameter -l may change something
from galig.
Anyway, I have no idea why you obtain an empty .mem file. The command you use seems the right one... If I manage to understand why, I'll let you know
from galig.
Sure I will check to change the parameters, thanks for the help. Such a prompt help is really great help for my work. Please check if you can find some issue and let me know.
from galig.
Can you please send me the command that you are using for SpliceAwareAligner, actually I want to see the command as well as the parameters I want ot change -l and -e but I cannot do it in my system still its not giving me any output
from galig.
The command I used is:
./bin/SpliceAwareAligner -l 15 -g /data/amrita1983_git/chr4.fa -a /data/amrita1983_git/annotation_NEK1.gtf -s /data/amrita1983_git/NEK1_NM1.fastq -o /data/amrita1983_git/NEK1_NM1.mem
from galig.
Related Issues (16)
- Asgal (from run in your provided docker): doesn't find Transcript file, looks for gtf instead. HOT 1
- Docker permissions issue HOT 15
- .mem and sam files generated as empty HOT 45
- logs files HOT 2
- Salmon installation not found after pip installation in Asgal HOT 3
- Another permissions question HOT 6
- Adds tmp directory to asgal script
- Novel intron retention events in genome-wide mode HOT 2
- ALK ATI event HOT 9
- Error: Reference genome not found Halting... HOT 1
- Computation of PSI (percent spliced in) HOT 1
- Installation Error: ModuleNotFoundError: No module named 'Bio' HOT 7
- FAILED dependencies of target libtbb; tbb-2018_U3.tgz HOT 15
- Issues with Dockerfile HOT 6
- Error when testing with example files HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from galig.