New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2331858	107 (0.860)	5 (0.040)	5/282	11.2	4.6%	pseudogene (2545/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
	?	NC_000913	= 2331859	106 (0.900)					pseudogene (2544/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Rejected: Position hash score below cutoff.

GCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331653‑2331858 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttgcgccTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCGTATTTCACGTTCCTGATGCCTG  <  NC_000913/2331859‑2331723                                                                                                                                                                                                                                                                                                                                                         GCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATATTTACTCGCCGCTCGCTGCGGTTAGT                                                                                                                                                 <  4:921732/200‑1                                                                                                                                           GTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATG                                                         >  2:227200/1‑150                                                                                                                                                CACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCGTATTTCACGTTCCTGATGCCTG  >  4:91084/1‑200                                                                                                                                                                             TTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAG                                >  3:494549/1‑141                                                                                                                                                                                              CTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCG                        >  1:496808/1‑132                                                                                                                                                                                                     TAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCGTATTTCACGTTCCTGA        >  3:815880/1‑141                                                                                                                                                                                                       GTTGCACCTGGCGCAACTAACCGAAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCGTATTACACG               <  1:175296/132‑1                                                                                                                                                                                                                                                                                                                                                         GCACGCCCCTGCCCGACCAGCCCGTCGCCGTTCATCCCACGCGCGGTGATACGCCAGCGGGTTAACGAATCAGGCATCAGGAACGTGAAATACGCTTTGCCTTGTTTATCGGTTGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2331653‑2331858 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttgcgccTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGTCACTCACAACCGATAAACAAGGCAAAGCGTATTTCACGTTCCTGATGCCTG  <  NC_000913/2331859‑2331723

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