New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 863172	31 (0.780)	3 (0.080)	3/316	NT	8.8%	coding (306/1017 nt)	hemH	Ferrochelatase
	?	AE015451	= 863193	33 (0.880)					coding (327/1017 nt)	hemH	Ferrochelatase

TACCTCAACCAGTTCCTCATGGACCCGTACGTGATCGACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGG  >  AE015451/862939‑863398                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   tACCTCAACCAGTTCCTCATGGACCCGTACGTGCTCGACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCgg                                                                                                                                                                                                                                   <  1:520421/235‑1 (MQ=255)      cAACCAGTTCCTCATGGACCCGTACGTGATCGACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGc                                                                                                                                                                                                                             >  1:578140/1‑236 (MQ=255)                   aTGGACCCGTACGGGATCGACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCt                                                                                                                                                                                                                           <  1:597789/225‑1 (MQ=255)                                     gACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGcg                                                                                                                                                                                                             >  2:303607‑M2/1‑219 (MQ=255)                                                        ggcgTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTg                                                                                                                                                                           >  1:352240/1‑236 (MQ=255)                                                                                    cAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGAc                                                                                                                                                >  1:737969/1‑235 (MQ=255)                                                                                    cAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGAc                                                                                                                                                <  1:426227/235‑1 (MQ=255)                                                                                     aaGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACa                                                                                                                                               >  2:457419/1‑235 (MQ=255)                                                                                               aGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGAcc                                                                                                                                    <  1:722351/236‑1 (MQ=255)                                                                                                       gCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGAt                                                                                                                            <  1:491651/236‑1 (MQ=255)                                                                                                                tACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCCTGGCCCGGGGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCac                                                                                                                                                                                                >  2:46361/1‑158 (MQ=255)                                                                                                                 aCGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGcggctggccgtagcgcatggcgatttccaccgggccatggggccagtgctcggccatcgcggcctg                                                                                                                                             <  1:259722‑M2/209‑67 (MQ=255)                                                                                                                 aCGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGc                                                                                                                   >  2:518812/1‑235 (MQ=255)                                                                                                                          tCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAg                                                                                                            >  2:487075/1‑233 (MQ=255)                                                                                                                                 ggaggaGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGc                                                                                                   <  1:487076/235‑1 (MQ=255)                                                                                                                                     gaggGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGAc                                                                                                          <  1:223912/224‑1 (MQ=255)                                                                                                                                       gggCTCGCCTCTGGTAGTGCTGGCCCGGCGCCTGCAGGCCGCGATGGCAGAGCACTGGCCCCGTGGCCCGGTGGAAATCCCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGccgctgaccgtagcgcatggcgatttccaccggg                                                                                                                                                                             >  2:465020‑M2/1‑121 (MQ=255)                                                                                                                                        ggCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCa                                                                                           <  1:518813/236‑1 (MQ=255)                                                                                                                                           tCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCa                                                                                                                                                                        >  2:445358/1‑156 (MQ=255)                                                                                                                                                tctGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         <  2:402611/230‑1 (MQ=255)                                                                                                                                                      aGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGa                                                                                                                                      <  1:445359/179‑1 (MQ=255)                                                                                                                                                          cTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGAt                                                                            <  2:578139/233‑1 (MQ=255)                                                                                                                                                            gACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGGGCTATGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGTCTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATgcg                                                                         >  2:406475/1‑234 (MQ=255)                                                                                                                                                              cccGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGcggctggccgtagcgcatggcgatttccaccgggccatggggccagtgctcgg                                                                                                                                                          >  2:259721‑M2/1‑98 (MQ=255)                                                                                                                                                               ccGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCa                                                                                                                                   <  2:130836/173‑1 (MQ=255)                                                                                                                                                                cGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCCAGGTGACTTTTGCACCGcc                                                                                                                                                                   >  2:59204/1‑139 (MQ=255)                                                                                                                                                                     ccTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCAcc                                                                                                                                                                      >  1:596899/1‑132 (MQ=255)                                                                                                                                                                     ccTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCAcc                                                                                                                                                                      <  2:596898/132‑1 (MQ=255)                                                                                                                                                                         cAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGcgccaacacatcaggcaaggccggctggccgtagcgcatggcgatttccaccgggccatggggccagtgctcggccatcgcggcctgcaggcgacgggtcagcactaccagaggcgagccctcctcccaccagatcgaagcgt                                                                                     <  1:259722‑M1/209‑144 (MQ=255)                                                                                                                                                                               gcgATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCa                                                                               >  2:565206/1‑209 (MQ=255)                                                                                                                                                                                   tGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGAGGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCg                                                                                                                          <  2:293341/162‑1 (MQ=255)                                                                                                                                                                                      ccGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGcgccaacacatcaggcaaggccggctggccgtagcgcatggcgatttccaccgggccatggggccagtgctcggccatcgcggcctgcaggcgacggg                                                                                                                                  >  2:259721‑M1/1‑53 (MQ=255)                                                                                                                                                                                                         cccGGTGGAAATCGCCATGCGCTACGGTCAGCGGcgccaacacatcaggcaaggccggctggccgtagcgcatggggatttccaccgggccacggggccagtgctctgccatcgcggcctgcaggcgccgggccagcactaccagaggcgagccc                                                                                                           >  2:465020‑M1/1‑34 (MQ=255)                                                                                                                                                                                                               ggAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGc                                                                                                                                                                         <  2:40857/87‑1 (MQ=255)                                                                                                                                                                                                               ggAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGc                                                                                                                                                                         >  1:40857/1‑87 (MQ=255)                                                                                                                                                                                                                    tCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTACCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTAc                                       <  2:626834/212‑1 (MQ=255)                                                                                                                                                                                                                       ccATGCGCGACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGGGCGCAAGGTGACTTTGGCAACGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCaacggtaac                                                        <  1:52546/192‑9 (MQ=255)                                                                                                                                                                                                                           gcgcTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCg                                             <  1:373299/199‑1 (MQ=255)                                                                                                                                                                                                                                  ggCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg  <  1:406476/235‑1 (MQ=255)                                                                                                                                                                                                                                         cgcgccaacacatcaggcaaggccggctggccgtagcgcatggcgatttccaccgggccatggggccagtgctcggccatcgcggcctgcaggcgacgggtcagcactaccagaggcgagccctcctcccaccagatcgaagcgtaggcatgtgccgactgctccgggcgcttgatcaggaccagcgacaccagcaaacgccgcaccggccacggcaggtc         >  2:303607‑M1/1‑2 (MQ=255)                                                                                                                                                                                                                                         cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         >  2:374092/1‑141 (MQ=255)                                                                                                                                                                                                                                         cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         <  1:374093/141‑1 (MQ=255)                                                                                                                                                                                                                                         cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCa                                                 >  1:306575/1‑181 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                   TACCTCAACCAGTTCCTCATGGACCCGTACGTGATCGACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGG  >  AE015451/862939‑863398

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