TopHat BAM file explanation
0
0
Entering edit mode
6.6 years ago
qudrat ▴ 100

When I see TopHat BAM header file it looks like

@SQ SN:chr10    LN:133797422
@SQ SN:chr10_GL383545v1_alt LN:179254
@SQ SN:chr10_GL383546v1_alt LN:309802
@SQ SN:chr10_KI270824v1_alt LN:181496
@SQ SN:chr10_KI270825v1_alt LN:188315

Here i have copied only chr10 header file and it has many lines. What are these multiple line indicates?

rna-seq Assembly • 2.0k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Thank you! I have one more question. I am doing Cufflinks analysis on hundreds of index BAM file. Can i do it chromosome wise to save the time as i am interested in only sex chromosome

ADD REPLY
0
Entering edit mode

no, you'll get false positives.

ADD REPLY
0
Entering edit mode

What if I split the BAM file by chromosome using following command samtools view -b input.bam chr1 chr2 chr3 chr4 > output.bam and then Cufflinks is done. How would it produce false positive? Can you help me understand or give me some reference? Actually it has been hampering my work for long and I would be extremely grateful of you.

ADD REPLY
0
Entering edit mode

sorry, i was wrong: I thought you wanted to map your reads on each distinct chromosome . I don't know what would happen for cufflink.

ADD REPLY

Login before adding your answer.

Traffic: 2503 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6