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.

GTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTC  >  NC_000913/3689160‑3689436                                                                                                                                               |                                                                                                                                       gtgAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACgcgctgg                                                                                                                                  >  2:60235/1‑149 (MQ=255)                  tGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTgcg                                                                                                                  >  2:41811/1‑148 (MQ=255)                   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                                                                                              <  2:153402/149‑1 (MQ=255)                                     aCTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAg                                                                                              <  1:41811/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                                                                                                            >  7:151403/1‑119 (MQ=255)                                                     gaAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATcgc                                                                                                            <  8:151403/119‑1 (MQ=255)                                                           gAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAgg                                                                                             <  7:91279/128‑1 (MQ=255)                                                           gAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAgg                                                                                             >  8:91279/1‑128 (MQ=255)                                                              aCGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGa                                                                     <  2:153536/149‑1 (MQ=255)                                                                      gATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCt                                                             >  5:46419/1‑149 (MQ=255)                                                                         cccAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAggg                                                                                            <  2:50356/115‑1 (MQ=255)                                                                         cccAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAggg                                                                                            >  1:50356/1‑115 (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                                                                     >  6:103492/1‑125 (MQ=255)                                                                                      aaaCAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGa                                                                     <  5:103492/125‑1 (MQ=255)                                                                                                aGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCcac                                   <  6:21199/149‑1 (MQ=255)                                                                                                 gCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCcaca                                  <  4:104617/149‑1 (MQ=255)                                                                                                         ttGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACttt                          <  1:80829/149‑1 (MQ=255)                                                                                                            tAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCt                                                      <  1:209254/118‑1 (MQ=255)                                                                                                            tAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCt                                                      >  2:209254/1‑118 (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  <  7:21650/149‑1 (MQ=255)                                                                                                                                 aaGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTtctc  <  5:13780/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                       GTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATCCCAGCCTTGTCCAAACAGGGTCAGCACATAGTTGTAATAGGCATCGCTGCCGGGAAAGTTATCGGCCACGCGCTGGCGCTGAACGGCCTGCGCATCGCGGTTTTGTAAAAAGGGCAGCATGGCGGCAGAAAAACCGACTGGTCCTTTACCCTGCGCTTTCCCCGTAGCCACATCCACTTTTTCCGGCGGATAACCGTTTTTCTC  >  NC_000913/3689160‑3689436

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.