New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 375811442 (0.650)4 (0.070) 4/248 NT 8.5% coding (427/1002 nt) hemBB porphobilinogen synthase
?NC_002947 = 3758136 47 (0.770)coding (449/1002 nt) hemBB porphobilinogen synthase

TCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCCGCGCT  >  NC_002947/3757966‑3758262
                                                                                                                                                    |                                                                                                                                                    
tCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACG                                                                                                                                                      >  5:37867‑M1/1‑149 (MQ=255)
   aGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCacgacg                                                                                                                                                   >  5:83026/1‑149 (MQ=255)
           gcAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGtgctg                                                                                                                                           <  5:249838/149‑1 (MQ=255)
             cgcaccAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTgcg                                                                                                                                         <  7:248501/144‑1 (MQ=255)
                       aGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGc                                                                                                                               >  7:207688‑M2/1‑148 (MQ=255)
                       aGGCCACGCAGGCCGTCGACTGGGACAACCTGATCTGCCGCGCGGCCAGGTTGCTCATGGCCACGTTCCCCAACCTTGGCTTGATCGGCGAAGGCGCACTGGACCCCTACACGGACCATGGCCACGACGGCGTGCTGGGCCAAGGCCGc                                                                                                                               >  6:343131‑M2/1‑148 (MQ=255)
                          gCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGcgtg                                                                                                                            >  1:263972‑M2/1‑145 (MQ=255)
                              gggtgccctCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGcgtggcc                                                                                                                         <  1:84753‑M2/144‑8 (MQ=255)
                              ccACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAAACTAGGCTTGATCGGCGCCGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGcgtggcc                                                                                                                         <  6:48282‑M2/148‑8 (MQ=255)
                                   cccTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGAc                                                                                                                   <  6:136177/149‑1 (MQ=255)
                                    ccTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGa                                                                                                                    <  2:106433/147‑1 (MQ=255)
                                    ccTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGAca                                                                                                                  >  1:106433/1‑148 (MQ=255)
                                        gACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACgag                                                                                                                <  3:322824/147‑1 (MQ=255)
                                             ggACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGcggccttggcgcagcacgccgtcgtgg                                                                                                                           <  6:285866‑M1/131‑28 (MQ=255)
                                             ggACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGcggccttggcgcagcacgccgtcgtgg                                                                                                                           >  5:285866‑M1/1‑104 (MQ=255)
                                             ggACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTg                                                                                                          >  1:269277/1‑148 (MQ=255)
                                              gACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCa                                                                                                        <  6:83026/149‑1 (MQ=255)
                                                  aCCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTa                                                                                                    >  7:70979/1‑149 (MQ=255)
                                                  aCCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTa                                                                                                    <  8:297503/149‑1 (MQ=255)
                                                   ccTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACATTGCAGGTAc                                                                                                   >  8:339355/1‑149 (MQ=255)
                                                          tGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTggcggc                                                                                            >  5:297005/1‑149 (MQ=255)
                                                           gccgcgccgcCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGcc                                                                                           >  1:183449/1‑149 (MQ=255)
                                                            ccgcgcctcCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGGAGGTACTggcgg                                                                                             <  8:70979/146‑1 (MQ=255)
                                                              gcgccgcCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAg                                                                                        >  6:109416/1‑149 (MQ=255)
                                                               cgccgcCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAgg                                                                                       >  5:133655/1‑149 (MQ=255)
                                                                    ccAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCtgtg                                                                                  <  8:1119/149‑1 (MQ=255)
                                                                    ccAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCtgtg                                                                                  >  6:8912/1‑149 (MQ=255)
                                                                        gTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTc                                                                               >  1:251167/1‑148 (MQ=255)
                                                                             tCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGcggccttggcgcagcacgccgtcgtggccatggtcggtgta                                                                                                             >  3:361478‑M1/1‑72 (MQ=255)
                                                                             tCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGcggccttggcgcagcacgccgtcgtggccatggtcggtgta                                                                                                             <  4:361478‑M1/113‑42 (MQ=255)
                                                                             tCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTccaggcc                                                                         >  2:216281/1‑149 (MQ=255)
                                                                               aaGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGTCCACGACGGc                                                                                                                                                 <  2:402641/75‑1 (MQ=255)
                                                                               aaGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGc                                                                                                                                                 >  1:402641/1‑75 (MQ=255)
                                                                                aGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCTAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTccaggcccag                                                                      <  2:76887/149‑1 (MQ=255)
                                                                                 ggCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTggcgg                                                                                             <  5:393644/125‑1 (MQ=255)
                                                                                 ggCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTggcgg                                                                                             >  6:393644/1‑125 (MQ=255)
                                                                                   ccACGTTCCCCACCCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCcaggcccaggcc                                                                   <  1:140404/149‑1 (MQ=255)
                                                                                   ccACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCcaggcccaggcc                                                                   >  2:140404/1‑149 (MQ=255)
                                                                                      cGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGc                                                                <  2:183449/149‑1 (MQ=255)
                                                                                             cAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGaca                                                         <  6:133655/149‑1 (MQ=255)
                                                                                              aaCCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGacac                                                        >  7:190610/1‑149 (MQ=255)
                                                                                                   tGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTgg                                                   >  5:245780/1‑149 (MQ=255)
                                                                                                      cTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCAc                                                >  7:290673/1‑149 (MQ=255)
                                                                                                         gATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCt                                              <  2:269277/148‑1 (MQ=255)
                                                                                                              gCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGAc                                        <  6:297005/149‑1 (MQ=255)
                                                                                                                 aCGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCAGGCCGGCATCGACACCTTGGCACCTTCGGACatga                                    <  6:37867/149‑1 (MQ=255)
                                                                                                                   gTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACatgat                                   <  2:251167/149‑1 (MQ=255)
                                                                                                                    tCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACatgatg                                  <  1:206853/149‑1 (MQ=255)
                                                                                                                      gCACTGGACCCCTACACCGTCCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGa                                >  7:49168/1‑149 (MQ=255)
                                                                                                                               cccTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCt                                              >  2:38502/1‑126 (MQ=255)
                                                                                                                               cccTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCt                                              <  1:38502/126‑1 (MQ=255)
                                                                                                                               cccTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCAGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCAt                       >  4:41397/1‑149 (MQ=255)
                                                                                                                                 cTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCg                     <  5:188556/149‑1 (MQ=255)
                                                                                                                                 cTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTACGGACATGATGGACGGGCGCATCg                     >  8:112501/1‑149 (MQ=255)
                                                                                                                                  tACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttga                                                        >  3:361478‑M2/1‑41 (MQ=255)
                                                                                                                                  tACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttga                                                        <  4:361478‑M2/113‑73 (MQ=255)
                                                                                                                                           cATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAg           <  8:190610/149‑1 (MQ=255)
                                                                                                                                              ggCCACGcggccttggcgcagcacgccgtcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtccgcgtcgagggcaccc         <  1:84753‑M1/148‑142 (MQ=255)
                                                                                                                                              ggCCACGcggccttggcgcagcacgccgtcgtggccatggtcggtgtaggggtccagtgcgacggcgccgatcaagcctagtttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtccgcgtcgagggcgtgg         <  6:48282‑M1/148‑142 (MQ=255)
                                                                                                                                                ccACGACGGCGTGCTGCGCCAAGGCCGcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtcc                        >  5:285866‑M2/1‑27 (MQ=255)
                                                                                                                                                ccACGACGGCGTGCTGCGCCAAGGCCGcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtcc                        <  6:285866‑M2/131‑105 (MQ=255)
                                                                                                                                                ccACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCcg      >  5:2651/1‑149 (MQ=255)
                                                                                                                                                ccACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCcg      >  7:362184/1‑149 (MQ=255)
                                                                                                                                                ccACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCcg      <  5:48282/149‑1 (MQ=255)
                                                                                                                                                 cACGcggccttggcgcagcacgccgtcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtccgcgtcgagggcgtgggtgc     >  1:263972‑M1/1‑4 (MQ=255)
                                                                                                                                                 cACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGc                                                  >  6:87375/1‑104 (MQ=255)
                                                                                                                                                 cACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGc                                                  <  5:87375/104‑1 (MQ=255)
                                                                                                                                                  acgacgGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCcgcg    >  3:147491/1‑149 (MQ=255)
                                                                                                                                                    gcggccttggcgcagcacgccgtcgtggccatggtcggtgtaggggtccagtgcgacgtcgccgatcaagccaaggttggggaacgtggccttgagcaacctggcggcgcggcagatcaggttgtccgcgtcgagggcgtgggtgccct  >  7:207688‑M1/1‑1 (MQ=255)
                                                                                                                                                    gcggccttggcccagcacgccgtcgtggccatggtccgtgtaggggtccagtgcgccttcgccgatcaagccaaggttggggaacgtggccatgagcaacctggccgcgcggcagatcaggttgtcccagtcgacggcctgcgtggcct  >  6:343131‑M1/1‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCCAGGTTGCTCAAGGCCACGTTCCCCAACCTTGGCTTGATCGGCGACGTCGCACTGGACCCCTACACCGACCATGGCCACGACGGCGTGCTGCGCCAAGGCCGCGTGGCCAATGACGAGACCTTGCAGGTACTGGCGGCCCAGGCTGTGGTCCAGGCCCAGGCCGGCATCGACACCTTGGCACCTTCGGACATGATGGACGGGCGCATCGGCGCGCTCAGGTCCGCGCT  >  NC_002947/3757966‑3758262

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.