Header Leaderboard Ad

Collapse

Error when bbmapping RNA-seq data

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Error when bbmapping RNA-seq data

    Hello

    I was using BBMap already for successful mapping of RNA-seq data (paired-end, Illumina) to a multi-FASTA transcriptome. Now I tried to map the very same data to the human GRCh38 genome (downloaded from NCBI human genome resources). I used the following command:
    bbmap.sh ref=$mapping_ref build="$BBmap_ref_ID" in=$BBduked_reads1 in2=$BBduked_reads2 \
    outu=unmapped_"$name1"_"$date".sam outm=mapped_"$name1"_"$date".sam \
    maxindel=200k mdtag=true sam=1.4 subfilter=3 pairedonly=true
    Unfortunately, I got the following error:
    [...]
    Retaining first best site only for ambiguous mappings.
    Writing reference.
    Executing dna.FastaToChromArrays2 [/scratch/hpc-prf-ptma2/ptma2001/bbmap/GRCh38_genome_16032022.fna, 6, writeinthread=false, genscaffoldinfo=true, retain, waitforwriting=false, gz=true, maxlen=536670912, writechroms=true, minscaf=1, midpad=300, startpad=8000, stoppad=8000, nodisk=false]

    Set genScaffoldInfo=true
    Writing chunk 1
    Writing chunk 2
    Exception in thread "main" java.lang.NumberFormatException: For input string: "ӌ�jHN�H��¬���l��,ΚK��/S��ͩm���T
    �!�l��X��o�ã’*s$�s^��m�h
    [email protected]�
    �m��$֧����@�H,���xI%_ Uz�?� S�B�2[*i�_����)i�NEր�ٜ�?Z�i£�� "
    at java.lang.NumberFormatException.forInputString(NumberFormatException.java:65)
    at java.lang.Integer.parseInt(Integer.java:580)
    at java.lang.Integer.parseInt(Integer.java:615)
    at dna.Data.setGenome2(Data.java:1016)
    at dna.Data.setGenome(Data.java:769)
    at align2.BBMap.loadIndex(BBMap.java:316)
    at align2.BBMap.main(BBMap.java:32)

    Would you have any idea, why this error comes up and how I could solve it?

    Thanks a lot in advance and have a great day,

    Elena

Latest Articles

Collapse

  • seqadmin
    Targeted Sequencing: Choosing Between Hybridization Capture and Amplicon Sequencing
    by seqadmin




    Targeted sequencing is an effective way to sequence and analyze specific genomic regions of interest. This method enables researchers to focus their efforts on their desired targets, as opposed to other methods like whole genome sequencing that involve the sequencing of total DNA. Utilizing targeted sequencing is an attractive option for many researchers because it is often faster, more cost-effective, and only generates applicable data. While there are many approaches...
    03-10-2023, 05:31 AM
  • seqadmin
    Expert Advice on Automating Your Library Preparations
    by seqadmin



    Using automation to prepare sequencing libraries isn’t a new concept, and most researchers are aware that there are numerous benefits to automating this process. However, many labs are still hesitant to switch to automation and often believe that it’s not suitable for their lab. To combat these concerns, we’ll cover some of the key advantages, review the most important considerations, and get real-world advice from automation experts to remove any lingering anxieties....
    02-21-2023, 02:14 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-17-2023, 12:32 PM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-15-2023, 12:42 PM
0 responses
17 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-09-2023, 10:17 AM
0 responses
66 views
1 like
Last Post seqadmin  
Started by seqadmin, 03-03-2023, 12:03 PM
0 responses
64 views
0 likes
Last Post seqadmin  
Working...
X