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.

ACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCAC  >  NC_000913/4465399‑4465690                                                                                                                                                   |                                                                                                                                                  accaccaGATAGAGGTAATCCGGCACGAGGCGTTTCAGGCGAGTTGCAATCACGCCCAGCGCCAGGCCGGCTAACAGTGCCGGGATCACCTGCGCCTGCTAGCCTACTTTGGCGATGCGGAACATGCCAAAGTCCCACACTTCCGGCAg                                                                                                                                                 <  7:51435/149‑1 (MQ=255)     ccaGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAgctgc                                                                                                                                             >  3:28489/1‑149 (MQ=255)         aTAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAgct                                                                                                                                               <  2:156631/143‑1 (MQ=255)         aTAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAgct                                                                                                                                               >  1:156631/1‑143 (MQ=255)             agGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCAGCTGCCGGaa                                                                                                                                     <  1:79581/149‑2 (MQ=255)                  aTCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGa                                                                                                                                 >  7:60897/1‑148 (MQ=255)                             tGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTtcatc                                                                                                                     <  8:60897/149‑1 (MQ=255)                                                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                          <  4:12566/143‑1 (MQ=255)                                                                                                                              gCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATaccgaccg                          >  3:12566/1‑143 (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)                                                                                                                                                   |                                                                                                                                                  ACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCAC  >  NC_000913/4465399‑4465690

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.