New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4465545 =	20 (0.590)	4 (0.120)	4/258	3.7	16.0%	coding (636/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
	?	NC_000913	4465581 =	23 (0.710)					coding (600/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAG  >  NC_000913/4465446‑4465729                                                                                                                                        |                                                                                                                                                     aaTAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAgtg                                                                                                                                         >  6:130890/1‑149 (MQ=255)              gcCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCAc                                                                                                                            <  8:10147/149‑1 (MQ=255)                                          ccTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtt                                                                                                >  3:67731/1‑149 (MQ=255)                                           cTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAttt                                                                                               >  4:20739/1‑149 (MQ=255)                                                    aTAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCAcc                                                                                       <  4:67731/148‑1 (MQ=255)                                                                          aaCATGCCAAAGTCCCACACTTACGGCAGCTTCTGCCCGCGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGACAAGCAAGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGAACgg                                                                >  8:30470/1‑149 (MQ=255)                                                                               gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgaccg                                                                 >  3:12566/1‑143 (MQ=255)                                                                               gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgaccg                                                                 <  4:12566/143‑1 (MQ=255)                                                                                   aaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAtt                                                                                                >  8:33488/1‑108 (MQ=255)                                                                                     gTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaaga                                                     <  3:20739/149‑1 (MQ=255)                                                                                     gTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaag                                                      >  6:44797/1‑148 (MQ=255)                                                                                        ccACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTagaag                                                      >  4:91396/1‑145 (MQ=255)                                                                                                 cGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCAc                                         <  3:91396/149‑1 (MQ=255)                                                                                                             cccGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACACCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACca                               >  3:12900/1‑147 (MQ=255)                                                                                                                    aGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATc                      >  6:146662/1‑149 (MQ=255)                                                                                                                      aTAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGt                    <  4:28489/149‑1 (MQ=255)                                                                                                                       tAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTa                   <  5:44797/149‑1 (MQ=255)                                                                                                                       tAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCAAAGCACGATACCAAGGATAGGCGTGCCGACCATTTTTTTAACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCTATCAACCACAGAAAATCGTa                   <  7:30470/149‑1 (MQ=255)                                                                                                                          gCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTa                   <  3:64329/146‑1 (MQ=255)                                                                                                                          gCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTa                   >  4:64329/1‑146 (MQ=255)                                                                                                                               catcaGCTGTGGAGAAACCAGTGTCACACCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGtt            <  4:12900/148‑1 (MQ=255)                                                                                                                                atcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTg          >  3:105709/1‑149 (MQ=255)                                                                                                                                 tcaGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGc         <  5:146662/149‑1 (MQ=255)                                                                                                                                   aGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCa        <  5:130890/148‑1 (MQ=255)                                                                                                                                   aGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAg       <  5:40785/149‑1 (MQ=255)                                                                                                                                        tgGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAg  >  5:91745/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                     AATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAG  >  NC_000913/4465446‑4465729

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.