New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4466003 =	33 (1.010)	5 (0.160)	5/230	3.2	14.5%	coding (178/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
	?	NC_000913	4466024 =	28 (0.920)					coding (157/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAA  >  NC_000913/4465856‑4466144                                                                                                                                                    |                                                                                                                                              agaTCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTg                                                                                                                                              <  7:45758/149‑1 (MQ=255)              aTGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCtt                                                                                                                                 <  1:119429/149‑1 (MQ=255)                                                aCCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGa                                                                                                                     <  1:158057/127‑1 (MQ=255)                                                aCCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGa                                                                                                                     >  2:158057/1‑127 (MQ=255)                                                   tGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTgg                                                                                             >  1:79283/1‑148 (MQ=255)                                                   tGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCAGGCATTg                                                                                                                                              >  6:43092/1‑99 (MQ=255)                                                              aaCCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGt                                                                                                                   <  8:89267/115‑1 (MQ=255)                                                              aaCCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGt                                                                                                                   >  7:89267/1‑115 (MQ=255)                                                                 cTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCggttggt                                                                              <  2:79283/149‑1 (MQ=255)                                                                  tGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCggttggtt                                                                             >  3:26581/1‑149 (MQ=255)                                                                    cgcCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCggttggtt                                                                             <  4:26581/147‑1 (MQ=255)                                                                     gcCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAg                                                                          >  1:99509/1‑149 (MQ=255)                                                                       cTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGa                                                                        >  8:140206/1‑149 (MQ=255)                                                                            cGGTTGACGCTATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACaaa                                                                    <  4:65896/148‑1 (MQ=255)                                                                                     cAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGa                                                          <  2:99509/149‑1 (MQ=255)                                                                                               ttGATGGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGt                                                >  4:35026/1‑149 (MQ=255)                                                                                               ttGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCAtt                                                                                         >  8:3539/1‑108 (MQ=255)                                                                                               ttGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCAtt                                                                                         <  7:3539/108‑1 (MQ=255)                                                                                                  aTAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTc                                                                                                                               >  8:121929/1‑67 (MQ=255)                                                                                                  aTAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTc                                                                                                                               <  7:121929/67‑1 (MQ=255)                                                                                                      tAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCAc                                         >  8:90920/1‑149 (MQ=255)                                                                                                                  ttGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATAtt                             >  4:52005/1‑149 (MQ=255)                                                                                                                            tCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGcgcccg                   >  4:95062/1‑149 (MQ=255)                                                                                                                                  ccTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAAtttct                                                                                                       <  7:134601/59‑1 (MQ=255)                                                                                                                                  ccTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAAtttct                                                                                                       >  8:134601/1‑59 (MQ=255)                                                                                                                                      aaaTTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGt          <  3:95062/148‑1 (MQ=255)                                                                                                                                        aTTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCa            >  7:4583/1‑144 (MQ=255)                                                                                                                                        aTTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCa            <  8:4583/144‑1 (MQ=255)                                                                                                                                             ccGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTtcaatcaa  >  8:40487/1‑149 (MQ=255)                                                                                                                                              cGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGa               <  3:120283/135‑1 (MQ=255)                                                                                                                                              cGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGa               >  4:120283/1‑135 (MQ=255)                                                                                                                                                   ttGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAg                                                                          <  2:67428/71‑1 (MQ=255)                                                                                                                                                   ttGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAg                                                                          >  1:67428/1‑71 (MQ=255)                                                                                                                                                    |                                                                                                                                              AGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAA  >  NC_000913/4465856‑4466144

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