search

milaboratory / mixcr

MiXCR is an ultimate software platform for analysis of Next-Generation Sequencing (NGS) data for immune profiling.
https://mixcr.com
Other
330 stars 79 forks source link

Optimize affine aligner scoring #180

Open PoslavskySV opened 7 years ago

PoslavskySV commented 7 years ago

Example:

»> Clone id: 1439
»> Abundance, reads (fraction): 1 (1.0621348911311737E-4)

                                                             FR1><CDR1              CDR1><FR2      
   Quality     74878878888858778878787878877686877878788777886787886786768686677668778776887777    
   Target0   0 GAGACTTGGTCCAGCCTGGAGGGTCCCTGAGACTCTCCTGTGCAGCCTCTGGATTCACCTTTAGTAGGTTTTTTATGAAC 79   Score (hit score)
IGHV3-7*00 182 gagGcttggtccagcctggGgggtccctgagactctcctgtgcagcctctggattcacctttagtagCtAttGGatgaGc 261  1885 (2517)

                                                        FR2><CDR2               CDR2><FR3          
   Quality     678867787888864777887887876768668678777766887878866787887776877866778868 6768887    
   Target0  80 TGGGTCCGCCAGGCTCCAGGGAAAGGACTGGAGTGGGTGGCCAACATAAGGCCAGATGGCAGTGGAGAAATA-TATATGG 158  Score (hit score)
IGHV3-7*00 262 tgggtccgccaggctccagggaaGggGctggagtgggtggccaacataaAgcAagatggAagtg-agaaataCtatGtgg 340  1885 (2517)

   Quality     67888867886678888888877788888877786887888887778888888878888887888777888768787788    
   Target0 159 ACTCTGTGAAGGGCCGATTCACCATCTCCAGAGACAACACCAAGAATTCGCTATTTTTACTAATGAACAGCCTGAGAGCC 238  Score (hit score)
IGHV3-7*00 341 actctgtgaagggccgattcaccatctccagagacaacGccaagaaCtcActGtAtCtGcAaatgaacagcctgagagcc 420  1885 (2517)

                                 FR3><CDR3   V>     <D D><J                  
    Quality     888887788888888888887888888788887787888888788888888887888    
    Target0 239 GACGACACGGCCGTCTATTATTGTGCGAGAGGAGAATCCAGTGCTTTTGATATGTGG 295  Score (hit score)
 IGHV3-7*00 421 gaGgacacggcTgtGtattaCtgtgcgagag                           451  1885 (2517)
IGHD1-20*00   5                                     tccag                 9    50 (50)
 IGHD1-7*00   5                                     tccag                 9    50 (50)
IGHD5-18*00  14                                    atcca                  18   50 (50)
IGHJ3*00  23                                          tgcttttgatatCtgg 38   131 (412)
PoslavskySV commented 6 years ago

Related to #213