New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AM260480	= 1755368	9 (0.320)	3 (0.110)	3/196	2.4	25.7%	coding (1181/2235 nt)	pulA	Type II secretory pathway, pullulanase PulA
	?	AM260480	= 1755387	9 (0.340)					coding (1200/2235 nt)	pulA	Type II secretory pathway, pullulanase PulA

GCGCGAGGGCAACGGCTTCGATGCGGGCTCGGGCTTTTTCGATGCACTGATGCAGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGGCGCGGCGACGCCGGCCACCGCGGCGAACTGGCGGCGCGGCTGGCGGCCTCGGCCGACCTCTTCGACCGCCGC  >  AM260480/1755255‑1755516                                                                                                                                     |                                                                                                                                  gcgcGAGGGCAACGGCTTCGATGCGGGCTCGGGCTTTTTCGATGCACTGATGCAGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACc                                                                                                                                  <  1:787596‑M2/134‑2 (MQ=255)        ggCAACGGCTTCGATGCGGGCTCGGGCTTTTTCGATGCACTGATGCAGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACccgcccgg                                                                                                                           <  2:137303‑M2/134‑9 (MQ=255)                                        cGATGCACTGATGCAGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGAc                                                                                            <  1:89600/133‑1 (MQ=255)                                                      aGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTc                                                                             >  2:225657/1‑134 (MQ=255)                                                       ggATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCt                                                                            >  2:464646/1‑134 (MQ=255)                                                          tCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGgc                                                                         <  2:461971/134‑1 (MQ=255)                                                                                              gTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCttt                                                                                                                       <  1:219351/52‑1 (MQ=255)                                                                                              gTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCttt                                                                                                                       >  2:219351/1‑52 (MQ=255)                                                                                                            cGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGGCGCGGCGACGCCGGCCACCGCGGCGAACTGGCGGCGCGGCTGGCGGCCTCg                       <  1:464646/134‑1 (MQ=255)                                                                                                             gggCGGCTACCAGCTTGGCAACCACccgcccgggcccaggtcccacggctcggagatcagcttgacgcggctcagcaccggatcctgcatcag                                                               >  2:51298‑M2/1‑25 (MQ=255)                                                                                                             gggCGGCTACCAGCTTGGCAACCAAccgcccgggcccaggtcccacggctcggagatcagcttgacgcggctcagcaccggatcctgcatcag                                                               <  1:51298‑M2/93‑69 (MQ=255)                                                                                                                             ggCAACCACccgcccgggcccaggtcccacggctcggagatcagcttgacgcggctcagcaccggatcctgcatcagtgcatcgaaaaagcccgagcccgcatcgaagccgttgccctcgcgcccgagggtggt      <  1:684464‑M2/134‑126 (MQ=255)                                                                                                                                aaCCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGGCGCGGCGACGCCGGCCACCGCGGCGAACTGGCGGCGCGGCTGGCGGCCTCGGCCGACCTCTTCGAccgccg   <  1:225657/134‑1 (MQ=255)                                                                                                                                  ccaccCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGGCGCGGCGACGCCGGCCACCGCGGCGAACTGGCGGCGCGGCTGGCGGCCTCGGCCGACCTCTTCGACcgccgc  >  2:607400/1‑133 (MQ=255)                                                                                                                                     |                                                                                                                                  GCGCGAGGGCAACGGCTTCGATGCGGGCTCGGGCTTTTTCGATGCACTGATGCAGGATCCGGTGCTGAGCCGCGTCAAGCTGATCTCCGAGCCGTGGGACCTGGGCCCGGGCGGCTACCAGCTTGGCAACCACCCGCCCGGCTTTGCCGAGTGGAACGACAAGTTCCGCGACGTGGCGCGCCGCTTCTGGCGCGGCGACGCCGGCCACCGCGGCGAACTGGCGGCGCGGCTGGCGGCCTCGGCCGACCTCTTCGACCGCCGC  >  AM260480/1755255‑1755516

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