New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1759424 =	30 (0.940)	6 (0.200)	5/242	2.9	17.6%	coding (1100/1221 nt)	sufS	cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
	?	NC_000913	1759462 =	28 (0.930)					coding (1062/1221 nt)	sufS	cysteine desulfurase, stimulated by SufE; selenocysteine lyase, PLP‑dependent
Rejected: Frequency below/above cutoff threshold.

AACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  >  NC_000913/1759314‑1759606                                                                                                                                                     |                                                                                                                                                 aatctgtaccggatctcactctctatggcccacaaaacaggctttgcgttattgcttttaatctcgggaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACTCACAGCAATGCCg                                                                                                                                                  >  6:106470‑M1/111‑149 (MQ=255)    ggTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCATaa                                                                                                                                               <  3:3/149‑1 (MQ=255)         ccggatctcactctctatggcccacaaaacaggcttggcgttattgcttttaatctcggtaaacaccacgcctatgatgttggcagttttctcgataattacGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATc                                                                                                                                          >  4:73841‑M1/103‑149 (MQ=255)         aTACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGCTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATc                                                                                                                                          <  5:81668/149‑1 (MQ=255)                              aGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAAcatca                                                                                                                       >  7:16765/1‑147 (MQ=255)                                                 gggAGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGAATTGCGCAGTGATGTCCGGTACGCACAGAAATGCAGTAATTATCGAGAAAAATGCCAACATCATATGCGTGGTGTTTAcc                                                                                                      >  8:24660/1‑145 (MQ=255)                                                      ttATACATGGCCAGCGACGCGCGACACTTCGCAGGGACGTTGTAATAGGATATCAATGGCATTTCGCAGTTATGTCCGGTACGCACAGCAATGGCGTAATTATCGAGAAAACTTCCAACATCATAGGCGTTGTGTTTACCGAGATTaa                                                                                              <  7:24660/148‑1 (MQ=255)                                                          aCATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGc                                                                                          >  7:6398/1‑148 (MQ=255)                                                           cATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTTCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCaa                                                                                        >  5:54744/1‑149 (MQ=255)                                                           cATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCaa                                                                                        <  8:16765/149‑1 (MQ=255)                                                              ggCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATa                                                                                                                     <  6:106467/117‑1 (MQ=255)                                                                ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCa                                                                                                                             <  7:95441/107‑1 (MQ=255)                                                                ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCa                                                                                                                             >  8:95441/1‑107 (MQ=255)                                                                ccAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACg                                                                                   >  3:84256/1‑149 (MQ=255)                                                                 cAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGc                                                                                  >  2:61903/1‑149 (MQ=255)                                                                         cccGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTTCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTg                                                                          <  6:54744/149‑1 (MQ=255)                                                                          ccGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGt                                                                         >  1:74661/1‑149 (MQ=255)                                                                                     cAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCa                                                              >  5:44321/1‑149 (MQ=255)                                                                                      aGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCa                                                              >  5:100202/1‑148 (MQ=255)                                                                                            gTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATagagag                                                       >  7:85182/1‑149 (MQ=255)                                                                                                tAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGt                                                      >  3:78787/1‑146 (MQ=255)                                                                                                tAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTgag                                                   <  1:61903/149‑1 (MQ=255)                                                                                                 aaTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAg                                                                                           <  1:35472/108‑1 (MQ=255)                                                                                                 aaTAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAg                                                                                           >  2:35472/1‑108 (MQ=255)                                                                                                      gataattacggcattgctgtgcgtaccggacatcactgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg                                             >  4:73841‑M2/48‑149 (MQ=255)                                                                                                      gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg                                             <  8:6398/149‑1 (MQ=255)                                                                                                      gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg                                             <  6:44321/149‑1 (MQ=255)                                                                                                      gCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCgg                                             <  2:74661/149‑1 (MQ=255)                                                                                                              cggcattgctgtgagtaccggacatcactgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTCCCGAGATTAAAAGCAATAACGCAAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGAtt                                     >  6:106470‑M2/40‑149 (MQ=255)                                                                                                               ggcattgctgtgcgtaccggacatcactgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGa                                       <  1:133888‑M2/108‑1 (MQ=255)                                                                                                               ggcattgctgtgcgtaccggacatcactgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGa                                       >  2:133888‑M2/39‑146 (MQ=255)                                                                                                                GCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTc                                    <  6:96352‑M1/148‑1 (MQ=255)                                                                                                                 CATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCa                                  >  4:84256‑M1/1‑149 (MQ=255)                                                                                                                                    gTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCTCATAATg               <  8:85182/149‑1 (MQ=255)                                                                                                                                        actgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAATCAATAACTCCAAGCCTGTTTTGTGGGACATAGAGAGTGATATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATc           <  5:106470‑M2/136‑1 (MQ=255)                                                                                                                                        actgcgcaatgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATc           <  3:73841‑M2/136‑1 (MQ=255)                                                                                                                                         cacaGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCa          <  4:54982/149‑1 (MQ=255)                                                                                                                                            aGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGAt       <  5:96352/149‑1 (MQ=255)                                                                                                                                                atgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTg   >  1:22401‑M2/6‑149 (MQ=255)                                                                                                                                                 tgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGt  <  6:100202‑M2/145‑1 (MQ=255)                                                                                                                                                 tgccGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGt  <  4:78787‑M2/145‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AACGGTGAATACGTTGCAGGCCGGTCACCAGACGATCCACTTCTTCATGGGTGTTATACATGGCCAGCGACGCCCGACACATCGCAGGGACGTTGTAATAGGCCATCAATGGCATTGCGCAGTGATGTCCGGTACGCACAGCAATGCCGTAATTATCGAGAAAACTGCCAACATCATAGGCGTGGTGTTTACCGAGATTAAAAGCAATAACGCCAAGCCTGTTTTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGT  >  NC_000913/1759314‑1759606

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