New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 581902	37 (0.720)	3 (0.060)	3/250	NT	7.2%	coding (1192/1428 nt)	selA	selenocysteine synthase
	?	NC_002947	= 581933	42 (0.850)					coding (1223/1428 nt)	selA	selenocysteine synthase

CTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACG  >  NC_002947/581755‑582049                                                                                                                                                    |                                                                                                                                                     cTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAg                                                                                                                                                     <  6:124932/149‑1 (MQ=255) cTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAg                                                                                                                                                     >  5:202404/1‑149 (MQ=255)   ggCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTg                                                                                                                                                   <  1:237266/149‑1 (MQ=255)    gCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTgg                                                                                                                                                  <  4:152152/149‑1 (MQ=255)                     cctgcgcctgcTGACCCGCAGCCCAGCCGAGATCCAGGCCAAAGCCTAGCGACTGGCGACTGAGATCCAGAAGCGTCTCAGAGATCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTgcgc                                                                                                                                 >  8:49770/1‑149 (MQ=255)                     cctgcgcctgcTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTgcgc                                                                                                                                 >  3:106781/1‑149 (MQ=255)                               ctgACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAg                                                                                                                       >  3:140555/1‑149 (MQ=255)                                tGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGt                                                                                                                      <  6:321661‑M2/149‑2 (MQ=255)                                  aCCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGt                                                                                                                      >  5:303933‑M2/1‑146 (MQ=255)                                        aGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCActgggcaagcgcgcaaccggctggctgccactgccgatca                                                                                                              <  1:282456‑M1/149‑41 (MQ=255)                                          ccAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCtgt                                                                                                            <  6:202404/149‑1 (MQ=255)                                            aGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGc                                                                                                          <  4:140555/149‑1 (MQ=255)                                              gCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTg                                                                                                                                                   >  4:312006/1‑106 (MQ=255)                                              gCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTg                                                                                                                                                   <  3:312006/106‑1 (MQ=255)                                                  agaTCAAGGCCCAGGCCTAGCGCCTGGCGCCTGAGCTCCAGGCTCGTCTCTGTGAGCAATGGGCAATCAGCGTAGAGCCCTCGTTTGGCATGATCGGCAGTGGCAGCAAGCCGTTTGCGCTCTTGCCCAGTGCCGCTCTGTGCCTgcgc                                                                                                    <  7:49770/149‑1 (MQ=255)                                                      ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCtt                                                                                                                             <  7:70515/120‑1 (MQ=255)                                                      ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCtt                                                                                                                             >  8:70515/1‑120 (MQ=255)                                                      ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGctc                                                                                                               >  7:166581/1‑134 (MQ=255)                                                      ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGctc                                                                                                               <  8:166581/134‑1 (MQ=255)                                                      ccaggcccaggcCGAGCGCCTGGCGCAGGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGGTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCg                                                                                                 >  1:183164/1‑148 (MQ=255)                                                        aggcccaggccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCaa                                                                                              <  4:106781/149‑1 (MQ=255)                                                          gcccaggccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCActgggcaagcgcgcaaccggctggctgccactgccgat                                                                                                                >  8:128830‑M1/1‑91 (MQ=255)                                                          gcccaggccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCActgggcaagcgcgcaaccggctggctgccactgccgat                                                                                                                <  7:128830‑M1/129‑39 (MQ=255)                                                              aggccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTc                                                                                         >  2:182477/1‑148 (MQ=255)                                                                gccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCg                                                                                        >  2:282286/1‑147 (MQ=255)                                                                       gcctggcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTg                                                                               >  2:68170/1‑149 (MQ=255)                                                                         ctggcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCgaagaa                                                                                   >  7:95914/1‑143 (MQ=255)                                                                         ctggcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCgaagaa                                                                                   <  8:95914/143‑1 (MQ=255)                                                                         ctggcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTgcg                                                                             >  5:159223/1‑149 (MQ=255)                                                                          tggcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTgcgc                                                                            >  1:290836/1‑149 (MQ=255)                                                                            gcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGTGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTgcgc                                                                            <  1:182477/147‑1 (MQ=255)                                                                            gcgcctgAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGcgcg                                                                           <  2:183164/148‑1 (MQ=255)                                                                                  gAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGc                                                                                                          <  6:296381/111‑1 (MQ=255)                                                                                     cTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCt                                                                 >  2:328273/1‑149 (MQ=255)                                                                                       ccAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTg                                                                >  8:147958/1‑148 (MQ=255)                                                                                        cAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGACGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCa                                                              >  6:152999/1‑149 (MQ=255)                                                                                            cgcgTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCAtgtg                                                          <  8:155842/149‑1 (MQ=255)                                                                                               gTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTg                                                       <  5:321661/149‑1 (MQ=255)                                                                                                          cAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCt                                           >  1:9949/1‑149 (MQ=255)                                                                                                            aTGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTgc                                         <  5:61907/149‑1 (MQ=255)                                                                                                                gCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGa                                     >  7:290061/1‑149 (MQ=255)                                                                                                                gCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCAATGTGTTCCAGCGGCCAGAAATGTCGAAGAAGTTTCCCCGGGCCAGCCTTCATGTGTTGGAGCC‑CGCCCAGCGCGa                                     >  1:267896/1‑149 (MQ=255)                                                                                                                      aGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCtgc                                <  4:136126/148‑1 (MQ=255)                                                                                                                              gCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCtgcgggtgccgg                       <  6:159223/149‑1 (MQ=255)                                                                                                                                ccGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCTAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGccggtgccggtg                     >  6:243650/1‑149 (MQ=255)                                                                                                                                          cATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCa           >  1:115293/1‑149 (MQ=255)                                                                                                                                            tGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGtgccgatcatgcccaacgcgggctctacgctgattgcccattgctcaccgagacgcgcctggagctcaggcgccaggcgctcggcctgggcctggatctcggcctggct          <  1:282456‑M2/149‑110 (MQ=255)                                                                                                                                             gATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCAtt         >  6:265564/1‑148 (MQ=255)                                                                                                                                              aTCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGtgccgatcatgcccaacgcgggctctacgctgattgcccattgctcaccgagacgcgcctggagctcaggcgccaggcgctcggcctgggc                            <  7:128830‑M2/129‑92 (MQ=255)                                                                                                                                              aTCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGtgccgatcatgcccaacgcgggctctacgctgattgcccattgctcaccgagacgcgcctggagctcaggcgccaggcgctcggcctgggc                            >  8:128830‑M2/1‑38 (MQ=255)                                                                                                                                              aTCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTga       <  1:282286/149‑1 (MQ=255)                                                                                                                                                  gCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCGTCGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTgacga    <  1:68170/149‑1 (MQ=255)                                                                                                                                                   cAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGc               <  7:195746/136‑1 (MQ=255)                                                                                                                                                   cAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGc               >  8:195746/1‑136 (MQ=255)                                                                                                                                                    actgggcaagcgcgcaaccggctggctgccactgccgatcatgcccaacgcgggctctacgctgattgcccattgctcaccgagacgcgcctggagctcaggcgccaggcgctcggcctgggcctggatctcggcctggctgcgggtca  <  6:321661‑M1/149‑149 (MQ=255)                                                                                                                                                    actgggcaagcgcgcaaccggctggctgccactgccgatcatgcccaacgcgggctctacgctgattgcccattgctcaccgagacgcgcctggagctcaggcgccaggcgctcggcctgggcctggatctcggcctggctgcgggt    >  5:303933‑M1/1‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACG  >  NC_002947/581755‑582049

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.