New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1084367 =	45 (0.550)	3 (0.040)	3/262	NT	6.3%	coding (992/1272 nt)	efeB	heme‑containing peroxidase/deferrochelatase
	?	NC_000913	1084406 =	46 (0.590)					coding (1031/1272 nt)	efeB	heme‑containing peroxidase/deferrochelatase

GGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCC  >  NC_000913/1084220‑1084513                                                                                                                                                    |                                                                                                                                                   cctccagcgcttcgccattgagccttttttgtactgtcaggaagcctttttccagatcgtgttggtagcagacaaacagcaaccccatatccagttgcccggagttggtgacgcccagtgaataactgtagccacgacgcagcatcaGGc                                                                                                                                                  >  2:362028‑M1/148‑150 (MQ=255)    ccagcgcttcgccattgagccttttttgtactgtcaggaagcctttttccagatcgtgttggtagcagacaaacagcaaccccatatccagttgcccggagttggtgacgcccagtgaataactctagccacgacgcagcatcaGGCTgg                                                                                                                                               <  1:496321‑M1/6‑1 (MQ=255)      aGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCg                                                                                                                                             <  1:638533/150‑1 (MQ=255)                tGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCa                                                                                                                                   >  2:496321/1‑150 (MQ=255)                tGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCa                                                                                                                                   <  1:1396811/150‑1 (MQ=255)                 gAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCAc                                                                                                                                  <  1:391726/150‑1 (MQ=255)                          gcagACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGagag                                                                                                                         >  2:126772/1‑150 (MQ=255)                          gcagACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGagag                                                                                                                         >  2:308968/1‑150 (MQ=255)                            agACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAgagtg                                                                                                                       >  2:1118776/1‑150 (MQ=255)                            agACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAgagtg                                                                                                                       <  2:880508/150‑1 (MQ=255)                            agACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAgagtg                                                                                                                       <  1:1165979/150‑1 (MQ=255)                            agACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAgagtg                                                                                                                       <  2:905460/150‑1 (MQ=255)                                 aTTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTcc                                                                                                                  <  2:478450/150‑1 (MQ=255)                                     ttGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGcc                                                                                                              <  1:1118776/150‑1 (MQ=255)                                         gaagcctttttccagatcgtgttggtagcagacaaacagcaaccccatatccagttgcccggagttggtgacgcccagtgaataactgtagccacgacgcagcatcaGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAt                                                                                                          >  1:401968‑M1/108‑150 (MQ=255)                                             gATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAtgctg                                                                                                      >  2:519073/1‑150 (MQ=255)                                                aaGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGcgt                                                                                                   >  2:1205736/1‑150 (MQ=255)                                                  gCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGcgtcg                                                                                                 >  2:1022489/1‑150 (MQ=255)                                                   cAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGcgtcgt                                                                                                >  2:11311/1‑150 (MQ=255)                                                    aaaCCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTg                                                                                               <  1:899135/150‑1 (MQ=255)                                                    aaaCCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTg                                                                                               <  2:330656/150‑1 (MQ=255)                                                       ccGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCt                                                                                            <  1:1185140/150‑1 (MQ=255)                                                       ccGGTGCGCCGCGGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCt                                                                                            <  1:710111/150‑1 (MQ=255)                                                             gttggtagcagacaaacagcaaccccatatccagttgcccggagttggtgacgcccagtgaataactgtagccacgacgcagcatcaGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAtgctg                                                                                                      <  2:840700‑M1/47‑1 (MQ=255)                                                             gttggtagcagacaaacagcaaccccatatccagttgcccggagttggtgacgcccagtgaataactgtagccacgacgcagcatcaGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAtgctg                                                                                                      >  1:840700‑M1/88‑134 (MQ=255)                                                              gccgcTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCAGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTa                                                                                     >  1:300726/1‑150 (MQ=255)                                                                   tGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGTCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCAc                                                                                <  2:946957/150‑1 (MQ=255)                                                                    gggAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACt                                                                               >  1:967856/1‑150 (MQ=255)                                                                     ggAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTg                                                                              >  1:577725/1‑150 (MQ=255)                                                                         tGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCg                                                                          <  1:37475/150‑1 (MQ=255)                                                                          gcagcaTGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGt                                                                         >  2:188989/1‑150 (MQ=255)                                                                              catgagcatgaTGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGTATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCAcc                                                                     >  1:67773/1‑150 (MQ=255)                                                                              catgagcatgaTGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCAcc                                                                     <  1:126772/150‑1 (MQ=255)                                                                                      tgatgTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCgg                                                             >  1:902742/1‑150 (MQ=255)                                                                                       gatgTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCCGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCAAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCggg                                                            >  2:1072453/1‑150 (MQ=255)                                                                                                    aCGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATgg                                               >  1:955071/1‑150 (MQ=255)                                                                                                              cccGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGtt                                     >  1:427725/1‑150 (MQ=255)                                                                                                                 ggaagggaagGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGt                                  >  1:1241059/1‑150 (MQ=255)                                                                                                                   aagggaaggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCt                                <  1:519073/150‑1 (MQ=255)                                                                                                                    agggaaggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGGTATGGGGTTGCTGTTTGTCTg                               <  2:109461/150‑1 (MQ=255)                                                                                                                     gggaaggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGc                              <  1:1205736/150‑1 (MQ=255)                                                                                                                      ggaaggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCt                             <  1:11311/150‑1 (MQ=255)                                                                                                                         aggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTAcc                          <  1:40273/150‑1 (MQ=255)                                                                                                                         aggTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTAcc                          >  2:63243/1‑150 (MQ=255)                                                                                                                             gATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAaca                      <  1:188989/150‑1 (MQ=255)                                                                                                                                         cAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGaaaa          >  2:771526/1‑150 (MQ=255)                                                                                                                                          aGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGaaaaa         <  1:1038381/150‑1 (MQ=255)                                                                                                                                           gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAtgctg                                                                                                      >  2:992204/1‑56 (MQ=255)                                                                                                                                           gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGAtgctg                                                                                                      <  1:992204/56‑1 (MQ=255)                                                                                                                                           gCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAg        >  1:1063844/1‑150 (MQ=255)                                                                                                                                                atCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCCGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTc   >  2:917806/1‑150 (MQ=255)                                                                                                                                                 tCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc  >  2:420832/1‑150 (MQ=255)                                                                                                                                                 tCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTcc  <  1:1407676/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GGGACAGAACGCCGCTGAAAGAACAGCAGACGATTTTTGGCCGTGATAAGCAAACCGGTGCGCCGCTGGGAATGCAGCATGAGCATGATGTGCCTGATTACGCCAGCGACCCGGAAGGGAAGGTGATCGCGCTGGACAGCCATATCCGGCTGGCGAATCCCCGCACGGCGGAGAGTGAGTCCAGCCTGATGCTGCGTCGTGGCTACAGTTATTCACTGGGCGTCACCAACTCCGGGCAACTGGATATGGGGTTGCTGTTTGTCTGCTACCAACACGATCTGGAAAAAGGCTTCC  >  NC_000913/1084220‑1084513

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