New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	969360 =	10 (0.530)	3 (0.180)	3/114	1.5	27.8%	intergenic (+8/‑29)	lpxK/ycaQ	lipid A 4'kinase/DUF1006 family protein with C‑terminal wHTH domain
	?	NC_000913	969402 =	7 (0.430)					coding (14/1233 nt)	ycaQ	DUF1006 family protein with C‑terminal wHTH domain

CTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACC  >  NC_000913/969335‑969462                                                                    |                                                             cTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCa                                                              >  2:230384/1‑68 (MQ=255)  tGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                                                             >  2:82365/1‑68 (MQ=255)   ccgcgcatcagcatgggagaggtgcggcagcgacattcatgactccatcactctaccgctgccgcAcc                                                            <  1:700589‑M2/3‑1 (MQ=255)      cgcatcagcaagggagaggtGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCGCctct                                                         <  2:541401‑M1/48‑1 (MQ=255)                    tACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCTGCACCTCTCCCTTGCTGATgcg                                           <  2:29034/68‑1 (MQ=255)                    agaggcgcggcagcgacattcatgactccatcaatcgaacgctgccgcACCTCTCCCTTGCTGATgcg                                           <  2:541401‑M2/20‑1 (MQ=255)                     aCGCCGCGGCAGCGTGCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATgcgc                                          <  1:194880/68‑1 (MQ=255)                       ggtGCGGCAGCGGTAGAGTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCATGCTGATgcgcgg                                        <  1:700589‑M1/65‑2 (MQ=255)                       gccgcGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGt                                        >  1:108836/1‑68 (MQ=255)                          GCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAt                                     >  1:351498‑M1/1‑68 (MQ=255)                                 CGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCAcc                              >  1:230384‑M1/1‑68 (MQ=255)                                        catgactccatcaatcgaacgctgccgcACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGc                       >  1:14541‑M2/29‑68 (MQ=255)                                               aGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                <  2:614091/68‑1 (MQ=255)                                                  cATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTg             >  2:295043/1‑68 (MQ=255)                                                        tgccgctgccgcACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTaaac       <  1:491328‑M2/56‑1 (MQ=255)                                                             ctgccgcACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAAcc  <  2:441461‑M2/61‑1 (MQ=255)                                                                    |                                                             CTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACC  >  NC_000913/969335‑969462

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