New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4507779 =	55 (1.020)	3 (0.060)	3/244	NT	5.2%	coding (1002/1392 nt)	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
	?	NC_002947	4507834 =	57 (1.120)					coding (947/1392 nt)	cysG	multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

GGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGCCTACA  >  NC_002947/4507689‑4507980                                                                                                                                                  |                                                                                                                                                   ggTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCtt                                                                                                                                                 >  4:302987/1‑149 (MQ=255)     ccGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGAt                                                                                                                                             >  8:204071/1‑149 (MQ=255)     ccGTCCTTGAGGTGCCCGGTGAAGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGa                                                                                                                                              >  7:215934/1‑148 (MQ=255)          cTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCtt                                                                                                                                        >  5:161637/1‑149 (MQ=255)                gcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCTAGct                                                                                                                                                    <  2:303135‑M1/56‑1 (MQ=255)                gcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCTAGct                                                                                                                                                    >  1:303135‑M1/76‑131 (MQ=255)                 tGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCATAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCg                                                                                                                                               <  8:241010/135‑1 (MQ=255)                 tGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCATAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCg                                                                                                                                               >  7:241010/1‑135 (MQ=255)                     cGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGTCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCgg                                                                                                                             >  5:57022/1‑149 (MQ=255)                      ggTGACGAAGCGCACCCACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCATCTCTTCGATCTCTTCGCCGCCCCGGc                                                                                                                            >  4:95374/1‑149 (MQ=255)                             aaGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGAt                                                                                                                     <  8:215934/149‑1 (MQ=255)                             aaGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGAt                                                                                                                     >  7:283703/1‑149 (MQ=255)                                ctcACCTACTGGGCATAGTCGAGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGTAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGaa                                                                                                                  <  8:318272/147‑1 (MQ=255)                                  cACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACg                                                                                                                >  4:324416/1‑149 (MQ=255)                                  cACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCACCGCCCCGGCCGAAGAt                                                                                                                     <  3:280681/144‑1 (MQ=255)                                  cACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCACCGCCCCGGCCGAAGAt                                                                                                                     >  4:280681/1‑144 (MQ=255)                                   aCCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGTCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACgg                                                                                                               <  6:57022/149‑1 (MQ=255)                                   aCCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACgg                                                                                                               >  8:169841/1‑149 (MQ=255)                                   aCCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCAGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACgg                                                                                                               >  3:59748/1‑149 (MQ=255)                                     cGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGt                                                                                                             <  1:58989/149‑1 (MQ=255)                                         tGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGcc                                                                                                         >  8:277530/1‑149 (MQ=255)                                         tGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGc                                                                                                          >  2:9593/1‑148 (MQ=255)                                          gggCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCa                                                                                                        >  4:219477/1‑149 (MQ=255)                                           ggCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCAc                                                                                                       <  1:136191/149‑1 (MQ=255)                                              aTAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCt                                                                                                    >  2:246668/1‑149 (MQ=255)                                                 gTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCa                                                                                                 <  6:161637/149‑1 (MQ=255)                                                   cgcgGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAgg                                                                                               >  8:270351/1‑149 (MQ=255)                                                       gTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCgca                                                                                           <  7:277530/149‑1 (MQ=255)                                                          aGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCgcagca                                                                                        <  3:302987/149‑1 (MQ=255)                                                           gTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCAc                                                                                       >  6:139408/1‑149 (MQ=255)                                                           gTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCAc                                                                                       <  8:283703/149‑1 (MQ=255)                                                            tCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCAcc                                                                                      >  5:352142/1‑149 (MQ=255)                                                              aGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGgcgc                                                                                            <  7:150228/141‑1 (MQ=255)                                                              aGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGgcgc                                                                                            >  8:150228/1‑141 (MQ=255)                                                              aGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCg                                                                                    >  2:30468/1‑149 (MQ=255)                                                                   aaTGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGc                                                                               >  1:184411/1‑149 (MQ=255)                                                                       ccgccATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCt                                                                            <  7:204071/148‑1 (MQ=255)                                                                       ccgccATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTg                                                                           >  2:78983/1‑149 (MQ=255)                                                                        cgccATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCAGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGc                                                                          <  4:59748/149‑1 (MQ=255)                                                                         gccATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGcc                                                                         <  7:169841/149‑1 (MQ=255)                                                                           cATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGa                                                                       <  1:246668/149‑1 (MQ=255)                                                                           cATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGa                                                                       <  1:30468/149‑1 (MQ=255)                                                                            aTACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAg                                                                      >  4:222214/1‑149 (MQ=255)                                                                             tACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGc                                                                     <  6:352142/149‑1 (MQ=255)                                                                               cgcgGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGcc                                                                    <  1:233895/148‑1 (MQ=255)                                                                                  ggAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGc                                                                     >  4:19368/1‑144 (MQ=255)                                                                                         ccGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCaa                                                         >  7:33341/1‑149 (MQ=255)                                                                                          cGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCaa                                                         >  3:165001/1‑148 (MQ=255)                                                                                          agctagccgagcatggcattccgttccaggtggtgccaggcattacggcggccagcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGc                                                                          >  1:303135‑M2/57‑131 (MQ=255)                                                                                          agctagccgagcatggcattccgttccaggtggtgccaggcattacggcggccagcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGc                                                                          <  2:303135‑M2/75‑1 (MQ=255)                                                                                            ctggccgagcatggcattccgttccaggtggtgccaggcattacggcggccagcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAg                                                      >  3:296907‑M2/55‑149 (MQ=255)                                                                                                 CGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGt                                                 <  1:45756‑M1/149‑1 (MQ=255)                                                                                                 CGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGt                                                 <  4:131887‑M1/149‑1 (MQ=255)                                                                                                   ccgTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTg                                               <  3:219477/149‑1 (MQ=255)                                                                                                    cgTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGa                                              <  3:95374/149‑1 (MQ=255)                                                                                                     gTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGAt                                             <  1:9593/149‑1 (MQ=255)                                                                                                      tAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTTATc                                            <  3:324416/149‑1 (MQ=255)                                                                                                            cTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGa                                                                 <  7:210632/122‑1 (MQ=255)                                                                                                            cTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGa                                                                 >  8:210632/1‑122 (MQ=255)                                                                                                               gCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGc                                   <  5:139408/149‑1 (MQ=255)                                                                                                                 accaccTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCgg                                 >  7:135523/1‑149 (MQ=255)                                                                                                                  ccaccTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCgg                                 >  6:321170/1‑148 (MQ=255)                                                                                                                   caccTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCa                               >  3:84080/1‑149 (MQ=255)                                                                                                                         gAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGaa                          <  4:84080/148‑1 (MQ=255)                                                                                                                         gAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAAt                         <  1:78983/149‑1 (MQ=255)                                                                                                                                      tGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGATCGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCgcggcg            <  2:184411/149‑1 (MQ=255)                                                                                                                                         tCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAgcgcatctatg         >  7:36875/1‑142 (MQ=255)                                                                                                                                           ggccagcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGc       <  8:33341‑M2/142‑1 (MQ=255)                                                                                                                                               agcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGCCTa    >  1:43971‑M2/4‑148 (MQ=255)                                                                                                                                               agcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGCCTAc   <  3:222214‑M2/146‑1 (MQ=255)                                                                                                                                                gcTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGCCTACa  <  3:19368‑M2/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAACAGGCGGTTGATCTGGTCCTGCGGCACTGAATGATCAGCGCGGCGCTTGCCTACA  >  NC_002947/4507689‑4507980

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