New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	996911 =	35 (1.180)	6 (0.220)	5/226	NT	15.0%	coding (688/2550 nt)	PP_0860	sulfite reductase flavoprotein subunit
	?	NC_002947	996946 =	36 (1.330)					coding (653/2550 nt)	PP_0860	sulfite reductase flavoprotein subunit

GGGCGGCAGGCGCAGGTTGTAGGTGGACAGGCCTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCAGTTGAAG  >  NC_002947/996762‑997059                                                                                                                                                      |                                                                                                                                                        gggcggCAGGCGCAGGTTGTAGGTGGACAGGCCTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTtg                                                                                                                                                         <  2:687946/150‑1 (MQ=255)      gaccttcgactgggccaagaaaggccgcgccttcaactgggacctgcacgccgttgccggcacctggtgcctgctgttctatcttctgttcgccctgaccgggctgttctggtcctatgagtggtaccgcgaaggcctgaacaagctgct                                                                                                                                                    <  2:114577‑M1/6‑1 (MQ=255)         cttcgactgggccaagaaaggccgcgccttcaactgggacctgcacgccgttgccggcacctggtgcctgctgttctatcttctgttcgccctgaccgggctgttctggtcctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGc                                                                                                                                                 >  2:519578‑M1/142‑150 (MQ=255)            gcAGGTTGTAGGTGGACAGGCCTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑Cg                                                                                                                                            >  1:66828/1‑150 (MQ=255)                gTTGTAGGTGGACAGGCCTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGc                                                                                                                                        >  2:286268/1‑150 (MQ=255)                            ggccgcgccttcaactgggacctgcacgccgttgccggcacctggtgcctgctgttctatcttctgttcgccctgaccgggctgttctggtcctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGc                                                                                                                           >  2:526798‑M1/123‑150 (MQ=255)                                 cTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCcagca                                                                                                                      >  1:219389/1‑150 (MQ=255)                                  tGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCca                                                                                                                         <  1:791381/146‑1 (MQ=255)                                  tGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCca                                                                                                                         >  2:791381/1‑146 (MQ=255)                                         gcggcgTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTc                                                                                                              <  2:219389/150‑1 (MQ=255)                                           ggcgTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAg                                                                                                            <  1:564198/150‑1 (MQ=255)                                                 cTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTc                                                                                                      >  1:238752/1‑150 (MQ=255)                                                      ggTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGt                                                                                                 >  1:400626/1‑150 (MQ=255)                                                        ttGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTAc                                                                                               <  2:66828/150‑1 (MQ=255)                                                             ccAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGc                                                                                                          >  2:132332/1‑134 (MQ=255)                                                             ccAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGc                                                                                                          <  1:132332/134‑1 (MQ=255)                                                                 aTGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAg                                                                                      <  1:484062/150‑1 (MQ=255)                                                                    gCGTCGTAGTCCACCACCAGAGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGAc                                                                                   <  2:668970/150‑1 (MQ=255)                                                                           aGTCCACCATCAGCGGCGGCGCATTCCTGTCGACAGTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACa                                                                            <  1:4338/150‑1 (MQ=255)                                                                               caccaccaGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGccc                                                                        >  1:501589/1‑150 (MQ=255)                                                                               caccaccaGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGccc                                                                        >  1:114577/1‑150 (MQ=255)                                                                                   accaGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTc                                                                    >  2:315604/1‑150 (MQ=255)                                                                                     caGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAg                                                                  >  2:91604/1‑150 (MQ=255)                                                                                     caGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAg                                                                  <  1:532055/150‑1 (MQ=255)                                                                                             gcgcATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACa                                                          <  2:501589/150‑1 (MQ=255)                                                                                                       gTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGaa                                                <  1:171044/150‑1 (MQ=255)                                                                                                         cgggctgttctggtcctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCcagc                                                                                                                       <  1:660365‑M1/32‑1 (MQ=255)                                                                                                         cgggctgttctggtcctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCcagc                                                                                                                       >  2:660365‑M1/46‑77 (MQ=255)                                                                                                          gACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAg                                                                                                            <  1:1534/87‑1 (MQ=255)                                                                                                          gACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAg                                                                                                            >  2:1534/1‑87 (MQ=255)                                                                                                             gTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGCACAGAAGATAGAAcagcag                                          >  1:29837/1‑150 (MQ=255)                                                                                                                gtggggccggggccggGGACGGCGTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTGAGGGCGAACAGAAGATAGAGCAGCAGGCa                                       <  2:29837/146‑1 (MQ=255)                                                                                                                cTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCa                                       <  1:286268/150‑1 (MQ=255)                                                                                                                        ctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAg                                                                                                            <  2:437054‑M1/43‑1 (MQ=255)                                                                                                                        ctatgagtggtaccgcgaaggcctgaacaaGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAg                                                                                                            >  1:437054‑M1/31‑73 (MQ=255)                                                                                                                        cGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGCACAGACGATAGAACAGAAGACACCatatca                               >  2:171100/1‑145 (MQ=255)                                                                                                                                 cACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACgg                      >  2:453749/1‑150 (MQ=255)                                                                                                                                    ggCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGt                   <  1:548691/150‑1 (MQ=255)                                                                                                                                             cccGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCAc                                      >  2:679644/1‑122 (MQ=255)                                                                                                                                             cccGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCAc                                      <  1:679644/122‑1 (MQ=255)                                                                                                                                              ctgaacaaGCTGCTGGCCGACGC‑GCCGGCTGCG‑‑GGCCAGCGGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCa         <  1:163489‑M1/142‑1 (MQ=255)                                                                                                                                              ccGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGt                                                                                                 >  2:23274/1‑62 (MQ=255)                                                                                                                                              ccGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGt                                                                                                 <  1:23274/62‑1 (MQ=255)                                                                                                                                              ccGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCa         >  1:299160/1‑150 (MQ=255)                                                                                                                                                    ttGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCAGTTGaa   <  2:238752/150‑1 (MQ=255)                                                                                                                                                     tgctgcTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCAGTTGAAg  <  1:91604/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                        GGGCGGCAGGCGCAGGTTGTAGGTGGACAGGCCTGGGCCCGCGGCGTCCTTGAGGTTGGCCCAGATGGCGTCGTAGTCCACCACCAGCGGCGGCGCATTCCTGTCGACGTTCTGCGGCCCGTGCCGGCCACGGCCTTCACCCCGTTTTTGCTGCTGGCC‑‑CGCAGCCGGC‑GCGTCGGCCAGCAGCTTGTTCAGGCCTTCGCGGTACCACTCATAGGACCAGAACAGCCCGGTCAGGGCGAACAGAAGATAGAACAGCAGGCACCAGGTGCCGGCAACGGCGTGCAGGTCCCAGTTGAAG  >  NC_002947/996762‑997059

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.