New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 800912	11 (0.370)	4 (0.140)	4/250	2.9	25.3%	coding (154/1053 nt)	ybhH	putative PrpF family isomerase
	?	NC_000913	= 800937	13 (0.460)					coding (179/1053 nt)	ybhH	putative PrpF family isomerase

GGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTGAAAATGGTT  >  NC_000913/800792‑801083                                                                                                                                                  |                                                                                                                                                   ggAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTcgc                                                                                                                                                 <  6:149424‑M2/149‑4 (MQ=255)         gAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCAttatta                                                                                                                                         <  5:149424/149‑1 (MQ=255)                            gCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCcgc                                                                                                                       <  7:90464/148‑1 (MQ=255)                                         ccGAAGATCAAACGCAGCGCCATAAAATATTGATGGCAATTATGGGGTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGCCCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGCTCAGCGTGCTGATGTCGAt                                                                                                         >  8:67435/1‑149 (MQ=255)                                           gAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTa                                                                                                       >  5:136331/1‑149 (MQ=255)                                                                              aaTTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTcgccaataccgtcaatttccagatcgttaccggaacccataattgccatcaatattttatcgcgctgcgtttgatcttcg                                                                    <  8:120318‑M2/149‑81 (MQ=255)                                                                                          cGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTcgccaataccgtcaatttccagatcgttaccggaacccataattgccatcaatattttatcgcgctgcgtttgatcttcgggtaaatgttc                                                         >  7:120318‑M2/1‑57 (MQ=255)                                                                                          cGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGc                                                        <  6:136331/149‑1 (MQ=255)                                                                                              aaCGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTaa                                                    >  1:139843/1‑149 (MQ=255)                                                                                                                  tATTGGCGGCGGTAATCCACTGACCAGTAAAGTcgccaataccgtcaatttccagatcgttaccggaacccataattgccatcaatattttatcgcgctgcgtttgatcttcgggtaaa                                                              >  5:55129‑M2/1‑33 (MQ=255)                                                                                                                  tATTGGCGGCGGTAATCCACTGACCAGTAAAGTcgccaataccgtcaatttccagatcgttaccggaacccataattgccatcaatattttatcgcgctgcgtttgatcttcgggtaaa                                                              <  6:55129‑M2/119‑87 (MQ=255)                                                                                                                            ggTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCCGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTACTCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGc                      <  7:67435/149‑1 (MQ=255)                                                                                                                                 tCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCg                 >  6:1644/1‑149 (MQ=255)                                                                                                                                    aCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGcc                                                       <  4:10721/108‑1 (MQ=255)                                                                                                                                    aCTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGcc                                                       >  3:10721/1‑108 (MQ=255)                                                                                                                                           aGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTGaaaa       >  4:5183/1‑149 (MQ=255)                                                                                                                                               aaGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTGAAAATGGt   >  5:30031/1‑149 (MQ=255)                                                                                                                                                aGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTGAAAATGGtt  <  5:1644/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCCATTATTAGCCGTTCCAGCGATCCGCGTGCTGATGTCGATTATCTGTTTGCGCAGGTAATCGTTCATGAGCAACGTGTCGATACCACGCCTAACTGCGGCAATATGCTGTCTGGTGTTGGGGCATTCGCCATTGAAAATGGTT  >  NC_000913/800792‑801083

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