New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4162839	121 (1.370)	3 (0.040)	3/418	NT	2.4%	coding (432/1101 nt)	trmA	tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
	?	NC_000913	= 4162849	126 (1.480)					coding (422/1101 nt)	trmA	tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Rejected: Frequency below/above cutoff threshold.

GAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4162597‑4162839 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcgcgATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTAC  <  NC_000913/4162849‑4162617                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              tAAAAGTTTCCGTTGCCGCAGTCCCGCTCCAGTAAATTGCGTTTTGAGCCTTTGGTTTGGTCCAGCCCCCAGTCCAGCAGTGGAACATTCATCGCCGTGTTCGCGTGGGTAAAGCTGTTTTCTTCCTGACGGTAGATCATCTCCTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTG                                                                                                                                                                                                                                               <  2:227950/238‑1  AAAAGTTCCCTTTGCCCGAATACAGCTCCAGTAAATCGACTTTTGAGCCTTGTGTTACGGCCAGCGCCCATTCCAGCATCTGAATATTCCTCGCCGGGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTGTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGC                                                                                                                                                                                                                                              <  2:831664/239‑1  AAAAGTTACCGTTGCCGCAGTACAGCTCCAGGAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGC                                                                                                                                                                                                                                              <  1:1163203/239‑1         ACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCTTCTATGTAATCCGTATCAAGCCCGTTTTTGGCTTTCGTTGCCCGACCAACCACATGCTCATTCAATTTCTGCccccaca                                                                                                                                                                                                                                       >  2:513453/1‑232                                                                                                                                                        CGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGAGGAGTTACCGTCAGGTAGAAAACAGCTTTACCC                                                                                        >  1:807306/1‑239                                                                                                                                                                                       ATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAAAAGCTTTACCCAGCCGAAGCCGGCGCGGAATATCCAGAg                                                            >  1:31778/1‑235                                                                                                                                                                                                       GTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATATTCAGATG                                                           <  1:79563/221‑1                                                                                                                                                                                                                                              TGCGCGATTCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTAC  <  1:1182066/239‑1                                                                                                                                                                                                                                              TGCGCGATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGCCGTAACCAAAGGCTCAAACGGCAATTTACTGGAGCTGa    >  1:1040191/1‑236                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              GAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4162597‑4162839 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcgcgATGCACTGCGCGCGCAGAATCTGAATGTGCATCTGATTGGTCGGGCAACGAAAACCAAAATCGAGCTGGATCAGGATTACATCGATGAACGTCTGCCGGTCGCAGGGAAAGAGATGATCTACCGTCAGGTAGAAAACAGCTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTAC  <  NC_000913/4162849‑4162617

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