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)77 (1.270) 53/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)58 (0.960) 47/144 0.2 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:887892/76‑1 (MQ=31)
 gAGGCTGCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                           <  2:69518/76‑1 (MQ=255)
 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                           <  2:1338196/76‑1 (MQ=32)
 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                           >  2:1300254/1‑76 (MQ=32)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTTAGTTAGCGCTTATGGGttt                                                                          <  1:275368/76‑1 (MQ=21)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                          <  2:1385945/76‑1 (MQ=34)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                          <  1:151038/76‑1 (MQ=34)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                         >  1:695260/1‑76 (MQ=35)
    gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                        <  1:1584350/76‑1 (MQ=38)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  1:1303530/1‑76 (MQ=40)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  1:1402023/1‑76 (MQ=40)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTTTAcc                                                                      >  1:1585512/1‑76 (MQ=17)
        cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    >  2:1283668/1‑76 (MQ=255)
        cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    >  1:877099/1‑76 (MQ=255)
        cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                    <  2:1068899/76‑1 (MQ=255)
          tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcg                                                                  >  2:1860838/1‑76 (MQ=255)
           cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgc                                                                 >  2:1079554/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                                >  2:677968/1‑76 (MQ=255)
            ttAACGAGGCAACGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                                >  1:66812/1‑76 (MQ=255)
             tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                               >  2:533305/1‑76 (MQ=255)
             tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                               >  2:146363/1‑76 (MQ=255)
              aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                              <  1:920881/76‑1 (MQ=255)
                cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggg                                                            >  2:156735/1‑75 (MQ=255)
                  aGGCACCGAGGCGTCACATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                          >  2:1107398/1‑76 (MQ=255)
                     cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTa                                                       <  2:1523885/76‑1 (MQ=255)
                      aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAc                                                      >  1:538615/1‑76 (MQ=255)
                       ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAct                                                     <  1:1181539/76‑1 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                    >  1:640390/1‑76 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                    >  2:1736223/1‑76 (MQ=255)
                          aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtct                                                  >  1:1851617/1‑76 (MQ=255)
                           ggCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtctt                                                 >  1:1313743/1‑76 (MQ=255)
                            gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACttcttc                                                <  1:827311/76‑1 (MQ=255)
                            gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACttcttc                                                <  1:1538373/76‑1 (MQ=255)
                              gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                              <  2:897729/76‑1 (MQ=255)
                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                         >  2:966926/1‑76 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                       >  2:1040048/1‑76 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                       <  2:1750857/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                       <  1:1905891/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCTATTCGTAAg                                      >  2:129444/1‑76 (MQ=255)
                                       ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGc                                     >  1:215999/1‑76 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                    <  2:1407655/76‑1 (MQ=255)
                                         cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                                   <  1:1283669/76‑1 (MQ=255)
                                          aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                                  <  1:808414/76‑1 (MQ=255)
                                           gATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgaga                                 <  1:1259908/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                                >  1:507835/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGCTTTACCCCCGCGCGGTGGGTACTTCTTCGCTTCGTAAGGGAGaa                                >  1:306854/1‑76 (MQ=21)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCTAGAACt                              <  2:935351/76‑1 (MQ=255)
                                               gTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTg                             >  2:961060/1‑76 (MQ=255)
                                                 tCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGNGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGcc                           >  1:1987743/1‑76 (MQ=255)
                                                  cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                          >  2:1121669/1‑76 (MQ=255)
                                                  cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                          >  2:397731/1‑76 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                      <  2:616812/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCCTAAAc                     <  1:129198/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTCTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                     >  2:1475497/1‑76 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                    <  2:975200/76‑1 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                    >  2:1030786/1‑76 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                    <  1:1608088/76‑1 (MQ=255)
                                                          aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                  >  1:1231038/1‑76 (MQ=255)
                                                          aaGTTAGCGCTTATGGGTTTTACCCCCGCGAGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                  >  1:1645908/1‑76 (MQ=255)
                                                             ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTNCTGCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa               <  1:1997707/76‑1 (MQ=255)
                                                                gcgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCa            <  2:1473448/76‑1 (MQ=255)
                                                                    ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACGACGtt        >  2:86725/1‑76 (MQ=255)
                                                                        gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct    <  1:695262/76‑1 (MQ=255)
                                                                          gTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGa  <  2:1851616/76‑1 (MQ=255)
                                                                          |                                                                           
CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGA  >  REL606/2774122‑2774271

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