New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3159416 =	57 (0.500)	3 (0.030)	3/262	NT	5.8%	coding (1731/2220 nt)	ygiQ	radical SAM superfamily protein YgiQ
	?	NC_000913	3159449 =	43 (0.400)					coding (1698/2220 nt)	ygiQ	radical SAM superfamily protein YgiQ

GAAGAAGGGCCGTTATCGAAGATGATGAAGCCGGGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3159565‑3159416 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatgaaATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATA  >  NC_000913/3159449‑3159592                                                                                                                                                                                                                                                                                                        GAAGAAGGGCCGTTATCGAAGATGATGAAGCCGGGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAA                                                                                                                                                  <  1:506186/150‑1           CGTTATCGAAGATGATGAAGCCGGGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGAT                                                                                                                                                     <  1:521866/137‑1           CGTTATCGAAGATGATGAAGCCGGGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGAT                                                                                                                                                     >  2:521866/1‑137                                  GGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAA                                                                                                                                                  <  1:139747/117‑1                                  GGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAA                                                                                                                                                  >  2:139747/1‑117                                                                                                               CCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCC                                    >  1:103896/1‑150                                                                                                                 CCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCT                                  >  1:300909/1‑150                                                                                                                  CTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTT                                 >  2:1766939/1‑150                                                                                                                                                 GATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATA  >  2:264698/1‑150                                                                                                                                                 GATGAAATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATA  >  2:80173/1‑150                                                                                                                                                                                                                                                                                                        GAAGAAGGGCCGTTATCGAAGATGATGAAGCCGGGCATGGGCAGCTATGACCGCTTTAAAGAGCTGTTCGATACTTACTCAAAACAGGCAGGTAAAGAGCAGTATCTGATCCCCTATTTCATCTCCGCGCACCCCGGTACGCGTGATGAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3159565‑3159416 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gatgaaATAGGGGATCAGATACTGCTCTTTACCTGCCTGTTTTGAGTAAGTATCGAACAGCTCTTTAAAGCGGTCATAGCTGCCCATGCCCGGCTTCATCATCTTCGATAACGGCCCTTCTTCGGTATGTTCCGGGGCAATCTTCAGATA  >  NC_000913/3159449‑3159592

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