Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 3,765,971 IS2 (+) +5 bp coding (3789‑3793/4134 nt) rhsA → Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3186096 =NA (NA)158 (1.010) 97/198 0.2 96.2% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 3765975 6 (0.040)coding (3793/4134 nt) rhsA Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
* ? NC_000913 = 3187426NA (NA)158 (1.010) 84/198 0.8 96.2% noncoding (1331/1331 nt) IS2 repeat region
?NC_000913 3765971 = 6 (0.040)coding (3789/4134 nt) rhsA Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

Only 100 of 203 total aligned reads displayed.

TTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATT  >  NC_000913/3185998‑3186192
                                                                                                  |                                                                                                
ttAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatta                                                                                                 >  1:24686/1‑100 (MQ=255)
ttAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGttatta                                                                                                 >  2:2991495/1‑100 (MQ=255)
 tAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAg                                                                                                >  1:2997102/1‑100 (MQ=255)
 naaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAg                                                                                                >  1:1157375/2‑100 (MQ=255)
  aaCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGa                                                                                               >  1:854923/1‑100 (MQ=255)
   aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              <  2:177518/100‑1 (MQ=255)
   aCGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGAc                                                                                              <  1:2522805/100‑1 (MQ=255)
    cGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACt                                                                                             <  1:2811860/100‑1 (MQ=255)
     ggATTAATGGTAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTg                                                                                            >  1:3984655/1‑100 (MQ=255)
     ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTg                                                                                            <  1:3903313/100‑1 (MQ=255)
     ggATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTg                                                                                            >  2:2369598/1‑100 (MQ=255)
       cgatctGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGc                                                                                          <  2:2021909/95‑1 (MQ=255)
        ttAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGcc                                                                                         <  1:2892627/100‑1 (MQ=255)
         tAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGccc                                                                                        <  2:1455641/100‑1 (MQ=255)
          aaTGATAAGTGGATCAGATGGATAAAAAAATAAAATAAACCACAAGAAGCGGAATAATAAAAAATATATAACGTTCATTATCGGGTAGTAGACTGGcccc                                                                                       >  2:3987496/1‑100 (MQ=255)
            tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCTCCCt                                                                                     >  2:463071/1‑100 (MQ=255)
            tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCt                                                                                     >  1:516437/1‑100 (MQ=255)
            tGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCt                                                                                     >  1:2147458/1‑100 (MQ=255)
             gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                    >  1:1829734/1‑100 (MQ=255)
             gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                    <  1:1398534/100‑1 (MQ=255)
             gATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTg                                                                                    <  1:3277735/100‑1 (MQ=255)
              aTAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGa                                                                                   <  2:552422/100‑1 (MQ=255)
               tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGaa                                                                                  <  2:1856427/100‑1 (MQ=255)
               tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGaa                                                                                  <  2:3453144/100‑1 (MQ=255)
               tAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGaa                                                                                  <  1:1481372/100‑1 (MQ=255)
                 aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                >  2:1573876/1‑100 (MQ=255)
                 aGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAAtc                                                                                >  2:688483/1‑100 (MQ=255)
                    gtaTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTTAAAATATTTATTGTTCATTATCTGTTCTTAGACTGGCCCCCTGAATCTcc                                                                             >  2:360258/3‑100 (MQ=255)
                    ggATCAGATGTATAAAAAATTAAAATTAACCATAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                             <  2:2811460/100‑1 (MQ=255)
                    ggATAAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTcc                                                                             <  1:3306797/100‑1 (MQ=255)
                    gcaTCAGATGTCTAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTTGCCCCCTGAATCTcc                                                                             <  1:188932/98‑1 (MQ=255)
                     gATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCa                                                                            <  2:941017/100‑1 (MQ=255)
                     gATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCa                                                                            <  2:35605/100‑1 (MQ=255)
                      aTCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAg                                                                           <  1:2609925/100‑1 (MQ=255)
                       tCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGa                                                                          <  1:1541570/100‑1 (MQ=255)
                       tCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGa                                                                          <  2:647360/100‑1 (MQ=255)
                        cAGATGTATAAAAAATTAAAATTAACCACAATAAGCGTATTGATAAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGAc                                                                         >  2:3655563/1‑100 (MQ=255)
                         aGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACa                                                                        >  2:2393506/1‑100 (MQ=255)
                           aTGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAAc                                                                      <  2:3612913/100‑1 (MQ=255)
                           aTGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCACCTGAATCTCCAGACAAc                                                                      <  2:451508/100‑1 (MQ=255)
                             gTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCa                                                                    >  2:3901231/1‑100 (MQ=255)
                               ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                  >  2:1477834/1‑100 (MQ=255)
                               ataAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAat                                                                  >  2:2970236/1‑100 (MQ=255)
                                taAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                 >  1:392883/1‑100 (MQ=255)
                                taAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                 >  2:633050/1‑100 (MQ=255)
                                taAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                 >  1:485951/1‑100 (MQ=255)
                                 aaaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                >  1:1319319/1‑100 (MQ=255)
                                 aaaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                >  1:2328749/1‑100 (MQ=255)
                                  aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                               >  1:1207315/1‑100 (MQ=255)
                                  aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                               <  2:2844179/100‑1 (MQ=255)
                                  aaaaaTTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATc                                                               >  1:3636379/1‑100 (MQ=255)
                                    aaattaaaattaACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAc                                                             <  2:373644/100‑1 (MQ=255)
                                      attaaaattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCCCCAGACAACCAATATCACtt                                                           <  2:2565922/100‑1 (MQ=255)
                                       ttaaaattaaCCACAATAAGCAAAATGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTa                                                          <  2:1509948/100‑1 (MQ=255)
                                          aaattgaCTAAAAAACGCGAATGTATTAAAAATATATATTGTTCACTATCCGTTATTAGACGAGCGCCCTGAATCTACAGTCAACCACTACCACACAAAt                                                       <  2:780778/100‑1 (MQ=255)
                                           aattaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATa                                                      >  1:2071472/1‑100 (MQ=255)
                                            attaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATaa                                                     >  1:3125363/1‑100 (MQ=255)
                                                aCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGa                                                 >  1:1467175/1‑100 (MQ=255)
                                                  cacaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATa                                               <  1:3782360/100‑1 (MQ=255)
                                                     gaTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCAGCCTGAATCTTCAGACAACGAATATCACTTATTTAAGTGATattg                                            <  2:2039244/99‑4 (MQ=255)
                                                       tAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattggt                                          <  2:2952436/100‑6 (MQ=25)
                                                       tAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTATTTAAGTGATattggt                                          <  2:22145/100‑6 (MQ=25)
                                                        aaGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTa                                         <  1:354427/100‑1 (MQ=255)
                                                          gCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                       <  2:836880/100‑1 (MQ=255)
                                                           cGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATa                                      >  2:2669909/1‑100 (MQ=255)
                                                             aaTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                    >  1:1037390/1‑100 (MQ=255)
                                                             aaTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACt                                    <  1:3386888/100‑1 (MQ=255)
                                                              aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTa                                   >  2:3513782/1‑100 (MQ=255)
                                                               ttGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                                  <  1:385133/100‑1 (MQ=255)
                                                               ttGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                                  <  2:3032458/100‑1 (MQ=255)
                                                                tGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                                 >  1:4017035/1‑100 (MQ=255)
                                                                 gATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                                >  2:2798194/1‑100 (MQ=255)
                                                                   ttAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                              >  1:536508/1‑100 (MQ=255)
                                                                   ttAAAAATATTTACTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                              <  2:3137926/100‑1 (MQ=255)
                                                                    tAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                             <  1:3706100/100‑1 (MQ=255)
                                                                     aaaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTa                            <  2:2831386/100‑1 (MQ=255)
                                                                      aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                           <  1:1432858/100‑1 (MQ=255)
                                                                        aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGAc                         <  2:1703084/100‑1 (MQ=255)
                                                                          tctttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTCATAGTCTTCATACTAGTTTTCAGACTa                       >  2:1221599/3‑100 (MQ=17)
                                                                             ttattGTTCATTATCCGTTATTAGACTTGCTCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTc                    <  2:385734/100‑1 (MQ=255)
                                                                              tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCa                   >  1:2252105/1‑100 (MQ=255)
                                                                                 tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                >  1:767730/1‑100 (MQ=255)
                                                                                  gtgcATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:1015791/97‑1 (MQ=255)
                                                                                  gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:918755/100‑1 (MQ=255)
                                                                                  gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               >  2:3219976/1‑100 (MQ=255)
                                                                                  cttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTgg               <  2:1678755/99‑1 (MQ=255)
                                                                                   ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGga              <  1:6564/100‑1 (MQ=255)
                                                                                    tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag             <  1:1057885/100‑1 (MQ=255)
                                                                                    tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgag             >  1:1706878/1‑100 (MQ=255)
                                                                                     cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga            <  2:1585260/100‑1 (MQ=255)
                                                                                      aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa           >  1:2761912/1‑100 (MQ=255)
                                                                                      aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGaa           >  1:2104410/1‑100 (MQ=255)
                                                                                       ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc          <  2:3968901/100‑1 (MQ=255)
                                                                                       ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAc          <  1:554499/100‑1 (MQ=255)
                                                                                        tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACa         <  1:2385219/100‑1 (MQ=40)
                                                                                          tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAga       <  2:2313130/100‑1 (MQ=38)
                                                                                          tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAga       >  2:1625380/1‑100 (MQ=38)
                                                                                           ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgat      <  1:2391247/100‑1 (MQ=37)
                                                                                            cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTACATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatg     <  2:2064744/100‑1 (MQ=25)
                                                                                              ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAgatgat   <  1:1931217/100‑1 (MQ=33)
                                                                                                  |                                                                                                
TTAACGGATTAATGATAAGTGGATCAGATGTATAAAAAATTAAAATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATT  >  NC_000913/3185998‑3186192

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