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

CATGTAGAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAA  >  NC_002947/4507632‑4507925                                                                                                                                                    |                                                                                                                                                   tgtaggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGc                                                                                                                                                   <  3:19368‑M1/2‑1 (MQ=255)  gtaggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCt                                                                                                                                                  <  3:222214‑M1/3‑1 (MQ=255)   taggcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCt                                                                                                                                                  >  1:43971‑M1/146‑148 (MQ=255)      gcaagcgccgcgctgatcattcagtgccgcaggaccagatcaaccgcctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGgcc                                                                                                                                              <  8:33341‑M1/7‑1 (MQ=255)      aGAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGgcc                                                                                                                                              >  7:318272/1‑149 (MQ=255)       gAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGgcc                                                                                                                                              <  7:310845/148‑1 (MQ=255)             cAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTaa                                                                                                                                       >  8:108783/1‑149 (MQ=255)             cAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTaa                                                                                                                                       >  1:56918/1‑149 (MQ=255)                  tCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCt                                                                                                                                  >  2:233895/1‑149 (MQ=255)                      cgcCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTgg                                                                                                                                <  7:108783/147‑1 (MQ=255)                      cgcCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTgg                                                                                                                                <  2:56918/147‑1 (MQ=255)                         cGGCGCCACCCGGTCATTCCACGGCAGGTTGCTGGGGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGGGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGcacca                                                                                                                           <  4:105512/149‑1 (MQ=255)                                            cAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTg                                                                                                                                 >  6:291154/1‑124 (MQ=255)                                            cAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTg                                                                                                                                 <  5:291154/124‑1 (MQ=255)                                               ggCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTc                                                                                                     <  6:69667/149‑1 (MQ=255)                                                     ctgttggtcgatctggctcggcagggcaagcgggtgctgcgcctgaagggtggcgacccgttcatcttcggccggggcggcgaagagatcgaagaGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAg                                                                                               >  3:296907‑M1/96‑149 (MQ=255)                                                          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                                                                                        >  7:241010/1‑135 (MQ=255)                                                                          tGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCATAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCg                                                                                        <  8:241010/135‑1 (MQ=255)                                                                              cGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGTCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCgg                                                                      >  5:57022/1‑149 (MQ=255)                                                                               ggTGACGAAGCGCACCCACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCATCTCTTCGATCTCTTCGCCGCCCCGGc                                                                     >  4:95374/1‑149 (MQ=255)                                                                                      aaGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGAt                                                              >  7:283703/1‑149 (MQ=255)                                                                                      aaGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGAt                                                              <  8:215934/149‑1 (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                                <  8:283703/149‑1 (MQ=255)                                                                                                                    gTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCAc                                >  6:139408/1‑149 (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)                                                                                                                                                    |                                                                                                                                                   CATGTAGAACACCAAGGTCTGCGCCGGCGCCACCAGGTCATTCCAAGGCAGGTTGCTGGTGCCGTCCTTGAGGTGCCCGGTGACGAAGCGCACCGACTGGGCATAGTCGCGGTGAGTCAGCGGAATGCCGCCATACGCGGAGCAGCCGCTGGCCGCCGTAATGCCTGGCACCACCTGGAACGGAATGCCATGCTCGGCCAGCTCTTCGATCTCTTCGCCGCCCCGGCCGAAGATGAACGGGTCGCCACCCTTCAGGCGCAGCACCCGCTTGCCCTGCCGAGCCAGATCGACCAA  >  NC_002947/4507632‑4507925

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.