New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 582276	37 (1.310)	4 (0.150)	4/232	NT	10.4%	coding (139/1923 nt)	selB	selenocysteyl‑tRNA‑specific translation elongation factor
	?	NC_002947	= 582287	34 (1.280)					coding (150/1923 nt)	selB	selenocysteyl‑tRNA‑specific translation elongation factor

CCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTGGCGGCCGACGACGGGGTCATGCCGCAGACCCGT  >  NC_002947/582130‑582419                                                                                                                                                   |                                                                                                                                                ccGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGa                                                                                                                                              <  2:389398/150‑1 (MQ=255)   ggTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGc                                                                                                                                            <  2:531247/150‑1 (MQ=255)       cAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAg                                                                                                                                        <  2:340276/150‑1 (MQ=255)                ttGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCgcgccgc                                                                                                                               <  1:510573‑M2/150‑8 (MQ=255)                        gCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGc                                                                                                                       <  1:180393/150‑1 (MQ=255)                         cTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCt                                                                                                                      >  2:162253/1‑150 (MQ=255)                          tGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCtt                                                                                                                     <  2:516210/150‑1 (MQ=255)                                     aCCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGc                                                                                                          <  2:647366/150‑1 (MQ=255)                                      ccACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGcc                                                                                                         >  2:255772/1‑150 (MQ=255)                                       cACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAGGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCg                                                                                                        >  2:486350/1‑150 (MQ=255)                                             aaGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCgg                                                                                                       <  2:704917/145‑1 (MQ=255)                                             aaGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCgg                                                                                                       >  1:704917/1‑145 (MQ=255)                                             aaGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCAc                                                                                                  <  2:319741/150‑1 (MQ=255)                                                      ctgctgctgGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCttt                                                                                         <  2:68890/150‑1 (MQ=255)                                                       tgctgctgGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTa                                                                                        >  2:171590/1‑150 (MQ=255)                                                                 gcTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACAt                                                                              <  1:99619/150‑1 (MQ=255)                                                                   tGACCGGCCAGGCCGGCGACCGGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGaa                                                                                                                                         >  2:627361/1‑89 (MQ=255)                                                                   tGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGaa                                                                                                                                         <  1:627361/89‑1 (MQ=255)                                                                      ccGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTgg                                                                         <  1:22298/150‑1 (MQ=255)                                                                       cGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCgccttcgctcagcgccgcatagcggtagccaaggtcga                                                                                                           <  2:453989‑M1/115‑39 (MQ=255)                                                                       cGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCgccttcgctcagcgccgcatagcggtagccaaggtcga                                                                                                           >  1:453989‑M1/1‑77 (MQ=255)                                                                       cGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGc                                                                        >  1:187145/1‑150 (MQ=255)                                                                           cAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCt                                                                           >  1:332542/1‑143 (MQ=255)                                                                           cAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCt                                                                           <  2:332542/143‑1 (MQ=255)                                                                            aGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCgccttcgctcagcgccgca                                                                                                                              >  1:410652‑M1/1‑72 (MQ=255)                                                                            aGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCgccttcgctcagcgccgca                                                                                                                              <  2:410652‑M1/91‑20 (MQ=255)                                                                               ccGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGccc                                                                <  2:510573/150‑1 (MQ=255)                                                                                cGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCa                                                               >  2:492572/1‑150 (MQ=255)                                                                                cGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCACGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGAGCGGCTTCATCGAAAAGCTTGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCa                                                               <  2:401932/150‑1 (MQ=255)                                                                                  gCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGc                                                                                                                       <  2:518097/92‑1 (MQ=255)                                                                                  gCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGc                                                                                                                       >  1:518097/1‑92 (MQ=255)                                                                                   cGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATgg                                                            <  1:255772/150‑1 (MQ=255)                                                                                          cGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGAc                                                     <  1:436523/150‑1 (MQ=255)                                                                                                  agagCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGTCCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCt                                             <  1:720604/150‑1 (MQ=255)                                                                                                    agCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGt                                           <  1:514147/150‑1 (MQ=255)                                                                                                     gcgcgcCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGCGCAACGAGCGCTTTATCCACAACATGCTGGCCGGCGCTCATGGCATCGACCTGGTGCTGtt                                          <  2:428527/150‑1 (MQ=255)                                                                                                      cgcgcCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTg                                         <  1:486350/150‑1 (MQ=255)                                                                                                           ccGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGAt                                                                                                              <  2:102938/76‑1 (MQ=255)                                                                                                           ccGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGAt                                                                                                              >  1:102938/1‑76 (MQ=255)                                                                                                                         tCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCgcgccgcatagcggtagccaaggtcgatggtcatgccgcgggcgcgctcttcttgccgctggtcgccggcctggccg                                                         <  2:453989‑M2/115‑78 (MQ=255)                                                                                                                         tCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCgcgccgcatagcggtagccaaggtcgatggtcatgccgcgggcgcgctcttcttgccgctggtcgccggcctggccg                                                         >  1:453989‑M2/1‑38 (MQ=255)                                                                                                                          cGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCGGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTGGCGGCCGACGACgg                     <  2:332248/150‑1 (MQ=255)                                                                                                                            aCCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTGGCGGCCGACGACgggg                   >  2:556998/1‑150 (MQ=255)                                                                                                                                    tACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTGGCGGCCGACGACGGGGTCATGCCg           >  1:604058/1‑150 (MQ=255)                                                                                                                                         cTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTggc                                 >  1:671332/1‑123 (MQ=255)                                                                                                                                         cTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTggc                                 <  2:671332/123‑1 (MQ=255)                                                                                                                                            tGCGGCGCTGAGCGAAGGCgcgccgcatagcggtagccaaggtcgatggtcatgccgcgggcgcgctcttcttgccgctggtcgccggcct                                                              <  2:410652‑M2/91‑73 (MQ=255)                                                                                                                                            tGCGGCGCTGAGCGAAGGCgcgccgcatagcggtagccaaggtcgatggtcatgccgcgggcgcgctcttcttgccgctggtcgccggcct                                                              >  1:410652‑M2/1‑19 (MQ=255)                                                                                                                                             gcggcgCgccttcgctcagcgccgcatagcggtagccaaggtcgatggtcatgccgcgggcgcgctcttcttgccgctggtcgccggcctggccggtcagcgccagcagcagtgcggtcttgccgtggtcgatgtggccagcggtaccaa  <  1:510573‑M1/150‑144 (MQ=255)                                                                                                                                                   |                                                                                                                                                CCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCTGCTGGCGCTGACCGGCCAGGCCGGCGACCAGCGGCAAGAAGAGCGCGCCCGCGGCATGACCATCGACCTTGGCTACCGCTATGCGGCGCTGAGCGAAGGCGCGCCGCTGACCGGCTTCATCGATGTGCCGGGCCACGAGCGCTTTATCCACAACATGCTGGCCGGCGCCCATGGCATCGACCTGGTGCTGTTGGTAGTGGCGGCCGACGACGGGGTCATGCCGCAGACCCGT  >  NC_002947/582130‑582419

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