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)72 (1.280) 53/144 0.0 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)52 (0.920) 40/144 0.3 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774122‑2774270
                                                                          |                                                                          
cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGt                                                                           >  1:1167549/1‑76 (MQ=31)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                        <  1:1858703/76‑1 (MQ=35)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                     >  2:1190555/1‑76 (MQ=40)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                     >  2:1476130/1‑76 (MQ=40)
           cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgc                                                                >  2:1278760/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                               >  2:80570/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                               >  1:352079/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                               <  2:1840101/76‑1 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                               >  2:186166/1‑76 (MQ=255)
             tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                              >  1:1850015/1‑76 (MQ=255)
               aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtgg                                                            <  1:1477998/76‑1 (MQ=255)
                 gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          >  1:864169/1‑76 (MQ=255)
                 gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          >  1:1062929/1‑76 (MQ=255)
                  aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                         <  1:119813/76‑1 (MQ=255)
                    gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                       <  1:130988/76‑1 (MQ=255)
                     cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTa                                                      >  2:572673/1‑76 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                   <  2:1453682/76‑1 (MQ=255)
                           ggCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtctt                                                <  1:1705778/76‑1 (MQ=255)
                              gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                             >  1:1347735/1‑76 (MQ=255)
                                    ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                       <  2:182193/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                      >  1:1750137/1‑76 (MQ=255)
                                       ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGc                                    >  1:713547/1‑76 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                   <  1:1798874/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               >  2:1501266/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               <  1:610401/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               >  1:204576/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               <  1:990520/76‑1 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAGGCGAGAAc                              <  1:1044979/76‑1 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                              <  1:1588362/76‑1 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                              <  2:1680077/76‑1 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGGGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                              <  1:95411/76‑1 (MQ=255)
                                                     cccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                      <  1:528420/76‑1 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                     <  1:1235265/76‑1 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                     <  1:546326/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                    >  1:1630067/1‑76 (MQ=255)
                                                         tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc                  >  2:469070/1‑76 (MQ=255)
                                                         tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc                  >  2:1263371/1‑76 (MQ=255)
                                                          aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                 <  1:1084856/76‑1 (MQ=255)
                                                          aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                 <  2:1821645/76‑1 (MQ=255)
                                                            gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt               >  2:1404034/1‑76 (MQ=255)
                                                             ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa              >  1:1500002/1‑76 (MQ=255)
                                                             ttAACGCTTTTGTGTTTTTCCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa              <  2:342714/76‑1 (MQ=25)
                                                              tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc             <  1:762822/76‑1 (MQ=255)
                                                               aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAcc            <  1:747294/76‑1 (MQ=255)
                                                               aGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAATCCATAcc            >  1:1126063/1‑76 (MQ=255)
                                                                 cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc          <  2:1625743/76‑1 (MQ=255)
                                                                 cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc          >  1:863357/1‑76 (MQ=255)
                                                                  gcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACg         >  2:1784104/1‑76 (MQ=255)
                                                                  gcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACg         <  1:969411/76‑1 (MQ=255)
                                                                      aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct     >  1:1031462/1‑76 (MQ=255)
                                                                      aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct     <  2:1278062/76‑1 (MQ=255)
                                                                       tGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtctc    >  2:1626129/1‑76 (MQ=255)
                                                                        gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct   >  1:602727/1‑76 (MQ=255)
                                                                         ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCCTCCCACGTTCTCTg  >  1:1406219/1‑76 (MQ=255)
                                                                         ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTg  <  1:1585244/76‑1 (MQ=255)
                                                                          |                                                                          
CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774122‑2774270

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