Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • meher
    Member
    • Jun 2011
    • 54

    #16
    Originally posted by Simon Anders View Post
    Just post the SAM file, so we can see what's wrong.
    Here are the first few lines of sam file:

    @HD VN:1.0 SO:coordinate
    @SQ SN:gi|17737322|ref|chrMT| LN:16727
    @SQ SN:gi|357579592|ref|chrX| LN:123869142
    @SQ SN:gi|357579593|ref|chr38| LN:23914537
    @SQ SN:gi|357579594|ref|chr37| LN:30902991
    @SQ SN:gi|357579595|ref|chr36| LN:30810995
    @SQ SN:gi|357579596|ref|chr35| LN:26524999
    @SQ SN:gi|357579597|ref|chr34| LN:42124431
    @SQ SN:gi|357579598|ref|chr33| LN:31377067
    @SQ SN:gi|357579599|ref|chr32| LN:38810281
    @SQ SN:gi|357579600|ref|chr31| LN:39895921
    @SQ SN:gi|357579601|ref|chr30| LN:40214260
    @SQ SN:gi|357579602|ref|chr29| LN:41845238
    @SQ SN:gi|357579603|ref|chr28| LN:41182112
    @SQ SN:gi|357579604|ref|chr27| LN:45876710
    @SQ SN:gi|357579605|ref|chr26| LN:38964690
    @SQ SN:gi|357579606|ref|chr25| LN:51628933
    @SQ SN:gi|357579607|ref|chr24| LN:47698779
    @SQ SN:gi|357579608|ref|chr23| LN:52294480
    @SQ SN:gi|357579609|ref|chr22| LN:61439934
    @SQ SN:gi|357579610|ref|chr21| LN:50858623
    @SQ SN:gi|357579611|ref|chr20| LN:58134056
    @SQ SN:gi|357579612|ref|chr19| LN:53741614
    @SQ SN:gi|357579613|ref|chr18| LN:55844845
    @SQ SN:gi|357579614|ref|chr17| LN:64289059
    @SQ SN:gi|357579615|ref|chr16| LN:59632846
    @SQ SN:gi|357579616|ref|chr15| LN:64190966
    @SQ SN:gi|357579617|ref|chr14| LN:60966679
    @SQ SN:gi|357579618|ref|chr13| LN:63241923
    @SQ SN:gi|357579619|ref|chr12| LN:72498081
    @SQ SN:gi|357579620|ref|chr11| LN:74389097
    @SQ SN:gi|357579621|ref|chr10| LN:69331447
    @SQ SN:gi|357579622|ref|chr9| LN:61074082
    @SQ SN:gi|357579623|ref|chr8| LN:74330416
    @SQ SN:gi|357579624|ref|chr7| LN:80974532
    @SQ SN:gi|357579625|ref|chr6| LN:77573801
    @SQ SN:gi|357579626|ref|chr5| LN:88915250
    @SQ SN:gi|357579627|ref|chr4| LN:88276631
    @SQ SN:gi|357579628|ref|chr3| LN:91889043
    @SQ SN:gi|357579629|ref|chr2| LN:85426708
    @SQ SN:gi|357579630|ref|chr1| LN:122678785
    @PG ID:TopHat VN:2.0.9 CL:/csc/lohi/dog_tools/tophat-2.0.9.Linux_x86_64/tophat -p 3 -o ./ /csc/lohi/canFam3_ref_dogData/ncbi_fasta_canFam3/canFam3 B00EJHW.t30l32.pair1.fastq.gz B00EJHW.t30l32.pair2.fastq.gz
    HISEQ4_0112:3:1101:1044:57612#GCCAAT 89 gi|357579618|ref|chr13| 621671950 101M * 0 0 TGCACTTGTTGAACTGCTGAAACACAAGCCCAAGGCAACAGATGAACAACTGAAAACTGTTATGGGAGATTTTGGAGCCTTTGTAGAGAAGTGCTGCGCNG @>3CDDC@3@DBC>6;6CC@EEBA?EEGBIGHCGIGFF@BC:JJJIIHCGD4BDGIFIGIGHEDGCJJIJJJIIHJJJJJIJJJJJHGHHHFFDDA4#C AS:i:-1 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:99A1 YT:Z:UUNH:i:1
    HISEQ4_0112:3:1101:1044:61660#GCCAAT 99 gi|357579610|ref|chr21| 328438950 6M1452N95M = 32845462 1652 CNCAAGAGTTTACTCTCCAAAGCTCGAGGAATCGATTCCAGCTCTGTTAAACTCCGAGGTGGTTCTTTATTCATGGATACAGAAAAATCAGGAAAAAGGGA C#1ADFFFHHHHHJJJJJJJJJJJJJJJJJIJJJJJJJJIIJJJJJJJJJJJJJJJHIJEHHHHHFFFFFFFFEEEEEEDDDDDDDDDDDDDDDDDDDDDD AS:i:-1 XN:i:0 XM:i:1 XO:i:0 XG:i:0 NM:i:1 MD:Z:1C99 YT:Z:UU XS:A:+ NH:i:1


    Should you require the whole sam file to be attached?

    Comment

    • Simon Anders
      Senior Member
      • Feb 2010
      • 995

      #17
      Your chomosome names don't match. HTSeq-count cannot guess that "gi|357579630|ref|chr1|" in the SAM file is the same as "chr1" in the GFF file.

      Comment

      • meher
        Member
        • Jun 2011
        • 54

        #18
        Originally posted by Simon Anders View Post
        Your chomosome names don't match. HTSeq-count cannot guess that "gi|357579630|ref|chr1|" in the SAM file is the same as "chr1" in the GFF file.

        Thank you. As i thought both have the same notation i.e. chr1,chr2.. they are in sync with each other but that seems to be costly.

        Comment

        Latest Articles

        Collapse

        • kinpict.com
          Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
          by kinpict.com
          The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
          Today, 06:37 PM
        • SEQadmin2
          16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
          by SEQadmin2





          NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
          10-30-2026, 06:46 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, 11-12-2026, 07:59 AM
        0 responses
        27 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-05-2026, 09:25 AM
        0 responses
        35 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-03-2026, 11:54 AM
        2 responses
        56 views
        0 reactions
        Last Post kinpict.com  
        Started by SEQadmin2, 10-28-2026, 09:57 AM
        0 responses
        175 views
        0 reactions
        Last Post SEQadmin2  
        Working...