New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3688674 =	0 (0.000)	21 (0.650)	18/268	0.9	56.0%	coding (500/3474 nt)	bcsC	cellulose synthase subunit
	?	NC_000913	= 4506221	33 (1.020)					coding (935/1119 nt)	yjhC	GFO/IDH/MOCA family putative oxidoreductase. NAD(P)‑dependent

CCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGCATCCGCATTGATGCGTTTTAACTGATTAATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCCTGCAATCGTGC  >  NC_000913/3688527‑3688819                                                                                                                                                    |                                                                                                                                                  cggtcagacaaagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggtaaggcaatatgacctcagatatggattgcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagAt                                                                                                                                                  >  4:114298‑M1/148‑149 (MQ=255)           aagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCg                                                                                                                                                <  1:96910‑M1/4‑1 (MQ=255)           aagcattttttggtccatgaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCg                                                                                                                                                >  2:96910‑M1/138‑141 (MQ=255)                             gaaacacaagaagaagatgatgatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTc                                                                                                                      <  2:6097‑M1/30‑1 (MQ=255)                                                  gatcgtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTc                                                                                                  <  7:81567‑M1/50‑1 (MQ=255)                                                      gtcggaaaggcaatatgacctcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGACACCGTACTCCAGTACTCGACAg                                                                                             >  3:107398‑M1/95‑149 (MQ=255)                                                                          tcagaaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCggcggc                                                                         <  8:45975‑M1/75‑1 (MQ=255)                                                                              aaatggatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCAc                                                                     <  6:7008‑M1/79‑1 (MQ=255)                                                                                   gatggcgctatagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGtt                                                                 <  1:15764‑M1/83‑1 (MQ=255)                                                                                             tagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTttg                                                       >  1:6097‑M1/56‑148 (MQ=255)                                                                                             tagcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTttg                                                       >  8:81567‑M1/56‑148 (MQ=255)                                                                                               gcatatggtcatccaggaaaaaaaacaccattatggcttgccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTttgttg                                                    >  7:45975‑M1/54‑149 (MQ=255)                                                                                                catatggtcatccaggaaaaaaaacaccattatggcttgccagattaataagATCGCTTCGCCACGGCGAGCCGCAATTTTCGCCACCGTACTCCAGTACTCGACAGCAAAGTCACCTTCCGGCGGCGCACCGTTGAACAGTttgttg                                                    <  1:60686‑M1/96‑1 (MQ=255)                                                                                                                                      gccagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGcc             <  3:19529‑M1/135‑1 (MQ=255)                                                                                                                                        cagtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCCTg           >  2:60686‑M1/13‑149 (MQ=255)                                                                                                                                          gtttaataagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCCTGc          <  1:38519‑M1/138‑1 (MQ=255)                                                                                                                                                 aagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCAt                           <  2:138236‑M1/121‑1 (MQ=255)                                                                                                                                                 aagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCAt                           >  7:37422‑M1/4‑124 (MQ=255)                                                                                                                                                 aagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCAt                           <  8:37422‑M1/121‑1 (MQ=255)                                                                                                                                                 aagATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCAt                           >  1:138236‑M1/4‑124 (MQ=255)                                                                                                                                                  agATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGAAATGTCACCTTCCGGCGGCGGACCGTTGAACAGTTTGTTGTAACTCGCCACATCTTCTTCTGCATGACCGGTCGTCGCCTGCAATCTTGc  >  4:19529‑M1/3‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCTTCGCGCCCGGCGTTCGATTTTGCCATCTGTTCCAGGACGGCAAAACCTTCGTCACGGCGATCGCTACTAAACAGCAATAGCGCCAGATTGTTTTGCAGGCCCGTATTGCCCGGTGCATCCGCATTGATGCGTTTTAACTGATTAATCGCTTCGCCACGGCGAGCCGGAATTTTCGCCACCGTACTCCAGTACTCGACAGCAATGTCACCTTCCGGCGGCGCACCGTTGAACAGTTTGTTGTAACTCGCCACAGCTTCTTCTGCATGACCGGTCGTCGCCTGCAATCGTGC  >  NC_000913/3688527‑3688819

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