New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	= 3870594	9 (0.720)	6 (0.510)	6/246	0.7	42.1%	coding (1137/1293 nt)	dgoT	D‑galactonate transporter
	?	NC_000913_3_pae_tpiA	= 3870607	8 (0.680)					coding (1124/1293 nt)	dgoT	D‑galactonate transporter

GTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGCACATAATCATCATCGGATCGTTAGTGT  >  NC_000913_3_pae_tpiA/3870524‑3870733                                                                                    |                                                                                                                               gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCcgccgg                                                                                                                          >  3:27339‑M2/1‑84 (MQ=255) gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCcgccgg                                                                                                                          <  4:27339‑M2/90‑7 (MQ=255) gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGc                                                                <  1:15550/148‑1 (MQ=255) gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGc                                                                >  2:15550/1‑148 (MQ=255) gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGACGCCGCCCAGACCACCGGCGAAGTTAAACACGCcgccggtggtctgg                                                                                                                  >  3:16010‑M2/1‑84 (MQ=255) gTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGACGCCGCCCAGACCACCGGCGAAGTTAAACACGCcgccggtggtctgg                                                                                                                  <  4:16010‑M2/98‑15 (MQ=255)                                caccaccaGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAg                               >  1:58760/1‑150 (MQ=255)                                caccaccaGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGAGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAg                               >  4:65333/1‑150 (MQ=255)                                   caccaGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGca                            <  2:58760/150‑1 (MQ=255)                                    accaGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGt                                                        >  1:39674/1‑121 (MQ=255)                                    accaGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGt                                                        <  2:39674/121‑1 (MQ=255)                                       aGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCcgccggtggtctgggcggcatcaccgttccgctggtggtggggtacctggcgcagggttacggtttcgcacctgcactggtttatatctccgccgtcgcgttga                        <  1:98142‑M2/150‑105 (MQ=255)                                           gAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGCACATAatca                    <  3:64576/150‑1 (MQ=255)                                                             gACCACCGGCGAAGTTAAACACGCcgccggtggtctgggcggcatcaccgttccgctggtggtggggtacctggcgcagggttacggtttcgcacctgcactggtttatatctccgccgtcgcgttgattggcgcgctctcttatatcct  <  4:13537‑M2/150‑127 (MQ=255)                                                                caccGGCGAAGTTAAACACGCcgccggtggtctgggcggcatcaccgttccgctggtggtggggtacctggcgcagggttacggttt                                                              <  1:88219‑M2/87‑67 (MQ=255)                                                                caccGGCGAAGTTAAACACGCcgccggtggtctgggcggcatcaccgttccgctggtggtggggtacctggcgcagggttacggttt                                                              >  2:88219‑M2/1‑21 (MQ=255)                                                                                    |                                                                                                                               GTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGCACATAATCATCATCGGATCGTTAGTGT  >  NC_000913_3_pae_tpiA/3870524‑3870733

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