small RNA creating .bam from .sam file
1
0
Entering edit mode
5 weeks ago

Hi Biostars community,

I need help to create .bam file from .sam.

I used the code: samtools view -bS x.sam > x.bam, and error report is:

[E::sam_hdr_create] Invalid header line: must start with @HD/@SQ/@RG/@PG/@CO.

I do not see a similar error report online.

The sequence is of barley small RNA unpaired reads that showed about 70% alignment rate.

Data was generated by one of the biggest sequencing companies, and they shared us cleaned and uncleaned data. From uncleaned data, x.bam file is created using the same code but not from the cleaned one.

Help highly appreciated!

Dem

RNA-Seq • 237 views
ADD COMMENT
0
Entering edit mode

How was the .sam file created? Provide program name and command line.

From uncleaned data, x.bam file is created using the same code but not from the cleaned one.

What does this mean. Normally smallRNA will need to be "cleaned" (adapter removed) before alignments are done.

ADD REPLY
0
Entering edit mode
5 weeks ago

the message is clear: there is not SAM header in your file x.sam or it is invalid.

ADD COMMENT

Login before adding your answer.

Traffic: 1925 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6