New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4129541 =	78 (0.640)	3 (0.030)	3/240	NT	3.9%	coding (870/1161 nt)	metB	cystathionine gamma‑synthase, PLP‑dependent
	?	NC_000913	4129565 =	72 (0.620)					coding (894/1161 nt)	metB	cystathionine gamma‑synthase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

GCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAA  >  NC_000913/4129394‑4129686                                                                                                                                                    |                                                                                                                                                  gctACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCcgcgc                                                                                                                                                  >  1:1111416/1‑149 (MQ=255)   tgccatgcctgcatgtgtcatggttgcggcgtgagagattaaactttccactccccctaatgattccgccagcgtaaacaacgacagcccgcccaggaaacgacgcagcgtctgctcatcgtcatccagttcaaaactcaacattgcgc                                                                                                                                                <  2:515105‑M1/4‑1 (MQ=255)    gccatgcctgcatgtgtcatggttgcggcgtgagagattaaactttccactccccctaatgattccgccagcgtaaacaacgacagcccgcccaggaaacgacgcagcgtctgctcatcgccatccagttcaaaactcaacattGCGcc                                                                                                                                               <  4:165438‑M1/5‑1 (MQ=255)     cGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCa                                                                                                                                              <  3:1036098/149‑1 (MQ=255)     cGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATGGCCGCGCGCCa                                                                                                                                              <  1:602379/149‑1 (MQ=255)             gCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAg                                                                                                                                      <  4:238443/149‑1 (MQ=255)               gAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAg                                                                                                                                      >  2:760513/1‑147 (MQ=255)               gAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAg                                                                                                                                      <  1:760513/147‑1 (MQ=255)               gAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGc                                                                                                                                    <  2:832644/149‑1 (MQ=255)                  cGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCttt                                                                                                                                 <  1:542651/149‑1 (MQ=255)                     tGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGgc                                                                                                                              >  2:792815/1‑149 (MQ=255)                         gCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCaa                                                                                                                          <  3:51735/149‑1 (MQ=255)                          ccGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAAt                                                                                                                         <  1:137477/149‑1 (MQ=255)                           cgcgTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATg                                                                                                                        >  1:542515/1‑149 (MQ=255)                                aTGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAg                                                                                                                   <  2:313078/149‑1 (MQ=255)                                   gagagattaaactttccactccccctaatgattccgccagcgtaaacaacgacagcccgcccaggaaacgacgcagcgtctgctcatcgccatccagttcaaaactcaacattGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGttt                                                                                                                >  1:1077725‑M1/114‑149 (MQ=255)                                   gAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGttt                                                                                                                <  2:1111416/149‑1 (MQ=255)                                     gCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTg                                                                                                              >  3:142158/1‑149 (MQ=255)                                      cTGGCGAAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCCCCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGa                                                                                                             <  1:139958/149‑1 (MQ=255)                                          cgcagcgcaACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTg                                                                                                         <  2:1109802/149‑1 (MQ=255)                                           gcagcgcaACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAg                                                                                                                                      <  3:1123298/119‑1 (MQ=255)                                           gcagcgcaACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAg                                                                                                                                      >  4:1123298/1‑119 (MQ=255)                                             agcgcaACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGAt                                                                                                      <  4:1070764/149‑1 (MQ=255)                                               cgcaACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATgg                                                                                                    <  2:321222/149‑1 (MQ=255)                                                    actccccctaatgattccgccagcgtaaacaacgacagcccgcccaggaaacgacgcagcgtctgctcatcgccatccagttcaaaactcaacattGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAg                                                                                                                   <  3:500448‑M1/33‑1 (MQ=255)                                                    actccccctaatgattccgccagcgtaaacaacgacagcccgcccaggaaacgacgcagcgtctgctcatcgccatccagttcaaaactcaacattGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAg                                                                                                                   >  4:500448‑M1/97‑129 (MQ=255)                                                     gcgcAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCGCCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGa                                                                                              <  3:44666/149‑1 (MQ=255)                                                      cgcAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAg                                                                                             <  1:158740/149‑1 (MQ=255)                                                             caTTGTGAAATACCTCCAAACCCAGCCGTTGGGGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGAGCATGAAATTGCCGCGCGCCAGCAAAAACGCTTTTGCGCAATGTGGAGTGTTTAACTGGGTGTCGAGGAGCAGACGc                                                                                      >  4:975030/2‑149 (MQ=255)                                                              aTTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATg                                                                                                     <  2:268221/133‑1 (MQ=255)                                                              aTTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATg                                                                                                     >  1:268221/1‑133 (MQ=255)                                                                  tgAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGcgt                                                                                 <  2:542515/149‑1 (MQ=255)                                                                           tGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTg                                                                        <  2:212164/149‑1 (MQ=255)                                                                            gCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGcgt                                                                                 <  2:875636/139‑1 (MQ=255)                                                                            gCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGcgt                                                                                 >  1:875636/1‑139 (MQ=255)                                                                            gCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgg                                                                       >  4:126575/1‑149 (MQ=255)                                                                              aaaCCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggc                                                                     >  4:1127768/1‑149 (MQ=255)                                                                                aCCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTTGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgg                                                                   >  3:987177/1‑149 (MQ=255)                                                                                aCCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgg                                                                   >  1:503953/1‑149 (MQ=255)                                                                                aCCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgg                                                                   >  3:293994/1‑149 (MQ=255)                                                                                  ccAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgggc                                                                 >  1:116956/1‑149 (MQ=255)                                                                                      ccGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTc                                                             <  2:503953/149‑1 (MQ=255)                                                                                        gTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCgt                                                           >  2:188888/1‑149 (MQ=255)                                                                                         ttGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCgtt                                                          >  2:67754/1‑149 (MQ=255)                                                                                           ggTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCgttgt                                                        <  1:792815/149‑1 (MQ=255)                                                                                            gTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgggc                                                                 >  3:693226/1‑139 (MQ=255)                                                                                            gTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTgggcgggc                                                                 <  4:693226/139‑1 (MQ=255)                                                                                             tGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGttt                                                      <  3:126575/149‑1 (MQ=255)                                                                                                 aaaaCTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGc                                                  >  1:74461/1‑149 (MQ=255)                                                                                                 aaaaCTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGc                                                  >  3:232827/1‑149 (MQ=255)                                                                                                      tGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTc                                                                           >  4:438542/1‑119 (MQ=255)                                                                                                      tGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTc                                                                           <  3:438542/119‑1 (MQ=255)                                                                                                      tGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCg                                             >  3:847205/1‑149 (MQ=255)                                                                                                          tCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAAt                                         >  2:628193/1‑149 (MQ=255)                                                                                                          tCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAAt                                         <  2:116956/149‑1 (MQ=255)                                                                                                          tCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAAt                                         >  1:167489/1‑149 (MQ=255)                                                                                                           cACCCGTCGTTGCCGGAAAATCAGGGGCATGAGATTGCCGCGCGCCAGCAAGGAGGGTTGGGTGCAATGTTGGGTGTTGAAGTGGATGGCGTGGAGCAGACGCTGCGTGGTTTCCTTGGCGGGCTGTCGTTGTTTACCCTGGCGGAATc                                        <  4:1036/149‑1 (MQ=255)                                                                                                           cACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATc                                        >  1:408520/1‑149 (MQ=255)                                                                                                           cACCCGTCGTTGCCGAAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATc                                        <  3:975030/149‑1 (MQ=255)                                                                                                             cccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCTGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGc                                              >  1:814571/1‑141 (MQ=255)                                                                                                             cccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGc                                              <  2:814571/141‑1 (MQ=255)                                                                                                              ccGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGACGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCAtt                                     >  1:515105/1‑149 (MQ=255)                                                                                                               cgtcgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTa                                    <  1:188888/149‑1 (MQ=255)                                                                                                                gtcgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAg                                   >  2:128184/1‑149 (MQ=255)                                                                                                                  cgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGGATCATTAggg                                 <  3:642142/149‑1 (MQ=255)                                                                                                                  cgtTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAggg                                 <  1:628193/149‑1 (MQ=255)                                                                                                                         ggAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGaa                                          >  3:472614/1‑133 (MQ=255)                                                                                                                         ggAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGaa                                          <  4:472614/133‑1 (MQ=255)                                                                                                                          gAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGCCGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGa                         <  4:446164/149‑1 (MQ=255)                                                                                                                           aaaaTCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGaa                        <  2:95074/149‑1 (MQ=255)                                                                                                                            aaaTCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGaaa                       <  1:67754/149‑1 (MQ=255)                                                                                                                               tCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGtt                    >  2:1148142/1‑149 (MQ=255)                                                                                                                               tCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGtt                    >  3:587493/1‑149 (MQ=255)                                                                                                                                  ggggCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCgt                                                           >  2:56730/1‑107 (MQ=255)                                                                                                                                  ggggCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCgt                                                           <  1:56730/107‑1 (MQ=255)                                                                                                                                  ggggCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGtttta                 <  2:74461/149‑3 (MQ=255)                                                                                                                                      cATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTCGCGGAATCATTTGGGGTAGTGGAAAGTTTTAcctc             <  4:232827/149‑1 (MQ=255)                                                                                                                                       aTGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAAtctct            >  3:1013826/1‑149 (MQ=255)                                                                                                                                        tGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAAtctctc           >  3:165438/1‑149 (MQ=255)                                                                                                                                         gAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAAt                                         >  3:928544/1‑118 (MQ=255)                                                                                                                                         gAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAAt                                         <  4:928544/118‑1 (MQ=255)                                                                                                                                           aaTTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGtt                    <  3:796444/137‑1 (MQ=255)                                                                                                                                           aaTTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGtt                    >  4:796444/1‑137 (MQ=255)                                                                                                                                           aaTTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCAcg        >  4:447118/1‑149 (MQ=255)                                                                                                                                                 cgcgcgCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCaa  >  2:775594/1‑149 (MQ=255)                                                                                                                                                 cgcgcgCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCaa  >  3:333508/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTACGTGGGTTGCGAACGCTGGTGCCGCGTATGGAGCTGGCGCAGCGCAACGCGCAGGCGATTGTGAAATACCTGCAAACCCAGCCGTTGGTGAAAAAACTGTATCACCCGTCGTTGCCGGAAAATCAGGGGCATGAAATTGCCGCGCGCCAGCAAAAAGGCTTTGGCGCAATGTTGAGTTTTGAACTGGATGGCGATGAGCAGACGCTGCGTCGTTTCCTGGGCGGGCTGTCGTTGTTTACGCTGGCGGAATCATTAGGGGGAGTGGAAAGTTTAATCTCTCACGCCGCAA  >  NC_000913/4129394‑4129686

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