Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • i19870503
    Junior Member
    • Aug 2014
    • 1

    #1

    some questions on cuffmerge and samtools flagstat output

    1.RNAseq is solid data mapped with bioscope, and I add XA tag manually to run cufflinks.
    when I use cuffmerge , it shows:Warning: couldn't find fasta record for 'chr1'! and so on.
    but I got merge.gtf, is there any different with normal output.

    2.I check the bam file with samtools flagstat, but I think there is something wrong
    54196276 + 0 in total (QC-passed reads + QC-failed reads)
    0 + 0 duplicates
    54196276 + 0 mapped (100.00%:-nan%)
    0 + 0 paired in sequencing
    0 + 0 read1
    0 + 0 read2
    0 + 0 properly paired (-nan%:-nan%)
    0 + 0 with itself and mate mapped
    0 + 0 singletons (-nan%:-nan%)
    0 + 0 with mate mapped to a different chr
    0 + 0 with mate mapped to a different chr (mapQ>=5)



    and also the samtools view output is here, did that mapped and sort correctly?
    @HD VN:1.0 SO:coordinate
    @SQ SN:chr1 LN:295529705 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:2cf1f3b70031dc473022a5a5893e8edb
    @SQ SN:chr2 LN:140133492 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:b709bf0683d1abdca4f6236a45524577
    @SQ SN:chr3 LN:123599780 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:2e896d4cc99e908405eb690a71467e37
    @SQ SN:chr4 LN:136254946 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:32ceede4917f4ea64a7ab93040211d21
    @SQ SN:chr5 LN:100516970 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:50bd953e6de61fb729390a122f830653
    @SQ SN:chr6 LN:123305171 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:2ed43939ee3490a30cbe40e96cf33289
    @SQ SN:chr7 LN:136409062 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:57ed74ddd703b537d8fece561887bcf4
    @SQ SN:chr8 LN:119985671 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:deadaedc161e58bf8c8d6b3ea1c25eba
    @SQ SN:chr9 LN:132468591 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:89eae5d26c774408170957bb006efbbd
    @SQ SN:chr10 LN:66736929 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:171e7a3f72ba195f457720406be3604b
    @SQ SN:chr11 LN:79814395 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:feeab3090c18a3cab9547f99d35052a8
    @SQ SN:chr12 LN:57431344 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:f3309f09ea225f9a191b5d0d9eee961d
    @SQ SN:chr13 LN:145235301 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:3a4b8ac5120541c225eb870fa4a2b501
    @SQ SN:chr14 LN:148510138 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:c607e774abadba7ed0b0185f18e3f9d5
    @SQ SN:chr15 LN:134541103 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:08f26a29f88b39bff614cd6fbfa7d2ef
    @SQ SN:chr16 LN:77435658 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:0f5873a922aa3327b4523def56cad5e6
    @SQ SN:chr17 LN:64395339 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:ac533e777ca34fe996fdf58598266b98
    @SQ SN:chr18 LN:54309914 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:1cc33aec4862e0d56e691e5dc5960501
    @SQ SN:chrX LN:125871292 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:d09f454e7e0e413bdd5e65af386c581d
    @SQ SN:chrM LN:16770 UR:file:/mnt/SolidPool/reference/genomes/susScr2/swine_susScr2_validated.fa M5:c32d2133c7327cc9e47776c5eed17977
    @RG ID:20121105185046758 LB:50F SM:SL002_251012_FC1_PG2_F3
    @PG ID:bioscope-genome-mapping VN:bioscope-1.2 CL:mapping.qual.pvalue=1; mapping.mismatch.penalty=-2.0; mapping.run.classic=false; mapping.qual.anchor.length=25; matching.max.hits=10; mapping.qual.filter.cutoff=0; mapping.classic.anchor.length=50; mapping.valid.adjacent=0; mapping.classic.mismatch=0; mapping.scheme.unmapped=25.2.0\:20; matching.use.iub.reference=false; mapping.qual.error.rate=0.2; mapping.qual.mismatch.level=3; mapping.qual.bvalue=1.0; read.length=50; mapping.scheme.repetitive=

Latest Articles

Collapse

  • SEQadmin2
    Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
    by SEQadmin2



    Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

    Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
    ...
    10-14-2026, 09:55 AM
  • 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
    ...
    09-28-2026, 10:25 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
31 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-15-2026, 10:52 AM
0 responses
39 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-14-2026, 05:39 AM
0 responses
45 views
0 reactions
Last Post SEQadmin2  
Working...