Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,873,031 IS1 (–) +9 bp coding (991‑999/1491 nt) dgcJ → putative diguanylate cyclase DgcJ

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 291401NA (NA)416 (1.160) 164/194 0.1 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 1873031 = 0 (0.000)coding (991/1491 nt) dgcJ putative diguanylate cyclase DgcJ
* ? NC_000913 = 18730390 (0.000)369 (1.020) 164/194 0.1 100% coding (999/1491 nt) dgcJ putative diguanylate cyclase DgcJ
?NC_000913 1978503 = NA (NA)noncoding (768/768 nt) IS1 repeat region

Only 100 of 783 total aligned reads displayed.

TCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAGCA  >  NC_000913/1872931‑1873131
                                                                                                    |                                                                                                    
 tcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCg                                                                                                     <  1:3558899‑M2/2‑1 (MQ=255)
   agtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGa                                                                                                   >  2:3373591‑M2/98‑101 (MQ=255)
    gtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGAt                                                                                                  >  1:3373499‑M2/97‑101 (MQ=255)
    gtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGAt                                                                                                  <  2:3558408‑M2/5‑1 (MQ=255)
    gtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGAt                                                                                                  <  2:3558575‑M2/5‑1 (MQ=255)
    gtCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGAt                                                                                                  >  1:1942757/1‑101 (MQ=255)
     tagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATg                                                                                                 <  1:3558021‑M2/6‑1 (MQ=255)
        ctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCAc                                                                                              >  1:3373736‑M2/93‑101 (MQ=255)
        cGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCAc                                                                                              >  2:1942330/1‑101 (MQ=255)
        cGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCAc                                                                                              >  2:3372982/1‑101 (MQ=255)
         tgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACt                                                                                             >  1:3373385‑M2/92‑101 (MQ=255)
         tgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACt                                                                                             >  2:3373771‑M2/92‑101 (MQ=255)
         gTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACg                                                                                             >  2:1942866‑M1/1‑100 (MQ=255)
          ttATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACgg                                                                                            >  2:1942445‑M1/1‑99 (MQ=255)
          gaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACtt                                                                                            >  1:3373915‑M2/91‑101 (MQ=255)
             tCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaa                                                                                         >  1:3372993‑M1/1‑96 (MQ=255)
              ccATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaat                                                                                        <  1:3372701‑M1/101‑7 (MQ=255)
               cATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaa                                                                                         <  2:3372417‑M1/99‑6 (MQ=255)
                gagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGttt                                                                                      <  2:3558253‑M2/17‑1 (MQ=255)
                gagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGttt                                                                                      <  1:3558678‑M2/17‑1 (MQ=255)
                aTTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatga                                                                                      <  2:3372637‑M1/101‑9 (MQ=255)
                    gacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATAc                                                                                  <  1:3558610‑M2/21‑1 (MQ=255)
                    atatTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactcc                                                                                  <  2:3372768‑M1/101‑13 (MQ=255)
                     tatTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactcca                                                                                 >  2:3373000‑M1/1‑88 (MQ=255)
                      cagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCa                                                                                <  2:3558801‑M2/23‑1 (MQ=255)
                      atTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaa                                                                                <  1:3372529‑M1/101‑15 (MQ=255)
                       ttCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaac                                                                               <  1:3372937‑M1/101‑16 (MQ=255)
                        gctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCaaa                                                                              >  2:3373684‑M2/77‑101 (MQ=255)
                        gctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCaaa                                                                              <  1:3558432‑M2/25‑1 (MQ=255)
                          tttactttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaactta                                                                            >  1:3373006‑M1/1‑83 (MQ=255)
                          tgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAAt                                                                            <  1:3558362‑M2/27‑1 (MQ=255)
                           ttactttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttat                                                                           <  1:3372760‑M1/101‑20 (MQ=255)
                           ttactttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttat                                                                           <  1:3372859‑M1/101‑20 (MQ=255)
                              ctttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattga                                                                        <  1:3372617‑M1/101‑23 (MQ=255)
                                  aTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagt                                                                    >  1:3373018‑M1/1‑75 (MQ=255)
                                   tCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtg                                                                   >  2:3373019‑M1/1‑74 (MQ=255)
                                     gaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgaga                                                                 >  2:3373434‑M2/64‑101 (MQ=255)
                                     aCGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgtt                                                                 <  2:3372594‑M1/101‑30 (MQ=255)
                                      aagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGaa                                                                >  1:3373625‑M2/63‑101 (MQ=255)
                                       agccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgag                                                                  >  2:3373934‑M2/62‑98 (MQ=255)
                                         cTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttat                                                             <  1:3372808‑M1/101‑34 (MQ=255)
                                         cTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttat                                                             >  2:3373029‑M1/1‑68 (MQ=255)
                                          ttCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatg                                                            <  2:3372857‑M1/101‑35 (MQ=255)
                                          ttCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatg                                                            >  1:1942239‑M1/1‑67 (MQ=255)
                                          cactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAata                                                            <  1:39719‑M2/43‑1 (MQ=255)
                                          cactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAata                                                            <  2:39818‑M2/43‑1 (MQ=255)
                                           tCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgt                                                           <  2:3372716‑M1/101‑36 (MQ=255)
                                           tCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgt                                                           <  1:3372741‑M1/101‑36 (MQ=255)
                                            ctggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGTGAAAATAtt                                                          <  1:39420‑M2/45‑1 (MQ=255)
                                                agcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTg                                                      <  1:39506‑M2/49‑1 (MQ=255)
                                                  aaaGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagata                                                    >  2:1942265‑M1/1‑59 (MQ=255)
                                                   acctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCaaa                                                                              >  1:1941394‑M2/50‑74 (MQ=255)
                                                   acctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATg                                                   <  1:39243‑M2/52‑1 (MQ=255)
                                                   aaGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataa                                                   >  1:3373042‑M1/1‑58 (MQ=255)
                                                    cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGc                                                  <  2:39566‑M2/53‑1 (MQ=255)
                                                    aGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataat                                                  <  2:3372533‑M1/101‑45 (MQ=255)
                                                    aGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataat                                                  <  2:3372791‑M1/101‑45 (MQ=255)
                                                    aGCGTTCTGTTCTGGATTTTGACGGCGGTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaact                                                                              >  2:3373044‑M1/1‑57 (MQ=255)
                                                      tcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGa                                                >  2:3373139‑M2/47‑101 (MQ=255)
                                                       gttctgttctGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcc                                               >  1:1942334‑M1/1‑54 (MQ=255)
                                                       caaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGAt                                               <  1:39655‑M2/56‑1 (MQ=255)
                                                         aaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATga                                             >  2:3373131‑M2/44‑101 (MQ=255)
                                                          aaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgac                                            >  1:3373126‑M2/43‑101 (MQ=255)
                                                           tgttctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgat                                           <  2:1943008‑M1/101‑52 (MQ=255)
                                                           tgttctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgat                                           <  1:1943010‑M1/101‑52 (MQ=255)
                                                            acaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactg                                          <  1:39849‑M2/61‑1 (MQ=255)
                                                             ttctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatga                                         >  1:523124‑M1/1‑48 (MQ=255)
                                                             caccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactgg                                         >  2:3373118‑M2/40‑101 (MQ=255)
                                                               ccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactggac                                       >  1:3373108‑M2/38‑101 (MQ=255)
                                                                catcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactggact                                      <  2:39876‑M2/65‑1 (MQ=255)
                                                                  gATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                    >  2:523115‑M1/1‑43 (MQ=255)
                                                                   catacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgaga                                                                 >  2:3373148‑M2/34‑71 (MQ=255)
                                                                   catacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGta                                   <  1:40163‑M2/68‑1 (MQ=255)
                                                                     tttGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                                 <  2:1942965‑M1/101‑62 (MQ=255)
                                                                      ttGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcat                                <  1:1942962‑M1/101‑63 (MQ=255)
                                                                      acactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATaa                                <  2:39623‑M2/71‑1 (MQ=255)
                                                                      acactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATaa                                >  2:3373100‑M2/31‑101 (MQ=255)
                                                                       tGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatg                               >  1:523109‑M1/1‑38 (MQ=255)
                                                                       cactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAAt                               <  1:39769‑M2/72‑1 (MQ=255)
                                                                        actaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATc                              <  1:40034‑M2/73‑1 (MQ=255)
                                                                        actaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATc                              >  1:3373091‑M2/29‑101 (MQ=255)
                                                                         ctaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCg                             >  2:3373085‑M2/28‑101 (MQ=255)
                                                                           ggCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagc                           <  1:1942934‑M1/101‑68 (MQ=255)
                                                                           aaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCa                           <  2:39935‑M2/76‑1 (MQ=255)
                                                                              tcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCaaaa                        <  2:40109‑M2/79‑1 (MQ=255)
                                                                              tcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCaaaa                        >  1:3373078‑M2/23‑101 (MQ=255)
                                                                                 aaTTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc                     <  1:1942914‑M1/101‑74 (MQ=255)
                                                                                  taagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAAtttt                    >  2:3373070‑M2/19‑101 (MQ=255)
                                                                                   aagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTa                   <  2:39671‑M2/84‑1 (MQ=255)
                                                                                    tttGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgat                  <  2:1942910‑M1/101‑77 (MQ=255)
                                                                                     ttGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt                 >  1:6544488‑M1/1‑24 (MQ=255)
                                                                                     gttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAAc                 >  1:3373064‑M2/16‑101 (MQ=255)
                                                                                      tGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattt                >  2:6543817‑M1/1‑23 (MQ=255)
                                                                                       gCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatttt               >  1:6544113‑M1/1‑22 (MQ=255)
                                                                                            gcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGaa          <  2:39753‑M2/93‑1 (MQ=255)
                                                                                            gcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGaa          >  2:3373054‑M2/9‑101 (MQ=255)
                                                                                            gcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGaa          <  1:39913‑M2/93‑1 (MQ=255)
                                                                                            gcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGaa          <  1:39560‑M2/93‑1 (MQ=255)
                                                                                             tatGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaac         >  2:6544169‑M1/1‑16 (MQ=255)
                                                                                                  ccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGa     <  2:39460‑M2/98‑1 (MQ=255)
                                                                                                    |                                                                                                    
TCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAGCA  >  NC_000913/1872931‑1873131

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