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.

ACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCC  >  NC_000913/3570931‑3571225                                                                                                                                                      |                                                                                                                                                  aCGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTcggcgg                                                                                                                                                  <  1:145916/151‑1 (MQ=255)              cAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGc                                                                                                                                      <  2:226200/150‑1 (MQ=255)              cAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGc                                                                                                                                      <  1:53081/150‑1 (MQ=255)                    cggcgTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGaa                                                                                                                                <  1:175297/150‑1 (MQ=255)                    cggcgTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGaa                                                                                                                                <  2:44513/150‑1 (MQ=255)                    cggcgTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGaa                                                                                                                                <  1:38128/150‑1 (MQ=255)                         ttGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGAt                                                                                                                          <  2:159645/151‑1 (MQ=255)                               cgATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCa                                                                                                                     <  1:286517/150‑1 (MQ=255)                                            cATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCgg                                                                                                       >  2:213880/1‑151 (MQ=255)                                                       gCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCg                                                                                            >  1:169235/1‑151 (MQ=255)                                                          tGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGc                                                                                          <  1:44317/150‑1 (MQ=255)                                                                          cTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGcgcattacgacggtcagggcgtcaacttcacacttttctccgctcatgccgagcgggtagaactgtgtgtcttt                                                                         >  2:87850‑M1/1‑77 (MQ=255)                                                                                              tGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCa                                                     <  1:263805/151‑1 (MQ=255)                                                                                                   cgcgcAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCa                                                <  2:289056/151‑1 (MQ=255)                                                                                                          ggATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGcgcattacgacggtcagggcgtcaacttcacacttttctccgctcatgccgagcgggtagaactgtgtgtctttgacgccaatggccaggaacatcgctatgactt                                         >  2:309515‑M1/1‑45 (MQ=255)                                                                                                                atcaaCAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGcgcattacgacggtcagggcgtcaacttcacacttttctccgctcatgccgagcgggtagaactgtgtgtctttgacgccaatggccaggaacatcgctatgacttgccagg                                   >  1:277706‑M1/1‑39 (MQ=255)                                                                                                                atcaaCAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTc                                     >  2:73560/1‑149 (MQ=255)                                                                                                                    acaacTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGcc                               <  2:198141/151‑1 (MQ=255)                                                                                                                               gggTTAAAGCGATGCCCCTCGGCGcgcattacgacggtcagggcgtcaacttcacacttttctccgctcatgccgagcgggtagaactgtgtgtctttgacgccaatggccaggaacatcgctatgacttgccagggcacagtggcgacat                    >  1:72730‑M1/1‑24 (MQ=255)                                                                                                                                        cGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGAcacaca           >  1:288662/1‑151 (MQ=255)                                                                                                                                                cTCGGCGcgcattacgacggtcagggcgtcaacttcacacttttctccgctcatgccgagcgggtagaactgtgtgtctttgacgccaatggccaggaacatcgctatgacttgccagggcacagtggcgacatttggcacggttatctgc   <  2:1120‑M1/151‑145 (MQ=255)                                                                                                                                                 tCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTAccc  <  1:132483/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGCCGCGCGCAAGGATCAATCAACAACTTCGCCGGGTTAAAGCGATGCCCCTCGGCGGGTTGCCAGGGGCCATGAACGCGATAACCATAACGCAAACCCGGGCGCGCATCCGGCAGATAACCGTGCCAAATGTCGCCACTGTGCCCTGGCAAGTCATAGCGATGTTCCTGGCCATTGGCGTCAAAGACACACAGTTCTACCC  >  NC_000913/3570931‑3571225

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