New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 863665	57 (1.060)	4 (0.080)	4/252	NT	6.1%	coding (799/1017 nt)	hemH	ferrochelatase
	?	NC_002947	= 863693	68 (1.320)					coding (827/1017 nt)	hemH	ferrochelatase

CAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATC  >  NC_002947/863551‑863836                                                                                                                                               |                                                                                                                                                cAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCcggtgt                                                                                                                                           >  2:84059‑M2/1‑143 (MQ=255)   aTGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGtgtg                                                                                                                                          >  4:303271/1‑148 (MQ=255)   aTGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGCCTGGATCAGCTGGCCAAGGCcggtctg                                                                                                                                          <  8:24476‑M2/148‑8 (MQ=255)    tGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGtgtg                                                                                                                                          >  1:339042/1‑147 (MQ=255)     gCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGtgtg                                                                                                                                          >  8:223564/1‑146 (MQ=255)      ccTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTgaag                                                                                                                                       <  3:296905/148‑1 (MQ=255)                   aaaGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGAt                                                                                                                         <  1:104752/149‑1 (MQ=255)                        cTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGcc                                                                                                                    <  8:153768/149‑1 (MQ=255)                             cgcAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCg                                                                                                               >  3:155242/1‑149 (MQ=255)                                 gagCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTtg                                                                                                           <  3:303271/149‑1 (MQ=255)                                 gagCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTtg                                                                                                           <  7:24476/149‑1 (MQ=255)                                   gCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTtgtg                                                                                                         <  7:32274/149‑1 (MQ=255)                                    cATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATCTGCACGGCGTTTGTgg                                                                                                        >  6:282707/1‑149 (MQ=255)                                     aTGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGc                                                                                                       <  1:136314/149‑1 (MQ=255)                                       gggCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCg                                                                                                      >  6:92119/1‑148 (MQ=255)                                       gggCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACAGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTTCTGGTCATGTGCCCGCCGTTTGTTGCg                                                                                                      >  6:243927/1‑148 (MQ=255)                                         gCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTgaagaa                                                                                                                                     >  2:211668/1‑115 (MQ=255)                                          cATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGAtt                                                                                                  <  2:64451/149‑1 (MQ=255)                                            tACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGc                                                                                                <  5:92119/149‑1 (MQ=255)                                              ccGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGATGGTGATGTGCCCGGCGTTTGTGGCGGATTGCa                                                                                               <  8:144193/148‑1 (MQ=255)                                               cGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATc                                                                                             <  6:185668/149‑1 (MQ=255)                                               cGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATc                                                                                             <  7:127852/149‑1 (MQ=255)                                               cGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATc                                                                                             <  2:106643/149‑1 (MQ=255)                                               cGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATc                                                                                             >  7:182119/1‑149 (MQ=255)                                                ggATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATc                                                                                             >  7:128834/1‑148 (MQ=255)                                                ggATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCg                                                                                            <  2:339042/149‑1 (MQ=255)                                                  aTGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTTGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCgag                                                                                          <  5:284816/149‑1 (MQ=255)                                                   tGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCgaga                                                                                         <  8:128834/149‑1 (MQ=255)                                                    ggCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGAc                                                                                        >  8:22386/1‑149 (MQ=255)                                                      cAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGc                                                                                      >  3:77402/1‑149 (MQ=255)                                                       aaaTGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCt                                                                                     <  6:223099/149‑1 (MQ=255)                                                       aaaTGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCt                                                                                     >  5:310285/1‑149 (MQ=255)                                                       aaaTGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCt                                                                                     >  1:330476/1‑149 (MQ=255)                                                                  cccggcccAGCCGCGACTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAg                                                                          >  2:138925/7‑149 (MQ=255)                                                                   tCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAg                                                                          <  7:22386/148‑1 (MQ=255)                                                                   tCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAg                                                                          <  7:223564/148‑1 (MQ=255)                                                                     gTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGTAGAGATAGGCa                                                                       >  8:214974/1‑149 (MQ=255)                                                                      ttCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCAt                                                                      >  1:258265/1‑149 (MQ=255)                                                                         cAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCg                                                                   >  2:153468/1‑149 (MQ=255)                                                                          aGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGt                                                                  >  3:201587/1‑149 (MQ=255)                                                                                   tGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAgaga                                                                             >  1:202061/1‑129 (MQ=255)                                                                                   tGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAgaga                                                                             <  2:202061/129‑1 (MQ=255)                                                                                       ccGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGt                                                     <  7:94956/149‑1 (MQ=255)                                                                                           gACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGACGGTCTGAACAAGCTGCTGGTGATGTGCCCGGAGATTGTGGCGGATTGCATCGAGACGCTGGAACAGATAGGCATGCGTGGCAGCGAGCAGttag                                                  >  8:175906/1‑146 (MQ=255)                                                                                                gTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAgg                                            >  4:181969/1‑149 (MQ=255)                                                                                                    aTCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCcggtgtagggttcgatccacttgtcccggcccagccgcgactggaacgataccgaccatttgccatccggtatgcccatgctctgcgcaaaggctttggccgtag                                        <  8:49604‑M2/149‑106 (MQ=255)                                                                                                       gAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAggcgggcg                                      <  1:138925/148‑1 (MQ=255)                                                                                                         aCCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCcggtgtagggttcgatccacttgtcccggcccagccgcgactggaacgataccgaccatttgccatccggtatgcccatgctctgc                                                           <  2:213469‑M2/125‑87 (MQ=255)                                                                                                         aCCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCcggtgtagggttcgatccacttgtcccggcccagccgcgactggaacgataccgaccatttgccatccggtatgcccatgctctgc                                                           >  1:213469‑M2/1‑39 (MQ=255)                                                                                                         aCCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGcc                                   <  2:258265/149‑1 (MQ=255)                                                                                                          cccTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGggcggc                                    >  5:306108/1‑147 (MQ=255)                                                                                                           ccTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTTCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAg                                 <  4:155242/149‑1 (MQ=255)                                                                                                           ccTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAg                                 >  3:298979/1‑149 (MQ=255)                                                                                                            cTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTTGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAgg                                <  3:181969/149‑1 (MQ=255)                                                                                                            cTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAg                                 >  5:65696/1‑148 (MQ=255)                                                                                                                acCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgct                            >  7:276652/1‑149 (MQ=255)                                                                                                                  cGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctgg                          >  2:94643/1‑149 (MQ=255)                                                                                                                  cGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctgg                          >  5:198308/1‑149 (MQ=255)                                                                                                                   gagaCGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctggt                         >  8:289120/1‑149 (MQ=255)                                                                                                                    agaCGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGTCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctggtg                        <  4:298979/149‑1 (MQ=255)                                                                                                                       cgcgGCTGGATGAGCTGGCCAAGGCcggtgtagggttcgatccacttgtcccggcccagccgcgactggaacgataccgaccatttgccatccggtatgcccatgctctgcgcaaaggctttggccgtagccaggcattgcccgcgata                     >  7:49604‑M2/1‑25 (MQ=255)                                                                                                                        gcgGCTGGATGAGCTGGCCAAGTCGGGGGTGAAGAAGCTGCTGGTGATGTGCGCGGCGTTTGTGGCGGTTTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGGTTGTCGAGGCGGGCGGCCAGGAgctggtgctgg                    <  5:294083/149‑1 (MQ=255)                                                                                                                        gcgGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctggtgctgg                    >  7:209774/1‑149 (MQ=255)                                                                                                                          ggCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGctggtgctggtg                  >  7:107768/1‑149 (MQ=255)                                                                                                                            cTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGgtgcc                >  4:237558/1‑149 (MQ=255)                                                                                                                              ggATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGgtgccgt              <  7:289120/149‑1 (MQ=255)                                                                                                                                aTGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGgtgccgtgc            <  1:153468/149‑1 (MQ=255)                                                                                                                                 tGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGgtgccgtgcc           <  8:107768/149‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTgg                                                                >  4:321163/1‑89 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTgg                                                                <  3:321163/89‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctg                           <  5:206832/126‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAgctg                           >  6:206832/1‑126 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAAt      <  1:94643/147‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAAt      <  6:306108/147‑1 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGa    >  8:114908/1‑149 (MQ=255)                                                                                                                                        gCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGa    <  6:198308/149‑1 (MQ=255)                                                                                                                                         ccAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCatgaatga    >  7:47834/1‑148 (MQ=255)                                                                                                                                          ccaGGCCGGTGTGTAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCATCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATc  <  2:235827/147‑1 (MQ=255)                                                                                                                                          cAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTATGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGAt   >  1:235827/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                CAATGCCTGGCTACGGCCAAAGCCTTTGCGCAGAGCATGGGCATACCGGATGGCAAATGGTCGGTATCGTTCCAGTCGCGGCTGGGCCGGGACAAGTGGATCGAACCCTACACCGAGACGCGGCTGGATGAGCTGGCCAAGGCCGGTGTGAAGAAGCTGCTGGTGATGTGCCCGGCGTTTGTGGCGGATTGCATCGAGACGCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATC  >  NC_002947/863551‑863836

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