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.

CCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTG  >  NC_000913/3686332‑3686600                                                                                                                            |                                                                                                                                                  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                                                                                             >  5:152001/1‑142 (MQ=255)                                     cGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGaa                                                                                             <  6:152001/142‑1 (MQ=255)                                                ggggAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTAcc                                                                           >  2:158172/1‑149 (MQ=255)                                                     gTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGatat                                                                      >  2:152942/1‑149 (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                                                                       >  5:148755/1‑124 (MQ=255)                                                                             cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGata                                                                       <  6:148755/124‑1 (MQ=255)                                                                             cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGcc                                              <  8:71555/149‑1 (MQ=255)                                                                             cAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGcc                                              <  2:173227/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)                                                                                                            ttAACGGCGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGa               <  8:111170/149‑1 (MQ=255)                                                                                                               aCGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACt             <  7:58155/148‑1 (MQ=255)                                                                                                                         ccAGCGGCCCGATATCATTGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATCTCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTg  >  5:117798/1‑149 (MQ=255)                                                                                                                         ccAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATc                  >  6:107945/1‑133 (MQ=255)                                                                                                                         ccAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATc                  <  5:107945/133‑1 (MQ=255)                                                                                                                         ccAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGc       >  4:17153/1‑144 (MQ=255)                                                                                                                         ccAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGc       <  3:17153/144‑1 (MQ=255)                                                                                                                            |                                                                                                                                                  CCACACCGTCGGCACGTACGCCACCCCATTTTTTCCCGGTATTGCTCGGGGAGTCTTTTTGCCCACCAAAGGCCAGCAAAGAACTGGAGATGGGCCGACGGTGGGCGTTAACGGTGTAACCCAGCGGCCCGATATCATCGCTGTAACTGATGCCGCCGACCACATCCACCACGTTGAAGCCCATCGGCGTGGTACCGATATCCCAGCTCCAGACGTCATTTCGCCAGCCGACCGCCACGCTGGCACCGGAATCCGACTGGCTGCGGTTG  >  NC_000913/3686332‑3686600

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.