Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • glacierbird
    Member
    • Dec 2009
    • 15

    samtools view error: CIGAR and sequence length are inconsistent (tophat/bowtie)

    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!
  • glacierbird
    Member
    • Dec 2009
    • 15

    #2
    I think the problem is caused by:

    36M92N36M473N536870911M

    But, I don't know how could it occur?

    Comment

    • glacierbird
      Member
      • Dec 2009
      • 15

      #3
      I think the solution could be:
      reformat the raw .fastq file.

      Maybe, somewhere is not Tab, but space.

      then, 536870911M
      will not occur.

      Comment

      Latest Articles

      Collapse

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, Today, 11:58 AM
      0 responses
      7 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-05-2026, 10:09 AM
      0 responses
      25 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-04-2026, 08:59 AM
      0 responses
      34 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 06-02-2026, 12:03 PM
      0 responses
      56 views
      0 reactions
      Last Post SEQadmin2  
      Working...