Miyakogusa Predicted Gene

Lj1g3v0129860.1
Show Alignment: 

BLASTN 2.2.25 [Feb-01-2011]


Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, 
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), 
"Gapped BLAST and PSI-BLAST: a new generation of protein database search
programs",  Nucleic Acids Res. 25:3389-3402.

Query= Lj1g3v0129860.1 tr|Q8S5C2|Q8S5C2_SOYBN 4-coumarate:CoA ligase
isoenzyme 3 OS=Glycine max GN=Gm4CL-3 PE=2 SV=1,84.06,0,Acetyl-CoA
synthetase-like,NULL; seg,NULL; no description,NULL;
AMP-binding,AMP-dependent
synthetase,NODE_26210_length_2390_cov_154.529709.path1.1
         (1746 letters)

Database: LJGI 
           47,486 sequences; 32,788,469 total letters

Searching..................................................done



                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gnl|LJGI|TC76024 similar to UniRef100_P31687 Cluster: 4-coumarat...  1283   0.0  
gnl|LJGI|TC70657 homologue to UniRef100_P31687 Cluster: 4-coumar...  1072   0.0  
gnl|LJGI|TC76208 homologue to UniRef100_Q8S5C2 Cluster: 4-coumar...  1045   0.0  
gnl|LJGI|TC76664 similar to UniRef100_Q8S5C2 Cluster: 4-coumarat...   854   0.0  
gnl|LJGI|FS350112 similar to UniRef100_P31687 Cluster: 4-coumara...   307   1e-82
gnl|LJGI|TC72583 homologue to UniRef100_P31687 Cluster: 4-coumar...   250   3e-65
gnl|LJGI|TC62381 homologue to UniRef100_Q8W558 Cluster: 4-coumar...    74   3e-12
gnl|LJGI|TC59755 homologue to UniRef100_Q8W558 Cluster: 4-coumar...    68   2e-10
gnl|LJGI|TC71364 similar to UniRef100_Q8W558 Cluster: 4-coumarat...    66   8e-10

>gnl|LJGI|TC76024 similar to UniRef100_P31687 Cluster: 4-coumarate--CoA ligase 2; n=1;
            Glycine max|Rep: 4-coumarate--CoA ligase 2 - Glycine max
            (Soybean), partial (36%)
          Length = 729

 Score = 1283 bits (647), Expect = 0.0
 Identities = 647/647 (100%)
 Strand = Plus / Plus

                                                                        
Query: 462  gaaactcgtcatcactcaggcaatgtacgtcgacaagctccggccactggaagacgttaa 521
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 1    gaaactcgtcatcactcaggcaatgtacgtcgacaagctccggccactggaagacgttaa 60

                                                                        
Query: 522  cggtaacgactacccgaagctcggggaggatttcaaggtggtcaccgtggacgaaccacc 581
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 61   cggtaacgactacccgaagctcggggaggatttcaaggtggtcaccgtggacgaaccacc 120

                                                                        
Query: 582  caaggattgcttggacttctccgtgatttcggaaggaaaagaggatgatctgccggaggt 641
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 121  caaggattgcttggacttctccgtgatttcggaaggaaaagaggatgatctgccggaggt 180

                                                                        
Query: 642  ggagatcaacccagaggacgcggtggcgctgccattctcctccggcacgacgggcctccc 701
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 181  ggagatcaacccagaggacgcggtggcgctgccattctcctccggcacgacgggcctccc 240

                                                                        
Query: 702  gaagggggtggttctgactcacaagagccttaccacgagcgtgggccagcaagtggatgg 761
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 241  gaagggggtggttctgactcacaagagccttaccacgagcgtgggccagcaagtggatgg 300

                                                                        
Query: 762  agagaatccgaacctttctctgggaacagaggatgttctgctttgcgtgcttcccctgtt 821
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 301  agagaatccgaacctttctctgggaacagaggatgttctgctttgcgtgcttcccctgtt 360

                                                                        
Query: 822  tcacattttctcgctgaacagcgttcttctctgcgcgctgagagcggggagcggtgtgct 881
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 361  tcacattttctcgctgaacagcgttcttctctgcgcgctgagagcggggagcggtgtgct 420

                                                                        
Query: 882  tctgatgccgaagtttgagatagggtcgttgctggagctgattcagaagcacagggtgtc 941
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 421  tctgatgccgaagtttgagatagggtcgttgctggagctgattcagaagcacagggtgtc 480

                                                                        
Query: 942  ggtggccatggtggtgccgccgctagtgctggctctggcgaagaatccgaaggtggcgga 1001
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 481  ggtggccatggtggtgccgccgctagtgctggctctggcgaagaatccgaaggtggcgga 540

                                                                        
Query: 1002 gtttgacctgagctccataaggttggtgctgtcaggggctgcgcctctcgggaaggagct 1061
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 541  gtttgacctgagctccataaggttggtgctgtcaggggctgcgcctctcgggaaggagct 600

                                                           
Query: 1062 cgaggaagctctccggagcagggtgcctcaagctgttttgggacagg 1108
            |||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 601  cgaggaagctctccggagcagggtgcctcaagctgttttgggacagg 647


