Hi,
When I try to use "samtools view" to creat .bam file from .sam, I get these error msg. The .sam file is created by mapping raw illumina reads with tophat.
samtools view -bS -t index.fa.fai -o accepted_hits.bam accepted_hits.sam
The tag [ID] required for [PG] not present.
[sam_header_read2] 7409 sequences loaded.
Parse error at line 8263: CIGAR and sequence length are inconsistent
This #line 8263 is like this:
head -n 8263 accepted_hits.sam |tail -n 1
HWI-EAS210R:5:55:1018:1930#0 0 SL1.00sc00002 428427 0 36M92N36M473N536870911M * 0 0 CTCATCTCTCATCTCAGAGAGGTCCTCCAGGAGCAGGAATATAAGTTGGATCAAGTACGCAATCATCTTCA XXQZXR^Ta\[a\ZVZ^ZYXWWYPSUMSSQXZTO^ZZPV_\`_U\VYY[W^G\SJSTTRTUTXFTTQTYVS NM:i:0 XS:A:+ NS:i:0
Anybody could explain me why Samtools gave such error msg?
Thanks!
When I try to use "samtools view" to creat .bam file from .sam, I get these error msg. The .sam file is created by mapping raw illumina reads with tophat.
samtools view -bS -t index.fa.fai -o accepted_hits.bam accepted_hits.sam
The tag [ID] required for [PG] not present.
[sam_header_read2] 7409 sequences loaded.
Parse error at line 8263: CIGAR and sequence length are inconsistent
This #line 8263 is like this:
head -n 8263 accepted_hits.sam |tail -n 1
HWI-EAS210R:5:55:1018:1930#0 0 SL1.00sc00002 428427 0 36M92N36M473N536870911M * 0 0 CTCATCTCTCATCTCAGAGAGGTCCTCCAGGAGCAGGAATATAAGTTGGATCAAGTACGCAATCATCTTCA XXQZXR^Ta\[a\ZVZ^ZYXWWYPSUMSSQXZTO^ZZPV_\`_U\VYY[W^G\SJSTTRTUTXFTTQTYVS NM:i:0 XS:A:+ NS:i:0
Anybody could explain me why Samtools gave such error msg?
Thanks!
Comment