Hi everybody,
I wannted to start working with MapSplice2, however, I run into the problem that it does not recognize my Bowtie1 indexes (build using v1.1.1).
I am using this command:
python ${PROGRAMDIR}MapSplice-v2.2.1/mapsplice.py -c ${GENOMEDIR}MapSplice -x ${GENOMEDIR}MapSplice/bt1_GRCh38_apr2019 -1 ${READ1}${SUFFIX2} -2 ${READ2}${SUFFIX2} --threads 6 --output ${align_DIR_MapSplice}${READ1:${#SOURCEDIR}:-${#SUFFIX1}-${#PAIR1}}MapSplice.sam
This code makes mapsplice2 run:
[Thu Jul 4 06:32:24 2019] Beginning Mapsplice run (MapSplice v2.2.1)
[Thu Jul 4 06:32:24 2019] Bin directory: /share/tools/MapSplice-v2.2.1/bin/
[Thu Jul 4 06:32:24 2019] Preparing output location /share/analysis/CEFIC/CEFIC_testdata/HeCaToS_APA/Trimmed_reads_trimmomatic/MapSplice2/APA_MapSplice.sam/
[Thu Jul 4 06:32:24 2019] Checking files or directory: /share/analysis/CEFIC/CEFIC_testdata/HeCaToS_APA/Trimmed_reads_trimmomatic/APA_The_000_2_R1.fastq
[Thu Jul 4 06:32:24 2019] Checking files or directory: /share/analysis/CEFIC/CEFIC_testdata/HeCaToS_APA/Trimmed_reads_trimmomatic/APA_The_000_2_R2.fastq
[Thu Jul 4 06:32:24 2019] Checking files or directory: /share/tools/data/genomes/GRCh38_r96_april2019/MapSplice/
[Thu Jul 4 06:32:24 2019] Checking Bowtie index files
[Thu Jul 4 06:32:24 2019] Building Bowtie index for reference sequence
But it can not find my bowtie1 indexes and tries to build them. Strangely, it actually starts overwriting the indexes I have already created. After about 6 hours, building the indexes fails, returning this error:
Error: Building Bowtie index files failed /share/tools/data/genomes/GRCh38_r96_april2019//MapSplice/bt1_GRCh38_apr2019.*
What is wrong with my code?
How do I get MapSplice to recognize the indezes I build myself?
Thank you already,
Marcha
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
Yesterday, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
58 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
53 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
45 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|