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

TAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGG  >  NC_000913_3_pae_tpiA/3570740‑3570998                                                                                                                                     |                                                                                                                               tAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCg                                                                                                                     <  1:189896/144‑1 (MQ=255) tAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCg                                                                                                                     >  2:189896/1‑144 (MQ=255)     gCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGcggcgg                                                                                                           <  4:2193/150‑1 (MQ=255)                  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  <  2:174111/141‑1 (MQ=255)                                                                                                                       aTAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCgg  >  1:174111/1‑141 (MQ=255)                                                                                                                                     |                                                                                                                               TAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTGCCCCACGGCGTGCGCGGCGGGGCATCATCTTCCCAGTCATAGTGATCAACCACCACTACGCATTTCGGCGCAATGGCGGCGTTGTCGCGATAGTCAGGTTCATTATGACCGGCGTGCAGCAGCGG  >  NC_000913_3_pae_tpiA/3570740‑3570998

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