New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3688674 =	0 (0.000)	29 (0.920)	26/268	0.3	60.4%	coding (500/3474 nt)	bcsC	cellulose synthase subunit
	?	NC_000913	= 4506221	38 (1.210)					coding (935/1119 nt)	yjhC	GFO/IDH/MOCA family putative oxidoreductase. NAD(P)‑dependent

TTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGCATCCGCATTGATGCGTTTTAACTGATTAATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCC  >  NC_000913/3688529‑3688808                                                                                                                                                  |                                                                                                                                       ttcagacaaagcattttttggtcgatgaaacacatcaagaagatgatgatcgtcggaaaggcaatatgccctcagaaatggatggcgatatagcatatggtgatccaggaaaaaaaacaccataatggcttgccagtttaataagATCg                                                                                                                                     <  4:10465‑M1/4‑1 (MQ=255)     gacaaagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTc                                                                                                                                 >  1:70062‑M1/142‑149 (MQ=255)        aaagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGc                                                                                                                               <  1:119131‑M1/10‑1 (MQ=255)        aaagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGc                                                                                                                               <  8:16314‑M1/10‑1 (MQ=255)          agcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCAc                                                                                                                            <  7:89884‑M1/13‑1 (MQ=255)                     gtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCgg                                                                                                                  >  8:89884‑M1/126‑148 (MQ=255)                                acaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCAc                                                                                                      >  7:16314‑M1/115‑149 (MQ=255)                                             gatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTAc                                                                                         <  2:114981‑M1/48‑1 (MQ=255)                                                gatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTc                                                                                       <  1:13522‑M1/50‑1 (MQ=255)                                                       ggaaaggaaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaccaccattatgtcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCaa                                                                               >  3:10465‑M1/92‑149 (MQ=255)                                                                      cctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCcggcg                                                                >  1:114981‑M1/77‑149 (MQ=255)                                                                       ctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTc                                                                           <  8:15208‑M1/62‑1 (MQ=255)                                                                       ctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTc                                                                           >  7:15208‑M1/76‑137 (MQ=255)                                                                          agaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTc                                                                                       <  6:31052‑M1/50‑1 (MQ=255)                                                                          agaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTc                                                                                       >  5:31052‑M1/73‑122 (MQ=255)                                                                             aatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCAcc                                                         <  4:81949‑M1/80‑1 (MQ=255)                                                                                 gatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGtt                                                      >  3:81949‑M1/66‑148 (MQ=255)                                                                                    ggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGaa                                                   <  8:99608‑M1/86‑1 (MQ=255)                                                                                          atagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTttg                                            <  4:33209‑M1/93‑1 (MQ=255)                                                                                            agcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTttgt                                           >  7:99608‑M1/55‑148 (MQ=255)                                                                                                     tcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGc                                 >  3:33209‑M1/46‑149 (MQ=255)                                                                                                           ccgatctaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGc                           >  1:3340‑M1/40‑149 (MQ=255)                                                                                                           aggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGc                           >  2:119131‑M1/40‑149 (MQ=255)                                                                                                            ggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACAGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGct                          <  3:96301‑M1/111‑1 (MQ=255)                                                                                                                   aaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGc                           <  2:3340‑M1/110‑1 (MQ=255)                                                                                                                         cattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGAc             >  8:88989‑M1/26‑149 (MQ=255)                                                                                                                           ttatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCg           >  2:13522‑M1/24‑149 (MQ=255)                                                                                                                                gcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGgtcgt      <  3:98005‑M1/131‑1 (MQ=255)                                                                                                                                    gccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGcc  <  2:136752‑M1/135‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                       TTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGCATCCGCATTGATGCGTTTTAACTGATTAATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCC  >  NC_000913/3688529‑3688808

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.