>gnl|LJGI|TC70657 homologue to UniRef100_P31687 Cluster: 4-coumarate--CoA ligase 2;
            n=1; Glycine max|Rep: 4-coumarate--CoA ligase 2 - Glycine
            max (Soybean), partial (31%)
          Length = 964

 Score = 1072 bits (541), Expect = 0.0
 Identities = 541/541 (100%)
 Strand = Plus / Plus

                                                                        
Query: 1206 gagaaatgctgagctcaaggttcttgaccctgaaactggttgttctcttggctataatca 1265
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 89   gagaaatgctgagctcaaggttcttgaccctgaaactggttgttctcttggctataatca 148

                                                                        
Query: 1266 acccggtgaaatttgcatccgaggccaacagatcatgaaaggatatttgaacgatgagaa 1325
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 149  acccggtgaaatttgcatccgaggccaacagatcatgaaaggatatttgaacgatgagaa 208

                                                                        
Query: 1326 ggccacagccacaacaattgatgcggaaggttggcttcacaccggtgatgttggttacat 1385
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 209  ggccacagccacaacaattgatgcggaaggttggcttcacaccggtgatgttggttacat 268

                                                                        
Query: 1386 tgatgatgatgatgagattttcatagttgacagggtgaaggaactcatcaaattcaaagg 1445
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 269  tgatgatgatgatgagattttcatagttgacagggtgaaggaactcatcaaattcaaagg 328

                                                                        
Query: 1446 cttccaagtgcccccagctgaactcgaaggccttctcgtaagccatccctcaattgcaga 1505
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 329  cttccaagtgcccccagctgaactcgaaggccttctcgtaagccatccctcaattgcaga 388

                                                                        
Query: 1506 tgcagctgtcgtcccacaaaaggatgctgctgctggtgaagttcccgttgcctttgttgt 1565
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 389  tgcagctgtcgtcccacaaaaggatgctgctgctggtgaagttcccgttgcctttgttgt 448

                                                                        
Query: 1566 gagatcaaatgggtttgatctaactgaagaggctgtaaaggagttcatatccaaacaggt 1625
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 449  gagatcaaatgggtttgatctaactgaagaggctgtaaaggagttcatatccaaacaggt 508

                                                                        
Query: 1626 ggtgttttataaaagactgcacaaggtttactttgttcatgcaattcccaagtctccatc 1685
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 509  ggtgttttataaaagactgcacaaggtttactttgttcatgcaattcccaagtctccatc 568

                                                                        
Query: 1686 aggaaagatattgaggaaagacctcagagcaaagctagaaaatggaaatcagaagccata 1745
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 569  aggaaagatattgaggaaagacctcagagcaaagctagaaaatggaaatcagaagccata 628

             
Query: 1746 g 1746
            |
Sbjct: 629  g 629


>gnl|LJGI|TC76208 homologue to UniRef100_Q8S5C2 Cluster: 4-coumarate:CoA ligase
            isoenzyme 3; n=1; Glycine max|Rep: 4-coumarate:CoA ligase
            isoenzyme 3 - Glycine max (Soybean), partial (31%)
          Length = 527

 Score = 1045 bits (527), Expect = 0.0
 Identities = 527/527 (100%)
 Strand = Plus / Plus

                                                                        
