Hello i'm trying to assembly a sequencing done 100% with sanger reads, using celera assembler
i have the fasta, fasta.qual, anc, clip and xml files from NCBI Trace
the xml files have the mate pair info, but when i run the assembler, skip this info
from the final .qc file
[Mates]
ReadsWithNoMate=691348(100.00%)
ReadsWithGoodMate=0(0.00%)
ReadsWithBadShortMate=0(0.00%)
ReadsWithBadLongMate=0(0.00%)
ReadsWithSameOrientMate=0(0.00%)
ReadsWithOuttieMate=0(0.00%)
ReadsWithBothChaffMate=0(0.00%)
ReadsWithChaffMate=0(0.00%)
ReadsWithBothDegenMate=0(0.00%)
ReadsWithDegenMate=0(0.00%)
ReadsWithBothSurrMate=0(0.00%)
ReadsWithSurrogateMate=0(0.00%)
ReadsWithDiffScafMate=0(0.00%)
ReadsWithUnassignedMate=0(0.00%)
TotalScaffoldLinks=0
MeanScaffoldLinkWeight=0.00
the xml file looks like
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>690</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>PACH2306</TEMPLATE_ID>
<TI>451441523</TI>
<TRACE_DIRECTION>FORWARD</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>PACH2306.x1</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>933</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>PACH2306</TEMPLATE_ID>
<TI>451443045</TI>
<TRACE_DIRECTION>REVERSE</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>PACH2306.y1</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>888</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>ACBG121468</TEMPLATE_ID>
<TI>452858124</TI>
<TRACE_DIRECTION>FORWARD</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>ACBG121468.b2</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>872</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>ACBG121468</TEMPLATE_ID>
<TI>452859852</TI>
<TRACE_DIRECTION>REVERSE</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>ACBG121468.g2</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
i need to extract this info in some way, and pass to celera assembler??
or i need a special flag so celera assembler use this??
thx in advance
Cristian
i have the fasta, fasta.qual, anc, clip and xml files from NCBI Trace
the xml files have the mate pair info, but when i run the assembler, skip this info
from the final .qc file
[Mates]
ReadsWithNoMate=691348(100.00%)
ReadsWithGoodMate=0(0.00%)
ReadsWithBadShortMate=0(0.00%)
ReadsWithBadLongMate=0(0.00%)
ReadsWithSameOrientMate=0(0.00%)
ReadsWithOuttieMate=0(0.00%)
ReadsWithBothChaffMate=0(0.00%)
ReadsWithChaffMate=0(0.00%)
ReadsWithBothDegenMate=0(0.00%)
ReadsWithDegenMate=0(0.00%)
ReadsWithBothSurrMate=0(0.00%)
ReadsWithSurrogateMate=0(0.00%)
ReadsWithDiffScafMate=0(0.00%)
ReadsWithUnassignedMate=0(0.00%)
TotalScaffoldLinks=0
MeanScaffoldLinkWeight=0.00
the xml file looks like
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>690</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>PACH2306</TEMPLATE_ID>
<TI>451441523</TI>
<TRACE_DIRECTION>FORWARD</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>PACH2306.x1</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>933</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>PACH2306</TEMPLATE_ID>
<TI>451443045</TI>
<TRACE_DIRECTION>REVERSE</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>PACH2306.y1</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>888</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>ACBG121468</TEMPLATE_ID>
<TI>452858124</TI>
<TRACE_DIRECTION>FORWARD</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>ACBG121468.b2</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
<trace>
<CENTER_NAME>JGI</CENTER_NAME>
<CLIP_LEFT>0</CLIP_LEFT>
<CLIP_RIGHT>872</CLIP_RIGHT>
<SOURCE_TYPE>GENOMIC</SOURCE_TYPE>
<SPECIES_CODE>XX</SPECIES_CODE>
<SUBMISSION_TYPE>NEW</SUBMISSION_TYPE>
<TEMPLATE_ID>ACBG121468</TEMPLATE_ID>
<TI>452859852</TI>
<TRACE_DIRECTION>REVERSE</TRACE_DIRECTION>
<TRACE_FORMAT>SCF</TRACE_FORMAT>
<TRACE_NAME>ACBG121468.g2</TRACE_NAME>
<TRACE_TYPE_CODE>WGS</TRACE_TYPE_CODE>
</trace>
i need to extract this info in some way, and pass to celera assembler??
or i need a special flag so celera assembler use this??
thx in advance
Cristian