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)139 (1.120) 90/178 NT 98.5% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 3411002 2 (0.020)coding (723/966 nt) dusB tRNA‑dihydrouridine synthase B
* ? NC_000913 = 3187426NA (NA)139 (1.120) 95/178 NT 98.5% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3410998 = 2 (0.020)coding (719/966 nt) dusB tRNA‑dihydrouridine synthase B

Only 100 of 148 total aligned reads displayed.

ATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATT  >  NC_000913/3185997‑3186192
                                                                                                   |                                                                                                
atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                  >  1:2536694/1‑100 (MQ=255)
atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                  <  1:1152512/100‑1 (MQ=255)
atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                  <  1:3055155/100‑1 (MQ=255)
atTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatt                                                                                                  >  2:1982940/1‑100 (MQ=255)
   aaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGa                                                                                               >  2:936831/1‑100 (MQ=255)
   aaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGa                                                                                               >  1:222020/1‑100 (MQ=255)
    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              <  1:3143130/100‑1 (MQ=255)
    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              >  2:2280863/1‑100 (MQ=255)
    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              <  1:1620339/100‑1 (MQ=255)
    aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              <  2:1897817/100‑1 (MQ=255)
     cGGATTAAGGATAAGTGGATCAGATGTATAAAAAAGTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACt                                                                                             >  1:2337555/1‑100 (MQ=255)
      ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTg                                                                                            <  2:2091495/100‑1 (MQ=255)
       gATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTgg                                                                                           >  1:2645684/1‑100 (MQ=255)
       gATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTgg                                                                                           <  2:1681719/100‑1 (MQ=255)
        aTTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGc                                                                                          <  2:352373/100‑1 (MQ=255)
           aaTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGcccc                                                                                       >  2:967815/1‑100 (MQ=255)
            aTGATAAGTGGATTAGATGTATAAAAAATTAAAATTAAGCACAATAAGAGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccccc                                                                                      >  1:2064100/1‑100 (MQ=255)
            aTGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccccc                                                                                      <  2:3059750/100‑1 (MQ=255)
             tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTAATTATCCGTTATTAGACTGGCCCCCt                                                                                     >  1:1318551/1‑100 (MQ=255)
              gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                    >  1:845498/1‑100 (MQ=255)
                tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGaa                                                                                  >  1:842563/1‑100 (MQ=255)
                 aaGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAt                                                                                 <  1:1625408/100‑1 (MQ=255)
                  aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                <  1:178589/100‑1 (MQ=255)
                  aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                <  2:145875/100‑1 (MQ=255)
                  aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                <  2:1365198/100‑1 (MQ=255)
                  aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                >  2:1162409/1‑100 (MQ=255)
                  aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                >  2:36152/1‑100 (MQ=255)
                    tGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtctc                                                                              >  1:3002057/1‑100 (MQ=255)
                     ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                             <  2:2171563/100‑1 (MQ=255)
                     ggATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                             >  2:849174/1‑100 (MQ=255)
                      gATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCa                                                                            <  1:949525/100‑1 (MQ=255)
                       aTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCan                                                                           <  1:1083112/100‑2 (MQ=255)
                             tGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAcc                                                                     >  1:1932060/1‑100 (MQ=255)
                              gTATAAAAAATTAAAATTAACCACAAGAAGCCAATTGATTAAAAAGATTTATTGTTGATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAc                                                                      <  1:1256702/98‑1 (MQ=38)
                                ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                  <  2:146393/100‑1 (MQ=255)
                                ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                  >  2:3206111/1‑100 (MQ=255)
                                 taAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                 >  2:1358217/1‑100 (MQ=255)
                                  aaaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                >  2:3078091/1‑100 (MQ=255)
                                   aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                               <  1:2723037/100‑1 (MQ=255)
                                   aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                               >  1:2959915/1‑100 (MQ=255)
                                   aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCACCCTGAATCTCCAGACAACCAATATc                                                               >  1:553294/1‑100 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                             <  1:550578/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                             <  1:2052586/100‑1 (MQ=255)
                                     aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                             <  2:814893/100‑1 (MQ=255)
                                       attaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                           <  1:950217/100‑1 (MQ=255)
                                        ttaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCACCTGAATCTCCAGACAACCAATATCACTTa                                                          >  1:2659702/1‑100 (MQ=255)
                                         taaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                           >  1:2877974/1‑98 (MQ=255)
                                           aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTattt                                                       <  1:1587273/100‑4 (MQ=255)
                                           aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                       >  1:3113469/1‑100 (MQ=255)
                                           aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                       <  1:2936827/100‑1 (MQ=255)
                                           aaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                       >  2:2665030/1‑100 (MQ=255)
                                               taaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGt                                                   <  1:358303/100‑1 (MQ=255)
                                                aaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCTCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                  >  1:1938414/1‑100 (MQ=255)
                                                aaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTg                                                  <  1:230031/100‑1 (MQ=255)
                                                  ccACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGAt                                                >  2:2200771/1‑100 (MQ=255)
                                                   cacaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                               >  2:1061229/1‑100 (MQ=255)
                                                    acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                              <  1:642952/100‑1 (MQ=255)
                                                    acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                              >  1:592900/1‑100 (MQ=255)
                                                    acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                              >  2:1278217/1‑100 (MQ=255)
                                                    acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCATGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                              <  1:2319212/100‑1 (MQ=255)
                                                     caATAATCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATatt                                             <  1:2184366/100‑3 (MQ=25)
                                                        tAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                          <  1:1369839/100‑1 (MQ=255)
                                                        tAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCtt                                          >  1:1780793/1‑100 (MQ=255)
                                                         aaGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                         <  1:1680929/100‑1 (MQ=255)
                                                          aGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                        <  2:2444265/100‑1 (MQ=255)
                                                          aGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTaa                                        >  1:1443519/1‑100 (MQ=255)
                                                           gCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                       <  1:2703543/100‑1 (MQ=255)
                                                            cGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                      <  1:450494/100‑1 (MQ=255)
                                                            cGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCATGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                       <  1:1055979/99‑1 (MQ=255)
                                                             gAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATGTACAGACAACAAATAACACTCAAATAAATGGTAGTCttg                                         >  2:3011529/1‑95 (MQ=25)
                                                              aaTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                    >  2:2110527/1‑100 (MQ=255)
                                                                  naTTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                >  2:2851850/2‑100 (MQ=255)
                                                                  gATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                >  2:2054387/1‑100 (MQ=255)
                                                                    ttAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAgat                                >  1:2308472/1‑96 (MQ=255)
                                                                      aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                            <  2:2342131/100‑1 (MQ=255)
                                                                      aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                            >  2:844604/1‑100 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                           <  1:238651/100‑1 (MQ=255)
                                                                       aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                           <  2:3200308/100‑1 (MQ=255)
                                                                         aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                         <  1:39537/100‑1 (MQ=255)
                                                                          atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACt                        <  2:1272755/100‑1 (MQ=255)
                                                                           tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTa                       >  2:2341826/1‑100 (MQ=255)
                                                                            atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAg                      <  1:2733250/100‑1 (MQ=255)
                                                                              ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt                     >  1:3302238/1‑99 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:592900/100‑1 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:2482903/100‑1 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               >  2:2962928/1‑100 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:1809465/100‑1 (MQ=255)
                                                                                   gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAACAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  1:2280846/100‑1 (MQ=255)
                                                                                       aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa           >  1:153545/1‑100 (MQ=255)
                                                                                       aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa           >  2:1926677/1‑100 (MQ=255)
                                                                                       aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa           <  1:2158457/100‑1 (MQ=255)
                                                                                         tATCCGTTATTAGACTGTCCCCCTGAATCTCCAGAAAACCAATATCTCTTAATTAAATGATAGACTTAATACTATTTTTTTGACTACTGATTTgatgaca         >  1:2633991/1‑95 (MQ=9)
                                                                                         tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga            <  1:2515932/97‑1 (MQ=40)
                                                                                          aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAg        >  1:1685793/1‑100 (MQ=38)
                                                                                            ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat      >  1:1437830/1‑100 (MQ=35)
                                                                                            ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat      >  1:2842274/1‑100 (MQ=35)
                                                                                             cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg     >  2:2748893/1‑100 (MQ=35)
                                                                                             cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg     >  1:795384/1‑100 (MQ=35)
                                                                                             cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg     >  1:1171706/1‑100 (MQ=35)
                                                                                                tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgatt  >  2:2883902/1‑100 (MQ=31)
                                                                                                   |                                                                                                
ATTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATT  >  NC_000913/3185997‑3186192

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.