Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,410,998 IS2 (+) +5 bp 100% coding (719‑723/966 nt) dusB → tRNA‑dihydrouridine synthase B

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3186096 =NA (NA)164 (1.250) 96/178 NT 99.4% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 3411002 1 (0.010)coding (723/966 nt) dusB tRNA‑dihydrouridine synthase B
* ? NC_000913 = 3187426NA (NA)129 (0.980) 81/178 NT 99.2% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3410998 = 1 (0.010)coding (719/966 nt) dusB tRNA‑dihydrouridine synthase B

Only 100 of 188 total aligned reads displayed.

ATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTG  >  NC_000913/3185997‑3186193
                                                                                                   |                                                                                                 
atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                   >  1:454537/1‑100 (MQ=255)
 ttAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatta                                                                                                  <  1:935831/100‑1 (MQ=255)
 ttAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatta                                                                                                  <  2:2453307/100‑1 (MQ=255)
  tAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAg                                                                                                 <  2:362080/100‑1 (MQ=255)
    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                               >  2:1259594/1‑100 (MQ=255)
     cGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACt                                                                                              >  1:585880/1‑100 (MQ=255)
     cGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCAGTATCCGTTATTAGACt                                                                                              >  1:2187086/1‑100 (MQ=255)
      ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACt                                                                                              >  1:3462002/1‑99 (MQ=255)
      ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACt                                                                                              <  1:2113054/99‑1 (MQ=255)
      ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTg                                                                                             >  1:2691710/1‑100 (MQ=255)
       gATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTgg                                                                                            <  1:1357374/100‑1 (MQ=255)
            aTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccccc                                                                                       <  1:3628586/100‑1 (MQ=255)
            aTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccccc                                                                                       <  2:550583/100‑1 (MQ=255)
            aTGATAAGTGGATCACATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccccc                                                                                       <  2:3670951/100‑1 (MQ=255)
              gATAAGTTCATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                     <  1:3498592/100‑1 (MQ=255)
              gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                     <  2:2804549/100‑1 (MQ=255)
              gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                     <  1:369982/100‑1 (MQ=255)
              gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                     <  2:3090572/100‑1 (MQ=255)
              gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                     <  1:1607331/100‑1 (MQ=255)
                    tGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtctc                                                                               >  1:161918/1‑100 (MQ=255)
                     ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                              <  2:1370268/100‑1 (MQ=255)
                     ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                              >  1:1046714/1‑100 (MQ=255)
                      gATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCa                                                                             >  1:3289397/1‑100 (MQ=255)
                       aTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAg                                                                            <  1:526198/100‑1 (MQ=255)
                       aTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAg                                                                            <  2:1463630/100‑1 (MQ=255)
                          aGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACa                                                                         <  2:2721328/100‑1 (MQ=255)
                          aGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATGGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACa                                                                         <  2:2652964/100‑1 (MQ=255)
                            atctataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAc                                                                       <  1:1715416/97‑1 (MQ=255)
                              gTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAGTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                     >  1:2317985/1‑100 (MQ=255)
                               tataAAAAATTAAAATTAACCACAATAAGCGAATTGAATAAAAATATTTATTGTTCATTATCCGTTATTAGACTGACCCCCTGAATTTCCATACaa                                                                        >  1:2434529/1‑96 (MQ=255)
                                ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAcagat                                                                   >  2:1715423/1‑96 (MQ=255)
                                ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                   >  1:3443895/1‑100 (MQ=255)
                                ataAAAAATTAAAATTAACCACAATAAGCGAATTGAATAAAAATATTTATTGTTCATTATCCATTATTAGACTGTCCCCCTGAATCTCCAGACAACCa                                                                     >  1:3351821/1‑98 (MQ=255)
                                 taAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAACATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                  <  1:2629362/100‑1 (MQ=255)
                                 taAAAAATTAAAATTAACCACAATAAGCGAATTAATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTAGCCCCCTGAATCTCCAGACATCCa                                                                     >  1:1214651/1‑97 (MQ=255)
                                   aaaaaTTAAAATTTACCACAATAAGCGAATTGATTAAACATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                                >  1:2301901/1‑100 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                              <  2:2812591/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                              <  1:886463/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                              <  1:813051/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                              <  2:2240827/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                              <  2:3417880/100‑1 (MQ=255)
                                      aattaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACt                                                             <  1:1948567/100‑1 (MQ=255)
                                        ttaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTa                                                           >  1:782237/1‑100 (MQ=255)
                                         taaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTaa                                                          >  1:2708739/1‑100 (MQ=255)
                                           aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                        <  2:3098285/100‑1 (MQ=255)
                                              ttaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAg                                                     <  2:3267305/100‑1 (MQ=255)
                                                agccACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                   <  1:2267598/98‑1 (MQ=255)
                                                aaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                   <  2:1922649/100‑1 (MQ=255)
                                                    acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                               >  2:635071/1‑100 (MQ=255)
                                                     caATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCACCCTGAATTTCCAGACAACCAATATCACTTAAATAAGTGATAGt                                              >  1:141082/1‑100 (MQ=255)
                                                      aaTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattg                                             >  2:934810/1‑97 (MQ=255)
                                                       aTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCt                                            >  1:967377/1‑100 (MQ=255)
                                                        tAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                           <  1:1117887/100‑1 (MQ=255)
                                                         aaGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattggtt                                          <  1:934805/100‑7 (MQ=25)
                                                          aGCGAATTGATTAAAAATATTTATTGTGCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattggttg                                         <  2:1040122/100‑8 (MQ=25)
                                                           gCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                        <  2:3185704/100‑1 (MQ=255)
                                                           gCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                        <  2:2795308/100‑1 (MQ=255)
                                                            cGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                       <  2:161920/100‑1 (MQ=255)
                                                             gAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                        <  1:353874/98‑1 (MQ=255)
                                                             gAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAc                                      >  2:1207691/1‑100 (MQ=255)
                                                                 tGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                                  <  2:2737401/100‑1 (MQ=255)
                                                                  gATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                 >  2:3580779/1‑100 (MQ=255)
                                                                   aTTAAAAATATTTATTGTTCATTATCCGTTATTAGACTAACATTCTTAATCTCTTGACAACCACTAACACTTCACTAAGTGTATTTCTTAATACTAgct                                 >  1:2050383/1‑97 (MQ=255)
                                                                    ttAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                               >  1:1901113/1‑100 (MQ=255)
                                                                     tAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                              >  2:2670201/1‑100 (MQ=255)
                                                                     tAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                              >  2:2059626/1‑100 (MQ=255)
                                                                     tAAAAATATTTATTGTTCATTATCCGTTATTAGACTGACCCCCTGAATCTCTAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                              >  1:3285799/1‑100 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            <  1:2174909/100‑1 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            >  2:1873170/1‑100 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            <  2:706541/100‑1 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            >  1:2839970/1‑100 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            >  1:729009/1‑100 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCACCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                            <  1:2558441/100‑1 (MQ=255)
                                                                            atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                       <  2:599537/100‑1 (MQ=255)
                                                                            atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                       <  2:352991/100‑1 (MQ=255)
                                                                             tttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTaat                      >  1:851284/1‑98 (MQ=255)
                                                                               tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCa                    <  1:3075419/100‑1 (MQ=255)
                                                                                  tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                 >  1:538567/1‑100 (MQ=255)
                                                                                  tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                 <  2:3443919/100‑1 (MQ=255)
                                                                                  tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                 >  2:1563995/1‑100 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                <  2:1994741/100‑1 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg                >  2:914413/1‑100 (MQ=255)
                                                                                    ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga               <  2:819648/100‑1 (MQ=255)
                                                                                    ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga               <  2:2187096/100‑1 (MQ=255)
                                                                                    ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga               <  1:1683512/100‑1 (MQ=255)
                                                                                     tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag              >  1:2504106/1‑100 (MQ=255)
                                                                                     tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag              >  2:2119432/1‑100 (MQ=255)
                                                                                     tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag              >  2:1421444/1‑100 (MQ=255)
                                                                                     tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag              <  1:861476/100‑1 (MQ=255)
                                                                                      cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga             <  1:1143318/100‑1 (MQ=255)
                                                                                        ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc           <  2:1158132/100‑1 (MQ=40)
                                                                                        ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc           <  1:3136224/100‑1 (MQ=40)
                                                                                        ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc           <  2:1754659/100‑1 (MQ=40)
                                                                                          aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAg         >  2:1271593/1‑100 (MQ=38)
                                                                                           tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAga        >  2:2725554/1‑100 (MQ=38)
                                                                                            ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat       >  2:2345339/1‑100 (MQ=35)
                                                                                             cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg      >  1:123301/1‑100 (MQ=35)
                                                                                               ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatgat    <  1:3522674/100‑1 (MQ=32)
                                                                                                tattaGACTGGCCCCCTTAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt   <  1:1055635/100‑1 (MQ=14)
                                                                                                 attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgattg  <  2:302623/100‑1 (MQ=11)
                                                                                                   |                                                                                                 
ATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTG  >  NC_000913/3185997‑3186193

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