New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_pgi	3570872 =	14 (0.860)	4 (0.260)	4/240	1.6	27.3%	coding (448/1974 nt)	glgX	glycogen debranching enzyme
	?	NC_000913_3_pae_pgi	3570898 =	8 (0.520)					coding (422/1974 nt)	glgX	glycogen debranching enzyme

TTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGC  >  NC_000913_3_pae_pgi/3570790‑3571026                                                                                                             |                                                                                                                                 ttATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTCCCCAGCCATAGTGATCAACCACCACTACGCATTTc                                                                                         >  2:132234/1‑150 (MQ=255)  taatgaacctgactatcgcgacaacgccgccattgcgccgaaatgcgtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                 >  3:9859‑M1/82‑109 (MQ=255)  taatgaacctgactatcgcgacaacgccgccattgcgccgaaatgcgtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                                                 <  4:9859‑M1/28‑1 (MQ=255)                   tcGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATACATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGTGGTTg                                                                       >  2:188449/1‑150 (MQ=255)                      aCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTcg                                                                    >  1:145283/1‑150 (MQ=255)                                 cGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGCGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAgg                                                         >  1:134848/1‑150 (MQ=255)                                                      tttGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAac                                                                                                         <  1:189322/81‑1 (MQ=255)                                                      tttGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAac                                                                                                         >  2:189322/1‑81 (MQ=255)                                                                       aaaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGgcg                                                                                     >  1:52702/1‑84 (MQ=255)                                                                       aaaTGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGgcg                                                                                     <  2:52702/84‑1 (MQ=255)                                                                                 gatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTa                                                  >  3:9859‑M2/29‑109 (MQ=255)                                                                                 gatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTa                                                  <  4:9859‑M2/81‑1 (MQ=255)                                                                                       ccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTAc                                                                                                <  1:23303‑M2/35‑1 (MQ=255)                                                                                       ccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTAc                                                                                                >  2:23303‑M2/23‑57 (MQ=255)                                                                                        CACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGc  <  3:14327‑M1/150‑1 (MQ=255)                                                                                                             |                                                                                                                                 TTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGGGTTATCTTTAAACTCCCCGTCAATTTGC  >  NC_000913_3_pae_pgi/3570790‑3571026

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