New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3989869 =	114 (0.950)	4 (0.030)	3/504	NT	3.6%	coding (898/942 nt)	hemC	hydroxymethylbilane synthase
	?	NC_000913	3989977 =	102 (0.880)					coding (790/942 nt)	hemC	hydroxymethylbilane synthase
Rejected: Frequency below/above cutoff threshold.

TATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3990161‑3989869 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgcgcCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA  >  NC_000913/3989977‑3990260                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   TATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGCAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACACATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGC                                                                                                                                                                                                                                                                                                  >  2:77572/1‑289  tTTGAGTGTCGCCGTGATGATGCACGGACTCGGCATCTGGTTGCCGAGCTGAATCACCACGAAACTGTACTGCGCGTTACCGCAGAACGAGCCATGAATACCCGTCTCGAAGGCGGATGTGAGGTGCCAATTGGGAGCTACGCCGAGCTTATTGATTGCGAAATCTGGCTGCGTGGGCTGGTCGGCGGGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGGGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  2:770470/288‑1  cTTGAATGCCGCCTTGATGATTCACGCACTCGCGGGCTGCTTGCCGCGCTGAATCACCACGAAACTGAACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  2:1040767/288‑1  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  1:222790/289‑1  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  1:249133/289‑1  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  1:458521/289‑1  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  1:980997/289‑1  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGAAGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACGGAATAACGGa                                                                                                                                                                                                                                                                                                  >  2:538814/1‑287  ATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGAAGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCG                                                                                                                                                                                                                                                                                                 <  1:300213/289‑1         GCCGCCTTGATGATTCACGCACTCGCGAGCTTCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGACATGAATACCCGTCGCGACGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGATTATTGATGGTGAACTCTGGCTGCGTGCGCTGGTCGGCGCGCCGCACGGTTCGACGATTATGCGCTGTGAACGCCACGGTGCCCCGCAAGATGCCGAACAACTGGGGATTTCGCTGCCCGATGGTGTACTGAATTACCGCGCGaccgt                                                                                                                                                                                                                                                                                          >  2:934667/1‑284                                                                                                                                                                                        gcGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTG                                                                                                                                                                         >  1:293441/3‑227                                                                                                                                                                                        gcGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGAGGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTG                                                                                                                                                                         <  2:293441/225‑1                                                                                                                                                                                            GCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGC                                                                                                                                     >  1:412967/1‑259                                                                                                                                                                                            GCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGC                                                                                                                                     <  2:412967/259‑1                                                                                                                                                                                                                                                                                              CGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACACACCGCACCTTGTACTACCGCCGGAAGAGAAATCTCGGGTGGCCACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG     >  2:8302/1‑289                                                                                                                                                                                                                                                                                                  gccCCGACCAGCCCAAGGCGCCAGATTTCGCCATAAGCAAGCGCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA  <  2:670676/285‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   TATTGAATGCCGCCTTGATGATTCACGCACTCGCGAGCTGCTTGCCGCGCTGAATCACCACGAAACTGCACTGCGCGTTACCGCAGAACGCGCCATGAATACCCGTCTCGAAGGCGGATGTCAGGTGCCAATTGGTAGCTACGCCGAGCTTATTGATGGCGAAATCTGGCTGCGTGCGCTGGTCGGCGCGCCGGACGGTTCGCAGATTATTCGCGGTGAACGCCGCGGTGCGCCGCAAGATGCCGAACAAATGGGGATTTCGCTGGCAGAAGAGCTACTGAATAACGGCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3990161‑3989869 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cggcgcgcCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA  >  NC_000913/3989977‑3990260

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