New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 581532	39 (0.920)	4 (0.100)	4/244	NT	9.0%	coding (822/1428 nt)	selA	selenocysteine synthase
	?	NC_002947	= 581553	44 (1.100)					coding (843/1428 nt)	selA	selenocysteine synthase

TTCACCACCCAAGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAA  >  NC_002947/581419‑581699                                                                                                                                        |                                                                                                                                                   ttCACCACCCAAGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGa                                                                                                                                      >  7:190393/1‑149 (MQ=255)   caccacCCAAGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATg                                                                                                                                    >  5:47539/1‑149 (MQ=255)          cAAGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGa                                                                                                                             <  8:67792/149‑1 (MQ=255)           aaGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGa                                                                                                                             >  8:190924/1‑148 (MQ=255)                     tGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCt                                                                                                                  >  2:128465/1‑149 (MQ=255)                               gCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGTCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGa                                                                                                        >  1:50197/1‑149 (MQ=255)                                 aCGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACa                                                                                                      <  6:47539/149‑1 (MQ=255)                                        gCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCGTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAActgct                                                                                               <  5:76895/149‑1 (MQ=255)                                            tccaACATGAATTTCCACTGCTTGAAGACCTGGGCAGAGTCATCCTTCTCTACCTTACCGGCTGGGGCCTGCCC‑GCCGAAGCCACGGTACGCCAGGCGCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTggg                                                                                            <  4:150064/147‑1 (MQ=255)                                            accaACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCATCCCA‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTACAGTGGCGACAAACTCATgggg                                                                                           >  3:177620/1‑149 (MQ=255)                                              caACATGAATTACCACTGCTTGAAGACCTGGGCAGCTGCAGACTGCTCGACCTTACACGCTGGGGCCTGCCC‑GCCGAACCCACGCTACGCCACGCCCTGGCCCATCGTGCCGATATCGTCACCTTCAGTGGCGACAAAcagctgcg                                                                                            >  3:150064/1‑144 (MQ=255)                                              caACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTggg                                                                                            >  3:11066/1‑146 (MQ=255)                                              caACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCACCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTggg                                                                                            <  4:11066/146‑1 (MQ=255)                                                    gAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGACTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGTGCCGc                                                                                    >  3:104647/1‑148 (MQ=255)                                                       ttACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGc                                                                                >  2:142238/1‑149 (MQ=255)                                                          ccACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCgg                                                                             >  1:129112/1‑149 (MQ=255)                                                           cACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCgg                                                                             >  4:173545/1‑148 (MQ=255)                                                              tGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCa                                                                                   >  6:265869/1‑139 (MQ=255)                                                              tGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCa                                                                                   <  5:265869/139‑1 (MQ=255)                                                                       aCCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGc                                                                <  3:173545/149‑1 (MQ=255)                                                                           gggCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCCGGGCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCca                                                              >  1:192170/1‑148 (MQ=255)                                                                                      gcCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGTCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATc                                                 <  2:50197/149‑1 (MQ=255)                                                                                       ccTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGa                                                                                                        >  7:141485/1‑93 (MQ=255)                                                                                       ccTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGa                                                                                                        <  8:141485/93‑1 (MQ=255)                                                                                         tGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGa                                                                           <  7:72448/120‑1 (MQ=255)                                                                                         tGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGa                                                                           >  8:72448/1‑120 (MQ=255)                                                                                            tCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGa                                                        >  4:115130/1‑136 (MQ=255)                                                                                            tCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGa                                                        <  3:115130/136‑1 (MQ=255)                                                                                                 cTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCaa                                      <  7:190924/149‑1 (MQ=255)                                                                                                  ttACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATc                                        >  6:96236/1‑146 (MQ=255)                                                                                                      ccGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGaaa                                  >  6:208695/1‑148 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGTTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATcc                             <  8:190393/147‑1 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGGTACGCcaggccccagcgggtaaggtcgagcaggctgccgctgcccaggtcttcaagcagtggtaattcatgttggtgggcgatacgtgccagctcggcagtaggcacttgggtggtg                                    >  1:152150‑M2/1‑28 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGGTACGCcaggccccagcgggtaaggtcgagcaggctgccgctgcccaggtcttcaagcagtggtaattcatgttggtgggcgatacgtgccagctcggcagtaggcacttgggtggtg                                    <  2:152150‑M2/140‑113 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATcc                             <  2:129112/147‑1 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGc                           >  5:220726/1‑149 (MQ=255)                                                                                                            gggCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGc                           <  5:96236/149‑1 (MQ=255)                                                                                                              gCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCt                          >  3:264390/1‑148 (MQ=255)                                                                                                               ccTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTgg                                                                                             >  1:128715/1‑80 (MQ=255)                                                                                                               ccTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTgg                                                                                             <  2:128715/80‑1 (MQ=255)                                                                                                               ccTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAACCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCta                         >  8:46985/1‑147 (MQ=255)                                                                                                                cTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTc                         >  8:34191/1‑147 (MQ=255)                                                                                                                  gCCC‑GCCGAACCCACGGTACGCcaggccccagcgggtaaggtcgagcaggctgccgctgcccaggtcttcaagc                                                                                                <  3:220184‑M2/74‑53 (MQ=255)                                                                                                                  gCCC‑GCCGAACCCACGGTACGCcaggccccagcgggtaaggtcgagcaggctgccgctgcccaggtcttcaagc                                                                                                >  4:220184‑M2/1‑22 (MQ=255)                                                                                                                     c‑gccgAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAgcgtg                  >  5:227637/1‑149 (MQ=255)                                                                                                                               cACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGcgcg         >  3:140650/1‑149 (MQ=255)                                                                                                                                aCGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGt        >  7:144141/1‑149 (MQ=255)                                                                                                                                  ggTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTc       >  3:187632/1‑148 (MQ=255)                                                                                                                                      cGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACaa  <  6:227637/149‑1 (MQ=255)                                                                                                                                       gCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAgcgt                   >  3:172112/1‑131 (MQ=255)                                                                                                                                       gCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAgcgt                   <  4:172112/131‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                   TTCACCACCCAAGTGCCTACTGCCGAGCTGGCACGTATCGCCCACCAACATGAATTACCACTGCTTGAAGACCTGGGCAGCGGCAGCCTGCTCGACCTTACCCGCTGGGGCCTGCCC‑GCCGAACCCACGGTACGCCAGGCCCTGGCCGATGGTGCCGATATCGTCACCTTCAGTGGCGACAAACTGCTGGGGGGGCCGCAAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAA  >  NC_002947/581419‑581699

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