New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3641820	48 (0.920)	3 (0.060)	3/254	NT	6.1%	coding (1385/1431 nt)	PP_3209	isoxanthopterin deaminase
	?	NC_002947	= 3641851	47 (0.950)					coding (1354/1431 nt)	PP_3209	isoxanthopterin deaminase

GAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAG  >  NC_002947/3641679‑3641963                                                                                                                                              |                                                                                                                                                gAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTg                                                                                                                                          >  7:126380/1‑149 (MQ=255)  aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                         <  8:152349/149‑1 (MQ=255)  aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                         >  2:38856/1‑149 (MQ=255)  aaGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTgg                                                                                                                                         >  4:88762/1‑149 (MQ=255)     ggCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTg                                                                                                                                      <  7:226452/149‑1 (MQ=255)       cAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGcc                                                                                                                                    <  4:281601/149‑1 (MQ=255)        aGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGccc                                                                                                                                   <  7:165815/149‑1 (MQ=255)          tggaactggCCCCGCCGGGGCCTGCTCCACACCCAGCACGACGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAg                                                                                                                                 >  2:73989/8‑149 (MQ=255)          gggAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGcc                                                                                                                                            >  7:250819/1‑138 (MQ=255)          gggAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGcc                                                                                                                                            <  8:250819/138‑1 (MQ=255)          gggAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAg                                                                                                                                 <  5:62126/149‑1 (MQ=255)                 ggCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCg                                                                                                                          <  3:88762/149‑1 (MQ=255)                   ccccGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAg                                                                                                                        >  1:18431/1‑149 (MQ=255)                    cccGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAgg                                                                                                                       >  5:227573/1‑149 (MQ=255)                               cAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCttggacctcgccgaactgggcagccaggcacgccagg                                                                                                            <  1:306193‑M1/149‑38 (MQ=255)                                    ccACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATc                                                                                                        <  2:18431/148‑1 (MQ=255)                                          ccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGt                                                                                                 >  5:117473/1‑149 (MQ=255)                                          ccAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGt                                                                                                 >  5:262196/1‑149 (MQ=255)                                           cAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGCCGAGGTCCAAGCCCGGAAACGCATCGtt                                                                                                >  3:44897/1‑149 (MQ=255)                                              cACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATTGGCGTACCCTTGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTc                                                                                             <  1:269252/149‑1 (MQ=255)                                               aCGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTcc                                                                                            <  2:191265/149‑1 (MQ=255)                                                     tcgGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCt                                                                                       >  4:233573/1‑148 (MQ=255)                                                        gCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAgg                                                                                   <  1:146334/149‑1 (MQ=255)                                                           tgcaGGTGCGGGCCCTGCGCCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAggcgg                                                                                >  3:284774/4‑149 (MQ=255)                                                             aaGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGttt                                                                                               >  4:151842/1‑132 (MQ=255)                                                             aaGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGttt                                                                                               <  3:151842/132‑1 (MQ=255)                                                             aaGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgcc                                                                              <  6:77389/149‑1 (MQ=255)                                                              aGGTGCGGGCCCTGCGCCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGgccg                                                                             <  4:284774/149‑1 (MQ=255)                                                              aGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCACGCGGCCTGTTGAAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGAGgccg                                                                             >  7:3209/1‑149 (MQ=255)                                                                     gggCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCa                                                                      <  1:38856/149‑1 (MQ=255)                                                                           gCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCttggacctcgccgaactgggcagccaggcacgccaggctg                                                                                                         <  7:295463‑M1/108‑41 (MQ=255)                                                                           gCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCttggacctcgccgaactgggcagccaggcacgccaggctg                                                                                                         >  8:295463‑M1/1‑68 (MQ=255)                                                                                 gCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCa                                                          >  6:87479/1‑149 (MQ=255)                                                                                            cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                               >  8:13321/1‑149 (MQ=255)                                                                                            cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                               >  8:280750/1‑149 (MQ=255)                                                                                            cAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACgg                                               <  8:118758/149‑1 (MQ=255)                                                                                                  aTGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGGCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCTGACGGCCa                                            <  6:262196/146‑1 (MQ=255)                                                                                                              ccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTggg                             <  1:58033/149‑1 (MQ=255)                                                                                                              ccGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTggg                             >  8:176517/1‑149 (MQ=255)                                                                                                                   gCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGAt                        <  7:280750/149‑1 (MQ=255)                                                                                                                     cTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGc                      >  1:185710/1‑149 (MQ=255)                                                                                                                        ttGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGcc                     <  1:169240/147‑1 (MQ=255)                                                                                                                        ttGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCa                   <  8:126380/149‑1 (MQ=255)                                                                                                                            aGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGAt                >  4:27661/1‑148 (MQ=255)                                                                                                                               cGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATg            <  2:185710/149‑1 (MQ=255)                                                                                                                               aactgggcAGCCAGGCACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATg            >  5:306145/8‑149 (MQ=255)                                                                                                                                gAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACTCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCAttc           <  8:3209/149‑3 (MQ=255)                                                                                                                                    ggcAGCCAGGCACGCCAGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGAc       >  3:304057/3‑149 (MQ=255)                                                                                                                                      cAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAgccaggctgtacttcgcctgcaacaggccgcgggggtacgcgcatgagctggacgcctctgcaggtgc                                              <  7:295463‑M2/108‑69 (MQ=255)                                                                                                                                      cAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAgccaggctgtacttcgcctgcaacaggccgcgggggtacgcgcatgagctggacgcctctgcaggtgc                                              >  8:295463‑M2/1‑40 (MQ=255)                                                                                                                                       aGCCTGGCGTGCCTGTCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCg     >  1:56192/1‑148 (MQ=255)                                                                                                                                       aGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCg     >  1:258642/1‑148 (MQ=255)                                                                                                                                        gCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCaa                                     >  7:49559/1‑115 (MQ=255)                                                                                                                                        gCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCaa                                     <  8:49559/115‑1 (MQ=255)                                                                                                                                        gCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGACGATGCCCAGATCATGCAGACCGaa   <  6:117473/149‑1 (MQ=255)                                                                                                                                         ccTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAgccaggctgtacttcgcctgcaacaggccgcgggggtacgcgcatgagctggacgcctctgcaggtgcgggccctgcgccttgaagccgacgacgtgctgggtgtggaactg  <  1:306193‑M2/149‑113 (MQ=255)                                                                                                                                              |                                                                                                                                                GAAGGGCAGGGGAGCAGGCCCCGGCGGGGCCAGTTCCACACCCAGCACGTCGTCGGCTTCAAGGCGCAGGGCCCGCACCTGCAGAGGCGTCCAGCTCATGCGCGTACCCCCGCGGCCTGTTGCAGGCGAAGTACAGCCTGGCGTGCCTGGCTGCCCAGTTCGGCGAGGTCCAAGCCAGGAATCGCATCGTTTTCCACGCTCAGGCGGCCGCCTACCAGCAATGCCCGCAGCGTCGGACGGCCACCGCCAACCACTGGGGCGATGCCCAGATCATGCAGACCGAAG  >  NC_002947/3641679‑3641963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.