SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
novoalign native output (strand+direction) lajoieb Bioinformatics 2 10-03-2011 06:28 PM
Use Illumina mate pair and paired ends with novoalign k2bhide Bioinformatics 1 07-23-2011 09:52 AM
Novoalign speed tips for 25bp reads cc472 Bioinformatics 0 05-10-2011 07:22 AM
Can novoalign use FASTA as a reference format? (Specifically for mapping paired-end) oiiio Bioinformatics 2 04-18-2011 10:47 AM
[Help!] How can I extract high quality reads from the output file of NovoAlign qc.share Illumina/Solexa 0 09-27-2010 10:41 AM

Reply
 
Thread Tools
Old 09-15-2010, 07:33 AM   #1
sandhya
Member
 
Location: India

Join Date: Sep 2010
Posts: 11
Default Novoalign native report - Paired reads

Hi,

If anyone has worked with a report (native format) generated using novoalign, please help me with these doubts. The datasets used are Illumina paired reads.

A) Below is the snapshot of a Novoalign report (native format) for Illumina paired reads.

Code:
@0:1:1:34:429 L GAAGNAAAAATAAAAGCATTAGNAGAAATTTGTACA IIII$IIIII&IIIIIIIIIII$IIIIIIIIIIIII U 14 91 >gi|9629357:1-9117 2177 F . 2308 R
@0:1:1:34:429 R TNCTTATTAAGCNCTCTGAAATNNANNNNTTTTCTC I$IIIIIIIIII$IIIIIIIII$$'$$$$IIIIIII U 126 91 >gi|9629357:1-9117 2308 R . 2177 F 25A>G 36G>A
With the help of the Novocraft Alignment suite pdf (section Output Formats, page 24), I was able to understand certain columns in the report but please help me identify what
1. Aligned Sequence
2. Aligned Offseet
3. Pair Sequence
4. Pair Offset
5. Mismatches
are, in the report.

B ) I was also looking for the aligned reads' start and end positions. Is that information available in this report?

C) At the end of the report are 3 columns given with data

# Fragment Length Distribution
# From To Count
# 27 29 4
# 30 32 30
# 33 35 141
# 36 38 696
# 39 41 1136 ..............etc


Does this mean that from positions 27 to 29, there are 4 reads and so on.

D) Finally, here were the report statistics.

# Paired Reads: 9686877
# Pairs Aligned: 6253455
# Read Sequences: 19373754
# Aligned: 14102273
# Unique Alignment: 14102068
# Gapped Alignment: 875179
# Quality Filter: 248607
#Homopolymer Filter: 1306


I understand that 2 times Paired Reads = Read Sequences. Please help me in understanding why 2 times Pairs Aligned < Aligned. Again if I add Gapped Alignment with Unique Alignment, I do not get Aligned.

Please advice.
sandhya 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 12:10 AM.


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