Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • bnfoguy
    Member
    • May 2011
    • 17

    Errors in mapping cancer exome data

    I am trying to map the reads from the ERP000703 study (http://www.ncbi.nlm.nih.gov/sra/ERX013480). I am using stampy to do the paired end mapping and when I execute the script with the following command:

    ./stampy.py -g hg18 -h hg18 -M /home/xxxxx/stampy-1.0.17/FastQ/CancerReads/ftp-trace.ncbi.nlm.nih.gov/sra/sra-instant/reads/ByStudy/sra/ERP/ERP000/ERP000703/ERR035488/ERR035488_1.fastq ERR035488_2.fastq > ERR035488_1_2.sam

    The script will run for about two hours and then abruptly crashes. When I try to conver the sam output to an unsorted bam file, I get an error saying "Parse error at line 1921134: CIGAR and sequence length are inconsistent" > aborted.

    Can someone advise me on troubleshooting this?

    Thanks,
  • dpryan
    Devon Ryan
    • Jul 2011
    • 3478

    #2
    Have you looked at the mentioned line in the SAM output to ensure that the CIGAR and sequence length agree? That would seem to be the obvious first step.

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    17 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    27 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    38 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    61 views
    0 reactions
    Last Post SEQadmin2  
    Working...