New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 219448	46 (0.710)	4 (0.070)	4/244	NT	8.2%	coding (147/708 nt)	tsaA	tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
	?	NC_000913	= 219458	47 (0.770)					coding (137/708 nt)	tsaA	tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent

TGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/219298‑219448 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaggccgACGCCGTTCGCGGCCTGGAAGCGTTCAGCCATTTATGGATCCTTTTCGTCTTTCATCAAACGATGGAAGGCGGCTGGCGTCCGACTGTGCGTCCACCGCGCCTCGGCGGTAACGCCAGAATGGGGGTTTTCGCAACACGCTCTAC  <  NC_000913/219458‑219314                                                                                                                                                                                                                                                                                                          TGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGC                                                                                                                                                     <  4:285545/149‑1                 GTAGAGCGTGTTGCGAAAACCCCCATCCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCCTTTGATGAACGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGACGCCGTTCGCGGC                                                                                                                                     <  2:241932/149‑1                    GAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGACGCCGTTCGCGGCCTG                                                                                                                                  <  2:530362/149‑1                                                                                             CTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGACGCCGTTCGCGGCCTGGAA                                                                                                                               >  1:387700/1‑79                                                                                             CTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGACGCCGTTCGCGGCCTGGAA                                                                                                                               <  2:387700/79‑1                                                                                                                                                CCAGGCCGACGCCGTTCGCGGCCTGGAAGCGTTCAGCCATTTATGGATCCTTTTCGTCTTTCATCAAACGATGGAAGGCGGCTGGCGTCCGACTGTGCGTCCACCGCGCCTCGGCGGTAACGCCAGAATGGGGGTTTTCGCAACACGCT      >  1:154917/1‑149                                                                                                                                                  AGGCCGACGCCGTTCGCGGCCTGGAAGCGTTCAGCCATTTATGGATCCTTTTCGTCTTTCATCAAACGATGGAAGGCGGCTGGCGTCCGACTGTGCGTCCACCGCGCCTCGGCGGTAACGCCAGAATGGGGGTTTTCGCAACACGCTCT    >  1:176673/1‑149                                                                                                                                                    GCCGACGCCGTTCGCGGCCTGGAAGCGTTCAGCCATTTATGGATCCTTTTCGTCTTTCATCAAACGATGGAAGGCGGCTGGCGTCCGACTGTGCGTCCACCGCGCCTCGGCGGTAACGCCAGAATGGGGGTTTTCGCAACACGCTCTAC  <  3:473708/149‑1                                                                                                                                                                                                                                                                                                          TGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/219298‑219448 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaggccgACGCCGTTCGCGGCCTGGAAGCGTTCAGCCATTTATGGATCCTTTTCGTCTTTCATCAAACGATGGAAGGCGGCTGGCGTCCGACTGTGCGTCCACCGCGCCTCGGCGGTAACGCCAGAATGGGGGTTTTCGCAACACGCTCTAC  <  NC_000913/219458‑219314

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