Query: 863  gagcggggagcggtgtgcttctgatgccgaagtttgagatagggtcgttgctggagctga 922
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 1    gagcggggagcggtgtgcttctgatgccgaagtttgagatagggtcgttgctggagctga 60

                                                                        
Query: 923  ttcagaagcacagggtgtcggtggccatggtggtgccgccgctagtgctggctctggcga 982
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 61   ttcagaagcacagggtgtcggtggccatggtggtgccgccgctagtgctggctctggcga 120

                                                                        
Query: 983  agaatccgaaggtggcggagtttgacctgagctccataaggttggtgctgtcaggggctg 1042
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 121  agaatccgaaggtggcggagtttgacctgagctccataaggttggtgctgtcaggggctg 180

                                                                        
Query: 1043 cgcctctcgggaaggagctcgaggaagctctccggagcagggtgcctcaagctgttttgg 1102
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 181  cgcctctcgggaaggagctcgaggaagctctccggagcagggtgcctcaagctgttttgg 240

                                                                        
Query: 1103 gacagggatacggaatgacggaggcagggccggttctgtccatgtccttaggatttgcaa 1162
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 241  gacagggatacggaatgacggaggcagggccggttctgtccatgtccttaggatttgcaa 300

                                                                        
Query: 1163 agcaaccgttccccacaaagtcaggttcatgtggtactgttgtgagaaatgctgagctca 1222
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 301  agcaaccgttccccacaaagtcaggttcatgtggtactgttgtgagaaatgctgagctca 360

                                                                        
Query: 1223 aggttcttgaccctgaaactggttgttctcttggctataatcaacccggtgaaatttgca 1282
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 361  aggttcttgaccctgaaactggttgttctcttggctataatcaacccggtgaaatttgca 420

                                                                        
Query: 1283 tccgaggccaacagatcatgaaaggatatttgaacgatgagaaggccacagccacaacaa 1342
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 421  tccgaggccaacagatcatgaaaggatatttgaacgatgagaaggccacagccacaacaa 480

                                                           
Query: 1343 ttgatgcggaaggttggcttcacaccggtgatgttggttacattgat 1389
            |||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 481  ttgatgcggaaggttggcttcacaccggtgatgttggttacattgat 527


>gnl|LJGI|TC76664 similar to UniRef100_Q8S5C2 Cluster: 4-coumarate:CoA ligase
           isoenzyme 3; n=1; Glycine max|Rep: 4-coumarate:CoA
           ligase isoenzyme 3 - Glycine max (Soybean), partial
           (19%)
          Length = 570

 Score =  854 bits (431), Expect = 0.0
 Identities = 438/439 (99%), Gaps = 1/439 (0%)
 Strand = Plus / Plus

                                                                       
Query: 1   atgttgtcagttgctccttctgttgatgccaaacagacctcaaccgcaccagaatcagtg 60
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 132 atgttgtcagttgctccttctgttgatgccaaacagacctcaaccgcaccagaatcagtg 191

                                                                       
Query: 61  tcaaatatttcagacaccaccccaacacaaaacgcatctaattcccaaactacccacata 120
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 192 tcaaatatttcagacaccaccccaacacaaaacgcatctaattcccaaactacccacata 251

                                                                       
Query: 121 ttcaaatctaaattaccagacataccaatctccaaccaccttcccctccacgcttactgc 180
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 252 ttcaaatctaaattaccagacataccaatctccaaccaccttcccctccacgcttactgc 311

                                                                       
Query: 181 ttcgagaatcttccacaattctccgaccgtccatgcctcatcctcggctccaccggaaaa 240
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 312 ttcgagaatcttccacaattctccgaccgtccatgcctcatcctcggctccaccggaaaa 371

                                                                       
Query: 241 atctacacctacgcc-gagactcacctcctctcccgccgcgtcgccgccggaatgtccaa 299
           ||||||||||||||| ||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 372 atctacacctacgccggagactcacctcctctcccgccgcgtcgccgccggaatgtccaa 431

                                                                       
Query: 300 gctcggaatcaccaaaggcgatgtcgtcatgatccttctccaaaactccgctgaattcgt 359
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 432 gctcggaatcaccaaaggcgatgtcgtcatgatccttctccaaaactccgctgaattcgt 491

                                                                       
Query: 360 cttctccttcttcgccgcctccatgatcggcgccgtcgccaccaccgccaacccgttcta 419
           ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 492 cttctccttcttcgccgcctccatgatcggcgccgtcgccaccaccgccaacccgttcta 551

                              
Query: 420 caccaccgctgagcttttc 438
           |||||||||||||||||||
Sbjct: 552 caccaccgctgagcttttc 570


