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.

ACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTG  >  NC_000913/4465903‑4466172                                                                                                                          |                                                                                                                                                     aCCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGa                                                                                                                                                 >  2:158057/1‑127 (MQ=255) aCCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGa                                                                                                                                                 <  1:158057/127‑1 (MQ=255)    tGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTgg                                                                                                                         >  1:79283/1‑148 (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)                                                       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                                           >  4:120283/1‑135 (MQ=255)                                                                                               cGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGa                                           <  3:120283/135‑1 (MQ=255)                                                                                                    ttGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAg                                                                                                      <  2:67428/71‑1 (MQ=255)                                                                                                    ttGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAg                                                                                                      >  1:67428/1‑71 (MQ=255)                                                                                                                        cATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAatat                                                          >  7:48262/1‑95 (MQ=255)                                                                                                                        cATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAatat                                                          <  8:48262/95‑1 (MQ=255)                                                                                                                          tAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTg  <  3:52005/149‑1 (MQ=255)                                                                                                                          |                                                                                                                                                     ACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTGACGCAATCAGTGCTTGATAGTAATCACCCACGTTGGTGCCAATCACCACCTGAAATTGCCCGGCATTGGTGAAACAGCCTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGATCGATATCCGTTTGGTTTATTTTG  >  NC_000913/4465903‑4466172

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.