New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_pae_tpiA	3570872 =	17 (1.090)	4 (0.270)	3/236	1.7	20.9%	coding (448/1974 nt)	glgX	glycogen debranching enzyme
	?	NC_000913_3_pae_tpiA	3570898 =	14 (0.950)					coding (422/1974 nt)	glgX	glycogen debranching enzyme

AAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGG  >  NC_000913_3_pae_tpiA/3570757‑3570998                                                                                                                                              |                                                                                                     aaaTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTcc                                                                                              >  1:176181/1‑150 (MQ=255)   aTAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCa                                                                                            <  3:30053/150‑1 (MQ=255)        tGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCAt                                                                                       <  1:125505/150‑1 (MQ=255)         gATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATa                                                                                      <  2:176181/150‑1 (MQ=255)                                tcataatgaacctgactatcgcgacaacgccgccattgcgccgaaatgcgtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGgcatcag                                                                                                   >  4:95132‑M1/85‑113 (MQ=255)                                tcataatgaacctgactatcgcgacaacgccgccattgcgccgaaatgcgtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGGcatcat                                                                                                   <  3:95132‑M1/30‑1 (MQ=255)                                                gATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATggcggc                                               <  2:75440/150‑1 (MQ=255)                                                    tcGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCGTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTg                                           >  4:64545/1‑150 (MQ=255)                                                      gACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATggcggc                                               >  2:98322/1‑144 (MQ=255)                                                      gACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATggcggc                                               <  1:98322/144‑1 (MQ=255)                                                              tctcCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGt                                 <  4:87132/150‑1 (MQ=255)                                                                             aCGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGAc                  <  3:64545/150‑1 (MQ=255)                                                                                 gtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGGCatcatc                                                                                                  >  3:80940‑M1/36‑66 (MQ=255)                                                                                 gtagtggtggttgatcactatgactgggaagatgaTGCCCCACGGCGTGCGCGGCGGGGCatcatc                                                                                                  <  4:80940‑M1/31‑1 (MQ=255)                                                                                  aaTCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                     <  2:187795/62‑1 (MQ=255)                                                                                  aaTCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGcatc                                                                                                     >  1:187795/1‑62 (MQ=255)                                                                                                      aTAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCgg  >  1:174111/1‑141 (MQ=255)                                                                                                      aTAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCgg  <  2:174111/141‑1 (MQ=255)                                                                                                               gatgatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTAc                                                                    <  4:80940‑M2/35‑1 (MQ=255)                                                                                                               gatgatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTAc                                                                    >  3:80940‑M2/32‑66 (MQ=255)                                                                                                                ctgatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGa                   >  4:95132‑M2/31‑114 (MQ=255)                                                                                                                atgatgccccgccgcgcacgccgtggggcaTCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGa                   <  3:95132‑M2/84‑1 (MQ=255)                                                                                                                                              |                                                                                                     AAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGG  >  NC_000913_3_pae_tpiA/3570757‑3570998

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