New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	16 (0.580)	14/468	0.7	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 family transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

GGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTT  >  NC_000913/258471‑258913                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      catttaacccggtgtatcggaattttgacaaaggtcttcaggaaatcaccgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCt                                                                                                                                                                                                  <  2:75863‑M2/38‑1 (MQ=255)                  cggaattttgacaaaggtcttcaggaaatcatcgactggatcgaaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagaggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCACCGTGa                                                                                                                                                                                  >  2:211977‑M2/197‑250 (MQ=255)                                                      tggatcgaaaaacgcttgtgttaaaaatttcataacttttgcccatctggtaaaacagctgcatcccaacaggagatagcaatgacgttccagagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAACCCGGCACCGGAAGTGCGTGCATTCTGGGGCtggtg                                                                                                                                             >  1:287975‑M2/161‑251 (MQ=255)                                                             caaaacgcttgttttaaaaatttgctaaattttgccaatttggaaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTAAACGTGAAACCGGCACCGGcagtgcgagcatt                                                                                                                                                         >  2:173120‑M2/154‑226 (MQ=255)                                                             aaaaacgcttgtgttaaaaatttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCGAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAAtt                                                                                                                                                         <  1:173120‑M2/79‑1 (MQ=255)                                                                   gcttgtgttaaaaatttgctaaattttgccaatgtggtaaaacagtttcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCGTTTGCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATGCTGGGGCTGGTGGATGGAGCTGGaa                                                                                                                                <  2:101537‑M2/104‑1 (MQ=255)                                                                                  ttgctaaattttgccaatttggtaaaacagttgcatcacaacaggagatagcaatgacgttaccgagtggacacccgaagcgcagattgatcaaaaaatttaccgcactaggcccgtatattcgagaaggtaAGTGCAAAGATACTCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCCGGAAGCGCAGGAATCCCGt                                                                                                                 >  2:270495‑M2/133‑251 (MQ=255)                                                                                   tgctaaatttagccaatttgggaaaacagtggcatcacaacaggagatataaatgacgttaccgagtggacacccgaaaagccgcttgatcaaaaaatttaccgcactaggcccgtatattcgtgacggtaAGTGCAAAGATAATAGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAGGTGCGTGCAGTCTGGGGCTGGGGGGAGGAGCTGGAAGCGAAGGACTCCCGtt                                                                                                                >  2:99848‑M2/132‑251 (MQ=255)                                                                                                                                                      ggacacccgaagagcagattgatcaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGa                                             <  1:248839‑M2/187‑1 (MQ=255)                                                                                                                                                             cgaagagcagattgatcaaaaaagttacctcactaggcccgtatatgcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCTTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCGGGAAGCGCAGGAATCCCGTTTTACCTACAGTGACCAGTGTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTATAATAcacttaa                                      <  2:187582‑M2/194‑4 (MQ=255)                                                                                                                                                               aagagcagattgataaaaaaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGt                                    >  1:236049‑M2/56‑251 (MQ=255)                                                                                                                                                                                 aaatttaccgcactaggcccgtatattcgtgaaggtaAGTGCCAAGATAAACGATTCTTTTTCGATTATCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGACTCCCGTTTTACCTACAGTTCACAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTCAAAGACACTGAAGTGGTTGACCGCCTGGAGc                   >  1:286515‑M2/38‑250 (MQ=255)                                                                                                                                                                                       accgcactagtgccggtgatttgggAAGGTAAGTGCAAAGAGAATCTATTCTTGTTCGAGTGTCTGGCTGGATGCGGCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAAGCCCGTGTTACCTACAGTTACCAGTTTGGGCTGTTCGATAAAGCAGTCTACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGTACGTCTGGAGCACACCCt            <  2:103276/227‑1 (MQ=255)                                                                                                                                                                                       accgcactaggcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCt            >  1:306242‑M2/32‑251 (MQ=255)                                                                                                                                                                                                 gcccgtatattcgtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGttt  >  2:251980‑M2/22‑251 (MQ=255)                                                                                                                                                                                                             gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGt                                                                                     >  1:311143‑M2/10‑156 (MQ=255)                                                                                                                                                                                                             gtgaaggtaAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGt                                                                                     <  2:311143‑M2/147‑1 (MQ=255)                                                                                                                                                                                                                     aaGTGCACAGCTAATCGATTCTTTTTCGCTTGTCTGGCTGTATGCGTCAACGTGCAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTaaa                                                >  2:30576/1‑185 (MQ=255)                                                                                                                                                                                                                     aaGTGCAAAGATAATCGATTCTTTGTCGATTGTCTGGCTGTATGCGGCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTaaa                                                <  1:30576/185‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      GGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCTGAAGGTAAGTGCAAAGATAATCGATTCTTTTTCGATTGTCTGGCTGTATGCGTCAACGTGAAACCGGCACCGGAAGTGCGTGAATTCTGGGGCTGGTGGATGGAGCTGGAAGCGCAGGAATCCCGTTTTACCTACAGTTACCAGTTTGGTCTGTTCGATAAAGCAGGCGACTGGAAGAGTGTTCCGGTAAAAGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTT  >  NC_000913/258471‑258913

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