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.

AATGCCAGACCGAAAGAAATTGAGCAACTCCCTATGGTGAAAGGCTGTTTCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4466057‑4466003
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttcaccaTAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGA  >  NC_000913/4466024‑4466148
                                                                                                                                                                                    
AATGCCAGACCGAAAGAAATTGAGCAACTCCCTATGGTGAAAGGCTGTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGT                                                            <  5:116160/122‑1
AATGCCAGACCGAAAGAAATTGAGCAACTCCCTATGGTGAAAGGCTGTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGT                                                            >  6:116160/1‑122
                               CTATGGTGAAAGGCTGTTTCACCATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGA  >  2:23330/1‑149
                                  TGGTGAAAGGCTGTTTCACCATAGGGAGTTGCTCAATTTCTTT                                                                                                         <  1:18360/43‑1
                                  TGGTGAAAGGCTGTTTCACCATAGGGAGTTGCTCAATTTCTTT                                                                                                         >  2:18360/1‑43
                                                     CATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATAT                                  >  7:48262/1‑95
                                                     CATAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATAT                                  <  8:48262/95‑1
                                                                                                                                                                                    
AATGCCAGACCGAAAGAAATTGAGCAACTCCCTATGGTGAAAGGCTGTTTCACCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4466057‑4466003
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tttcaccaTAGGGAGTTGCTCAATTTCTTTCGGTCTGGCATTGGCCGGTTGGTTGAGGACAAAGCGTAGGCGAGTAATACAGTGGCTCACCGTCGCAATATTGCCGCGCCCGCCGACCAGTTCAATCAACCGA  >  NC_000913/4466024‑4466148

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