Header Leaderboard Ad

Collapse

MapSplice2 does not recognize my Bowtie indexes

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Marcha
    started a topic MapSplice2 does not recognize my Bowtie indexes

    MapSplice2 does not recognize my Bowtie indexes

    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

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 06-01-2023, 08:56 PM
0 responses
11 views
0 likes
Last Post seqadmin  
Started by seqadmin, 06-01-2023, 07:33 AM
0 responses
14 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-31-2023, 07:50 AM
0 responses
52 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-26-2023, 09:22 AM
0 responses
56 views
0 likes
Last Post seqadmin  
Working...
X