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)

AACCTGCGCGTAGCCAAATTCACGTTCGTTAGAGGCTTCAACGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCACGC  >  NC_000913/2632197‑2632329                                                                   |                                                                   aaCCAGCGCGTAGCCAAATTCACGTTCGTTAGAGGCTTCAACGTGACCGCCCAACTGCATTGCCATgg                                                                   <  1:308648/68‑1 (MQ=255)    cTGCGCGTAGCCAAATTCACGTTCGTTAGAGGCTTCAACGTGACCGCCCAACTGCATTGCCATGGTCt                                                                <  2:658792/68‑1 (MQ=255)                  cgtccgacttcatcaccgtagccagcaccgaaagctgcccgtttgccatGGTc                                                                 <  1:417767‑M2/4‑1 (MQ=255)                  cgtccgacttcatcaccgtagccagcaccgaaagctgcccgtttgccatGGTc                                                                 >  2:417767‑M2/50‑53 (MQ=255)                  cgtccgacttcatcaccgtagccagcaccgaaagctgcccgtttgccatGGTCTGCATGCCATAGCaa                                                  >  1:229566‑M2/50‑68 (MQ=255)                  cgtccgacttcatcaccgtagccagcaccgaaagctgcccgtttgccatGGTCTGCATGCCATAGCaa                                                  >  2:369706‑M2/50‑68 (MQ=255)                  cgtccgacttcatcaccgtagccagcaccgaaagctgcccgtttgccatGGTCTGCATGCCATAGCaa                                                  >  2:656925‑M2/50‑68 (MQ=255)                      aCGTTCGTTAGAGGCTTCAACGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGc                                              <  1:809247/68‑1 (MQ=255)                             atcaccgtagccagcaccgaaagctgcccgtttgccatGGTCTGCATGCCATAGCAAACGCCGAATAc                                       <  1:656925‑M2/30‑1 (MQ=255)                                        aaCGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTg                            <  1:389035/68‑1 (MQ=255)                                         aCGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGc                           <  1:347937/68‑1 (MQ=255)                                         aCGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGc                           >  1:740097/1‑68 (MQ=255)                                          cGTGAAAGCCCAACTGCATTGCCATGGTCTGCATGAAATAGAAAACTCAGAATACCGGTACGCCTGCt                          >  1:334939/1‑68 (MQ=255)                                                  cccAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGAc                  <  2:968692/68‑1 (MQ=255)                                                    ctgcccgtttgccatGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACAt                <  2:229566‑M2/53‑1 (MQ=255)                                                        tGCATTGCCGTGGTCTGCGTGCCATAGCAAACGCCGAATACCGGTACGCCTGCGTCAAAGCCATAcag            >  2:196444/1‑66 (MQ=255)                                                             tgccatGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTgcggcg       >  2:96775‑M2/7‑68 (MQ=255)                                                                catGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGCCATACTGCGGCGCAc    >  2:99941‑M2/4‑68 (MQ=255)                                                                 atGGGCTGCATGCCATAGCAAACGCCGATTACCGGTACGCCTGCTTCAAAGCCATACTGCGGCGCAcg   <  2:864587‑M2/66‑1 (MQ=255)                                                                  tGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCAcgc  <  1:641234‑M2/67‑1 (MQ=255)                                                                   |                                                                   AACCTGCGCGTAGCCAAATTCACGTTCGTTAGAGGCTTCAACGTGACCGCCCAACTGCATTGCCATGGTCTGCATGCCATAGCAAACGCCGAATACCGGTACGCCTGCTTCAAAGACATACTGCGGCGCACGC  >  NC_000913/2632197‑2632329

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