diff options
author | 2011-01-06 11:43:55 +0100 | |
---|---|---|
committer | 2011-01-06 11:43:55 +0100 | |
commit | 88432c8385bd33e16c980fdc97cc2343204a5634 (patch) | |
tree | e6fbb2dee07db181020c1da0af28c8a7f1eba8fb /sci-biology/align_to_scf/files/README | |
parent | sci-biology/roche454ace2caf: improved description and polished deps (diff) | |
download | sci-88432c8385bd33e16c980fdc97cc2343204a5634.tar.gz sci-88432c8385bd33e16c980fdc97cc2343204a5634.tar.bz2 sci-88432c8385bd33e16c980fdc97cc2343204a5634.zip |
sci-biology/align_to_scf: QA fix, changed pkgdescr and created a README file
(Portage version: 2.1.9.26/git/Linux i686, unsigned Manifest commit)
Diffstat (limited to 'sci-biology/align_to_scf/files/README')
-rw-r--r-- | sci-biology/align_to_scf/files/README | 37 |
1 files changed, 37 insertions, 0 deletions
diff --git a/sci-biology/align_to_scf/files/README b/sci-biology/align_to_scf/files/README new file mode 100644 index 000000000..be51f6fdd --- /dev/null +++ b/sci-biology/align_to_scf/files/README @@ -0,0 +1,37 @@ +* What's about "align_to_scf" ? + + This is important if you want to see the real traces within the contig editor from STADEN at a given position. + + Note: there are no information given within the ACE assembly to recalculate the trace position backward from assembly. + STADEN "trev" will access/open traces as given from "sequencer machine" and Roche assembly give you only the complete assembly + information. This is the problem ! + + Are you firm with gap4 ? + The program "trev" is automatically called within "gap4". + And you start this with: + -> open gap4 database + --> open contig editor, select contig + ---> double-click on one letter will open this trace at this position and set cursor to exact this position ! + double-click on the consensus line will open all traces at this position ! + + You can also see the "slack" by enable/disable "show cutoff"; + + Example ? + Try this very simple test data and verify this letter by letter: + + + % align_to_scf -f test.qry -o test.aln + + --- test.qry --- + >F0FH8BE01EUM6Y <- name of read from sff archive + ATGCATCATG-ATGC <- fasta sequence of this read from ace file + ATGCATGCATGCATGC <- fasta sequence of this read from sff trace file + + + And you should get this output "test.aln" ( see also CAF format description ) : + + >F0FH8BE01EUM6Y + Align_to_SCF 1 6 1 6 + Align_to_SCF 7 10 8 11 + Align_to_SCF 12 15 13 16 + |