New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5196440	21 (0.710)	4 (0.150)	4/226	NT	17.3%	coding (123/690 nt)	ybgJ	allophanate hydrolase subunit 1
	?	NC_002947	= 5196465	19 (0.700)					coding (98/690 nt)	ybgJ	allophanate hydrolase subunit 1

TGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAAT  >  NC_002947/5196320‑5196591                                                                                                                                                  |                                                                                                                               tgtgTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggc                                                                                                                            >  2:401007‑M2/1‑146 (MQ=255)       aGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcat                                                                                                                      <  2:517311‑M2/150‑11 (MQ=255)           tGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCAt                                                                                                                  <  1:46000/150‑1 (MQ=255)               aGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCa                                                                                                              >  1:32299/1‑150 (MQ=255)                         cgcCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGGCTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATACTCGCCGAACGCCGCGCTCAAGCGCAGGCTGGCGGCGAGCATCCATGGCATGTTGACCTc                                                                                                    >  2:107139/1‑150 (MQ=255)                           ccTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCAt                                                                                                  >  2:384316/1‑150 (MQ=255)                               aGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGAt                                                                                              <  1:517311/150‑1 (MQ=255)                                gCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATa                                                                                             <  1:250295/150‑1 (MQ=255)                                           ccTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAg                                                                                  <  1:249045/150‑1 (MQ=255)                                                ccTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGNGAGCATNCATGGCATGTTGNCCTCATCGATACGGTCGAACAGGCGTa                                                                             <  1:768145/150‑1 (MQ=255)                                                      cTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCa                                                                       <  1:384316/150‑1 (MQ=255)                                                        ggCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAgg                                                                     >  1:329367/1‑150 (MQ=255)                                                        ggCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAgg                                                                     <  2:32299/150‑1 (MQ=255)                                                                 gCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATg                                                            <  2:580309/150‑1 (MQ=255)                                                                  cAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggc                                                                                                                            <  1:225855‑M2/85‑5 (MQ=255)                                                                  cAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggc                                                                                                                            >  2:225855‑M2/1‑81 (MQ=255)                                                                     aCTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGa                                                        >  2:460836/1‑150 (MQ=255)                                                                           ccATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACtt                                                  <  1:460836/150‑1 (MQ=255)                                                                                       tataAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGACTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCt                                      >  1:577194/1‑150 (MQ=255)                                                                                             aCGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtcatca                                >  2:614843/1‑150 (MQ=255)                                                                                                   ccAGGTCCAGCAAATGCTCGCCGAACGCCGCGCCAAATCCAGGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTc                          <  1:439565/150‑1 (MQ=255)                                                                                                          cAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctggacggcttg                   <  2:176965‑M2/150‑110 (MQ=255)                                                                                                          cAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCa                   >  2:738817/1‑150 (MQ=255)                                                                                                               aaTGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAg              <  2:524358/150‑1 (MQ=255)                                                                                                                    tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcagg                                 <  1:484671‑M2/126‑96 (MQ=255)                                                                                                                    tCGCCGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcagg                                 >  2:484671‑M2/1‑31 (MQ=255)                                                                                                                        cGAACGCCGCGCTCAAGCGCCGGCTGGcggcgagcatttgctggacctggtgccgtcttatactacgctgatggtgcagtctgccctggcgccaggcgaggccagggcgcgaatcgctcaggcgctggacggcttgcagcctgacacagg     <  2:381629‑M2/150‑124 (MQ=255)                                                                                                                           aCGCCGCGCTCAAGCGCCAGCGGGAGGCGAGTATCCATGTCATGTTGCCCTCATCGATACGGTCGAACAGCCGTACCAGCAGGCTGTCGATGGCGACCCCTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAAt  <  1:107139/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                               TGTGTCAGGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAAT  >  NC_002947/5196320‑5196591

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