New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4349069 =	29 (1.250)	7 (0.310)	7/242	NT	24.1%	intergenic (‑390/+84)	PP_3820/galU	group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
	?	NC_002947	4393427 =	16 (0.720)					intergenic (‑369/+101)	PP_3868/PP_3869	group II intron‑encoding maturase/sheath protein

TTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAT  >  NC_002947/4348920‑4349209                                                                                                                                                      |                                                                                                                                             gtgcttgccctgcaataccagccccgcctgcagtacgcggcatcgtaactgaatcatccaccaaggcccgccaagtgcgggcctttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccgcttggttgtggtggccaagtg                                                                                                                                              <  1:118595‑M1/1‑1 (MQ=255)   gcttgccctgcaataccagccccgcctgcagtacgcggcatcgtaactgaatcatccaccaaggcccgccaagtgcgggcctttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccgcttggttgtggtggccaagtggg                                                                                                                                            >  2:118595‑M1/148‑150 (MQ=255)     tGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGtt                                                                                                                                          >  2:327418/1‑150 (MQ=32)          gTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCt                                                                                                                                     >  2:168155/1‑150 (MQ=32)               gCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCTGGTTGCGCTTAcgc                                                                                                                                <  1:178205/150‑1 (MQ=17)                ccAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTAcgcg                                                                                                                               >  2:354295/1‑150 (MQ=32)                  aGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCa                                                                                                                             >  1:49726/1‑150 (MQ=32)                        ttGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACgcgc                                                                                                                       <  2:306642/150‑1 (MQ=32)                                  cttcattcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGgc                                                                                                             >  2:57833/1‑150 (MQ=17)                                       ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGcaca                                                                                                        >  1:531474/1‑150 (MQ=32)                                          aGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCCCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACaa                                                                                                     >  1:296451/1‑150 (MQ=17)                                              tCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACaaaaaa                                                                                                 <  1:57833/150‑1 (MQ=32)                                                        ccGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGCTTCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGa                                                                                       >  2:575881/1‑150 (MQ=255)                                                          gcgACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAg                                                                                     <  1:422145/150‑1 (MQ=40)                                                                cccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCCCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGc                                                                               <  2:296451/150‑1 (MQ=255)                                                                   ttGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCttt                                                                            <  2:136920/150‑1 (MQ=255)                                                                     tccTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCCCCTCCAAGTCACCCACTTGGCCACCACACCCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                          <  2:114805/149‑1 (MQ=255)                                                                     gCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                          <  2:93727/150‑1 (MQ=255)                                                                       cTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAgg                                                                        >  2:158757/1‑150 (MQ=255)                                                                       cTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAgg                                                                        >  1:377714/1‑150 (MQ=255)                                                                                    ttCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTa                                                           <  1:306234/150‑1 (MQ=255)                                                                                     tCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTAt                                                          <  1:158757/150‑1 (MQ=255)                                                                                      ccccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATc                                                         <  2:488686/150‑1 (MQ=255)                                                                                        ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCcagcag                                                                 <  1:134153/140‑1 (MQ=255)                                                                                        ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCcagcag                                                                 >  2:134153/1‑140 (MQ=255)                                                                                                ggTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGACACCACCAACACGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCtt                                               <  2:377714/150‑1 (MQ=255)                                                                                                 gtgtgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                          >  2:535353/1‑122 (MQ=255)                                                                                                 gtgtgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                          <  1:535353/122‑1 (MQ=255)                                                                                                           gACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTGGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgt                                    <  2:282580/150‑1 (MQ=255)                                                                                                            ccTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATCCCTCGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgta                                   <  2:425571/149‑1 (MQ=255)                                                                                                            aCTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgta                                   <  1:354295/150‑1 (MQ=255)                                                                                                                        aaGTCACCCACTTGGCCACCACAACCAAGCGGGTTTCGCGGACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                       >  1:50051/1‑150 (MQ=255)                                                                                                                            aacccgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                  <  1:491393‑M1/61‑1 (MQ=255)                                                                                                                            aacccgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                  >  2:491393‑M1/27‑87 (MQ=255)                                                                                                                              cccgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACACCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATGAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGaa                 <  1:348967‑M1/126‑1 (MQ=255)                                                                                                                               ccgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGACTTATCAGTAAGCCTTTCCAGTCTTGTAGTAGTTCTCGAAGCAgaaa                <  2:101717‑M1/127‑2 (MQ=255)                                                                                                                                cgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCAt                                                                                                                    <  2:144429‑M1/27‑1 (MQ=255)                                                                                                                                cgcttggttgtggtggccaagtGGGTTGCGCTTACGCGCAACGCGCCAt                                                                                                                    >  1:144429‑M1/23‑49 (MQ=255)                                                                                                                                             tggccaagtGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAt  >  1:101717‑M1/10‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                             TTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAT  >  NC_002947/4348920‑4349209

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