New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	1629047 =	39 (0.850)	4 (0.090)	4/272	6.3	10.3%	intergenic (+9/‑168)	rspR/ydfZ	transcriptional repressor for rspAB/selenoprotein, function unknown
	?	NC_000913_3_bme_pgi	1629071 =	32 (0.730)					intergenic (+33/‑144)	rspR/ydfZ	transcriptional repressor for rspAB/selenoprotein, function unknown
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTGATATTTTCTCTGCCCTGCAAAAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGT  >  NC_000913_3_bme_pgi/1628905‑1629191                                                                                                                                               |                                                                                                                                                 ttGATATTTTCTCTGCCCTGCAAAAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTcat                                                                                                                                            <  1:241355/149‑1 (MQ=255)  tGATATTTTCTCTGCCCTGCAAAAACGTGATGGCGATGCGGGAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTcatc                                                                                                                                           <  2:465195/149‑1 (MQ=255)  tGATATTTTCTCGGCCCTGCAAAAACGGGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTcatc                                                                                                                                           <  1:589109/149‑1 (MQ=255)    atatTTTCTCTGCCCTGCAAAAACGTGAGGGCGAAGCGGGAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTcatccc                                                                                                                                         <  1:69112/149‑1 (MQ=255)           ctctTCCCTGCAAAAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCggg                                                                                                                                  <  1:102969/149‑1 (MQ=255)           ctctGCCCTGGAAAAACGGGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCggg                                                                                                                                  <  2:261856/149‑1 (MQ=255)                 ccTGCAAAAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTgaga                                                                                                                            <  1:7906/149‑1 (MQ=255)                        aaaCGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAAAGTCtttt                                                                                                                     >  1:567844/1‑149 (MQ=255)                            gTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCCTTTCCTCTCATCCCATCCGGGGGGGGAGTCTTTTcccc                                                                                                                 >  1:310273/1‑149 (MQ=255)                             tGATAGAGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTccccc                                                                                                                >  2:75297/1‑149 (MQ=255)                                              aaCGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTcatgca                                                                                               >  2:149353/1‑149 (MQ=255)                                                cGTGGAATGACGGAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTcatgcatg                                                                                             <  2:254164/149‑1 (MQ=255)                                                     aaTGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCa                                                                                        >  1:46910/1‑149 (MQ=255)                                                           gCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCaaaaagg                                                                                  >  1:257601/1‑147 (MQ=255)                                                               cATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGGTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGAtgt                                                                              <  1:622692/149‑1 (MQ=255)                                                               cATTTGCAGGAAATCAGCCAATGCGTGCGCTAGACGCGTGAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATTGCTCATGCATGCATCAAAAAAGAtgt                                                                              <  2:448485/149‑1 (MQ=255)                                                                       ggAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAATGATGTGAGGTTGa                                                                      >  1:443043/1‑149 (MQ=255)                                                                        gAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCctc                                                                                                                                                <  1:186841/74‑1 (MQ=255)                                                                        gAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCctc                                                                                                                                                >  2:186841/1‑74 (MQ=255)                                                                            tCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGACGAGTAATTCATTTCCTCTCATCCCATCCGGGCAGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCCTCAAAAAACGTCTGTGCTTCATccac                                                                 >  2:235396/1‑146 (MQ=255)                                                                                   aTCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATaaaaaacaaaa                                                          >  1:583316/1‑149 (MQ=255)                                                                                     ccGTGCGCCAGAGCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACaaaaaa                                                        <  2:39754/149‑1 (MQ=255)                                                                                       gTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGCCTCATGCATTCATCAAAAAAGATGTGAGCTTGCTCAACAACAAAAAata                                                      >  2:140389/1‑149 (MQ=255)                                                                                                aTCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCg                                             >  1:364629/1‑149 (MQ=255)                                                                                                      gaGGAAAGCGGCGAGTGGGTTGGGGAAGAGGAATTAAGTTCCGCTCAGGCGATCCGGGGTGAGAGTGTTTTCCCCCGCATTAGGGCTCATGCATGCATAAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGa                                       <  2:123827/148‑1 (MQ=255)                                                                                                      cAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGa                                       <  1:417539/149‑1 (MQ=255)                                                                                                               aGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTtata                              >  2:284746/1‑149 (MQ=255)                                                                                                                   aCTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTgc                          >  1:772289/1‑149 (MQ=255)                                                                                                                     tGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTgcgc                        <  1:462817/149‑1 (MQ=255)                                                                                                                       gTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGGATTTatatgactccc                      >  2:702944/1‑142 (MQ=255)                                                                                                                              gAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGt               <  1:284746/149‑1 (MQ=255)                                                                                                                                agagTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTgg             >  2:430989/1‑149 (MQ=255)                                                                                                                                   gTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATC‑AAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTc                                              <  2:1600/112‑1 (MQ=255)                                                                                                                                   gTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATC‑AAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTc                                              >  1:1600/1‑112 (MQ=255)                                                                                                                                    tAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGACTTATGGCTCATGCATGCATCAAAAAAGATGTAAGCTTGATCAAAAACAAAAAATATTTCACTCGACCGGAGTAGTTATATTGCGCCCGTTACGTGGGCtt         >  2:252028/1‑149 (MQ=255)                                                                                                                                       ttCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGa      <  1:739199/149‑1 (MQ=255)                                                                                                                                       tgcATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGa      <  1:149353/147‑1 (MQ=255)                                                                                                                                         cATTTCCCCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTAGGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACt    <  2:601264/149‑1 (MQ=255)                                                                                                                                           tttCCTCTCATCCCATCCTGGGTGAGAGTCTTTTCCCCCGCCTTATCGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGCAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGt  <  1:690822/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                 TTGATATTTTCTCTGCCCTGCAAAAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGT  >  NC_000913_3_bme_pgi/1628905‑1629191

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