New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4037399 =	NA (NA)	8 (0.210)	6/252	3.9	12.7%	intergenic (+64/‑120)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
	?	NC_000913	4168481 =	56 (1.440)					intergenic (+34/‑160)	gltT/rrlB	tRNA‑Glu/23S ribosomal RNA of rrnB operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGG  >  NC_000913/4037340‑4037542                                                            |                                                                                                                                                tGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatga                                                        >  4:147318/1‑149 (MQ=34)   tAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgatgat                                                       >  7:13689/1‑148 (MQ=34)        aTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATc                                                  >  4:111746/1‑148 (MQ=35)            aataaaAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaag                                             <  2:125921/149‑1 (MQ=37)               aaaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa                                          >  1:97604/1‑149 (MQ=37)                        ttCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg                                 >  4:128073/1‑149 (MQ=37)                          cAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg                                 <  3:128073/147‑1 (MQ=37)                          cAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg                                 <  6:40768/147‑1 (MQ=37)                            gagTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtg                             <  4:20622/149‑1 (MQ=37)                            gagTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtg                             <  4:22818/149‑1 (MQ=37)                             agTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtgt                            >  3:143231/1‑149 (MQ=37)                             agTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtgt                            <  8:13689/149‑1 (MQ=37)                                 tACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAgg                        <  5:20308/149‑1 (MQ=37)                                     tGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGtt                      >  4:148309/1‑147 (MQ=37)                                      gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAg                   >  3:54378/1‑149 (MQ=37)                                       cAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGc                  <  3:148309/149‑1 (MQ=37)                                            gTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt              >  1:187765/1‑148 (MQ=37)                                                aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc          <  3:111746/148‑1 (MQ=37)                                                 cTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGt        <  3:147318/149‑1 (MQ=37)                                                       aGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg  <  8:62172/149‑1 (MQ=32)                                                            |                                                                                                                                                TGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGG  >  NC_000913/4037340‑4037542

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