New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3989606 =	172 (1.200)	4 (0.030)	4/318	NT	2.7%	coding (223/741 nt)	hemD	uroporphyrinogen III synthase
	?	NC_000913	3989703 =	131 (0.970)					coding (126/741 nt)	hemD	uroporphyrinogen III synthase
Rejected: Frequency below/above cutoff threshold.

GAGACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3989841‑3989606 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcagcaagTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCG  >  NC_000913/3989703‑3989928                                                                                                                                                                                                                                                                                                                                                                                                                                                                                GAGACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCA                                                                                                                                                                                                                                            >  2:503347/1‑228  AGACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCAT                                                                                                                                                                                                                                           >  2:309098/1‑228    ACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCA                                                                                                                                                                                                                                         >  2:1515496/1‑228    ACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCA                                                                                                                                                                                                                                         >  2:923207/1‑228      GCCCCGGCATGAGTATCCTTGTCGCCCGCCCGTCGCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTGTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGC                                                                                                                                                                                                                                       <  1:659771/228‑1      GCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGC                                                                                                                                                                                                                                       <  1:309625/228‑1       CCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCA                                                                                                                                                                                                                                      <  2:1653894/228‑1                                CGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGCGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCA                                                                                                                                                                                                                                            >  1:1767966/1‑197                                 GCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCAT                                                                                                                                                                                                                                           <  2:1764842/197‑1                                                                                 GGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGC                                                                                                                                                                                                                                       >  2:1879451/1‑153                                                                                  GCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCA                                                                                                                                                                                                                                      <  1:1882781/153‑1                                                                                                                              AACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGC                                                                                                                                                                                                                                       >  2:1524188/1‑108                                                                                                                               ACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCA                                                                                                                                                                                                                                      <  1:1526911/108‑1                                                                                                                                           ACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCAT                                                                                                                                                                                                                                           >  1:1771867/1‑91                                                                                                                                            CTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATC                                                                                                                                                                                                                                          <  2:1768732/91‑1                                                                                                                                                                                       TTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTC                                                                                                                                             >  1:1866439/1‑141                                                                                                                                                                                        TGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCT                                                                                                                                            <  2:1863143/141‑1                                                                                                                                                                                               TCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGC                                                                                                                                                                                                                                       >  1:209828/1‑43                                                                                                                                                                                                CGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCA                                                                                                                                                                                                                                      <  2:209467/43‑1                                                                                                                                                                                                               CTTTTGCCCAATCACAGCTGCATCAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAA                                                                                                                                                                                           >  2:615453/1‑71                                                                                                                                                                                                                TTTTGCCCAATCACAGCTGCATCAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAA                                                                                                                                                                                          <  1:616461/71‑1                                                                                                                                                                                                                                    ATCAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATCCTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCG            <  1:1003064/225‑1                                                                                                                                                                                                                                    ATCAGCAAGTTGCGGTAATTGTTGACCCGCACAAACCTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCA         <  2:1010769/228‑1                                                                                                                                                                                                                                      CAGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTC             >  1:839146/1‑222                                                                                                                                                                                                                                       AGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCT      >  2:139559/1‑228                                                                                                                                                                                                                                       AGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCT      >  2:489846/1‑228                                                                                                                                                                                                                                       AGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCG            <  2:837688/222‑1                                                                                                                                                                                                                                       AGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGGGACAAGGCTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGGGCCGTTATTCAGTAGCTCTTCTGCCCGCGAAATCCCCATTTGTTCGGCATCT      >  1:1848937/1‑228                                                                                                                                                                                                                                       AGCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACG                                                                                                                                                         >  2:822496/1‑81                                                                                                                                                                                                                                        GCAAGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGG                                                                                                                                                        <  1:823918/81‑1                                                                                                                                                                                                                                           AGTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCG  <  1:121783/228‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                GAGACGCCCCGGCATGAGTATCCTTGTCACCCGCCCGTCTCCCGCTGGAGAAGAGTTAGTGAGCCGTCTGCGCACACTGGGGCAGGTGGCCTGGCATTTTCCGCTGATTGAGTTTTCTCCGGGTCAACAATTACCGCAACTTGCTGATCAACTGGCAGCGCTGGGGGAGAGCGATCTGTTGTTTGCCCTCTCGCAACACGCGGTTGCTTTTGCCCAATCACAGCTGCATCAGCAAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3989841‑3989606 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcagcaagTTGCGGTAATTGTTGACCCGGAGAAAACTCAATCAGCGGAAAATGCCAGGCCACCTGCCCCAGTGTGCGCAGACGGCTCACTAACTCTTCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCG  >  NC_000913/3989703‑3989928

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