New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1821027	0 (0.000)	4 (0.040)	4/234	NT	100%	coding (748/1422 nt)	PP_1624	group II intron‑encoding maturase
	?	NC_002947	= 2070376	NA (NA)					coding (707/1422 nt)	PP_1846	group II intron‑encoding maturase

GGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAATATCCTGCTCGACGAACTCGACCGCGAGCTGGAGCGGCGGGGCCATCGCTTCGTACGTTATGCCGATGATGCGAACATCTATGTACGCAGTCCCCGGGCCGGCGAACGGGTGTTGGTCAGCGTCGAGCGCTTCCTGAGA  >  NC_002947/1820996‑1821176                                |                                                                                                                                                      ggggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgca                                                           <  1:2197093‑M1/123‑93 (MQ=255) ggggACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgca                                                           >  2:2197093‑M1/1‑31 (MQ=255)    gACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtc                                            <  1:211019‑M1/135‑108 (MQ=255)    gACGCCGCAAGGCGGCCC‑CTCCTGCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtc                                            >  2:211019‑M1/1‑28 (MQ=255)                            gCCGTgggctggcgtcactgctcgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagtt      <  2:857726‑M1/150‑146 (MQ=255)                            gCCGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtccgcacttgtctttgatccgacgttcgatgcaggccatgtggatgtcgttatttacccgatcgaaaaacttctccagatcgagtt      <  2:348947‑M1/150‑146 (MQ=255)                             ccGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttc     <  1:2245637‑M1/150‑147 (MQ=255)                             ccGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaa                                                                                                      <  1:1988107‑M1/53‑50 (MQ=255)                             ccGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaa                                                                                                      >  2:1988107‑M1/1‑4 (MQ=255)                              cGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttcc    <  2:150634‑M1/150‑148 (MQ=255)                              cGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttcc    >  1:1539595‑M1/1‑3 (MQ=255)                              cGTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttcc    <  2:2045928‑M1/150‑148 (MQ=255)                               gTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgtcgtttatgatgatgtcgtgattgacccgatcgaaaaacttctccagatcgagttcca   <  2:3049756‑M1/150‑149 (MQ=255)                               gTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgtgctaaaaacttctccagatcgagttcca   <  2:1886166‑M1/150‑149 (MQ=255)                               gTgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttcca   <  1:726133‑M1/150‑149 (MQ=255)                                tgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgt                                          >  2:1674113‑M1/1‑1 (MQ=255)                                tgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtacgcacttgtctttgatccgacgttcgatgcaggccatgaggatgtcgt                                          <  1:1674113‑M1/110‑110 (MQ=255)                                tgggctgacgacaccgcccgacatgatcccggcttcaaggtagcggcgaataagcctgagtaagctggttggtttgatccgacgttcgatgcaggccatgaggatgtcgtgattgacccgatcgaaaaacttctccagatcgagttccac  <  1:968597‑M1/150‑150 (MQ=255)                                |                                                                                                                                                      GGGGACGCCGCAAGGCGGCCCGCTCTCGCCGTTGCTGTCGAATATCCTGCTCGACGAACTCGACCGCGAGCTGGAGCGGCGGGGCCATCGCTTCGTACGTTATGCCGATGATGCGAACATCTATGTACGCAGTCCCCGGGCCGGCGAACGGGTGTTGGTCAGCGTCGAGCGCTTCCTGAGA  >  NC_002947/1820996‑1821176

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