New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_tpiA	4465486 =	28 (0.970)	5 (0.180)	4/264	3.4	16.8%	coding (677/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
	?	NC_000913_3_hsa_tpiA	4465514 =	23 (0.850)					coding (649/1422 nt)	treB	trehalose‑specific PTS enzyme: IIB and IIC component
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGACGGGTACCACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGA  >  NC_000913_3_hsa_tpiA/4465369‑4465657                                                                                                                                                  |                                                                                                                                                cAGACGGCTACCACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCca                                                                                                                                              <  2:478151/149‑1 (MQ=255)   gACGGGTACCACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCcaca                                                                                                                                            <  1:400717/149‑1 (MQ=255)                         agGTAATCCGGGACGATGGGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAg                                                                                                                      <  2:488767/149‑1 (MQ=255)                               tCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGGTa                                                                                                                >  2:138780/1‑149 (MQ=255)                                  ggCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGAGAGCCTACTTGGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGc                                                                                                             <  1:227349/149‑1 (MQ=255)                                             tttagtgCGGGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCt                                                                                                  <  2:206879/143‑1 (MQ=255)                                             tttAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGAGAATGGTAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCt                                                                                                  >  1:498873/1‑149 (MQ=255)                                                  ggCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGTCCGAGCAGATAAGCGTTCATCAGCTGTTga                                                                                             >  2:109391/1‑149 (MQ=255)                                                  ggCGAGTTTAAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGga                                                                                             <  2:48859/149‑1 (MQ=255)                                                   gCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTACCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACAATTCCGGCAGCTGCTGCCCGAGCAGATAACCGTTCATCAGCTGTGgag                                                                                            >  1:100497/1‑149 (MQ=255)                                                           cAATAACGCCCAGCGCCCGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAgt                                                                                    <  1:423131/149‑1 (MQ=255)                                                             aTAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGTAGATAAGCGTTCATCAGCTGTGGAGAAACCAgtgt                                                                                  <  1:273009/149‑1 (MQ=255)                                                                            gttCCGGCTAACAGTGCCGGGATCACCCGCGAAAAAACGTCAACTTTTTCTATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGTAGATAAGCGTTCATCTGCTGTGGAGAAACCAGTGTCACGCCAAGCACGAt                                                                   <  2:236104/147‑1 (MQ=255)                                                                            aGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCTGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGAt                                                                   <  1:441164/149‑1 (MQ=255)                                                                                 ggCGAACAGTGCCGGGATCACCTGCGCCTGATAGGCTACGTTGGGGATGCTGAACATGCCAAAGGCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCGGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCaa                                                              <  2:245389/149‑1 (MQ=255)                                                                                         gtccctggATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGc                                                      <  1:173293/143‑1 (MQ=255)                                                                                                          gcCTGATAGTCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCAttttttt                                     <  1:254851/149‑1 (MQ=255)                                                                                                              gATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCAcc                                 >  1:59208/1‑149 (MQ=255)                                                                                                               aTAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCg                                <  2:27543/149‑1 (MQ=255)                                                                                                                 aGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCt                              <  1:135151/149‑1 (MQ=255)                                                                                                                   ccTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCtaa                            >  2:19029/1‑147 (MQ=255)                                                                                                                   ccTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGa                            >  2:474787/1‑149 (MQ=255)                                                                                                                               aTGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTTGCCCGCAAATccc                >  2:186358/1‑149 (MQ=255)                                                                                                                                         cgCCAAAGTCCCACACCTCCGGCATCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAgg      <  1:138780/148‑1 (MQ=255)                                                                                                                                             aaaGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTaga  >  2:51470/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                CAGACGGGTACCACCACCAGATAGAGGTAATCCGGCACGATGCGTTTAAGGCGAGTTTCAATAACGCCCAGCGCCAGTCCGGCTAACAGTGCCGGGATCACCTGCGCCTGATAGCCTACTTTGGCGATGCTGAACATGCCAAAGTCCCACACTTCCGGCAGCTGCTGCCCGAGCAGATAAGCGTTCATCAGCTGTGGAGAAACCAGTGTCACGCCAAGCACGATACCAAGGATCGGCGTGCCGCCCATTTTTTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGA  >  NC_000913_3_hsa_tpiA/4465369‑4465657

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