New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	572405 =	683 (16.650)	13 (0.420) +AGCGTTTGACCTCTGCGGAGGGGAAGCGTGA	10/202	NT	2.4%	coding (38/102 nt)	ylcH	DLP12 prophage; uncharacterized protein YlcH
	?	NC_000913	= 1095504	NA (NA)					intergenic (‑67/‑40)	insE‑4/ymdE	IS3 element protein InsE/protein YmdE
Rejected: Frequency below/above cutoff threshold.

AGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGTCAATTATCAA  >  NC_000913/1095367‑1095641                                                                                                                                          |                                                                                                                                          gtggCTTCACTTCGGAATTCTGGCGGATGCTGGTTATGTGGTTTTTTTCTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTg                                                                                                                                      <  2:722799/141‑1 (MQ=11)   ggCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTTATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTCCTGGGTAAGATCAGTGCCTGGc                                                                                                                                    <  2:271900/143‑1 (MQ=25)      ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCCCTCTGAGTGTTTCCTCACTTAGGCGCGTTTCCACTATTGCTGGGTAAGATCAGTGCCTGGcaga                                                                                                                                 >  1:110515/1‑141 (MQ=11)         aCTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                                                              <  1:120544/143‑1 (MQ=40)           tGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGt                                                                                                                            >  1:681276/1‑143 (MQ=255)                    ccggCGAATGATGTTTACGGGGTTTTTTACTGGTGGATACTGTTTTTGGCATGTGAGTCACCTCTGACTGAGAGTTTCCTCACGTAGCCGAGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGaa                                                                                                                   <  1:81596/141‑1 (MQ=17)                         gAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACgg                                                                                                              >  1:636612/1‑143 (MQ=255)                           aTGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACGTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCg                                                                                                            <  1:321682/143‑1 (MQ=255)                                     cGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                                                  <  1:693644/143‑1 (MQ=255)                                        tgTTTTTTGCTGGTTGATGCTGTTTTTGCCATGTGCGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTCGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATgcg                                                                                               <  2:351776/142‑1 (MQ=21)                                                agtggtggGTACTGTGTTTGTCAGGTGGGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCGACTATTGCGGGGTAAGATCAGTGCCTGGCACTCGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                                                       <  2:37104/141‑1 (MQ=21)                                                   ccttGGTACGGTTTTTGGCGTGTGAGTCACCTCTGTGGGCGAGTTTACTCACTGAGCCGCGTGTCCACTATTGCTGGGGAAGAGCAGTGCCTCGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTg                                                                                    <  2:129270/141‑1 (MQ=21)                                                          aCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTATGATGGCAGACGGCGCGCAGGCAATGCGCGATTATCTTGCCGGACt                                                                             >  2:386823/1‑143 (MQ=255)                                                           cTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTg                                                                            >  1:242400/1‑143 (MQ=255)                                                            tGTTTTTGTGATGTGTGTGACCTCTGACTGAGAGCTTATTCTCTTAGTCGCGTGTTCTCTATTGCTGTGTAAGATCAGTGCCTGGCACACGCGTTTGATGGGAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTgg                                                                           <  2:374367/143‑1 (MQ=255)                                                            tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCCATGCGCGATTATCTTGCCCGACTgg                                                                           >  2:470951/1‑143 (MQ=255)                                                            tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAGTGCGCGATTATCTTGCCGGACTgg                                                                           >  2:470733/1‑143 (MQ=255)                                                                  tGTCATGTGAGTCACCTCTGACTGAGAGGTTGCTGACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCg                                                                     <  1:518274/143‑1 (MQ=255)                                                                                  ttcgactgaGTGTTTGCTTTTGTAGCGGGGTGTTCACTATTGTTGGGTAAGTGCAGTGCCTGGCACACGCGTTTGCGGGAATACGGCGTTTAGTGAATGCGCGATTATCTTGTCGTACTCGATATCGCTTCGCCAGAGTAtga                                                     <  2:476028/140‑3 (MQ=255)                                                                                   ctgactgaGTGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAg                                                    <  1:537047/143‑1 (MQ=255)                                                                                    ggactgaGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                                                   <  1:39908/142‑1 (MQ=255)                                                                                                   tCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTaa                                    >  1:667122/1‑143 (MQ=255)                                                                                                     aCTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTACCa                                  >  2:596853/1‑143 (MQ=255)                                                                                                                           gCTGGGTAAGGTGAGTGGCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGt            <  2:117632/143‑1 (MQ=255)                                                                                                                           gCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGt            >  1:535456/1‑143 (MQ=255)                                                                                                                                 tAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGTCAATTa      >  1:19816/1‑143 (MQ=255)                                                                                                                                     aTCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGTCAATTATCaa  <  1:449645/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                          AGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGGTCGCTCCGTCAATTATCAA  >  NC_000913/1095367‑1095641

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