Question: "Error in splitted_reads[[1]] : subscript out of bounds"
0
gravatar for ignizamutdinov
3 months ago by
ignizamutdinov0 wrote:

I downloaded SRA file and converted it to SAM file:

sam-dump -1 SRR6040553 > SRR6040553.sam

and after to BAM file:

samtools view -S -b -T GRCh38_latest_genomic.fna SRR6040553.sam > SRR6040553.bam

But NIPTeR cannot read it file:

NIPTsample <- bin_bam_sample("SRR6040553.bam", do_sort = FALSE, separate_strands = F, custom_name = NULL)

The error is:

Error in splitted_reads[[1]] : subscript out of bounds

Thank you in advance!

bam nipter • 156 views
ADD COMMENTlink modified 3 months ago by ATpoint29k • written 3 months ago by ignizamutdinov0
0
gravatar for ATpoint
3 months ago by
ATpoint29k
Germany
ATpoint29k wrote:

What you downloaded are raw data in SAM format and contains no positional information. You have to align them first prior to any downstream analysis. Get the fastq files e.g. with this tutorial Fast download of FASTQ files from the European Nucleotide Archive (ENA) or use samtools fastq to get fastq from that SAM file and then use an aligner of choice, e.g. bwa, bowtie2, hisat2 for alignment. I guess the error here is that the script tries to access the first of the bins via [[1]] but as there are no positional records (and by this no bin counts at all) there is this out-of-bounds error.

ADD COMMENTlink modified 3 months ago • written 3 months ago by ATpoint29k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1858 users visited in the last hour