New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	414338 =	13 (0.760)	3 (0.180)	3/268	1.9	20.3%	coding (153/1050 nt)	waaF	ADP‑heptose:LPS heptosyltransferase II
	?	NC_002947	414374 =	11 (0.660)					coding (189/1050 nt)	waaF	ADP‑heptose:LPS heptosyltransferase II

ATGGCGCAGACCCTGTTCCAGTGCCTGAAACAGCAGCACCCCGACTGCGTGATCGACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTCCCAAACGCACCGGCTGGCGCGGTGAAATGC  >  NC_002947/414231‑414519                                                                                                                                                |                                                                                                                                                  aTGGCGCAGACCCTGTTCCAGTGCCTGAAACAGCAGCACCCCGACTGCGTGATCGACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCATCGTGGCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAg                                                                                                                                             >  1:127263/1‑150 (MQ=255)              cagcgccgacttgagcgagttgggcagcaggatggcctggtcgtactggccggccagggacttgccgatgcgccgccgcgtggccagctccagcGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTg                                                                                                                                          >  1:402256‑M1/95‑140 (MQ=255)              cagcgccgacttgagcgagttgggcagcaggatggcctggtcgtactggccggccagggacttgccgatgcgccgccgcgtggccagctccagcGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTg                                                                                                                                          <  2:402256‑M1/46‑1 (MQ=255)                                   agcaCCCCGACTGCGTGATCGACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTgg                                                                                                           <  2:390694/150‑1 (MQ=255)                                                    aTCGACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAg                                                                                          <  2:198481/150‑1 (MQ=255)                                                       gACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCaa                                                                                                                   <  1:396782/122‑1 (MQ=255)                                                       gACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCCAGCGCATGCCGGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCaa                                                                                                                   >  2:396782/1‑122 (MQ=255)                                                            gCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCct                                                                                  >  2:260093/1‑150 (MQ=255)                                                                   ccTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTa                                                                                                 >  1:58427/1‑128 (MQ=255)                                                                   ccTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTAGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTa                                                                                                 <  2:58427/128‑1 (MQ=255)                                                                         tGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCATGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCCACTCg                                                                     >  2:82916/1‑150 (MQ=255)                                                                                                  cagctccagcgcgccgtggccgagcgggaagctcaaggcctggcgcGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTg                                                      >  1:402256‑M2/47‑140 (MQ=255)                                                                                                  cagctccagcgcgccgtggccgagcgggaagctcaaggcctggcgcGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTg                                                      <  2:402256‑M2/94‑1 (MQ=255)                                                                                                            gcgccgtggccgagcgggaagctcaaggcctggcgcGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCAt                                  >  1:183245‑M2/37‑150 (MQ=255)                                                                                                                CAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTccc                              <  2:376409‑M1/150‑1 (MQ=255)                                                                                                                                cGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTCCCAAACGCACCGGCTGgc              <  1:260093/150‑1 (MQ=255)                                                                                                                                         aCGGCGCGCTGGAGCTGGCCACGCGGCGGGGCCTCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTCCCAAACGCACCGGCTGGCGCGGTGaaa     <  2:127263/150‑1 (MQ=255)                                                                                                                                            gcacGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTCCCAAACGCTCCGGCTGGCGCGGTGAAATgc  >  2:213766‑M2/5‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                  ATGGCGCAGACCCTGTTCCAGTGCCTGAAACAGCAGCACCCCGACTGCGTGATCGACGTGCTGGCCCCTGAATGGAGCCGGCCGATCCTCGAGCGCATGCCTGAAGTGCGCCAGGCCTTGAGCTTCCCGCTCGGCCACGGCGCGCTGGAGCTGGCCACGCGGCGGCGCATCGGCAAGTCCCTGGCCGGCCAGTACGACCAGGCCATCCTGCTGCCCAACTCGCTCAAGTCGGCGCTGGTGCCCTTTTTTGCCGGCATTCCCAAACGCACCGGCTGGCGCGGTGAAATGC  >  NC_002947/414231‑414519

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