SOAP output format - description?

[Aengus]

Hi,

I have searched the SOAP website and read the paper but I cant find a description for all the columns in the ouput.

Col1 = Read Seq_ID
Col2 = Read Sequence
Col3 = Read Sequence Quality
Col4 = ? I have the digit 1 here
Col5 = ? I have the letter 'a' here
Col6 = Read Length
Col7 = Strand
Col8 = Reference Seq_ID
Col9 = Position of start of read alignment
Col10= 0/1/2 - I think this is the number of mismatches?
Col11= base change - cant work out this out eg G->10C7


Can anyone fill in the blanks for me?

Regards
Aengus

[swbarnes2]

column 4 is uniqueness.

if you gave it paired end data, a and b signify the two halves of the same paried end read. With single end data, its always a.

Yes, column 10 is the number of mismatches. The base change the reference letter and its position, and the read's letter and its quality score. The output is a littel trickier with indels.

[AXW]

I realize this is an old thread, but I had trouble finding the specifics as well. Ended up going to the SOAPv1 page, found the output file description there.

http://soap.genomics.org.cn/soap1/#Formatofoutput