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.

AAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTAT  >  NC_000913_3_bme_pgi/1628928‑1629216                                                                                                                                                |                                                                                                                                                  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)                                                     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)                                                                                                                                   tccGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCAGAAAAAAAGATGTGAGCTTGCTCAAAAACAAAAAATATTTCACTCGACAGGCTGATTTATATTGCGCCCGTTACGTAGGCTGAGACTGTAAATGAAGAAGgaga            >  1:368320/1‑149 (MQ=255)                                                                                                                                         gTGACAGGCTTTCCCCTCGCATTATGGCTAATACATACATCTAAACAGATGTTAGCTTGATCAAAAGCCAAAAATATTTCATCTGACAGGAGGATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAacac      <  2:62554/149‑1 (MQ=255)                                                                                                                                             gagTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGGGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATTAGTAAGGGGAGAACACCTAt  >  1:147239/1‑149 (MQ=255)                                                                                                                                             gagTCTATTCCGCCGCCTTATGGCTAATGCAGGCATCAAAAAAGATGTGAGCTTAAACAAAAACAAAAAATATTTCACTCGCCAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTTGACTGTAAATCAGAAAGGAGAAAACACCTAt  <  2:630832/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  AAACGTGATGGCGATGCGGTAGAACGTGCAATGACGCAACATTTGCAGGAAATCAGCGAATCCGTGCGCCAGATCCGCCAGGAAAACAGCGACTGGTTTAGCGAAGAGTAATTCATTTCCTCTCATCCCATCCGGGGTGAGAGTCTTTTCCCCCGCCTTATGGCTCATGCATGCATCAAAAAAGATGTGAGCTTGATCAAAAACAAAAAATATTTCACTCGACAGGAGTATTTATATTGCGCCCGTTACGTGGGCTTCGACTGTAAATCAGAAAGGAGAAAACACCTAT  >  NC_000913_3_bme_pgi/1628928‑1629216

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