SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools flagstat output nguyendofx Bioinformatics 23 01-22-2014 04:45 AM
samtools flagstat output manore Bioinformatics 11 07-23-2013 11:19 PM
SAMtools flagstat output interpretation a2z@blr Bioinformatics 2 10-20-2011 01:23 PM
Samtools flagstat Anelda Bioinformatics 0 09-26-2011 03:55 AM
samtools flagstat bair Bioinformatics 3 05-28-2010 06:15 AM

Reply
 
Thread Tools
Old 09-11-2014, 06:41 PM   #1
i19870503
Junior Member
 
Location: china

Join Date: Aug 2014
Posts: 1
Default 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=
i19870503 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:27 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO