New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2290113	NA (NA)	17 (0.500)	16/486	0.9	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	4581029 =	0 (0.000)					coding (434/1395 nt)	hsdS	type I restriction enzyme EcoKI specificity subunit

ATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGG  >  NC_000913/4580796‑4581276                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        ataaggtcgccggccattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCAccagcc                                                                                                                                                                                                                                        >  2:32697‑M2/234‑250 (MQ=255)               cattaccagccttggggtaaaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAAt                                                                                                                                                                                                                          <  1:32697‑M2/32‑1 (MQ=255)                                  aaactgctctattacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaatgtaagttgatgactcatgatgaaccctgttctatggctccagattacaaacatgatctcatatcagggacttgttcgcaccttccTTAATAGTTGCGCCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGagaa                                                                                                                                                                                                       <  1:38060‑M2/51‑1 (MQ=255)                                  aaacggctcgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAgaga                                                                                                                                                                                                        >  1:252842‑M2/201‑250 (MQ=255)                                          cgatgacttccaccatgttttgccatggcagaatctgctccatgcgggacaagcaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTtgt                                                                                                                                                                                               >  2:176791‑M2/193‑251 (MQ=255)                                                                              gctccatgcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTc                                                                                                                                                           >  2:176803‑M2/157‑251 (MQ=255)                                                                                 ccatgcgggacaagaaaatctcttgtctggtctgacggcgcttactgctgaattcacgggcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATTAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAatataagttgaacatg                                                                                                                                                        <  2:192166‑M2/98‑8 (MQ=255)                                                                                     gcgggacaagaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGt                                                                                                                                                    <  1:176791‑M2/102‑1 (MQ=255)                                                                                             agaaaatctcttttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACAcc                                                                                                                                            <  2:5858‑M2/110‑1 (MQ=255)                                                                                                        tttctggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATgcgc                                                                                                                                 >  1:341424‑M2/131‑251 (MQ=255)                                                                                                             ggtctgacggcgcttactgctgaattcactgtcggcgaaggtaagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCg                                                                                                                               <  1:260932‑M2/123‑1 (MQ=255)                                                                                                                          ttaccgctgaattcactgtaggcgaaggtaagttggtgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTCCGATAAAGAGAAGATTTTGTGAAATGAGCGATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATg                                                                                                               >  2:325227‑M2/113‑251 (MQ=255)                                                                                                                                                       aagttgatgactcatgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACt                                                                                  <  1:176803‑M2/168‑1 (MQ=255)                                                                                                                                                                   catgatgaaccctgttctatggctccagatgacaaacatgatctcatatcagggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTTATGTGCGGATTCACCAACTTCGGATTTGCTc                                                                      <  2:262849‑M2/180‑1 (MQ=255)                                                                                                                                                                                                                      gggacttgttcgcaccttccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGActttct                   >  1:352788‑M2/21‑251 (MQ=255)                                                                                                                                                                                                                                    ccttccTTAATATGTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCTGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAAGATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGAttt       <  1:202952‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                       tccTTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGCAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAgg  >  1:215572‑M2/4‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        ATACAGAATGTTGCGGCTCAAAATTACGCCATTTTTCTGTCAATTTTCCATTAACTGCGCCCCCCAATACCGCTTGACGAAAACGTTTCAGGATTTGTGGGATTTGCTCAAAACGTGCTTTGGTGCTGTCTACCTGCGCCAGCAGCGTATCGAGTTTTTCAGCGATGATTTTTTGTTCGGCAAGTGGTGGGATTGGTATATTTATCAAATCAAAGCTTGCCGGCTTAATATTATTAATATTTGCACCAGCAGAAAGTGATGAAATTTTGTTTCGATAAAGAGAAGATTTTGTGAAATGAGCAATAAAACCAGAAAATATAAGTTTTTCAGGACGTAATACACCGCAAAATGCGCCGAAACTACATTCAAATGGTAGATGCTGATGTGCGGATTTACCAACTACGGATTTGCTCCCTGATGACATTGCAATAACAATATCTTCAGGAGATATTTTTTGACTTTCTTTAACAAGATTTTTAGG  >  NC_000913/4580796‑4581276

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