New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	1950990 =	8 (0.290)	3 (0.120)	3/196	2.3	29.5%	coding (1259/1563 nt)	pntA4	NAD(P) transhydrogenase subunit alpha
	?	AM260480	1951044 =	7 (0.280)					coding (1205/1563 nt)	pntA4	NAD(P) transhydrogenase subunit alpha

CAAAGATCGCCAGCACGAACACCGTGAAATGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGCTGCGCGGCGGCTCGGCCTCGGCCTGCACCGGCGGCGGCGCGGCCGGCGCCTGCTGCTGCGGGATCGCCACCGTCAGCGGCGGCGGTGGCCAGGTCACCGTGCCCTGGTGCAGCACCGTGGCACCGCGGATCATCTCGTCGTC  >  AM260480/1950913‑1951174                                                                                                                                    |                                                                                                                                   cAAAGATCGCCAGCACGAACACCGTGAAATGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGCTGCGCGGCGGctcg                                                                                                                                  <  1:500933/134‑1 (MQ=255)                      gagctcccgcagcagcaggcgccggccgcgccgccgccggtgcgggccgaggccgagccgccgcggagccgcaccgacgcggcgcggatggcgctggccaccgcggcggcTCGGCCTCGGCCTGCACcggcggc                                                                                                             <  1:573537‑M2/24‑1 (MQ=255)                           aaaTGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGCTGCGCGGCGGCTcggcctcggcct                                                                                                                        >  1:380011/1‑118 (MQ=255)                           aaaTGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGCTGCGCGGCGGCTcggcctcggcct                                                                                                                        <  2:380011/118‑1 (MQ=255)                            aaTGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGGTGCGCGGCGGCTCGGCCTCGGCCTGCACCGGCGgcggcgcg                                                                                                       <  2:37259/134‑1 (MQ=255)                                                      gccgccggtgcaggccgaggccgaGCCGCCGCGGTGGCCAGCGCCATCCGCGCCGCGTCGGTGCGGCTCCGCGGCGGCTCGGCCTCGGCCCGCACCGGCGGCGGCGCGGCCGGCGCCTGCTGCTGcgggagctc                                                                             <  1:573537‑M1/110‑5 (MQ=255)                                                                                   cgcgcagccgcaccggcgcgacgctgatcgcgctggccatcgcggcggcTCGGCCTCGGCCTGCACCGGCGgcggcgcg                                                                                                       <  1:57059‑M2/30‑1 (MQ=255)                                                                                   cgcgcagccgcaccggcgcgacgctgatcgcgctggccatcgcggcggcTCGGCCTCGGCCTGCACCGGCGgcggcgcg                                                                                                       >  2:57059‑M2/50‑79 (MQ=255)                                                                                       ctGGCCAGCGCTATCCGCGTCGCGCCGGTGCCGGTGCGCGGGGGCTCGGCGTCGGCGTGCACCGGCGGCGGCGCGGCCGGCGCCTGCTGCTGCGGGATCGCCACCGTCAGCGGCGGCGGTGGCCAGGTCACCGt                                            <  2:573537/133‑1 (MQ=255)                                                                                                                 ggTGCGGCTGCGCGGCGGCTCGGCCTCGGCCTGCACCGGCGGCGGCGCGGCCGGCGCctgctgct                                                                                       <  1:970451/65‑1 (MQ=255)                                                                                                                 ggTGCGGCTGCGCGGCGGCTCGGCCTCGGCCTGCACCGGCGGCGGCGCGGCCGGCGCctgctgct                                                                                       >  2:970451/1‑65 (MQ=255)                                                                                                                                 ggcTCGGCCTCGGCCTGCACCGGCGGCGGCGCGGCCGGCGCCTGCTGCTGCGGGATCGCCACCGTCAGCGGCGGCGGTGGCCAGGTCACCGTGCCCTGGTGCAGCACCGTGGCACCGCGGATCATCTCgtcgtc  >  1:313706‑M2/4‑134 (MQ=255)                                                                                                                                    |                                                                                                                                   CAAAGATCGCCAGCACGAACACCGTGAAATGCGCCATGAAGGCCGGCGGCGCGACCGCGCCCAGGCCCAGCAGCGCCGCCGCCGCGATGGCCAGCGCGATCAGCGTCGCGCCGGTGCGGCTGCGCGGCGGCTCGGCCTCGGCCTGCACCGGCGGCGGCGCGGCCGGCGCCTGCTGCTGCGGGATCGCCACCGTCAGCGGCGGCGGTGGCCAGGTCACCGTGCCCTGGTGCAGCACCGTGGCACCGCGGATCATCTCGTCGTC  >  AM260480/1950913‑1951174

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