Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2772854 =NA (NA)49 (1.150) 39/144 0.1 100% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)31 (0.730) 26/144 0.7 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGA  >  REL606/2774122‑2774271
                                                                          |                                                                           
cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGt                                                                            >  2:989795/1‑76 (MQ=31)
 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                           >  2:659237/1‑76 (MQ=32)
 gAGGCATCGTCTTAACGAGGCACCGAGGAGTCGCGTTATTCAGATGGNTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                           >  1:125621/1‑76 (MQ=11)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                          <  2:501425/76‑1 (MQ=34)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                         <  2:324853/76‑1 (MQ=35)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                         <  2:524016/76‑1 (MQ=35)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGCtttt                                                                         <  2:1386113/76‑1 (MQ=17)
    gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                        <  2:759301/76‑1 (MQ=38)
    gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                        <  1:866240/76‑1 (MQ=38)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  2:1233071/1‑76 (MQ=40)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  2:923770/1‑76 (MQ=40)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  2:1244370/1‑76 (MQ=40)
        cGTCTTAANGTGGCACCGANGCGTCGCATTCTTCAGATGGTTNAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    <  2:1177780/76‑1 (MQ=255)
        cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    <  2:1083595/76‑1 (MQ=255)
        cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    >  1:582770/1‑76 (MQ=255)
          tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGATAGCGCTTATNGGTTTTAccctca                                                                  >  1:180605/1‑73 (MQ=255)
           cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgc                                                                 >  1:373020/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                                >  2:35547/1‑76 (MQ=255)
             naaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGNTTTTACCCCCGCgt                                                               >  2:201060/2‑76 (MQ=255)
              aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                              >  1:780924/1‑76 (MQ=255)
               aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtgn                                                             <  1:290246/76‑2 (MQ=255)
               aCGAAGCACCGAGGCGTGGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACACCGGCgngg                                                             >  2:1321645/1‑76 (MQ=14)
                  aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                          >  2:617470/1‑76 (MQ=255)
                     cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTa                                                       <  2:526179/76‑1 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                    <  1:608467/76‑1 (MQ=255)
                         gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtc                                                   >  2:394055/1‑76 (MQ=255)
                           ggCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtctt                                                 <  1:1191275/76‑1 (MQ=255)
                            gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACttcttc                                                <  2:266983/76‑1 (MQ=255)
                             cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCg                                               <  1:394055/76‑1 (MQ=255)
                              gTCGCATTCTTCTGATGGTTCAACCCTTAAGTNAGAGCTTANGGGTTGTACCCCCGCGTGGNGGCTACTTCTTCGa                                              >  2:1154034/1‑76 (MQ=21)
                              gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                              >  2:967567/1‑76 (MQ=255)
                               tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTACTTCTTCGAt                                             >  1:564710/1‑76 (MQ=255)
                                cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                            >  1:1009501/1‑76 (MQ=255)
                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                          >  1:1154493/1‑76 (MQ=255)
                                  cATTCTTCAGATGGGTCAACCCTTGAGTCAGNGCTTATGGTTTTTANCCCCGNGTGGTGGCTACTTGTTCGATTCg                                          <  2:178418/76‑1 (MQ=25)
                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                         >  1:1100291/1‑76 (MQ=255)
                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                         >  1:376273/1‑76 (MQ=255)
                                    ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                        <  2:1100291/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                       <  1:991481/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                       <  1:567863/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTGGATTCGTAAg                                      <  1:1292169/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                      >  1:1100510/1‑76 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                      <  1:326647/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                      >  1:597958/1‑76 (MQ=255)
                                       ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGc                                     <  1:914419/76‑1 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                    >  1:652320/1‑76 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                    <  2:480284/76‑1 (MQ=255)
                                          aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                                  <  1:762816/76‑1 (MQ=255)
                                          aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                                  >  1:851815/1‑76 (MQ=255)
                                           gATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgaga                                 <  2:641266/76‑1 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                              <  1:967567/76‑1 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                              >  2:1234448/1‑76 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                              <  1:217469/76‑1 (MQ=255)
                                                ttCAACCCTTAAGTTAGCNCTTATGGGTTTTACCCCCGCTTGGTGGGTACTTCTTCGATTCGTAAGCGAGAACTGc                            >  1:1336500/1‑76 (MQ=25)
                                                  cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                          <  1:648178/76‑1 (MQ=255)
                                                     cccTNAAGTTAGCGCTTATGGGTCTTACACCGGCGTGGTNGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                       <  1:173922/76‑1 (MQ=25)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                      <  2:317398/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                     >  2:1119516/1‑76 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                     <  2:843361/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                     <  1:438505/76‑1 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                    >  1:800675/1‑76 (MQ=255)
                                                         tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc                   >  2:903858/1‑76 (MQ=255)
                                                            gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt                >  2:794992/1‑76 (MQ=255)
                                                            gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt                >  2:795008/1‑76 (MQ=255)
                                                            gTTAGCGCTTATGGGTTTTACCCCCGCGTAGTGGCTACTTCTTCGATTCGTAAGCGANAACTGCCTTAAACTCCAt                <  1:1333126/76‑1 (MQ=25)
                                                              tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc              >  2:373757/1‑76 (MQ=255)
                                                                 cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATTCCAc           <  2:1386856/76‑1 (MQ=255)
                                                                  gcgtATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACg          <  2:326487/73‑1 (MQ=255)
                                                                     tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTCAACTCCATACCACGTtc       >  1:156127/1‑76 (MQ=255)
                                                                      aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct      <  2:164664/76‑1 (MQ=255)
                                                                       tGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCGCGTtctc     >  2:609248/1‑76 (MQ=255)
                                                                       tGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtctc     >  2:609230/1‑76 (MQ=255)
                                                                        gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct    <  2:649780/76‑1 (MQ=255)
                                                                          gTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGa  <  1:325114/76‑1 (MQ=255)
                                                                          |                                                                           
CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGA  >  REL606/2774122‑2774271

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