Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

error in runing map_genes.sh #127

Open
rdf1993 opened this issue Oct 25, 2023 · 1 comment
Open

error in runing map_genes.sh #127

rdf1993 opened this issue Oct 25, 2023 · 1 comment

Comments

@rdf1993
Copy link

rdf1993 commented Oct 25, 2023

bash map_genes.sh --tr1 /mnt/volume3/rdf/01_refData/GRCh38_fa/Homo_sapiens.GRCh38.dna.primary_assembly.fa --t1 nucl --n1 hs --tr2 /mnt/volume3/rdf/01_refData/GRCm38_fa/Mus_musculus.GRCm38.dna.primary_assembly.fa --t2 nucl --n2 mm --threads 8
Running tblastx in both directions
map_genes.sh: line 124: 308696 Segmentation fault (core dumped) tblastx -query $tr1 -db $tr2 -outfmt 6 -out "maps/${n1}${n2}/${n1}to${n2}.txt" -num_threads ${n_threads} -max_hsps 1 -evalue 1e-6
map_genes.sh: line 124: 310032 Segmentation fault (core dumped) tblastx -query $tr2 -db $tr1 -outfmt 6 -out "maps/${n1}${n2}/${n2}to${n1}.txt" -num_threads ${n_threads} -max_hsps 1 -evalue 1e-6

@jylee43
Copy link

jylee43 commented Dec 12, 2023

Hi, you need to use sequences from mRNAs or proteins not chromosomes. "*.dna.primary_assembly.fa" are for genome sequences, so they are not appropriate for SAMap.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants