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)

CTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACT  >  NC_000913/2631969‑2632098                                                                    |                                                               gcgtgcgccgcagtatgtctttgaagcaggcgtaccggtattcggcgtttgctatggcatgcagacca                                                                <  1:641234‑M1/1‑1 (MQ=255)  cgtgcgccgcagtatggctttgaagcaggcgtaccggtaatcggcgtttgctatggcatgcagcccAt                                                               <  2:864587‑M1/2‑1 (MQ=255)   gtgcgccgcagtatggctttgaagcaggcgtaccggtattcggcgtttgctatggcatgcagaccATg                                                              >  2:99941‑M1/66‑68 (MQ=255)      cgccgcagtatgtctttgaagcaggcgtaccggtattcggcgtttgctatggcatgcagaccATGGCa                                                           >  2:96775‑M1/63‑68 (MQ=255)          aTGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACg                                                       >  2:20715/1‑68 (MQ=255)              gTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGcg                                                   >  1:621845/1‑67 (MQ=255)               atgtctttgaagcaggcgtaccggtattcggcgtttgctatggcatgcagaccATGGCAAACGGGCAg                                                  <  2:229566‑M1/15‑1 (MQ=255)                   ttCCGGGTGGAACTGTACGCAATAGAAGCGATTTTCTTCGTTAGCCATAATGGCAAACGGGAAGCttt                                              >  1:484343/1‑68 (MQ=255)                       gggTGGAACTGTACGCCATAGAAGCGCCTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCggc                                          <  2:489033/68‑2 (MQ=255)                         gTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAACCGGGCAGCATTCGGTGc                                        >  2:617288/1‑68 (MQ=255)                           ggAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTg                                      <  1:50164/68‑1 (MQ=255)                           ggAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTg                                      <  1:950395/68‑1 (MQ=255)                            gAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTgg                                     >  2:303847/1‑68 (MQ=255)                            gAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTgg                                     <  2:387340/68‑1 (MQ=255)                                      gtattcggcgtttgctatggcatgcagaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGtgat                           <  1:656925‑M1/38‑1 (MQ=255)                                        aTAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGtgatga                         >  1:854542/1‑68 (MQ=255)                                             aGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcg                    >  1:493607/1‑68 (MQ=255)                                             aGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcg                    >  1:386585/1‑68 (MQ=255)                                              gCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcgg                   <  1:303847/68‑1 (MQ=255)                                               cGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcgga                  >  1:44184/1‑68 (MQ=255)                                                 ttgctatggcatgcagaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                >  2:656925‑M1/20‑68 (MQ=255)                                                 ttgctatggcatgcagaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                >  1:229566‑M1/20‑68 (MQ=255)                                                 ttgctatggcatgcagaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                >  2:369706‑M1/20‑68 (MQ=255)                                                       ttcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAg          >  1:759350/1‑68 (MQ=255)                                                       ttcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAg          <  2:654070/68‑1 (MQ=255)                                                        tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc         <  2:386585/68‑1 (MQ=255)                                                        tcGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGc         <  1:517881/68‑1 (MQ=255)                                                               ccATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACt  <  1:881334/68‑1 (MQ=255)                                                               ccATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACt  >  1:389079/1‑68 (MQ=255)                                                               ccATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACt  >  1:18456/1‑68 (MQ=255)                                                                gaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                >  2:417767‑M1/5‑53 (MQ=255)                                                                gaccATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTcggacg                <  1:417767‑M1/49‑1 (MQ=255)                                                                    |                                                               CTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACT  >  NC_000913/2631969‑2632098

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.