New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3943623	NA (NA)	7 (0.150)	6/268	5.7	17.0%	intergenic (+114/‑80)	gltU/rrlC	tRNA‑Glu/23S ribosomal RNA of rrnC operon
	?	NC_000913	= 4037456	36 (0.730)					intergenic (+122/‑62)	alaT/rrlA	tRNA‑Ala/23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGC  >  NC_000913/3943540‑3943718                                                                                    |                                                                                                cTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                                      <  1:691918/143‑1 (MQ=32) cTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa                                      <  2:423020/143‑1 (MQ=32)     atatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa                                  <  1:42012/143‑1 (MQ=32)     atatCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaaa                                  <  1:771619/143‑1 (MQ=32)          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt                             >  1:618944/1‑143 (MQ=32)          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACAGCt                             >  2:652502/1‑143 (MQ=17)          tcAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCGCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCt                             >  1:30127/1‑143 (MQ=17)           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTGAAAATCTGGATCAAGCTGAAAATTGAACAACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCAAACGAACCATCtt                            >  1:582383/1‑143 (MQ=11)           cAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTATTTAAAAATCCGGATCAAGTTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATTATGAATCGAAAGAAACATGtt                            >  2:726307/1‑143 (MQ=11)              caCTCATCTTCGGCTGATGTTTGAGATACTTGCTCTTTAAAAATCTGGATCACGCTGAAAATTGAAACACCGAACACCGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg                         <  2:61629/142‑1 (MQ=11)                 tCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGtt                      >  2:267670/1‑143 (MQ=32)                   aTCTTCGGGGGATGTTTGAGATAGTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAATACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtgt                    <  2:432605/143‑1 (MQ=12)                     gttCTGGTGATGTTTGAGATATGTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGa                  <  1:507372/142‑1 (MQ=17)                     cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGTATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGa                  <  1:28990/143‑1 (MQ=17)                     cTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAGCATCTTCAGGTTGCGa                  >  1:286424/1‑143 (MQ=12)                             gATGTGTGAGGTATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGTAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGACTCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGc          <  1:259846/143‑1 (MQ=255)                                 tttGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGGTAAGCGACt      >  2:569187/1‑143 (MQ=17)                                    gagaTATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTACGCGACTAAg   >  2:207514/1‑143 (MQ=17)                                     agaTATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc  <  2:789834/143‑1 (MQ=32)                                                                                    |                                                                                                CTTAATATCTCAAAACTCATCTTCGGGTGATGTTTGAGATATTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGC  >  NC_000913/3943540‑3943718

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