>gnl|LJGI|FS350112 similar to UniRef100_P31687 Cluster: 4-coumarate--CoA ligase 2; n=1;
            Glycine max|Rep: 4-coumarate--CoA ligase 2 - Glycine max
            (Soybean), partial (36%)
          Length = 630

 Score =  307 bits (155), Expect = 1e-82
 Identities = 422/511 (82%)
 Strand = Plus / Minus

                                                                        
Query: 941  cggtggccatggtggtgccgccgctagtgctggctctggcgaagaatccgaaggtggcgg 1000
            ||||||| ||||| ||||||||| | ||| |||| ||||||||| | |||| ||||| ||
Sbjct: 542  cggtggcgatggtagtgccgccgatggtggtggcgctggcgaaggacccgagggtggagg 483

                                                                        
Query: 1001 agtttgacctgagctccataaggttggtgctgtcaggggctgcgcctctcgggaaggagc 1060
            ||| |||||| ||||||||| |||||||| |||||||||||||||| || |||||||| |
Sbjct: 482  agtatgacctcagctccatacggttggtgatgtcaggggctgcgccgctggggaaggacc 423

                                                                        
Query: 1061 tcgaggaagctctccggagcagggtgcctcaagctgttttgggacagggatacggaatga 1120
            | ||||| |||||||||||| || ||||||  || |||||||||||||| || || ||||
Sbjct: 422  tggaggatgctctccggagccggttgcctcgcgcagttttgggacagggttatgggatga 363

                                                                        
Query: 1121 cggaggcagggccggttctgtccatgtccttaggatttgcaaagcaaccgttccccacaa 1180
            | ||||| ||||||||| |||| ||||  || || ||||| ||||| ||||| || || |
Sbjct: 362  cagaggcggggccggttttgtcaatgtgtttggggtttgcgaagcagccgtttccaacga 303

                                                                        
Query: 1181 agtcaggttcatgtggtactgttgtgagaaatgctgagctcaaggttcttgaccctgaaa 1240
            || |||| || || || || || ||||| ||||||||| | |||||  |||||||||| |
Sbjct: 302  agccagggtcgtgcggcacggtggtgaggaatgctgagttgaaggtgattgaccctgata 243

                                                                        
Query: 1241 ctggttgttctcttggctataatcaacccggtgaaatttgcatccgaggccaacagatca 1300
             |||||| ||  ||     ||| |   |||| || ||||||||||| |||||||| ||||
Sbjct: 242  gtggttgctccgttcctcgtaaccgctccggcgagatttgcatccgtggccaacaaatca 183

                                                                        
Query: 1301 tgaaaggatatttgaacgatgagaaggccacagccacaacaattgatgcggaaggttggc 1360
            ||||||||||  |||| |||||||||||||| ||  |||| ||||||||  |||||||||
Sbjct: 182  tgaaaggatacctgaatgatgagaaggccactgcagcaactattgatgcacaaggttggc 123

                                                                        
Query: 1361 ttcacaccggtgatgttggttacattgatgatgatgatgagattttcatagttgacaggg 1420
            ||||||| |||||  |||||||  | ||||||||||||||||||||||| ||||| ||||
Sbjct: 122  ttcacactggtgacattggttatgtagatgatgatgatgagattttcattgttgataggg 63

                                           
Query: 1421 tgaaggaactcatcaaattcaaaggcttcca 1451
            | ||||||||||||||||||||||| |||||
Sbjct: 62   ttaaggaactcatcaaattcaaaggtttcca 32


>gnl|LJGI|TC72583 homologue to UniRef100_P31687 Cluster: 4-coumarate--CoA ligase 2;
            n=1; Glycine max|Rep: 4-coumarate--CoA ligase 2 - Glycine
            max (Soybean), partial (16%)
          Length = 1000

 Score =  250 bits (126), Expect = 3e-65
 Identities = 126/126 (100%)
 Strand = Plus / Plus

                                                                        
Query: 1621 caggtggtgttttataaaagactgcacaaggtttactttgttcatgcaattcccaagtct 1680
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 641  caggtggtgttttataaaagactgcacaaggtttactttgttcatgcaattcccaagtct 700

                                                                        
Query: 1681 ccatcaggaaagatattgaggaaagacctcagagcaaagctagaaaatggaaatcagaag 1740
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 701  ccatcaggaaagatattgaggaaagacctcagagcaaagctagaaaatggaaatcagaag 760

                  
Query: 1741 ccatag 1746
            ||||||
