New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 5578021	59 (1.100)	3 (0.060)	3/252	NT	5.3%	coding (163/810 nt)	rhdA	thiosulfate sulfurtransferase
	?	NC_002947	= 5578058	51 (0.990)					coding (200/810 nt)	rhdA	thiosulfate sulfurtransferase

GCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCAACCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5577875‑5578021 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctcgCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGATCAGTTGCGGCGAATCCAGGCGT  <  NC_002947/5578058‑5577912                                                                                                                                                                                                                                                                                                        GCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCAACCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCG                                                                                                                                                     <  6:145453/147‑1                                          CTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCATCCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCGCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCT                                                                                                          <  8:316101/149‑1                                             GATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCAACCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCGCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGG                                                                                                        <  2:102701/148‑1                                              ATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCAACCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCGCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGT                                                                                                      <  3:103749/149‑1                                                                                                                                                 TCGCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGATCAGTTGCGGCGAATCCAGGCG   >  2:121576/1‑149                                                                                                                                                 TCGCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGATCAGTTGCGGCGAATCCAGGCG   >  7:248435/1‑149                                                                                                                                                   GCCGAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGTTCTGTTGCGGCGAATCCAGGCGT  <  5:94170/148‑1                                                                                                                                                   GCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGATCAGTTGCGGCGAATCCAGGCGT  <  6:180954/148‑1                                                                                                                                                   GCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTG                                               >  6:113555/1‑103                                                                                                                                                                                                                                                                                                        GCCTGCCCCTGGTGATCGAAGCTGCAGACCTGCTGCCACGCCTGGATTCGCCGCAACTGATCCTGGTCGACCTGAGCAGCGCCAACCGCTATGTCAGCGGGCATATCCCCGGCGCGCGCTTTGTCGAAGGCAAGCGCACCCAGCTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5577875‑5578021 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑agctcgCCCAGGGCTGGCAGCAGGCCGGGCGCCGGGGGCTGGCCGAGCTGGGTGCGCTTGCCTTCGACAAAGCGCGCGCCGGGGATATGCCCGCTGACATAGCGGTTGGCGCTGCTCAGGTCGACCAGGATCAGTTGCGGCGAATCCAGGCGT  <  NC_002947/5578058‑5577912

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.