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

CGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAGGGCAAGCAC  >  NC_002947/4393286‑4393575                                                                                                                                              |                                                                                                                                                     atcgacgcgaccaacttctgcttcgagaactactacaagactggcaaggcttactgataagccctgctggcctcaaagccatccttcggggtggcttttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccACTTGGcca                                                                                                                                              >  1:101717‑M2/142‑150 (MQ=255)               tttctgcttcgagaactactacaagactggaaaggcttactgataagtcctgctggcctcaaagccatccttcggggtggcttttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCgg                                                                                                                                <  2:101717‑M2/23‑1 (MQ=255)                ttctgcttcgagaactactacaagactggcaaggcttactcataagccctgctggcctcaaagccatccttcggtgtggcttttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCggg                                                                                                                               <  1:348967‑M2/24‑1 (MQ=255)                                                             gcgGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGc                                                                                  <  2:441178/150‑1 (MQ=31)                                                                    cccGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACt                                                                           >  1:355542/1‑150 (MQ=33)                                                                            acCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGc                                                                   <  2:167903/150‑1 (MQ=39)                                                                              ccTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCt                                                                 >  2:142543/1‑150 (MQ=40)                                                                                 atccttcggggtggcttttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCGGGtt                                                                                                                             >  2:491393‑M2/62‑87 (MQ=255)                                                                                 atccttcggggtggcttttttgtgtgcgccaggcatggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCGGGtt                                                                                                                             <  1:491393‑M2/26‑1 (MQ=255)                                                                                  gCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTtggt                                                             <  2:463525/150‑1 (MQ=255)                                                                                   ccTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTtggtg                                                            >  2:222347/1‑150 (MQ=255)                                                                                    cTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTtggtgg                                                           >  1:249248/1‑150 (MQ=255)                                                                                        ggCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGgatga                                                       <  2:249248/150‑1 (MQ=255)                                                                                          cTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGAtt                                                     <  1:530343/150‑1 (MQ=255)                                                                                               acTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGtt                                                <  1:222347/150‑1 (MQ=255)                                                                                                cTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTa                                               >  1:88052/1‑150 (MQ=255)                                                                                                                   atggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCg                                                                                                                                 <  2:144429‑M2/22‑1 (MQ=255)                                                                                                                   atggcgcgttgcgcgtaagcgcaacccACTTGGCCACCACAACCAAGCg                                                                                                                                 >  1:144429‑M2/28‑49 (MQ=255)                                                                                                                     gagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGcc                                        >  2:120856/1‑136 (MQ=255)                                                                                                                     gagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGcc                                        <  1:120856/136‑1 (MQ=255)                                                                                                                     gagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCg                          >  2:43287/1‑150 (MQ=255)                                                                                                                           tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTg                    <  2:88052/150‑1 (MQ=255)                                                                                                                                  tCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGc             <  2:137516/150‑1 (MQ=255)                                                                                                                                     aaGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTCACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAgg          <  2:355542/150‑1 (MQ=255)                                                                                                                                         cacacaCGTGGTCACCACTCCCAACCGGTTTGCGCTTCCGCCCAATGCGCCATGCCTGTCACACACAAAAAGTCCCGCACTTGTCGGTGCTGGGTTGATGATTCACTTACGATCCCGCGTACTGCATGCGGGGCTGGTATTGCAGGGCaa      <  2:202819/150‑1 (MQ=255)                                                                                                                                         cACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAGGGCaa      <  1:360573/150‑1 (MQ=255)                                                                                                                                           cccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAGGGCAAGc    >  2:118595‑M2/4‑150 (MQ=255)                                                                                                                                             cACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAGGGCAAGCAc  <  1:118595‑M2/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                     CGCCTCTATCCGCTTCCTGTTCGTCAGGCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAACACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAACGCGCCATGCCTGGCGCACACAAAAAGGCCCGCACTTGGCGGGCCTTGGTGGATGATTCAGTTACGATGCCGCGTACTGCAGGCGGGGCTGGTATTGCAGGGCAAGCAC  >  NC_002947/4393286‑4393575

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.