- comme on ne sait pas si la banque contient nos deux vecteurs on peut se faire une banque qui ne contiendrait que nos vecteurs
créez une arborescence ably3
nedit mes_vecteurs&
copier/coller des sequences des deux vecteurs au format fasta
enregistrez
cross_match traces.fasta mes_vecteurs -minmatch 12 -penalty -2 -minscore 20 -screen
cp traces.fasta.qual traces.fasta.screen.qual
more traces.fasta.screen
- on va voir essayer d'assembler des contigs pour voir s'il partagent des relations de séquences
cap3 traces.fasta.screen > traces.fasta.screen.capout
more traces.fasta.screen.capout
- ou encore mieux utiliser le programme cross_match pour faire une comparaison 2 a 2 des sequences consensus
cross_match traces.fasta.screen.contigs | more