New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3795313	18 (1.110)	3 (0.190)	3/248	NT	16.9%	coding (1184/1884 nt)	fcs	feruloyl‑CoA‑synthetase
	?	NC_002947	= 3795375	12 (0.770)					coding (1122/1884 nt)	fcs	feruloyl‑CoA‑synthetase

GTCGCCCACCGGCACCAGCTTCACTTCGCAGCCAGGTGCCGGCAGCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGACCAGTTCCTCCCA  >  NC_002947/3795234‑3795521                                                                                                                                              |                                                                                                                                                   gTCGCCCACCGGCACCAGCTTCACTTCGCAGCCAGGTGCCGGCAGCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAgtg                                                                                                                                            >  1:273192/1‑150 (MQ=255)                cAGCTTCACTTCGCAGCCAGGTGCCGGCAGCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCg                                                                                                                             >  2:132576/1‑150 (MQ=255)                       aCTTCCCACCCAGGTGCCGCCAACCCGACATAGCCGGTCACCGACCAAAGCCCGGTGGTGAATGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGc                                                                                                                      <  1:4324/150‑1 (MQ=255)                                    gTGCCGGCAGCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGGGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGc                                                                                                         >  2:39143/1‑150 (MQ=255)                                             gCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACa                                                                                                <  1:39143/150‑1 (MQ=255)                                              cccGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGt                                                                                                                           >  2:79827/1‑122 (MQ=255)                                              cccGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGt                                                                                                                           <  1:79827/122‑1 (MQ=255)                                                                           cGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGccaccgggcctttgtcgatggccggctatgtcgggctgccggcacctggctgcgaagtgaagctggtgccggtgggcgacaa                                                                  <  2:225946‑M2/150‑83 (MQ=255)                                                                                 tgAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCggg                                                            <  1:164167/150‑1 (MQ=255)                                                                                                                 ccAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGc                            >  2:266822/1‑150 (MQ=255)                                                                                                                 ccAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGAAACGCGGTCCAGCCCGTTCCAGACGCTGTGCTACAGGCCTGCGGCGTCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGc                            <  2:273192/150‑1 (MQ=255)                                                                                                                    aGGCCGGCCATCATGCGGATGCGTTCGccaccgggcctttgtcgatggccggctatgtcgggctgccggcacctgg                                                                                                   >  2:98562‑M2/1‑27 (MQ=255)                                                                                                                    aGGCCGGCCATCATGCGGATGCGTTCGccaccgggcctttgtcgatggccggctatgtcgggctgccggcacctgg                                                                                                   <  1:98562‑M2/76‑50 (MQ=255)                                                                                                                               catGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGAc              <  2:364875/150‑1 (MQ=255)                                                                                                                                 tGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGACCa            >  1:384/1‑150 (MQ=255)                                                                                                                                         cGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGACCAGTtcctcc    <  2:350410/150‑1 (MQ=255)                                                                                                                                           ttCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGACCAGTTCCTCCCa  <  1:185338/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   GTCGCCCACCGGCACCAGCTTCACTTCGCAGCCAGGTGCCGGCAGCCCGACATAGCCGGCCATCGACAAAGGCCCGGTGGTGAAGGTGCACGATGGCGAGGCTTCGGTCATGCCAAGGCCGGCCATCATGCGGATGCGTTCGCCACAGTGTTGCTCGGCAATGCGGTCCAGCCGGTCCCAGACGCTTTGCGACAGGCCTGCGGCGGCAAAGAAGAACAGCTTGATGCGGGCAAAGAACACCTCGCGTAGCGCGGGGTCCTGCTCCAGTGCCTTGACCAGTTCCTCCCA  >  NC_002947/3795234‑3795521

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