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

CTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGA  >  NC_002947/581785‑582080                                                                                                                                                     |                                                                                                                                                     ctgACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAg                                                                                                                                                      >  3:140555/1‑149 (MQ=255)  tGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGt                                                                                                                                                     <  6:321661‑M2/149‑2 (MQ=255)    aCCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGt                                                                                                                                                     >  5:303933‑M2/1‑146 (MQ=255)            ccAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCtgt                                                                                                                                           <  6:202404/149‑1 (MQ=255)              aGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGc                                                                                                                                         <  4:140555/149‑1 (MQ=255)                    agaTCAAGGCCCAGGCCTAGCGCCTGGCGCCTGAGCTCCAGGCTCGTCTCTGTGAGCAATGGGCAATCAGCGTAGAGCCCTCGTTTGGCATGATCGGCAGTGGCAGCAAGCCGTTTGCGCTCTTGCCCAGTGCCGCTCTGTGCCTgcgc                                                                                                                                   <  7:49770/149‑1 (MQ=255)                        ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGctc                                                                                                                                              <  8:166581/134‑1 (MQ=255)                        ccaggcccaggcCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGctc                                                                                                                                              >  7:166581/1‑134 (MQ=255)                        ccaggcccaggcCGAGCGCCTGGCGCAGGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGGTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCg                                                                                                                                >  1:183164/1‑148 (MQ=255)                          aggcccaggccGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCaa                                                                                                                             <  4:106781/149‑1 (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                                              >  8:195746/1‑136 (MQ=255)                                                                                                                     cAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGc                                              <  7:195746/136‑1 (MQ=255)                                                                                                                        tGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATg                              <  2:115293/149‑1 (MQ=255)                                                                                                                          gCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGcc                            >  2:228822/1‑149 (MQ=255)                                                                                                                               cagccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGCCGATGCCCtgtg                       >  8:290906/1‑149 (MQ=255)                                                                                                                                agccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTg                                                 <  7:17361/122‑1 (MQ=255)                                                                                                                                agccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTg                                                 >  8:17361/1‑122 (MQ=255)                                                                                                                                agccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTgg                      <  2:9949/149‑1 (MQ=255)                                                                                                                                  ccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGc                                                                                                  <  2:301095/72‑1 (MQ=255)                                                                                                                                  ccGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGc                                                                                                  >  1:301095/1‑72 (MQ=255)                                                                                                                                   cGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTg                   <  2:267896/149‑1 (MQ=255)                                                                                                                                                 cccAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGc      <  1:328273/148‑1 (MQ=255)                                                                                                                                                  ccAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTg    >  5:266760/1‑149 (MQ=255)                                                                                                                                                    aGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGa  <  5:296751/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTGGCAGCCAGCCGGTTGCGCGCTTGCCCAGTGCCGCTCTGTGCCTGCGCCCGCAAGTGTCGAAGAAGTTGCGCGGGCGCAGCCTGCATGTGTTGGAGCG‑CGCCCTGCGCGACCTGCCGGTGCCGGTGCTGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGA  >  NC_002947/581785‑582080

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.