New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3901949	92 (0.790)	3 (0.030)	3/244	NT	3.0%	coding (272/1617 nt)	bglH	carbohydrate‑specific outer membrane porin, cryptic
	?	NC_000913	= 3901970	110 (1.010)					coding (251/1617 nt)	bglH	carbohydrate‑specific outer membrane porin, cryptic
Rejected: Frequency below/above cutoff threshold.

TGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3901795‑3901949 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaccagtaGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTATTACCGTT  <  NC_000913/3901970‑3901824                                                                                                                                                                                                                                                                                                                TGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTG                                                                                                                                                            <  4:594140/148‑1  GAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGAC                                                                                                                                                          >  2:332204/1‑149  GAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGA                                                                                                                                                           <  4:426287/148‑1    ACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACC                                                                                                                                                         >  4:975263/1‑148     CCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAG                                                                                                                                                       <  2:538715/149‑1                                                            TTTCATCTTTCACAAATTTACTGAATTCTGTCAGAGACATGGAAGAATAAATCGAGCCTCTCTGACTGGCATTTATCTCTTCATTTTTGACCAGTAGTGCCGCAAAACCTGATGCTGTCCTGGTCAAAAATGAGGAGAAAAATGCCAG                                                                                                 <  4:886137/148‑1                                                                 TCTTTCACACATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTC                                                                                                                        >  2:345701/1‑120                                                                 TCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTC                                                                                                                        <  1:345701/120‑1                                                                                                                                                    GACCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTATT        >  1:722389/1‑149                                                                                                                                                    GACCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTAT         <  4:640181/148‑1                                                                                                                                                    GACCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAACGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTAT         >  4:630175/1‑148                                                                                                                                                      CCAGTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTATTAC      >  1:65973/1‑149                                                                                                                                                         GTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTATTACCGT   >  1:770908/1‑149                                                                                                                                                         GTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAA                <  2:1059945/136‑1                                                                                                                                                         GTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAA                >  1:1059945/1‑136                                                                                                                                                         GTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAA                >  1:831227/1‑136                                                                                                                                                         GTAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAA                <  2:831227/136‑1                                                                                                                                                          TAGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGTAATTGGCTTTAGTTACAACGGCTATTACCGTT  >  2:191021/1‑149                                                                                                                                                                                                                                                                                                                TGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGAAGAATAAATCGAGCCTGTCTCACTGGCATTTTTCTCTTCATTTTTGACCAGTA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3901795‑3901949 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gaccagtaGTGCCGCAAAACCTGATGCTGTACTGGTCAAAAATGAAGAGAAAAATGCCAGTGAGACAGGCTCGATTTATTCTTCCATGACTCTGAAAGATTTCAGTAAATTTGTGAAAGATGAAATTGGCTTTAGTTACAACGGCTATTACCGTT  <  NC_000913/3901970‑3901824

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