New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5791892 =	44 (1.020)	3 (0.080)	3/226	NT	6.7%	intergenic (‑119/+76)	hemE/gltD	uroporphyrinogen decarboxylase/glutamate synthase subunit beta
	?	NC_002947	5791957 =	45 (1.170)					intergenic (‑184/+11)	hemE/gltD	uroporphyrinogen decarboxylase/glutamate synthase subunit beta

GCAGGAAACGATCGTTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCG  >  NC_002947/5791743‑5792041                                                                                                                                                      |                                                                                                                                                      ggcaagttcaagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgt                                                                                                                                                       <  2:520358‑M1/1‑1 (MQ=255)   caagttcaagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTc                                                                                                                                                     <  1:388186‑M1/3‑1 (MQ=255)     gAAACGATCGTTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTcgc                                                                                                                                                   <  2:494440/150‑1 (MQ=255)     gAAACGATCGTTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTcgc                                                                                                                                                   >  1:494440/1‑150 (MQ=255)      aaaCGATCGTTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGAGCAAATGGCACGGCCATGGAAGCGCGCTTTGGCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTcgcg                                                                                                                                                  <  1:690309/150‑1 (MQ=255)       ttcaagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCgg                                                                                                                                                 >  2:483062‑M1/144‑150 (MQ=255)       ttcaagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCgg                                                                                                                                                 <  1:320289‑M1/7‑1 (MQ=255)       aaCGATCGGTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCgg                                                                                                                                                 <  1:501340/150‑1 (MQ=255)        tcaagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCggg                                                                                                                                                >  2:685999‑M1/143‑150 (MQ=255)          aagcaccagaccagcaacccgaaagtgtttgccggtggcgacatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCGGGCa                                                                                                                                              >  1:943236‑M1/141‑150 (MQ=255)             caccagaccagcaacccgaaagtgtttgccggtggcgacctggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCGGGCAtgc                                                                                                                                           >  1:520358‑M1/138‑150 (MQ=255)                tcttcAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCAtgcctg                                                                                                                                        >  1:841145/1‑150 (MQ=255)                    cAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGcccc                                                                                                                                    >  2:1103399/1‑150 (MQ=255)                              aTGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCt                                                                                                                          <  1:687810/150‑1 (MQ=255)                               tGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTg                                                                                                                         <  2:662807/150‑1 (MQ=255)                                 tCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCAtgc                                                                                                                                           <  1:339168/130‑1 (MQ=255)                                 tCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCAtgc                                                                                                                                           >  2:339168/1‑130 (MQ=255)                                      aTCCGGAAAAATCGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGCTGCCACACATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTaa                                                                                                                  <  2:617401/150‑1 (MQ=255)                                            aaaaaTAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGc                                                                                                            <  2:841145/150‑1 (MQ=255)                                                  acatggtgcgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGc                                                                                                      <  2:989479‑M1/50‑1 (MQ=255)                                                       gggCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCa                                                                                                 >  1:53361/1‑150 (MQ=255)                                                          cgcggttcggacctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGc                                                                                              >  2:435158‑M1/93‑150 (MQ=255)                                                          cATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGCGCAGGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGc                                                                                              >  1:910258/1‑150 (MQ=255)                                                           aTTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGcc                                                                                             >  2:464495/1‑150 (MQ=255)                                                                 tcAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGaa                                                                                       >  2:650096/1‑150 (MQ=255)                                                                 tcAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGaa                                                                                       >  2:709416/1‑150 (MQ=255)                                                                   aGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACa                                                                                     <  1:464495/150‑1 (MQ=255)                                                                    acctggtggtgactgcgatcttcgaagggcgtacggctgccgagggtatcctggactacctcgaagtctgagcttcgacctgTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAg                                                                                    >  1:989479‑M1/83‑150 (MQ=255)                                                                          aGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGaa                                                                              <  2:285190/150‑1 (MQ=255)                                                                             gcAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCt                                                                           >  2:941657/1‑150 (MQ=255)                                                                               aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                         >  1:792332/1‑150 (MQ=255)                                                                               aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                         >  2:272793/1‑150 (MQ=255)                                                                               aaaaGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTACACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                         <  2:762564/150‑1 (MQ=255)                                                                                                     tgtTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATa                                                   <  2:162136/150‑1 (MQ=255)                                                                                                       ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                 >  1:778542/1‑134 (MQ=255)                                                                                                       ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                 <  2:778542/134‑1 (MQ=255)                                                                                                       ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                 <  2:775609/134‑1 (MQ=255)                                                                                                       ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGa                                                                 >  1:775609/1‑134 (MQ=255)                                                                                                       ttttGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATAcc                                                 >  2:388186/1‑150 (MQ=255)                                                                                                              tATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCa                                          <  1:1103399/150‑1 (MQ=255)                                                                                                                 tGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGcc                                       >  1:1085487/1‑150 (MQ=255)                                                                                                                        cacaGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCACAGTGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCTAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGcc                                >  1:29008/1‑150 (MQ=255)                                                                                                                         acaGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                         >  2:947185/1‑108 (MQ=255)                                                                                                                         acaGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCa                                                                         <  1:947185/108‑1 (MQ=255)                                                                                                                               cAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGaa                         >  1:453294/1‑150 (MQ=255)                                                                                                                                    gTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGccc                               <  2:272719/139‑1 (MQ=255)                                                                                                                                    gTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGccc                               >  1:272719/1‑139 (MQ=255)                                                                                                                                           tataCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGCGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTcac             >  2:24117/1‑150 (MQ=255)                                                                                                                                           tataCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTcac             <  1:709416/150‑1 (MQ=255)                                                                                                                                           tataCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTcac             <  1:404339/150‑1 (MQ=255)                                                                                                                                              aCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGa                                                                       >  1:865538/1‑89 (MQ=255)                                                                                                                                              aCAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGa                                                                       <  2:865538/89‑1 (MQ=255)                                                                                                                                               cAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCt                              <  2:729146/129‑1 (MQ=255)                                                                                                                                               cAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGGTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCt                              >  1:729146/1‑129 (MQ=255)                                                                                                                                                aaaaGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCAccacca        <  1:272793/150‑1 (MQ=255)                                                                                                                                                aaaaGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCAccacca        >  2:1075323/1‑150 (MQ=255)                                                                                                                                                    gCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTc    <  1:313630/150‑1 (MQ=255)                                                                                                                                                      TTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCg  <  2:320289‑M1/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCAGGAAACGATCGTTCTTCAGGGCAGTCATGTCGGCATCCGGAAAAATAGTGCGGGCATTTTCTCAGACACCAGCGCAAAAGGCACGGCCATGGCCATGTGTTTTGTCTATTGGATGCCACAGATCAAGCGTTTTTGTATACAAAAGCTTCGCGGGCATGCCTGCCCCCACGGGTTCTGCGCTTAACCCTGCGGGAGCGGGCATGCCCGCGAACAGGTCGAAGCTCAGACTTCGAGGTAGTCCAGGATACCCTCGGCAGCCGTACGCCCTTCGAAGATCGCAGTCACCACCAGGTCCG  >  NC_002947/5791743‑5792041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ 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.