Cufflinks invalid BAM binary header
I get this warning when running cufflinks on my RNA-seq. Does not prevent me from doing my analysis, I' am just curious what it is. In short, I have paired end RNA-seq mapped with STAR (v2.4.1d) against GRCm38.
STAR \
--outSAMtype BAM SortedByCoordinate \
--runThreadN 8 \
--outSAMstrandField intronMotif \
--genomeDir ./GRCm38/star \
--readFilesIn [fastqfiles] \
--readFilesCommand zcat \
--outFileNamePrefix Sample_het1.
Then quantify with cufflinks (v2.2.1).
cufflinks \
-q \
-p 8 \
-o ./ \
-b Mus_musculus.GRCm38.dna.primary_assembly.fa \
-G Mus_musculus.GRCm38.79.gtf Sample_het1.Aligned.sortedByCoord.out.bam
I get the following message
[bam_header_read] invalid BAM binary header (this is not a BAM file)
• 3,932 views
•
link
updated
by
Ram
•
written
by
murphy.charlesj
•
written
by
murphy.charlesj
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How I should set `STAR` parameters to keep all the intronic reads?
written by Dan •
-
unmapped output of the STAR
written by Wakala •
-
Aligning multiple fastq files with genome in one script/one line with STAR
written by j_eag •
-
STAR alignment EXITING because of FATAL ERROR: could not create output file
written by gnomee •
-
ERROR In aligmnet of RNA_seq data using STAR aligner_updated
written by akhan •
-
STAR+Cufflinks for RNA-Seq analysis
written by wangdp123
-
Alignment loop using STAR
written by lamia_203
-
Strange Error while running STAR.
written by jaidev53ster •
-
cufflinks large bundles warning
written by Sam
-
Run STAR with 4 fastq files, ERROR
written by shangqianwang •
Try
samtools viewto see if the header is present.You can try to rehead the file using
reheaderCheck the syntax here or on the samtools page Bam Header Edit
Indeed there is, here are the first few lines:
Run
If it doesn't say it's a "gzip" file, then it's NOT bam.
Ran the command and its says it is a gzip file