New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3498808 =	47 (1.180)	4 (0.100)	3/254	NT	9.0%	coding (981/1374 nt)	cysG	siroheme synthase
	?	NC_000913	3498902 =	36 (0.940)					coding (1075/1374 nt)	cysG	siroheme synthase

AAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGAAAACCTGGC  >  NC_000913/3498666‑3498942                                                                                                                                               |                                                                                                                                       ccaccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtcgcggcctgattcaaccccatatagaacaccagcgtctgtttttctgccgccgggatttcccagtccagctcgccaccggtttttacgtgt                                                                                                                                        <  1:743177‑M1/1‑1 (MQ=255)  caccggcatttcgcctggcattccgtggtcaatcagcgtttgctgaatagtcgcggcctgattcaacctcttatagaacaccagcgtctgtttttctgccgccaggttttcccagtccagctcgccaccggtttttaagtgTc                                                                                                                                       >  1:295053‑M1/142‑143 (MQ=255)    ccggcatttcgcctggcattccgtgttcaatcagcttttgctgaatagtcgcggcctgattcaaccccatatagaacaccagcgtctgtttttctgccgccaggttttcccagtccagctcgccaccggtttttaagtgTCCg                                                                                                                                     >  1:605688‑M1/140‑143 (MQ=255)    ccggcactttgcctcgcattccgtgttcaatcagcttttgctgaatagtcgcggcctgattcaaccccatatagaacaccagcgtctgtttttctgccgccaggttttcccagtccagctcgccaccggtttttaagtgTCCg                                                                                                                                     >  2:487293‑M1/140‑143 (MQ=255)    agaTTAACCAGGTCCTGCGGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCg                                                                                                                                     <  1:59931/143‑1 (MQ=255)      aTTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGTGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTGCggg                                                                                                                                   >  1:188073/1‑143 (MQ=255)           ccAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTCTTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTGCCggtgatga                                                                                                                              >  2:593086/1‑137 (MQ=255)               tctctgctgCGGGAAGCGGAAAAAGGCAAACGCGTGGTGCGGCGGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGc                                                                                                                          <  1:726235/140‑1 (MQ=255)                tCCGGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCa                                                                                                                         <  1:787882/143‑1 (MQ=255)                  ctgctgCGGGAAGCGCAAAAAGGCAAACGCGTGTTGCGGCTGATACCTTTGTTTTCGTTTATTGTTGGCCCCGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGc                                                                                                                       <  1:147063/143‑1 (MQ=255)                         ggctaGCGCAAAAGGGCAAACGCGTGCTGCGGCTGCAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACCCTGGGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGt                                                                                                                <  1:101206/139‑1 (MQ=255)                          gggaGCGCAGAAAGGGAAACGCGTGGGGCGGCTGAAGGGTGGCGATCCGTGTATTGTTGGCCGTGGTGGGGAAGAGCTGGAAACGCTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGtt                                                                                                               <  1:755707/140‑1 (MQ=255)                           gcctcccAAACACGCAAACGCGTCGCGCCCCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTTGAACCACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTg                                                                                                              <  1:515623/137‑1 (MQ=255)                           gAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCTATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAATCACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTg                                                                                                              <  1:337549/143‑1 (MQ=255)                             aGCGCAAAAAGGCAAACGCGTGGAACGGCTTCGTTGGGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGct                                                                                                            <  1:25524/143‑1 (MQ=255)                               cgcAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGctct                                                                                                          >  1:99624/1‑143 (MQ=255)                                gcAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTCCCGCAGCTTCTGATTGctggg                                                                                                         >  1:568580/1‑140 (MQ=255)                                 cAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGc                                                                                                        >  1:597549/1‑143 (MQ=255)                                  aaaaaGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGcc                                                                                                       <  1:322219/143‑1 (MQ=255)                                      aGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTCGCGGGTACAACCGCAGCTCATGGCTGCTCTGCCtctg                                                                                                   >  2:618255/1‑140 (MQ=255)                                       ttttgctgaatagtcgcggcctgattcaaccccatatagaacaccagcgtctgtttttctgccgccaggttttcccagtccagctcgccaccggtttttaagtgTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTc                                                                                                  <  2:366706‑M1/39‑1 (MQ=255)                                          aaaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                                                                               >  1:14705/1‑143 (MQ=255)                                          aaaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCggg                                                                                               >  2:91995/1‑143 (MQ=255)                                          aaaCGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTAGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTCCCGGGTTCTACCGCAGCTTCTGGTTGCTCTGCGTAGTCggg                                                                                               >  2:279390/1‑143 (MQ=255)                                              gcgTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTGGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTAtt                                                                                           <  1:604840/143‑1 (MQ=255)                                                   gtgCGGCGGAAAGGGGGCGATCCGTTTAGTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACt                                                                                      <  1:135099/143‑1 (MQ=255)                                                     gCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACACGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCAGGTCTCCCACTCa                                                                                    >  1:711729/1‑143 (MQ=255)                                                             aaGGCGGCGACCCGTGTATTTTGGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcgc                                                                            <  1:647893/143‑1 (MQ=255)                                                              aGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATcggg                                                                           >  1:774945/1‑141 (MQ=255)                                                                  ggCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTa                                                                       <  1:91995/143‑1 (MQ=255)                                                                     gATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCTGGGTATTACCGCAACGTGTCGATGCTCTGACTATTCAGGTAGCTCACTCCCGCAACGCGATTATgg                                                                    >  2:310012/1‑142 (MQ=255)                                                                         cGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCCGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAGTATGccgcg                                                                >  2:331201/1‑140 (MQ=255)                                                                               tttttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGAGTATGCCCAGAGCGTa                                                          >  1:94511/1‑143 (MQ=255)                                                                                 tttGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACg                                                        <  1:449954/143‑1 (MQ=255)                                                                                   tGGCCGGGGGGGCGAAGAGCTGCAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCt                                                      <  1:378269/143‑1 (MQ=255)                                                                                   cagcgtctgtttttctgccgccaggttttcccagtccagctcgccaccggtttttaagtgTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCt                                                      >  1:366706‑M1/61‑143 (MQ=255)                                                                                      ccGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTCGCGGGGATAACCGAAGCATCTGGTTGCTCTGCCCATTCGGGTATTCCACTAACGCACCAAGCATATGCGCCGAGAGTCCGCttga                                                   >  2:576495/1‑141 (MQ=255)                                                                                        gtggtgGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGCATTAGCGGTGCCCCACCCCGCTCTGCCTATTCGGGTATTCCCCTCACGCATCCCGCTTATGCCCAGAGCGCACGCTTAAtt                                                 >  1:784864/1‑143 (MQ=255)                                                                                           gtggCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTAcc                                              <  1:600028/143‑1 (MQ=255)                                                                                            gggCGAAGAGCTGGAAACGCTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGTCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCg                                             <  1:615233/142‑1 (MQ=255)                                                                                               cGAAGAGCCGCAAACACTGCGCACCGCGGGTATCCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGac                                          <  1:136050/143‑1 (MQ=255)                                                                                                gAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTGCGGGAACTATAACAAATACTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTCTGCCCAGAGCGTACGCTTAATTACCGGcca                                         >  2:8549/1‑143 (MQ=255)                                                                                                 aaGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGCATTACCGCAGATTGGGGGTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGacac                                        <  1:364025/143‑1 (MQ=255)                                                                                                  agagCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTCGCGGGTACAACCGAAACATCTGCTTGCTCTGCCTCTTCGGGTATTCCACTCACGCATCGCGATTACGCCCAGAGCGTACGCCTAATTACCGGACACt                                       >  2:140172/1‑143 (MQ=255)                                                                                                          aaaCACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAAcc                               <  1:512870/143‑1 (MQ=255)                                                                                                                tgtgCAACGCGGGTATTCCGTTCTCGGTGGTCGCCGGGACATCCACAACAAAAGGTTGCTCCGCCTATTCGAGTCTTCCACTCACGAATCACGATTATGACAACAGCGGACACCTAATTACCGGGCACTTAAAAACCGCTGGc                         >  2:434254/1‑143 (MQ=255)                                                                                                                   gCAACGCGGGTCTTCCGTTCTCGGTGGTTGCGGGAACAACCGCAGCAACTGGTTGCTCTGCCTATTCGGGCCTTCCCCTCACGCATCACGCTCATGCCCAGCGCGTCCGCTTACTTACCGGACACGTCAAAACCGGTAGccgg                      >  2:353948/1‑140 (MQ=255)                                                                                                                     aaCGCGGGTAGTCCGTTGTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACCCACGCAGCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCt                    <  2:142860/143‑1 (MQ=255)                                                                                                                        gcgGGTATTCCGTTCTCGGTGGTTCCGGGTACTCCCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGa                 <  1:212263/143‑1 (MQ=255)                                                                                                                             tagcccGTTCTGGGTCGATGGGGGGAGTACCGGAGCCCCTGGGTGCGCTGCCTATTCGGGTATTCCACTCGCGCGGCGCGATGGTGCCGAGAGCGTACGCTTAAGTACCGGACACGAAAAGACCGGTGGCGAGCGGGACTggg            <  2:605542/139‑1 (MQ=255)                                                                                                                               ttCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCGTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGaa          <  1:310313/143‑1 (MQ=255)                                                                                                                                 ccgggcgcGGTGGTGCCGGGTAGGACCGCAGCGTCTGGTTGCTCTGCCTATTCGGGTATGCCACTCACGCGTCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGAGACTTAAAAACCGGTGGCGAGCTGGACTGGGaaaa        <  1:197599/136‑1 (MQ=255)                                                                                                                                 ccGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGaaaa        >  1:717969/1‑143 (MQ=255)                                                                                                                                      ctcGGTGGTGGCGCGTATTGCCGCAGCTGCTGGGTGCCCTGCCTGGGCGGGGCTTCCACGCACGCACCGCGTGTAGGCCCAGGGCGTACGCTTAATTACCGGACACTGAAAAACCGGTGGCGAGCTGGACTGGGAAAACCTgg   <  2:611282/143‑1 (MQ=255)                                                                                                                                       tcGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGCTTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGAAAACCTggc  >  1:554880/1‑143 (MQ=255)                                                                                                                                         taagtgTCCGGGTCTTACCGCATCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGAAAACCtga   >  1:324018‑M1/7‑139 (MQ=255)                                                                                                                                         taagtgTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGAAAACCtga   <  2:324018‑M1/134‑2 (MQ=255)                                                                                                                                               |                                                                                                                                       AAGAGATTAACCAGATCCTGCTGCGGGAAGCGCAAAAAGGCAAACGCGTGGTGCGGCTGAAAGGTGGCGATCCGTTTATTTTTGGCCGTGGTGGCGAAGAGCTGGAAACACTGTGCAACGCGGGTATTCCGTTCTCGGTGGTTCCGGGTATTACCGCAGCTTCTGGTTGCTCTGCCTATTCGGGTATTCCACTCACGCATCGCGATTATGCCCAGAGCGTACGCTTAATTACCGGACACTTAAAAACCGGTGGCGAGCTGGACTGGGAAAACCTGGC  >  NC_000913/3498666‑3498942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ 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.