New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4436762 =	25 (0.980)	4 (0.160)	4/486	2.4	13.8%	coding (8/741 nt)	cysQ	3'(2'),5'‑bisphosphate nucleotidase
	?	NC_000913	= 4499511	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region
Rejected: Frequency below/above cutoff threshold.

CAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGA  >  NC_000913/4436526‑4437011                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          cAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGc                                                                                                                                                                                                                                             >  1:98383/1‑251 (MQ=255)     cAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGATAATAAACGTCTTACTTATAGAACAGGGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACATTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGccctcttg                                                                                                                                                                                                                                          <  2:47174/250‑7 (MQ=255)             atagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATCAAGTATGCCAGCTTGCACGGAATg                                                                                                                                                                                                                                 <  1:37823‑M1/27‑1 (MQ=255)                           cTAAACTTAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCAtt                                                                                                                                                                                                                   >  1:17014/1‑251 (MQ=255)                                  tAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAg                                                                                                                                                                                                             <  1:208724/250‑1 (MQ=255)                                        aagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTacg                                                                                                                                                                                                      <  1:90235‑M1/54‑1 (MQ=255)                                                              gggAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGAcgt                                                                                                                                                                                >  1:46490/1‑251 (MQ=255)                                                                                   aaaTAAACGTCTTACTTATAGAACAGTGAAGAATTCCACAATTTTACGCTTTGAAAATGATGTCACTAGCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGTTGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATtc                                                                                                                                                           <  2:184821/250‑1 (MQ=255)                                                                                         aCGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTCTGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACACTTCTCCGGt                                                                                                                                                     >  1:237868/1‑251 (MQ=255)                                                                                                   tataGAACAGTGAAGAATGCCACAATTTGACGCTTTGAAAATGATGACACTATCACATTTGGCGCATGCATTAACGATAGGGTATAAGTAAAACAATAAGTTATCACCGCTCACAGAGACGAGGTGGAGAAATGGTAGATCAAGTCTGCCAGCTTGCACGGAATGCAGGCGATGCCACTATGCAGGTCTACGACGGGTCGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTATCCGTCAACGGCac                                                                                                                                             <  2:190566/249‑2 (MQ=255)                                                                                                    ataGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCgg                                                                                                                                          <  1:292624/251‑1 (MQ=255)                                                                                                      aGAACAGTGAAGACTGCCACAATTTTACGCTTTGAAAATACTGACACTATCACAGTTGGAGCATTCATTAACGTTAGGGTATAAGTAAAACAATAAGTTCACACCGCTCACAGAGAAGAGGAGGAGACATGTTAGATCAAGTATGCCAGCTTGCACGGACTGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGTCAGCggct                                                                                                                                        >  2:273138/1‑249 (MQ=255)                                                                                                                     gCCACAATTTTACGCTTTGAAAATGATGACCATATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGCCGACAGCAAAGCGGCCAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACAcc                                                                                                                         >  2:272305/1‑251 (MQ=255)                                                                                                                              atggcatccgcatagggcgcggggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAGTTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATg                                                                                                                >  2:37823‑M1/112‑251 (MQ=255)                                                                                                                                          aaTGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTtacgg                                                                                                    >  1:138187/1‑250 (MQ=255)                                                                                                                                                       tatcgctcgccacgggaatagctgcggcagcgggcttggaatgggttaagtgataacagatgtctggaaatataggggcaaatccaATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTCCGCTGACACCg                                                                                        >  2:90235‑M1/87‑250 (MQ=255)                                                                                                                                                         acaGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGAAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCgg                                                                                       >  2:57932/1‑249 (MQ=255)                                                                                                                                                               ggCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTcc                                                                                <  1:9648/250‑1 (MQ=255)                                                                                                                                                                                   ggTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTACTTTCTGAAGAAGAtcc                                                           >  1:250695/1‑251 (MQ=255)                                                                                                                                                                                   ggTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGAAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGAGGATATTGCCGCTCACACCGTTATCATGGTCGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGCAGAtcc                                                           >  1:178014/1‑251 (MQ=255)                                                                                                                                                                                        aaGTAAAACAATAAGTTAACCCCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCg                                                      <  2:65303/251‑1 (MQ=255)                                                                                                                                                                                          gTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGGATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCCCCCGGt                                                    <  2:46490/251‑1 (MQ=255)                                                                                                                                                                                                cAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGaa                                              <  1:272305/251‑1 (MQ=255)                                                                                                                                                                                                    aaGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTcc                                          >  1:166501/1‑251 (MQ=255)                                                                                                                                                                                                             ccGCTCACAGAGAAGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGGCGTCAGCAAAGCGGACAATTCTCCGGGAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTGCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACt                                 <  2:17014/251‑1 (MQ=255)                                                                                                                                                                                                                  cacaGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGc                              >  1:57304/1‑249 (MQ=255)                                                                                                                                                                                                                        gaCGAGGGGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGGTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGccattgcagctttac                      >  2:113412/1‑237 (MQ=255)                                                                                                                                                                                                                                          tAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCt     <  1:195640/250‑1 (MQ=255)                                                                                                                                                                                                                                            gATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGGTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGa  <  1:261354/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          CAATCAGCGTGGCCAGGAGCGTTGCGCTAAACTTAATCATCAGGGACATCCTTTTATCATCGGGAATACGAAAGAAAAGGGAGAATAAACGTCTTACTTATAGAACAGTGAAGAATGCCACAATTTTACGCTTTGAAAATGATGACACTATCACAGTTGGCGCATTCATTAACGATAGGGTATAAGTAAAACAATAAGTTAACACCGCTCACAGAGACGAGGTGGAGAAATGTTAGATCAAGTATGCCAGCTTGCACGGAATGCAGGCGATGCCATTATGCAGGTCTACGACGGGACGAAACCGATGGACGTCGTCAGCAAAGCGGACAATTCTCCGGTAACGGCAGCGGATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGA  >  NC_000913/4436526‑4437011

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