New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1812849 =	57 (0.880)	6 (0.090)	6/258	NT	9.7%	coding (328/1116 nt)	frmA	glutathione‑dependent formaldehyde dehydrogenase
	?	NC_002947	1812891 =	56 (0.880)					coding (370/1116 nt)	frmA	glutathione‑dependent formaldehyde dehydrogenase

TGGTGCCGTCCGGCATCAGGCCCTTGCCCTGAGTGGCGCGGATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/1812894‑1812849 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCTGCCTGCTGGGTTGCGGC  >  NC_002947/1812891‑1813037                                                                                                                                                                                                   TGGTGCCGTCCGGCATCAGGCCCTTGCCCTGAGTGGCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATC                                               >  2:380249/1‑148  GGTGCCGTCCGGCATCAGGCCCTTGCCCTGAGTGGCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATC                                               >  4:116031/1‑147    TGCCGTCCGGCATCAGGCCCTTGCCCTGAGTGGCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAG                                            >  6:311722/1‑148                  AGGCCCTTGCCCTGAGTGGCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCT                              <  3:116031/148‑1                            CCTGAGTGGCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCT                   <  5:311722/149‑1                                    GCGCGGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCTGCCTGCTG           <  1:380249/149‑1                                        GGATGGCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCTGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCTGCCTGCTGGGTT       <  4:357211/149‑1                                             GCACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCTGCCTGCTGGGTTGCGGC  <  8:54176/149‑1                                                                                                                                                                                                   TGGTGCCGTCCGGCATCAGGCCCTTGCCCTGAGTGGCGCGGATGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/1812894‑1812849 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcACCACGCGCTTCTCGTACAAGGGCCAGCAGCTGTTCCACTACATGGGCACCTCGACGTTCTCCGAATACACCGTGCTGCCGGAAATTTCCGTGGCCAAGATCCAGAAAGAAGCACCCCTGGAAAAGGTCTGCCTGCTGGGTTGCGGC  >  NC_002947/1812891‑1813037

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.