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

GATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAT  >  NC_000913/969299‑969427                                                              |                                                                    ggtttgtttaacaggccttgtgcggcaaggtgaagattacgcgcatcagcaagggagaggtGCGGCAg                                                               <  2:441461‑M1/7‑1 (MQ=255) gatgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAg                                                               >  1:461345/1‑68 (MQ=255)      gtttaacaggccttgtgcggcaaggtgaagattacgcgcatcagcaagggagaggtGCGGCagcggca                                                          <  1:491328‑M1/12‑4 (MQ=255)      aCCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTc                                                          >  2:199032/1‑68 (MQ=255)                    ttACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCa                                            <  2:523949/68‑1 (MQ=255)                      gcggcaaggtgaagattacgcgcatcagcaagggagaggtGCGGCAGCGTTCGATTGATGGAGTCATg                                          >  1:14541‑M1/41‑68 (MQ=255)                             tAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCg                                   <  1:311784/68‑1 (MQ=255)                                ccTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTg                                >  1:29034/1‑68 (MQ=255)                                     cTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCa                           >  2:230384/1‑68 (MQ=255)                                      tGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                          >  2:82365/1‑68 (MQ=255)                                          cgcatcagcaagggagaggtGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCGCctct                      <  2:541401‑M1/48‑1 (MQ=255)                                                        tACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCTGCACCTCTCCCTTGCTGATgcg        <  2:29034/68‑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)                                                              |                                                                    GATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAAT  >  NC_000913/969299‑969427

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