New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 86499054 (0.920)3 (0.050) 3/288 7.8 5.4% coding (140/750 nt) moeB molybdopterin synthase sulfurylase
?NC_000913_3_bme_pgi = 865073 52 (0.920)coding (57/750 nt) moeB molybdopterin synthase sulfurylase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCT  >  NC_000913_3_bme_pgi/864841‑865138
                                                                                                                                                     |                                                                                                                                                    
tCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGCTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCg                                                                                                                                                      <  1:92787/150‑1 (MQ=255)
     gTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCt                                                                                                                                                 <  2:254227/150‑1 (MQ=255)
      gcAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTg                                                                                                                                                <  1:459495/149‑1 (MQ=255)
        cgggCGGCACGGGCGGATTCCACGTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGGCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGca                                                                                                                                              <  2:302765/149‑1 (MQ=255)
         gggCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGcac                                                                                                                                             <  1:38822/150‑1 (MQ=255)
          ggCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGcaca                                                                                                                                            <  1:602521/150‑1 (MQ=255)
                   ggCGGATTCCACCTTCGGTTGCGCGCCCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGGCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCAGAGCCGAGgcc                                                                                                                                   <  2:210146/150‑1 (MQ=255)
                       gATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGCGGCCGccc                                                                                                                               >  2:360551/1‑150 (MQ=255)
                        aatCCACCTTCGGTTTCCCGAACGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTTCCAGATACTGAGAGGCTGCACAGCCGAGGCCGCCCa                                                                                                                              <  2:631226/148‑1 (MQ=255)
                            cACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc                                                                                                                          <  2:126153/150‑1 (MQ=255)
                              ccTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCa                                                                                                                        >  1:330662/1‑150 (MQ=255)
                                ttCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACt                                                                                                                      <  1:595107/150‑1 (MQ=255)
                                   ggTTGCCCGACCGTGGCATCACTGTGGAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTGCCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGACGCCGCCCAGGCCCCCTATc                                                                                                                   >  2:233906/1‑150 (MQ=255)
                                    gTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTACACAGCCGAGGCCGCCCAGGCCCACTATCa                                                                                                                  <  1:379040/150‑1 (MQ=255)
                                       gCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCa                                                                                                               <  1:38020/150‑1 (MQ=255)
                                        cccGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCa                                                                                                                        <  2:814440/140‑1 (MQ=255)
                                        cccGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCa                                                                                                                        >  1:814440/1‑140 (MQ=255)
                                         ccGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTAACATCAcg                                                                                                             >  2:393731/1‑150 (MQ=255)
                                            aCCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCga                                                                                                          <  1:264964/150‑1 (MQ=255)
                                             ccGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGTCCACTATCAGCACGCgag                                                                                                         >  1:196702/1‑150 (MQ=255)
                                             ccGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCgag                                                                                                         <  2:927865/150‑1 (MQ=255)
                                                tGGCATCACTGTGCAGTGGCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAAt                                                                                                      <  2:758039/150‑1 (MQ=255)
                                                gggCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGGCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAAt                                                                                                      <  2:762822/149‑1 (MQ=255)
                                                  gCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGATCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCt                                                                                                    <  2:793818/150‑1 (MQ=255)
                                                     gcACTGTGCAGTGGCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGGCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTc                                                                                                 <  1:51600/149‑1 (MQ=255)
                                                         tgtgCAGTGTCTGGCGTTGCAGATTCGCGAGCGAAACCGTGTCGAAGGTGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAgcg                                                                                             <  1:309334/150‑1 (MQ=255)
                                                            gCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGcct                                                                                          >  2:137707/1‑150 (MQ=255)
                                                              aGTGTCTGGCGGTGCAGATTCGAGAGCGAAACCGTGTCGAAGGCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGTGcctcc                                                                                        <  2:222288/150‑1 (MQ=255)
                                                                    tGGCGTTGCCGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTAGCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCg                                                                                  <  1:180877/150‑1 (MQ=255)
                                                                      gCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                                                                                <  2:330662/150‑1 (MQ=255)
                                                                           gCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAAt                                                                           <  2:785789/150‑1 (MQ=255)
                                                                           gCAGATTCGAGAGCGAAACCGTGGCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAAt                                                                           <  1:877904/150‑1 (MQ=255)
                                                                             aGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCa                                                                         >  2:350942/1‑150 (MQ=255)
                                                                              gATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGCGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCaa                                                                        >  1:217918/1‑150 (MQ=255)
                                                                                ttCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAg                                                                      <  1:679763/150‑1 (MQ=255)
                                                                                  cGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCCCAGCCGAGGCCGCCCAGGCCCACTATCAGCGCGCGAGAATCTTTCAGCGCCTCCTGGCCGTCACAATCAACGcc                                                                    >  1:528704/1‑150 (MQ=255)
                                                                                     gagCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCcgcg                                                                 >  1:3363/1‑150 (MQ=255)
                                                                                       ccGAAACCGTGTCGAAGTCGAGCAGCGGCAGGTTACCGACACCGGCGCTTGCCTGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGCCCCACTAGCAGCACGCGAGAATCTATCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCa                                                               <  2:932318/149‑1 (MQ=255)
                                                                                          aaaCCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGaa                                                            >  1:877842/1‑150 (MQ=255)
                                                                                                 gtCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGTCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTg                                                     <  2:196702/150‑1 (MQ=255)
                                                                                                    gAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCCGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGCATGATTTGCCg                                                  >  2:332889/1‑150 (MQ=255)
                                                                                                     aaGTAGAGCAGCGTCAGGTTACCGACACCGTCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAAGCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCgg                                                 <  1:184788/150‑1 (MQ=255)
                                                                                                              agcGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAgcgc                                        <  1:739747/150‑1 (MQ=255)
                                                                                                                cggcAGGTTCCCGACACCGGCGCTTGCGAGATCCTGCGAGGCGGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCCCGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCgcag                                      <  1:694333/147‑1 (MQ=255)
                                                                                                                cGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCgcag                                      <  1:824707/150‑1 (MQ=255)
                                                                                                                   cAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGTGAGAATCTTTCAGCGCCTCCTCGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCAt                                   >  1:742493/1‑150 (MQ=255)
                                                                                                                       ttACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGCGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTcc                               <  2:217918/150‑1 (MQ=255)
                                                                                                                         aCCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTg                             <  2:828388/150‑1 (MQ=255)
                                                                                                                          ccGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGTGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGa                            >  2:890546/1‑150 (MQ=255)
                                                                                                                                  ggCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAg                    <  1:350942/150‑1 (MQ=255)
                                                                                                                                      gttGCCAGATACTGCGAGGGTGTACAGCCGAGGCCGCCCAGGGCCACTAGCAGGACGCGAGAATCTTTCAGGGCCTCCTGGACTTCAAATTCAAAGCCGCGAAGAATGTTTTGCCTTTTGTAGCGCAGCTTCTCCTGATCGCTTAGttgg                <  2:657147/149‑3 (MQ=255)
                                                                                                                                         gCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGcc             <  2:858095/150‑1 (MQ=255)
                                                                                                                                                 cTGCGAGGCTGCTCAGCCGAGGCCGCCCAGGCCGACTATCTGCAGGCGAGATTCTTTCAGCGCCTCCTGGCCGTCAAAATCAAATCCGCGCAGAATCATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAgg     <  1:233906/150‑1 (MQ=255)
                                                                                                                                                   gCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGcc   <  2:3363/150‑1 (MQ=255)
                                                                                                                                                    cGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCACAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCCGCATCTCCTGATCGCGGAGTTCCGCCATTACAGGCCt  >  1:880167/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCT  >  NC_000913_3_bme_pgi/864841‑865138

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