How to extract all sequences mapped to a transcript from Kallisto output
I ran Kallisto with the --pseudobam option. How do I extract all the short reads that are mapped to a single transcript (e.g. ENST00000367969.8)? As a person without any previous SAM/BAM experience, I tried the following things without success. Really appreciate some help here.
I first tried:
samtools view -S pseudoalignments.bam ENST00000367969.8
but got the following error:
[E::idx_find_and_load] Could not retrieve index file for 'pseudoalignments.bam' [main_samview] random alignment retrieval only works for indexed BAM or CRAM files.
I tried to index the file using:
samtools index pseudoalignments.bam
but get the following error:
[E::hts_idx_push] NO_COOR reads not in a single block at the end 79729 -1
[E::sam_index] Read 'A00742:184:HLTF2DSXY:4:1101:4417:1000' with ref_name='ENST00000556431.1', ref_length=4142, flags=83, pos=2098 cannot be indexed samtools index: failed to create index for "pseudoalignments.bam"
1 answer
Log in to answer this question.
More posts like this
-
samtools index error - "Unsorted positions on sequence" even though file is sorted
written by Geoffrey •
-
samtools index: failed to create index. ([E::hts_idx_push] NO_COOR reads not in a single block )
written by DareDevil
-
BRAKER for large genome?
written by Dave •
-
bam index chromosome size limit
written by николай •
-
Is my output .bam file from STAR ok?
written by elee228 •
-
samtools: Numerical result out of range
written by Gregor Rot
-
MethylDackel Error running on HPC server
written by magicridge •
-
samtools index error
written by pt.taklifi
-
cant make index file for my bam file that sorted by -n
written by meisam.radfar •
-
cant make index file for Wheat accepted_hits.bam file
written by meisam.radfar •
Cross-posted on bioinfo SE: https://bioinformatics.stackexchange.com/q/17685/650
Don't post your question on mutliple sites. That's inconsiderate, like asking multiple people to help you with something just to ensure you get help faster.