New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5196431 =	34 (1.470)	4 (0.190)	4/230	NT	12.1%	coding (132/690 nt)	ybgJ	allophanate hydrolase subunit 1
	?	NC_002947	5196476 =	27 (1.270)					coding (87/690 nt)	ybgJ	allophanate hydrolase subunit 1

GGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCC  >  NC_002947/5196327‑5196604                                                                                                                                                      |                                                                                                                                 ggCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATc                                                                                                                                  >  1:286745/1‑150 (MQ=255)         gCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCAt                                                                                                                          >  1:546098/1‑150 (MQ=255)               cAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAc                                                                                                                    <  2:564171/150‑1 (MQ=255)                aGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATg                                                                                                                              <  2:280525/139‑1 (MQ=255)                aGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATg                                                                                                                              >  1:280525/1‑139 (MQ=255)                     cTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATc                                                                                                              <  2:293651/150‑1 (MQ=255)                      tGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCg                                                                                                             >  1:571851/1‑150 (MQ=255)                         gCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATa                                                                                                          <  1:450474/150‑1 (MQ=255)                         gCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATa                                                                                                          >  1:529495/1‑150 (MQ=255)                                     cTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAgg                                                                                              >  2:489393/1‑150 (MQ=255)                                        gCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATg                                                                                                                         >  2:296446/1‑120 (MQ=255)                                        gCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATg                                                                                                                         <  1:296446/120‑1 (MQ=255)                                                tGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAg                                                                                   <  1:145923/150‑1 (MQ=255)                                                           cAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATgg                                                                        >  1:246280/1‑150 (MQ=255)                                                            aGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGc                                                                       <  2:144553/150‑1 (MQ=255)                                                              aCTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGa                                                                     >  2:510814/1‑150 (MQ=255)                                                                tacgcctgttcgaccgtatcgatgaggtcaacatgccatggATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCa                                                                                                                                >  2:603784‑M1/42‑89 (MQ=255)                                                                tacgcctgttcgaccgtatcgatgaggtcaacatgccatggATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCa                                                                                                                                <  1:603784‑M1/48‑1 (MQ=255)                                                                   aCCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACt                                                                >  1:456646/1‑150 (MQ=255)                                                                    ccATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACtt                                                               <  2:546098/150‑1 (MQ=255)                                                                    ccATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGCCGACCACtt                                                               >  2:62765/1‑150 (MQ=255)                                                                     cATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTc                                                              >  1:375955/1‑150 (MQ=255)                                                                      aTCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCg                                                             >  2:269223/1‑150 (MQ=255)                                                                                      aCGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtcatca                                             <  2:490763/150‑1 (MQ=255)                                                                                      aCGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtcatca                                             <  1:269223/150‑1 (MQ=255)                                                                                      aCGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCACGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTtcatca                                             <  1:208160/150‑1 (MQ=255)                                                                                                     gCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATc                              <  1:510814/150‑1 (MQ=255)                                                                                                      tggatgctcgccgccagccggcgcttgagcgcggcgttcggcgagcatCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTa                                                                                          <  1:603784‑M2/41‑1 (MQ=255)                                                                                                      tggatgctcgccgccagccggcgcttgagcgcggcgttcggcgagcatCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTa                                                                                          >  2:603784‑M2/49‑89 (MQ=255)                                                                                                        aaTGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAg                           <  2:246280/150‑1 (MQ=255)                                                                                                             TCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTAc                                                                                         >  1:612518‑M1/1‑83 (MQ=255)                                                                                                             TCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTAc                                                                                         <  2:612518‑M1/83‑1 (MQ=255)                                                                                                               GCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTgg                    >  2:571851‑M1/1‑150 (MQ=255)                                                                                                                   ccagccggcgcttgagcgcggcgttcggcgagcatCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACtt                                                               <  2:391618‑M2/68‑1 (MQ=255)                                                                                                                   ccagccggcgcttgagcgcggcgttcggcgagcatCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACtt                                                               >  1:391618‑M2/36‑103 (MQ=255)                                                                                                                           gcTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGAACTCATCGATACGGTCGAAGAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGGTTTGCGGCTACATAACACCTCCAGCCCATCGAGGGCCTGGCGGATCTTcccc        >  2:27705/1‑148 (MQ=255)                                                                                                                             tCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCAc      >  1:398849/1‑150 (MQ=255)                                                                                                                                 gcgcCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCAccgcc  <  2:398849/150‑1 (MQ=255)                                                                                                                                 gcgcCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCAGCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGTCGAATCTGCCGCAccgcc  <  1:489393/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                 GGCTGCAAGCCGTCCAGCGCCTGAGCGATTCGCGCCCTGGCCTCGCCTGGCGCCAGGGCAGACTGCACCATCAGCGTAGTATAAGACGGCACCAGGTCCAGCAAATGCTCGCCGAACGCCGCGCTCAAGCGCCGGCTGGCGGCGAGCATCCATGGCATGTTGACCTCATCGATACGGTCGAACAGGCGTACCATCAGGCTGTCGATGGCGACCACTTCGATTCGCGGCTTCATCACACCTCCAGCCCATCGAGGGCCTGGCGAATCTGCCGCACCGCC  >  NC_002947/5196327‑5196604

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