New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3571080	17 (0.840)	4 (0.210)	4/278	1.7	17.8%	coding (240/1974 nt)	glgX	glycogen debranching enzyme
	?	NC_000913	= 3571279	21 (1.110)					coding (41/1974 nt)	glgX	glycogen debranching enzyme
Rejected: Frequency below cutoff threshold.

GCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGTGCCGCCA  >  NC_000913/3571136‑3571424                                                                                                                                                |                                                                                                                                                  gCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGcgccgag                                                                                                                                            <  2:1120‑M2/151‑8 (MQ=255)       aaCCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGc                                                                                                                                      <  2:7646/151‑1 (MQ=255)          cGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCggg                                                                                                                                   <  2:169234/151‑1 (MQ=255)                  aTGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGcgccgaggggcatcgctttaaccc                                                                                                                           >  1:72730‑M2/1‑127 (MQ=255)                          aCTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGt                                                                                                                   <  2:64672/151‑1 (MQ=255)                              tgCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCCAggaac                                                                                                                <  1:73560/150‑5 (MQ=255)                               gCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTgtgt                                                                                                              >  2:228710/1‑151 (MQ=255)                                 ccTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGcgccgaggggcatcgctttaacccggcgaagttgttgat                                                                                                            >  1:277706‑M2/1‑112 (MQ=255)                                       aaGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGcgccgaggggcatcgctttaacccggcgaagttgttgattgatcc                                                                                                      >  2:309515‑M2/1‑106 (MQ=255)                                        aGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGc                                                                                                      >  2:239819/1‑150 (MQ=255)                                                cGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAc                                                                                             <  2:141019/151‑1 (MQ=255)                                                     ttCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAAccagc                                                                                         >  1:211046/1‑150 (MQ=255)                                                                       aaaGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGcgccgaggggcatcgctttaacccggcgaagttgttgattgatccttgcgcgcggcaaattgacggggagtttaaag                                                                      >  2:87850‑M2/1‑74 (MQ=255)                                                                           acacacAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCggtggt                                                                    <  2:92502/149‑1 (MQ=255)                                                                           acacacAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAg                                                                  <  1:184893/151‑1 (MQ=255)                                                                                         ccGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTg                                                    <  2:286259/151‑1 (MQ=255)                                                                                                            aaaaGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGtggc                                  <  2:250298/150‑1 (MQ=255)                                                                                                                       gTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCAt                      >  2:59553/1‑151 (MQ=255)                                                                                                                              cccTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCtgtg               >  2:66782/1‑151 (MQ=255)                                                                                                                                 tGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTAc            <  2:299490/151‑1 (MQ=255)                                                                                                                                 tGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTAc            <  2:146927/151‑1 (MQ=255)                                                                                                                                 tGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTAc            <  1:189979/151‑1 (MQ=255)                                                                                                                                    ccGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGt         <  2:264793/151‑1 (MQ=255)                                                                                                                                    ccGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGt         <  2:304687/151‑1 (MQ=255)                                                                                                                                       tcgtAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGTgc       <  2:165735/150‑1 (MQ=255)                                                                                                                                           aaTGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGTGCCGCCa  <  2:214493/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  GCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCCGCTCGGCATGAGCGGAGAAAAGTGTGAAGTTGACGCCCTGACCGTCGTAATGCGCGCCGAGGGGAGCGGGTTTGCCAATGGCGAGTTGTGTCATTCTGCCTCCCGAACCAGCCAGATAGTGGCCAGCGGTGGTAGCGTCAGGCTTAGTGAATGCTGACGACCGTGGCTGGCAATCTCATCGCTGTGTACCGTGCCGCCA  >  NC_000913/3571136‑3571424

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