New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2632036 =	21 (0.810)	7 (0.300)	4/118	1.9	33.5%	coding (500/1578 nt)	guaA	GMP synthetase (glutamine aminotransferase)
	?	NC_000913	2632263 =	9 (0.390)					coding (273/1578 nt)	guaA	GMP synthetase (glutamine aminotransferase)

GCGATAAAGTTACCGCTATTCCGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2632105‑2632036 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgccatGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCACGC  >  NC_000913/2632263‑2632329                                                                                                                                           GCGATAAAGTTACCGCTATTCCGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCC                                                                       >  1:606733/1‑68                      CGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTCTGCATGCCATAGCAA                                                  >  1:229566/1‑68                      CGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTCTGCATGCCATAGCAA                                                  >  2:369706/1‑68                      CGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTCTGCATGCCATAGCAA                                                  >  2:656925/1‑68                      CGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTC                                                                 <  1:417767/53‑1                      CGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTC                                                                 >  2:417767/1‑53                                 ATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCATGGTCTGCATGCCATAGCAAACGCCGAATAC                                       <  1:656925/68‑1                                                        CTGCCCGTTTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACAT                <  2:229566/68‑1                                                                 TGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCG       >  2:96775/1‑68                                                                    CATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGCCATACTGCGGCGCAC    >  2:99941/1‑68                                                                     ATGGGCTGCATGCCATAGCAAACGCCGATTACCGGTACGCCTGCTTCAAAGCCATACTGCGGCGCACG   <  2:864587/68‑1                                                                      TGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCACGC  <  1:641234/68‑1                                                                                                                                           GCGATAAAGTTACCGCTATTCCGTCCGACTTCATCACCGTAGCCAGCACCGAAAGCTGCCCGTTTGCCAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/2632105‑2632036 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttgccatGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCACGC  >  NC_000913/2632263‑2632329

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