New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1941980 =	8 (0.590)	4 (0.320)	4/248	NT	29.1%	coding (604/933 nt)	znuA	Zn(2(+)) ABC transporter periplasmic binding protein
	?	NC_000913	1942149 =	12 (0.950)					coding (435/933 nt)	znuA	Zn(2(+)) ABC transporter periplasmic binding protein

GCATCGTGAAAAACGAAATAACCTTTCCCCTTGAGCGGCGCGAGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTCACATCTTCAAGCTGC  >  NC_000913/1941958‑1942269                                                                                                                                                                                                |                                                                                                                         gCATCGTGAAAAACGAAATAACCTTTCCCCTTGAGCGGCGCGAGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATc                                                                                       >  1:232978/1‑227 (MQ=255)         aaaaaCGAAATAACCTTTCCCCTTGAGCGGCGCGAGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTc                                                               <  1:264085/243‑1 (MQ=255)                                        gcgAGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTg                                                                                                           >  1:186553/1‑168 (MQ=255)                                           aGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGaaa                                                                                                        <  2:186552/168‑1 (MQ=255)                                                                                    cTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACttttt                                                                          >  1:211708/1‑157 (MQ=255)                                                                                    cTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTCACATCTTCAAGc     >  1:226825/1‑226 (MQ=255)                                                                                       aaaaTCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTcc                                                                       <  2:211707/157‑1 (MQ=255)                                                                                       aaaaTCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTCACATCTTCAAGCTgc  <  2:226824/226‑1 (MQ=255)                                                                                                                                               tcaaccgaaacgcaggttggtaacgagctcgcgccgctcaaggggaaagCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGa                                                           >  1:239195‑M2/50‑113 (MQ=255)                                                                                                                                                  accgaaacgcaggttggtaacgagctcgcgccgctcaaggggaaagCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGatca                                                        <  2:239194‑M2/67‑1 (MQ=255)                                                                                                                                                                             gcgccgctcaaggggaaagCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCcg                                                                      >  1:251126‑M2/20‑72 (MQ=255)                                                                                                                                                                               gcTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATc                                                                                       >  2:91590/1‑53 (MQ=255)                                                                                                                                                                               gcTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGt                    >  2:101717/1‑120 (MQ=255)                                                                                                                                                                                ccgctcaaggggaaagCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGt                                                                   <  2:251125‑M2/56‑1 (MQ=255)                                                                                                                                                                                  aTCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTc                                                                                    <  1:91591/53‑1 (MQ=255)                                                                                                                                                                                  aTCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTc                 <  1:101718/120‑1 (MQ=255)                                                                                                                                                                                              agCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTCACATCTTCAAGc     >  1:263538‑M2/3‑120 (MQ=255)                                                                                                                                                                                                |                                                                                                                         GCATCGTGAAAAACGAAATAACCTTTCCCCTTGAGCGGCGCGAGCTCGTTACCAACCTGCGTTTCGGTTGAGGCTAATTGTGCCTCAAAATCCTTCAGGTTGGCGTCAAGTTTGGCTCGACTTTGCGGCATAAGTTCCACTAATTTTCCATGGATTGCAACCGCTGTAGCCCGCGCTATCTCTGGGGAAAGCCAAAGATGCATGTTGAAATCGCCGTGATGGTGATCTTCGTCACTTTTTTCCGCGTGGTCGTGATCATCATCATCGCCGTGAATACTTTTCATCAGTAGCGGTTTCACATCTTCAAGCTGC  >  NC_000913/1941958‑1942269

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