Unconfigured Ad

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

    #1

    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

      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        Yesterday, 10:25 AM
      • SEQadmin2
        How Immunogenomics Decodes Immunity’s Genetic Blueprint
        by SEQadmin2




        The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.

        This convergence of genetics, immunology, and computation...
        09-01-2026, 05:41 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 09-25-2026, 09:06 AM
      0 responses
      28 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 09-23-2026, 11:05 AM
      0 responses
      24 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 09-18-2026, 11:37 AM
      1 response
      46 views
      0 reactions
      Last Post pekgio
      by pekgio
       
      Started by SEQadmin2, 09-16-2026, 10:23 AM
      1 response
      55 views
      0 reactions
      Last Post pekgio
      by pekgio
       
      Working...