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)268 (1.040) 142/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)243 (0.940) 129/194 0.6 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 510 total aligned reads displayed.

TCGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAGCA  >  NC_000913/1872931‑1873131
                                                                                                    |                                                                                                    
gtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccc                                                                                                      >  2:2578040‑M2/101‑101 (MQ=255)
gtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccc                                                                                                      <  2:2577488‑M2/1‑1 (MQ=255)
 cGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCg                                                                                                     >  2:2442612/1‑101 (MQ=255)
 cGTGTCTCGTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCg                                                                                                     <  1:2442430/101‑1 (MQ=255)
   agtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGa                                                                                                   >  1:2442955‑M2/98‑101 (MQ=255)
     tagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATg                                                                                                 <  1:2577779‑M2/6‑1 (MQ=255)
       gctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCa                                                                                               <  1:2577529‑M2/8‑1 (MQ=255)
         gTTATCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACg                                                                                             >  1:2442617‑M1/1‑100 (MQ=255)
            acaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTcc                                                                                          <  1:2577575‑M2/13‑1 (MQ=255)
             tCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaa                                                                                         >  2:2442622‑M1/1‑96 (MQ=255)
             tCCATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaa                                                                                         <  2:2442422‑M1/101‑6 (MQ=255)
               ggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGtt                                                                                       <  2:2577632‑M2/16‑1 (MQ=255)
               ggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGtt                                                                                       <  2:2577584‑M2/16‑1 (MQ=255)
               cATTGATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatg                                                                                       >  1:1374198‑M1/1‑94 (MQ=255)
                   gATATTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactc                                                                                   <  2:2442366‑M1/101‑12 (MQ=255)
                    atatTCTTTACTTTATCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactcc                                                                                  >  1:2442624‑M1/1‑89 (MQ=255)
                       agctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCaa                                                                               >  1:2443179‑M2/78‑101 (MQ=255)
                          tttactttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaactta                                                                            <  1:2442412‑M1/101‑19 (MQ=255)
                          tgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAAt                                                                            >  1:2442876‑M2/75‑101 (MQ=255)
                           ttactttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttat                                                                           >  2:1374134‑M1/1‑82 (MQ=255)
                            atagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAAtgt                                                                          >  2:2443178‑M2/73‑101 (MQ=255)
                            atagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAAtgt                                                                          <  1:2577706‑M2/29‑1 (MQ=255)
                             actttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattg                                                                         <  1:2442492‑M1/101‑22 (MQ=255)
                                ttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgata                                                                      <  2:2442238‑M1/101‑25 (MQ=255)
                                ttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgata                                                                      >  1:2442634‑M1/1‑77 (MQ=255)
                                ttaTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgata                                                                      >  2:1374353‑M1/1‑77 (MQ=255)
                                aaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTcg                                                                      <  2:27479‑M2/33‑1 (MQ=255)
                                 taTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatag                                                                     <  2:2442446‑M1/101‑26 (MQ=255)
                                 aacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTcgc                                                                     <  1:27260‑M2/34‑1 (MQ=255)
                                  aTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagt                                                                    <  1:2442332‑M1/101‑27 (MQ=255)
                                  aTCACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagt                                                                    <  2:2442466‑M1/101‑27 (MQ=255)
                                    cACGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgt                                                                  >  2:2442639‑M1/1‑73 (MQ=255)
                                    agaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgag                                                                  >  1:2443067‑M2/65‑101 (MQ=255)
                                    agaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgag                                                                  >  2:2443050‑M2/65‑101 (MQ=255)
                                     gaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCgaga                                                                 <  1:27657‑M2/38‑1 (MQ=255)
                                      cGTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttt                                                                >  2:1374210‑M1/1‑71 (MQ=255)
                                      aagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGaa                                                                <  2:27175‑M2/39‑1 (MQ=255)
                                      aagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTCTATACCAAAATGTGTCGCGAGaa                                                                <  1:27440‑M2/39‑1 (MQ=255)
                                       gTCTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgtttt                                                               >  1:1374282‑M1/1‑70 (MQ=255)
                                         cTTCGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttat                                                             >  1:2442642‑M1/1‑68 (MQ=255)
                                          cactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAata                                                            >  1:2443011‑M2/59‑101 (MQ=255)
                                            cGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgtt                                                          <  2:2442328‑M1/101‑37 (MQ=255)
                                            cGTGGAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgtt                                                          <  1:2442454‑M1/101‑37 (MQ=255)
                                               ggAAAAGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcag                                                       >  2:2442650‑M1/1‑62 (MQ=255)
                                               gagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGt                                                       >  2:2442985‑M2/54‑101 (MQ=255)
                                                agcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTg                                                      <  1:27474‑M2/49‑1 (MQ=255)
                                                agcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTg                                                      <  1:27560‑M2/49‑1 (MQ=255)
                                                 aaaaGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagat                                                     >  1:2442652‑M1/1‑60 (MQ=255)
                                                  cacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGAt                                                    >  1:2442736‑M2/51‑101 (MQ=255)
                                                    cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGc                                                  <  2:27543‑M2/53‑1 (MQ=255)
                                                    aGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataat                                                  <  2:2442518‑M1/101‑45 (MQ=255)
                                                    aGCGTTCTGTTCTGGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataat                                                  <  2:2442490‑M1/101‑45 (MQ=255)
                                                     ctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCg                                                 >  2:2442729‑M2/48‑101 (MQ=255)
                                                     ctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGAGGCg                                                 <  2:27367‑M2/54‑1 (MQ=255)
                                                      tcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGa                                                >  1:2442726‑M2/47‑101 (MQ=255)
                                                        aaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATg                                              <  1:27096‑M2/57‑1 (MQ=255)
                                                         aaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATga                                             >  2:2442717‑M2/44‑101 (MQ=255)
                                                         aaaacaccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATga                                             >  1:2442718‑M2/44‑101 (MQ=255)
                                                             ttctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatga                                         >  1:371296‑M1/1‑48 (MQ=255)
                                                             caccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactgg                                         >  2:2442710‑M2/40‑101 (MQ=255)
                                                              tctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgac                                        <  2:1374662‑M1/101‑55 (MQ=255)
                                                               ctgGATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgact                                       >  2:371293‑M1/1‑46 (MQ=255)
                                                               ccatcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactggac                                       <  2:27261‑M2/64‑1 (MQ=255)
                                                                catcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactggact                                      <  1:27213‑M2/65‑1 (MQ=255)
                                                                 aacatacactaaatcagtaagttggcagcatcaccCGGAt                                                                                                  <  2:1373479‑M2/5‑1 (MQ=255)
                                                                  tcatacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGt                                    >  2:2442703‑M2/35‑101 (MQ=255)
                                                                  gATTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                                    <  1:1374651‑M1/101‑59 (MQ=255)
                                                                   catacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGta                                   >  1:2442701‑M2/34‑101 (MQ=255)
                                                                   catacactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGta                                   <  2:27414‑M2/68‑1 (MQ=255)
                                                                   aTTTTGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                   >  1:371288‑M1/1‑42 (MQ=255)
                                                                     tttGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                                 <  2:1374645‑M1/101‑62 (MQ=255)
                                                                      ttGACGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcat                                >  1:371283‑M1/1‑39 (MQ=255)
                                                                      acactaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATaa                                >  2:2442696‑M2/31‑101 (MQ=255)
                                                                        actaaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATc                              <  1:27594‑M2/73‑1 (MQ=255)
                                                                         aCGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgca                             <  2:1374636‑M1/101‑66 (MQ=255)
                                                                          taaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGc                            >  1:2442687‑M2/27‑101 (MQ=255)
                                                                          cGGCGTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcag                            <  1:1374630‑M1/101‑67 (MQ=255)
                                                                           aaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGaa                                                                <  1:1373483‑M2/39‑1 (MQ=255)
                                                                           aaatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCa                           >  2:2442686‑M2/26‑101 (MQ=255)
                                                                            aatcagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCaa                          <  1:27368‑M2/77‑1 (MQ=255)
                                                                              gTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcc                        >  1:4683529‑M1/1‑31 (MQ=255)
                                                                              gTTAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcc                        >  2:1374159‑M1/1‑31 (MQ=255)
                                                                               ttAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcca                       >  2:4683519‑M1/1‑30 (MQ=255)
                                                                               cagtaagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAAt                       >  2:2442679‑M2/22‑101 (MQ=255)
                                                                                tAATTTTGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccac                      >  1:4683606‑M1/1‑29 (MQ=255)
                                                                                  taagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAAtttt                    <  2:27617‑M2/83‑1 (MQ=255)
                                                                                   aagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTa                   <  2:27652‑M2/84‑1 (MQ=255)
                                                                                   aagttggcagcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTa                   >  1:2442675‑M2/18‑101 (MQ=255)
                                                                                    tttGCTGAATATGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgat                  <  2:1374609‑M1/101‑77 (MQ=255)
                                                                                    tttGCTGAATATGGAGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgc                              >  2:4683275‑M1/1‑25 (MQ=255)
                                                                                      tGCTGAATATGGTGCGGATGCACgggaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattt                >  1:2576978‑M1/1‑23 (MQ=255)
                                                                                            gcatcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGaa          >  2:2442667‑M2/9‑101 (MQ=255)
                                                                                            atatGGTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaa          <  1:1374592‑M1/101‑85 (MQ=255)
                                                                                             catcaccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAAc         >  1:2442666‑M2/8‑101 (MQ=255)
                                                                                                ggTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcag                                                       <  1:1374692‑M1/52‑40 (MQ=255)
                                                                                                 gTGCGGATGCACggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgaca     <  2:1374579‑M1/101‑90 (MQ=255)
                                                                                                  ccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATgactgga                                        <  2:27317‑M2/63‑1 (MQ=255)
                                                                                                  ccCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAg    <  2:27837‑M2/99‑1 (MQ=255)
                                                                                                   cCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAGc   >  2:2442657‑M2/2‑101 (MQ=255)
                                                                                                   cCGGATGCACTTCCGTTTATACCAAAATGTGTCGCGAGAAAATATTAGTGATGCGATGACTGGACTGTATAATCGCAAAATTTTAACCCCTGAACTGGAGc   <  1:27769‑M2/100‑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.