New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 863375	32 (0.730)	4 (0.100)	4/324	NT	11.4%	coding (509/1017 nt)	hemH	Ferrochelatase
	?	AE015451	= 863376	32 (0.770)					coding (510/1017 nt)	hemH	Ferrochelatase

CCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTT  >  AE015451/863166‑863577                                                                                                                                                                                                                   |                                                                                                                                                                                                          ccAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAg                                                                                                                                                                                   >  1:721226/1‑235 (MQ=255)                tGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgc                                                                                                                                                                                                         >  2:635474‑M2/1‑196 (MQ=255)                  aTGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTg                                                                                                                                                                 <  1:2470/236‑1 (MQ=255)                      gtTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                                                                                              <  2:491009/235‑1 (MQ=255)                       ttGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgccagc                                                                                                                                                                                                     >  2:728931‑M2/1‑189 (MQ=255)                        tGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                                                                                              <  1:795921/233‑1 (MQ=255)                                  ggCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgccagcgcctcgatgtacgccgggtgctcatagaacggctgcagcacgcgcatctgca                                                                                                                                                 <  1:120950‑M2/236‑59 (MQ=255)                                      gccCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                                                                                             >  1:152141/1‑236 (MQ=255)                                                                         ttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgccagcgcctcgatgtacgc                                                                                                                                                                                       >  2:357893‑M2/1‑139 (MQ=255)                                                                           aTCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgcca                                                                                                                                                                                                       >  2:518742‑M2/1‑137 (MQ=255)                                                                                               aCGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcgc                                                                                    <  1:514702/236‑1 (MQ=255)                                                                                                cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCa                                                                                                >  2:446049/1‑223 (MQ=255)                                                                                                       cccTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                                                                                      <  2:34933/128‑1 (MQ=255)                                                                                                                  aGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAggggg                                                                                                       <  1:118814/198‑1 (MQ=255)                                                                                                                  aGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAggggg                                                                                                       >  2:118812/1‑198 (MQ=255)                                                                                                                    cAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTg                                                                <  1:694743/235‑1 (MQ=255)                                                                                                                      ggCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGc                                                                                                                                  <  2:741674/167‑1 (MQ=255)                                                                                                                       gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACgg                                                                  >  2:34154/1‑230 (MQ=255)                                                                                                                          aaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGa                                                                                             <  2:49491/200‑1 (MQ=255)                                                                                                                          aaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc                                                         <  1:562761/236‑1 (MQ=255)                                                                                                                              aGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAgccagcgcctcgatgtacgccgggtgctcatagaacggctgcagcacg                                                                                                                                                           >  2:78188‑M2/1‑86 (MQ=255)                                                                                                                                    tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcgcc                                                                                   >  1:120703/1‑199 (MQ=255)                                                                                                                                            cACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc                                                            <  1:549579/215‑1 (MQ=255)                                                                                                                                            cACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGt                                          <  2:551309/233‑1 (MQ=255)                                                                                                                                                 gTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaagt                                                   >  2:607832/1‑219 (MQ=255)                                                                                                                                                    aCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc                                                            <  1:783180/207‑1 (MQ=255)                                                                                                                                                    aCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATcgcg                               <  1:446051/236‑1 (MQ=255)                                                                                                                                                       gTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgcc                                                        >  2:681204/1‑208 (MQ=255)                                                                                                                                                              aTGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg                                                                                                                                                                                     <  2:348985/76‑1 (MQ=255)                                                                                                                                                              aTGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg                                                                                                                                                                                     >  1:348987/1‑76 (MQ=255)                                                                                                                                                                 cgcgTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaagtaagt                                               <  2:699379/207‑1 (MQ=255)                                                                                                                                                                      gctgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACgg             >  1:128155/1‑236 (MQ=255)                                                                                                                                                                        tgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaa                                                     >  2:620617/1‑194 (MQ=255)                                                                                                                                                                             ccGTTCTAGTAGCACCCGGCGTACATCGAGTCGCTGGCAGCAAGCGCCCGGCCGCACCCGGAGCAGGGTTACGACCACCTGCTACTGACCTTCCATGGTTTGCTGGAACGGCACCTGAAGAAACTg                                                                                                                    <  2:783089/126‑1 (MQ=255)                                                                                                                                                                                 tCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCttt  <  1:607834/236‑1 (MQ=255)                                                                                                                                                                                          aCCCGGCGTACATCGAGGCGCTGGCAgccagcgcctcgatgtacgccgggtgctcatagaacggctgcagcacgcgcatctgcaacggtaactggtgctcgctgaccgtctgttttgcctgctcgatcac                                                                                                   >  2:120948‑M2/1‑26 (MQ=255)                                                                                                                                                                                             cGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTgg                                                                                                                                                                                     <  2:293143/45‑1 (MQ=255)                                                                                                                                                                                                  tACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGa                                                                                                             >  1:152550/1‑112 (MQ=255)                                                                                                                                                                                                     aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                                                                          <  2:152548/96‑1 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                          CCAGCCGGCCTTGCCTGATGTGTTGGCGCGGCTGGCCGCCCAAGGTGTGCGCAAGGTGACTTTGGCACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTAAGTTCGGTTTGCTATCGCGGGCAATGCCTGGCTACGGCCAAAGCCTTT  >  AE015451/863166‑863577

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