New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3976787 =	NA (NA)	20 (0.980)	14/208	NT	100%	intergenic (+177/‑203)	PP_3501/PP_3502	transposase/transposase
	?	NC_002947	4578248 =	0 (0.000)					intergenic (+63/‑74)	PP_4054/glgX	hypothetical protein/glycogen debranching enzyme

AAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCC  >  NC_002947/3976656‑3976888                                                                                                                                    |                                                                                                      aacgcagcTCGACAGGGTGGGCGAGTCCATACATCGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTAt                                                                                                   <  2:316318/133‑1 (MQ=255)      agcTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAg                                                                                              <  2:19417/136‑1 (MQ=255)           gACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACtt                                                                                         <  1:166224/136‑1 (MQ=255)                   gggCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag                                                                                 >  1:282058/1‑136 (MQ=255)                   gggCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag                                                                                 >  1:3259/1‑136 (MQ=255)                   gggCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggag                                                                                 >  1:506820/1‑136 (MQ=255)                   gggCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTAGAAGggat                                                                                 >  2:654740/1‑135 (MQ=255)                    gtcGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagg                                                                                <  2:310220/134‑1 (MQ=255)                    ggCTAGTCATTACTTAGGAGGGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagg                                                                                <  1:342519/136‑1 (MQ=21)                    ggCGAGTCCATACATAGGTGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGggagg                                                                                <  2:592084/136‑1 (MQ=255)                             atacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGcc                                                                       >  2:319724/1‑136 (MQ=255)                             atacataGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGcc                                                                       >  2:375824/1‑136 (MQ=255)                               acataGGTCCGGGTCAACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACTTAGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGAGTTGAAGGGAGGAGGAGGGCCAt                                                                     <  2:580631/136‑1 (MQ=25)                                      aGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGc                                                              >  2:30383/1‑136 (MQ=255)                                           gTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTg                                                         >  1:223267/1‑136 (MQ=255)                                             ttACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTAt                                                                                                   <  1:56331/92‑1 (MQ=255)                                             ttACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTAt                                                                                                   >  2:56331/1‑92 (MQ=255)                                                  ggcgAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTc                                                  <  1:636670/135‑1 (MQ=255)                                                  cgcgAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTc                                                  <  1:30383/136‑1 (MQ=255)                                                       aTACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCa                                             <  2:238370/136‑1 (MQ=255)                                                         aCGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGt                                           <  1:460279/136‑1 (MQ=255)                                                         aCGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGt                                           <  1:319724/136‑1 (MQ=255)                                                                ccGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGCCTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGccc                                    >  2:553546/1‑136 (MQ=255)                                                                        aTGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAggtta                            <  2:384123/136‑4 (MQ=255)                                                                        aTGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGgaga                            <  1:177230/136‑1 (MQ=255)                                                                                cGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCAccc                    <  1:667984/136‑1 (MQ=255)                                                                                cGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCAccc                    <  1:283525/136‑1 (MQ=255)                                                                                          aCATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCttt          <  2:7746/136‑1 (MQ=255)                                                                                                  tagtagGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACAccc  <  2:493001/136‑1 (MQ=40)                                                                                                  tagtagGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACAccc  <  1:351197/136‑1 (MQ=40)                                                                                                                                    |                                                                                                      AAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGTTTACCCGCGAATACGGTAGCCGCGACAATGGTGAACGGCAGGTGGACATTGGCTAGTAGGCCCGGCTCTTTCGCGGGTAAACCCGCTTGTATGGTAGGACTTGAAGGGAGGAGGAGGGCCATAACTCGCTTCTGACTGTTCGCGCAGTACAGACCGTAGGAGAAATCACCCCTCCTCCTTTCAACACCC  >  NC_002947/3976656‑3976888

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