New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3689302 =	45 (1.150)	5 (0.130)	5/252	3.8	9.7%	coding (960/1107 nt)	bcsZ	endo‑1,4‑D‑glucanase
	?	NC_000913	3689326 =	49 (1.290)					coding (936/1107 nt)	bcsZ	endo‑1,4‑D‑glucanase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTAAACCGGTTGAGCAT  >  NC_000913/3689178‑3689472                                                                                                                                                     |                                                                                                                                                   ggCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCa                                                                                                                                                    >  3:146548/1‑149 (MQ=255)        gTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGt                                                                                                                                             >  7:14943/1‑149 (MQ=255)         tCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGtt                                                                                                                                            >  3:104617/1‑149 (MQ=255)           aGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGTGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGtttt                                                                                                                                          <  8:14943/149‑1 (MQ=255)           aGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGtttt                                                                                                                                          <  1:60235/149‑1 (MQ=255)             gTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTg                                                                                                                                         >  8:67657/1‑148 (MQ=255)               aataaCTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTa                                                                                                                                       >  1:153402/1‑148 (MQ=255)                   aCTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAg                                                                                                                                  <  1:41811/149‑1 (MQ=255)                   aCTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAg                                                                                                                                  <  2:153402/149‑1 (MQ=255)                    cTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAgg                                                                                                                                 >  3:164389/1‑149 (MQ=255)                     tCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAggg                                                                                                                                >  5:21199/1‑149 (MQ=255)                        ccTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGgcag                                                                                                                             <  4:146548/149‑1 (MQ=255)                         cTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGgcagc                                                                                                                            <  4:164389/149‑1 (MQ=255)                          ttttGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGgcagca                                                                                                                           >  2:80829/1‑149 (MQ=255)                                   gtaGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCGGGGAAAGTTATCGGCCACGCGCTGGCGCGGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAg                                                                                                                  <  1:87404/147‑1 (MQ=255)                                   gaAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATcgc                                                                                                                                                <  8:151403/119‑1 (MQ=255)                                   gaAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATcgc                                                                                                                                                >  7:151403/1‑119 (MQ=255)                                         gAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAgg                                                                                                                                 >  8:91279/1‑128 (MQ=255)                                         gAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAgg                                                                                                                                 <  7:91279/128‑1 (MQ=255)                                            aCGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGa                                                                                                         <  2:153536/149‑1 (MQ=255)                                                    gATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCt                                                                                                 >  5:46419/1‑149 (MQ=255)                                                       cccAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAggg                                                                                                                                >  1:50356/1‑115 (MQ=255)                                                       cccAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAggg                                                                                                                                <  2:50356/115‑1 (MQ=255)                                                        ccAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTAc                                                                                             >  8:21650/1‑149 (MQ=255)                                                            ccTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCGGCATGGCGGCAGAAAAACCGACTgg                                                                                                     >  1:93857/1‑137 (MQ=255)                                                            ccTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCGGCATGGCGGCAGAAAAACCGACTgg                                                                                                     <  2:93857/137‑1 (MQ=255)                                                            ccTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTg                                                                                         >  2:149267/1‑149 (MQ=255)                                                              ttGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTgcg                                                                                       >  6:13780/1‑149 (MQ=255)                                                                    aaaCAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGt                                                                                                                                        >  2:79405/1‑94 (MQ=255)                                                                    aaaCAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGt                                                                                                                                        <  1:79405/94‑1 (MQ=255)                                                                    aaaCAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGa                                                                                                         <  5:103492/125‑1 (MQ=255)                                                                    aaaCAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGa                                                                                                         >  6:103492/1‑125 (MQ=255)                                                                              aGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCcac                                                                       <  6:21199/149‑1 (MQ=255)                                                                               gCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCcaca                                                                      <  4:104617/149‑1 (MQ=255)                                                                                       ttGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACttt                                                              <  1:80829/149‑1 (MQ=255)                                                                                          tAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCt                                                                                          >  2:209254/1‑118 (MQ=255)                                                                                          tAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCt                                                                                          <  1:209254/118‑1 (MQ=255)                                                                                                     gCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAAc                                                >  5:179021/1‑149 (MQ=255)                                                                                                          cGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGttt                                           <  1:149267/149‑1 (MQ=255)                                                                                                           gggAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTCCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGtttt                                          >  6:123169/1‑149 (MQ=255)                                                                                                            ggAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGtttt                                          >  1:199200/1‑148 (MQ=255)                                                                                                               aaGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTtctc                                      <  5:13780/149‑1 (MQ=255)                                                                                                               aaGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTtctc                                      <  7:21650/149‑1 (MQ=255)                                                                                                                                   gcgctgAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGtt                  <  6:46419/149‑1 (MQ=255)                                                                                                                                      ctgAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTaa               <  2:199200/149‑1 (MQ=255)                                                                                                                                            ggCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTAAACCGGt         <  7:67657/149‑1 (MQ=255)                                                                                                                                              ccTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGtt                                            >  7:14685/1‑112 (MQ=255)                                                                                                                                              ccTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGtt                                            <  8:14685/112‑1 (MQ=255)                                                                                                                                                tGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTAAACCGGTTGAg     >  2:16934/1‑149 (MQ=255)                                                                                                                                                   gcATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTAAACCGGTTGAGCAt  <  5:123169/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTCAGTGAATGTCGCCATCGGTTTAAACCGGTTGAGCAT  >  NC_000913/3689178‑3689472

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