New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3685751	43 (1.260)	6 (0.180)	6/260	3.2	14.9%	coding (3423/3474 nt)	bcsC	cellulose synthase subunit
	?	NC_000913	= 3685795	27 (0.820)					coding (3379/3474 nt)	bcsC	cellulose synthase subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATT  >  NC_000913/3685647‑3685923                                                                                                                                                     |                                                                                                                                 ccgccCGATTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCt                                                                                                                                  <  8:132260/149‑1 (MQ=255)  cgccCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTg                                                                                                                                 <  1:131373/149‑1 (MQ=255)  cgccCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTg                                                                                                                                 >  8:92495/1‑149 (MQ=255)            tagccGATCGCAGCGGCATTAAGTGAGGCGATATCTGAAAACTTACCAGTCGGCGTTAGGTATCAGCGGCTGCGGCGGTAAATCAATGTGACCCGGCCATCCGTCGGCGGAATAAAGTTCGTAGAGCAGGAAATGGCTGGGTGCGTaa                                                                                                                        <  4:99776/145‑1 (MQ=255)            taCCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAAt                                                                                                                       <  1:165780/149‑1 (MQ=255)             aCCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATc                                                                                                                      >  6:7475/1‑149 (MQ=255)             aCCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATc                                                                                                                      >  6:9502/1‑149 (MQ=255)               ccGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCtt                                                                                                                    >  2:58859/1‑149 (MQ=255)                   tCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGc                                                                                                                <  6:94865/149‑1 (MQ=255)                             aTTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGata                                                                                                      <  5:132847/149‑1 (MQ=255)                             aTTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAAGCTTTCGCCTGCTGGata                                                                                                      >  3:28440/1‑149 (MQ=255)                                     agGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCt                                                                                                                                  >  1:150357/1‑113 (MQ=255)                                     agGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCt                                                                                                                                  <  2:150357/113‑1 (MQ=255)                                     agGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGc                                                                                              >  3:103175/1‑149 (MQ=255)                                        cgcTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGt                                                                                           <  4:103175/149‑1 (MQ=255)                                               tGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCaa                                                                                                  <  3:27707/135‑1 (MQ=255)                                               tGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCaa                                                                                                  >  4:27707/1‑135 (MQ=255)                                                  aaaCTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAAc                                                                                 <  5:7475/149‑1 (MQ=255)                                                  aaaCTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAAc                                                                                 <  5:9502/149‑1 (MQ=255)                                                    aCTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCa                                                                               <  1:58859/149‑1 (MQ=255)                                                        aCCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGGGCCAACAAACCAGtt                                                                            <  3:111572/148‑1 (MQ=255)                                                         ccAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTgg                                                                          <  1:133714/149‑1 (MQ=255)                                                                 cGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAaca                                                                  <  7:92495/149‑1 (MQ=255)                                                                      ggTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACAcgtcg                                                             >  6:89381/1‑149 (MQ=255)                                                                            aGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAg                                                       <  5:143751/149‑1 (MQ=255)                                                                                   gcggcgGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCcaac                                                                                     >  2:41985/1‑112 (MQ=255)                                                                                   gcggcgGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCcaac                                                                                     <  1:41985/112‑1 (MQ=255)                                                                                            aaTCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCgt                                        <  5:89381/148‑1 (MQ=255)                                                                                                              aTCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTc                                                                                                                  >  7:37889/1‑56 (MQ=255)                                                                                                              aTCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTc                                                                                                                  <  8:37889/56‑1 (MQ=255)                                                                                                                 cggcggcggAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAAGCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGActgct                   <  4:28440/148‑1 (MQ=255)                                                                                                                                 tacgtaGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATAGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCACGCGACGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGAtt  >  4:33985/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                 CCGCCCGAGTATACCCGATCGCAGCGGCATTAAGAGAGGCGCTATCTGAAAACTTACCAGTCGGCGTAAGGTATCAGCGGCTGCGGCGGTAAATCCATGTCACCCTGCCATCCGGCGGCGGAATAACGTACGTAGAGCAGGAAATGGCTGGGTGCGTAATCTTTCGCCTGCTGGATATCAATTGCCGTGCCAACAAACCAGTTGGAAGTAACACGTCGTTCAAGTAATGCCCGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATT  >  NC_000913/3685647‑3685923

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.