New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3989869 =	80 (0.990)	3 (0.040)	3/294	9.9	3.6%	coding (898/942 nt)	hemC	hydroxymethylbilane synthase
	?	NC_000913	3989977 =	87 (1.140)					coding (790/942 nt)	hemC	hydroxymethylbilane synthase
Rejected: Position hash score below cutoff.

CGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3990073‑3989869 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgcgcCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGG  >  NC_000913/3989977‑3990104                                                                                                                                                                                                                                                                                                                                               CGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCTGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGG                                                                                                                                       <  2:442457/200‑1   CGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                     <  2:348584/200‑1     CCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCG                                                                                                                                   >  2:514755/1‑200                                                 AGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAAC                                                                                                                                         <  3:389355/150‑1                                                 AGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAAC                                                                                                                                         <  3:556944/150‑1                                                   CTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGG                                                                                                                                       <  3:502145/150‑1                                                   CTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGAGGATTCGCGGTGAACGCCGCGGTGCGCCGAAGGATGCCGAACAAATGAGGATTTCGCTGGCAGAAGAGCTACTGAATAACGG                                                                                                                                       <  3:211603/150‑1                                                   CTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTGCGCAGAGTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGG                                                                                                                                       <  3:266616/150‑1                                                            GCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCGGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGG                                                                                                                                       <  4:331217/141‑1                                                                      GGCGAAATCTGGCTGCCGGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGC                                                                                                                    <  3:350652/150‑1                                                                                                        CGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGGGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCT                                                                                           <  4:150008/141‑1                                                                                                                        TATTGGCGGTGAACGACGGGGTGGGCCGCAAGATGCCGAACAAATGGGGATTTCGCAGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAT                                                                           <  4:242136/141‑1                                                                                                                                                                                                          GCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGG  <  1:228533/132‑1                                                                                                                                                                                                                                                                                                                                               CGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3990073‑3989869 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgcgcCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGG  >  NC_000913/3989977‑3990104

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