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

GACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGG  >  AE015451/862975‑863398                                                                                                                                                                                                                           |                                                                                                                                                                                                              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                                                                                                                                                                      <  2:596898/132‑1 (MQ=255)                                                                                                                                 ccTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCAcc                                                                                                                                                                      >  1:596899/1‑132 (MQ=255)                                                                                                                                           gcgATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCa                                                                               >  2:565206/1‑209 (MQ=255)                                                                                                                                               tGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGAGGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCg                                                                                                                          <  2:293341/162‑1 (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)                                                                                                                                                                                                     cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         <  1:374093/141‑1 (MQ=255)                                                                                                                                                                                                     cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         >  2:374092/1‑141 (MQ=255)                                                                                                                                                                                                     cGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCa                                                 >  1:306575/1‑181 (MQ=255)                                                                                                                                                                                                           tGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGt                                                                                         <  2:600841/135‑1 (MQ=255)                                                                                                                                                                                                                     gtTGGCGCGGCTGGCCGCCCAAGGCGCGCGCAAGGTGCCTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGCTCGAGCAGGCac                                                                                                                 >  2:120246/1‑100 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                              GACCTGCCGTGGCCGGTGCGGCGTTTGCTGGTGTCGCTGGTCCTGATCAAGCGCCCGGAGCAGTCGGCACATGCCTACGCTTCGATCTGGTGGGAGGAGGGCTCGCCTCTGGTAGTGCTGACCCGTCGCCTGCAGGCCGCGATGGCCGAGCACTGGCCCCATGGCCCGGTGGAAATCGCCATGCGCTACGGCCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGG  >  AE015451/862975‑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.