Sbjct: 761  ccatag 766



 Score =  208 bits (105), Expect = 9e-53
 Identities = 105/105 (100%)
 Strand = Plus / Plus

                                                                        
Query: 1521 acaaaaggatgctgctgctggtgaagttcccgttgcctttgttgtgagatcaaatgggtt 1580
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 400  acaaaaggatgctgctgctggtgaagttcccgttgcctttgttgtgagatcaaatgggtt 459

                                                         
Query: 1581 tgatctaactgaagaggctgtaaaggagttcatatccaaacaggt 1625
            |||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 460  tgatctaactgaagaggctgtaaaggagttcatatccaaacaggt 504



 Score =  135 bits (68), Expect = 1e-30
 Identities = 68/68 (100%)
 Strand = Plus / Plus

                                                                        
Query: 1453 gtgcccccagctgaactcgaaggccttctcgtaagccatccctcaattgcagatgcagct 1512
            ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Sbjct: 107  gtgcccccagctgaactcgaaggccttctcgtaagccatccctcaattgcagatgcagct 166

                    
Query: 1513 gtcgtccc 1520
            ||||||||
Sbjct: 167  gtcgtccc 174


>gnl|LJGI|TC62381 homologue to UniRef100_Q8W558 Cluster: 4-coumarate:CoA ligase; n=1;
            Amorpha fruticosa|Rep: 4-coumarate:CoA ligase - Amorpha
            fruticosa, partial (59%)
          Length = 1215

 Score = 73.8 bits (37), Expect = 3e-12
 Identities = 106/129 (82%)
 Strand = Plus / Plus

                                                                        
Query: 1351 gaaggttggcttcacaccggtgatgttggttacattgatgatgatgatgagattttcata 1410
            ||||||||| | || || |||||| |||| ||||| || || ||||||||| | ||||| 
Sbjct: 567  gaaggttggttgcatacaggtgatattggatacatcgacgaagatgatgagttgttcatc 626

                                                                        
Query: 1411 gttgacagggtgaaggaactcatcaaattcaaaggcttccaagtgcccccagctgaactc 1470
            ||||| ||| |||||||| | |||||||  ||||| || |||||| | ||||||||||||
Sbjct: 627  gttgataggttgaaggaattgatcaaatataaaggttttcaagtggctccagctgaactc 686

                     
Query: 1471 gaaggcctt 1479
            |||| ||||
Sbjct: 687  gaagccctt 695


>gnl|LJGI|TC59755 homologue to UniRef100_Q8W558 Cluster: 4-coumarate:CoA ligase; n=1;
            Amorpha fruticosa|Rep: 4-coumarate:CoA ligase - Amorpha
            fruticosa, partial (34%)
          Length = 833

 Score = 67.9 bits (34), Expect = 2e-10
 Identities = 88/106 (83%)
 Strand = Plus / Plus

                                                                        
Query: 1376 ttggttacattgatgatgatgatgagattttcatagttgacagggtgaaggaactcatca 1435
            ||||||| ||||| |||||||||||| | ||||| ||||| ||| |||||||| | ||||
Sbjct: 195  ttggttatattgacgatgatgatgagttgttcatcgttgataggttgaaggaattgatca 254

                                                          
Query: 1436 aattcaaaggcttccaagtgcccccagctgaactcgaaggccttct 1481
            ||| |||||| || || ||  | ||||||||||| |||| ||||||
Sbjct: 255  aatacaaaggatttcaggttgctccagctgaacttgaagcccttct 300


>gnl|LJGI|TC71364 similar to UniRef100_Q8W558 Cluster: 4-coumarate:CoA ligase; n=1;
           Amorpha fruticosa|Rep: 4-coumarate:CoA ligase - Amorpha
           fruticosa, partial (25%)
          Length = 528

 Score = 65.9 bits (33), Expect = 8e-10
 Identities = 57/65 (87%)
 Strand = Plus / Plus

                                                                       
Query: 313 aaaggcgatgtcgtcatgatccttctccaaaactccgctgaattcgtcttctccttcttc 372
           ||||| || ||| |||||||||| ||||  |||||| |||||||||||||||||||| ||
Sbjct: 317 aaaggagacgtcatcatgatcctcctccccaactcccctgaattcgtcttctccttcctc 376

                
Query: 373 gccgc 377
           |||||
Sbjct: 377 gccgc 381