New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 3145831	25 (0.620)	4 (0.110)	4/294	NT	13.6%	coding (237/1188 nt)	metC	cystathionine beta‑lyase, PLP‑dependent
	?	CP009273	= 3146357	27 (0.710)					coding (763/1188 nt)	metC	cystathionine beta‑lyase, PLP‑dependent

AGACATCCAGACGTATAAAAACAGGAATCCCGACATGGCGGACAAAAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGT  >  CP009273/3145561‑3146091                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                     agaaaaccacanatataAAAAAAGGAATCCCGACATGGCAGACAAAAAGCTTGATACTCAATTGGTGAATGCAGGACGCAGCAAAAAATCCACGCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTGCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCCCGCCAATGGAGAGTTGT‑TCGATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCg                                                                                                                                                                                                                                                                  <  2:67389/262‑1 (MQ=255)   aCATCCAAACGTAGAACAAAAGGAATCACTGCATGCCGGACAAACAGCTTGTGAGTCAACTGTTGAATGCAGGACGCAGCAAAAAATACACGCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCGGGTCTTTGACATTGTAGAAGACAAAAAACCCGCCACACGTAATCGCGCTAATTGACAGTTGT‑TCTACGTACGGCGTGGAACGTTAACCCATCCCTCCTTACACAAACCTATGTGTGAACTGGAAGCTTGCGCAGGCTGCgcg                                                                                                                                                                                                                                                                <  2:308033/275‑3 (MQ=255)                           aTCCAGACATGGAGGACAAAAAGCTTGATACGCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGGAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCg                                                                                                                                                                                                                                        <  2:250818/275‑1 (MQ=255)                            tCCCGACATGGCGGACAAAAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGggcggc                                                                                                                                                                                                                                            >  1:250821/1‑270 (MQ=255)                                                      aTACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGATGGTCTTTGACAGTGTAGAAGCCAAAAAACAAGAGACACGTCATCGCGCCACTGGCGAGTTGT‑TCTATGGAAGACGCGGCACGTTAACCCATGTCACCTGACAACCAGCGAGGTGTGAACTGGCACGTGGCGCATGATGAGTTCTAGTTCCCGGCGGGGCGGCTGAGGTGGCTAATGCCATTCTTGCTTTTa                                                                                                                                                                                                               >  2:324602/1‑273 (MQ=255)                                                                 gTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCTAACAGGGCg                                                                                                                                                                                                   >  1:297866/1‑274 (MQ=255)                                                                           gACGCACAAAAAAATCCACCCGCGGGGAGTTAAATAGCGTGATTCAGCGCGCTTCTGCGGTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TATATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtg                                                                                                                                                                                        <  2:297861/275‑1 (MQ=255)                                                                                agcaAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGATGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTAGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGCGAACTGGAAGGTNGCGCAGGCNGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCg                                                                                                                                                                                                            >  1:493946/1‑250 (MQ=255)                                                                                                                                    cGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAAACCATTTCTCCTTACAACAAACGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTCTGAACCGAGTCAGGATTTTTGTAACAAAATCCCCAGCaa                                                                                                                                >  1:198405/1‑274 (MQ=255)                                                                                                                                                   gtgtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTcc                                                                                                                                                                                                                            >  1:81568/1‑167 (MQ=255)                                                                                                                                                   gtgtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                                                                                                >  1:221174/1‑275 (MQ=255)                                                                                                                                                    tgtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAATGGGGGTAACggaa                                                                                                                >  2:197172/1‑271 (MQ=255)                                                                                                                                                       aGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTc                                                                                                                                                                                                                        <  2:81567/167‑1 (MQ=255)                                                                                                                                                                aaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCAtgt                                                                                                                                                                                            >  2:161475/1‑186 (MQ=255)                                                                                                                                                                    cccccgacacGTAACCGCGCCTCTGGATA‑TTGTGTCTATGGACGTCGCGGAACGTTCACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgcacacggcagtgccaatcatcgcatctgaatgcccaaccagatatttggtggcggcttgaatagaaacatcgatgcca                                                                                               <  2:339252‑M1/271‑168 (MQ=255)                                                                                                                                                                    cacGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtg                                                                                                                                                                                        <  1:161477/186‑1 (MQ=255)                                                                                                                                                                     acGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGGTAACCCATTTCTCCTTACAACAAGCGATGTGTGCACTGGAACGGGGCGCAGGCcacggctggttatataggcggtatcggaatcgaccatctggcccatc                                                                                                                                                                                                                       >  1:440448‑M1/1‑107 (MQ=255)                                                                                                                                                                      cgcgACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgcacacggcagtgccaatcatcgcatctgaatgcccaaccagatatttggtggcggcttgaatagaaacatcgatgccaa                                                                                              >  1:339258‑M1/1‑106 (MQ=255)                                                                                                                                                                         gACACGTAATCGCGCCATTGGAGAGTTGT‑TCTCTGGACGGCGCGGAACGTTAACCCATTTCTCCTTCCACCAAGCGATGTGTGAACTAGAAGGTGGCGCAGGCcacggctggttatataggcggtatcggcatcgtccatctggcccatcagat                                                                                                                                                                                                                   <  2:440441‑M1/154‑52 (MQ=255)                                                                                                                                                                                     cgcgAATGTAGAGTGGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTTTGTACAACAAGCGATGTGTGAACTGGATGGTGGCGCAGTCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCg                                                                                    <  2:198403/269‑1 (MQ=255)                                                                                                                                                                                     cgcCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTc                                                                                                                                                                                                                             >  1:146044/1‑132 (MQ=255)                                                                                                                                                                                       ccAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCt                                                                                                                                                                                                                                  >  1:483534/1‑125 (MQ=255)                                                                                                                                                                                         aaTGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGGGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCAtt                                                                                                                                                                                                                         <  2:146043/132‑1 (MQ=255)                                                                                                                                                                                           tGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAAtt                                                                                                                                                                                                                              <  2:483525/125‑1 (MQ=255)                                                                                                                                                                                                  ttgt‑tCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGa                                                                                                                  >  1:397926/1‑226 (MQ=255)                                                                                                                                                                                                       tCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTg                                                                                                                                                                                                                                                                    <  1:240109/76‑1 (MQ=255)                                                                                                                                                                                                       tCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                                                                                              <  2:397919/226‑1 (MQ=255)                                                                                                                                                                                                                  cgcgGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTc                                                                                                                                                                                                                        >  1:135364/1‑109 (MQ=255)                                                                                                                                                                                                                      gAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGc                                                                                                                                                                                                                    <  2:135363/109‑1 (MQ=255)                                                                                                                                                                                                                                                                  gAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                                                                                                     >  1:531856/1‑176 (MQ=255)                                                                                                                                                                                                                                                                   aaGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTCCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATAGCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCAAGACGt  >  1:690122/1‑274 (MQ=255)                                                                                                                                                                                                                                                                      gTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCg                                                                                                 <  2:531847/176‑1 (MQ=255)                                                                                                                                                                                                                                                                            cAGGCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgcacacggcagtgccaatcatcgcatctgaatgcccaaccagatattt                                                                                                                              >  1:108241‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                                            cAGGCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcat                                                                                                                                                                               >  1:387668‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                                            cAGGCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccag                                                                                                                                                                                         >  1:243707‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                                              ggCcacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgcacacggcagtgccaatcatcttatctgactgccctaccagatatttggtg                                                                                                                          <  1:509044‑M1/143‑141 (MQ=255)                                                                                                                                                                                                                                                                                ccacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgcacacggcagtgccaatcatcgcatctgaatgcccaaccagatatttggtg                                                                                                                          <  2:108240‑M1/141‑141 (MQ=255)                                                                                                                                                                                                                                                                                ccacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaacgggcattgca                                                                                                                                                                           <  2:387661‑M1/92‑92 (MQ=255)                                                                                                                                                                                                                                                                                ccacggctggttatataggcggtatcggcatcgaccatctggcccatcagataggcattttcccgtagctgctcccagcaac                                                                                                                                                                                     <  2:243704‑M1/82‑82 (MQ=255)                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                     AGACATCCAGACGTATAAAAACAGGAATCCCGACATGGCGGACAAAAAGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGT‑TCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGT  >  CP009273/3145561‑3146091

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ 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.