New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3686428	32 (0.970)	4 (0.130)	4/246	3.5	10.9%	coding (2746/3474 nt)	bcsC	cellulose synthase subunit
	?	NC_000913	= 3686455	35 (1.110)					coding (2719/3474 nt)	bcsC	cellulose synthase subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGG  >  NC_000913/3686303‑3686574                                                                                                                              |                                                                                                                                                   ttCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTgg                                                                                                                                             >  8:57366/1‑133 (MQ=255) ttCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTgg                                                                                                                                             <  7:57366/133‑1 (MQ=255)                              ccACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAACGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGCGTAACCCAGCGGCCCGATATCATCGCTGTAACt                                                                                                 >  7:111170/1‑148 (MQ=255)                              ccACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACt                                                                                                 >  4:136341/1‑148 (MQ=255)                              ccACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACt                                                                                                 >  7:71555/1‑148 (MQ=255)                              ccACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACt                                                                                                 <  3:136341/148‑1 (MQ=255)                                                      cccATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGat                                                                                                                <  1:13734/109‑1 (MQ=255)                                                      cccATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGat                                                                                                                >  2:13734/1‑109 (MQ=255)                                                      cccATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACg                                                                        >  1:173227/1‑149 (MQ=255)                                                      cccATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACg                                                                        >  8:120271/1‑149 (MQ=255)                                                        cATTTTTTCGCGGTATTGCTCGGGGAGTCGTTTTGCCCCCCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGGTACCGGTGTACCCCAGCGGCCCGAGATCATCGCTGTACCTGATGCCGCCGACCACATCCACCCCGtg                                                                      >  8:58155/1‑148 (MQ=255)                                                           tttttCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAATGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCCGCGGCCCGATGTCATCGCTGTAAGTGATGCCGCCGCCCCCATCCACCACGTTGac                                                                   >  7:136248/1‑148 (MQ=255)                                                                  cGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGaa                                                                   <  6:152001/142‑1 (MQ=255)                                                                  cGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGaa                                                                   >  5:152001/1‑142 (MQ=255)                                                                             ggggAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTAcc                                                 >  2:158172/1‑149 (MQ=255)                                                                                  gTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGatat                                            >  2:152942/1‑149 (MQ=255)                                                                                     tttttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACg                                                                                                                                    <  7:82276/58‑1 (MQ=255)                                                                                     tttttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACg                                                                                                                                    >  8:82276/1‑58 (MQ=255)                                                                                     tttttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGCCCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATcc                                          >  1:176904/1‑148 (MQ=255)                                                                                     tttttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCg                                                           >  6:89230/1‑131 (MQ=255)                                                                                     tttttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATcc                                          >  2:7989/1‑148 (MQ=255)                                                                                       tttGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATccc                                         <  7:120271/147‑1 (MQ=255)                                                                                        ttGCCCACCAAAGGCCAGCAAAGAACGGGAGATGGGCCGACGGTGGGCGGTAACGGGGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAg                                       <  8:136248/148‑1 (MQ=255)                                                                                                      ccAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCAttt                        <  1:152942/149‑1 (MQ=255)                                                                                                          cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGata                                             <  6:148755/124‑1 (MQ=255)                                                                                                          cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGata                                             >  5:148755/1‑124 (MQ=255)                                                                                                          cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGcc                    <  2:173227/149‑1 (MQ=255)                                                                                                          cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGcc                    <  8:71555/149‑1 (MQ=255)                                                                                                                    gagaTGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCAGCGTGGTACCGATATCCCAGCTCCAGAc                               >  2:76438/1‑128 (MQ=255)                                                                                                                    gagaTGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCAGCGTGGTACCGATATCCCAGCTCCAGAc                               <  1:76438/128‑1 (MQ=255)                                                                                                                           gCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTg   <  1:7989/149‑1 (MQ=255)                                                                                                                            ccGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTgg  >  4:126698/1‑149 (MQ=255)                                                                                                                              |                                                                                                                                                   TTCACCTTTATCGTAGCTCAGACTTAGCCCCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGG  >  NC_000913/3686303‑3686574

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