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)54 (0.900) 46/144 0.2 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)50 (0.840) 42/144 0.3 97.9% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.020)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTC  >  REL606/2774122‑2774266
                                                                          |                                                                      
cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGt                                                                       >  1:1352629/1‑76 (MQ=31)
 gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                      <  1:1149710/76‑1 (MQ=32)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                     <  2:101163/76‑1 (MQ=34)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                    <  2:704139/76‑1 (MQ=35)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                    <  2:1881408/76‑1 (MQ=35)
    gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                   <  2:449189/76‑1 (MQ=38)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                 <  2:1201706/76‑1 (MQ=40)
        cGTNTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                               >  1:1834267/1‑76 (MQ=255)
         gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                              >  2:950394/1‑76 (MQ=255)
             tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                          >  2:1262307/1‑76 (MQ=255)
                   ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGc                                                    >  2:738848/1‑76 (MQ=255)
                   ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGc                                                    >  2:261153/1‑76 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                               <  1:566825/76‑1 (MQ=255)
                         gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtc                                              >  2:1950384/1‑76 (MQ=255)
                         gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtc                                              <  2:1945412/76‑1 (MQ=255)
                                cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                       >  1:924553/1‑76 (MQ=255)
                                 gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTc                                      <  1:1518614/76‑1 (MQ=255)
                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                     <  1:984151/76‑1 (MQ=255)
                                   aTTCTTCAGATGGTTCANCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCTTGGTGGCTACCTCATCGATTCgc                                    >  1:252406/1‑75 (MQ=21)
                                   aTTCTTCAGATGGTTCANCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACATCCTCGCTTCGt                                    >  1:252384/1‑76 (MQ=21)
                                    ttcttcAGNTGGCTCAACCCTTAAGTTAGCGCTGCGGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                   <  2:346975/76‑1 (MQ=37)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                  <  1:1903907/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTACTTCTTCGATTCGtta                                  >  1:630304/1‑74 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                               <  2:1205599/76‑1 (MQ=255)
                                         cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                              >  1:265152/1‑76 (MQ=255)
                                          aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                             >  2:681780/1‑76 (MQ=255)
                                          aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                             >  2:1246178/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                           <  1:1074440/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                           >  2:434163/1‑76 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                          >  1:80442/1‑76 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGNACt                         <  1:1830882/76‑1 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAcn                         <  1:1254556/76‑2 (MQ=255)
                                              ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                         <  2:1177140/76‑1 (MQ=255)
                                               gTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTg                        <  2:458123/76‑1 (MQ=255)
                                                ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                       <  1:830676/76‑1 (MQ=255)
                                                ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                       <  1:926474/76‑1 (MQ=255)
                                                ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                       >  2:1314581/1‑76 (MQ=255)
                                                ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                       >  1:417196/1‑76 (MQ=255)
                                                 tCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGcc                      >  2:1574642/1‑76 (MQ=255)
                                                  cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                     >  2:220139/1‑76 (MQ=255)
                                                   aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                    <  2:1971385/76‑1 (MQ=255)
                                                    aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTa                   <  1:361426/76‑1 (MQ=255)
                                                     cccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                  <  1:1176785/76‑1 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                >  2:1218134/1‑76 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGGTCCTTTTTTTGTTCTTAAACGGGAACAGCCTTtaaa                >  1:979617/1‑75 (MQ=12)
                                                         tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc              >  1:343385/1‑76 (MQ=255)
                                                           aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa            >  2:1324502/1‑76 (MQ=255)
                                                             ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa          <  1:1265050/76‑1 (MQ=255)
                                                              tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc         <  1:717648/76‑1 (MQ=255)
                                                              tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc         >  1:448591/1‑76 (MQ=255)
                                                              tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGAAAGCGAGAACTGCCTTAAAATCCATAc         >  1:197851/1‑76 (MQ=255)
                                                                 cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc      <  1:969801/76‑1 (MQ=255)
                                                                   cTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGt    >  2:1259354/1‑76 (MQ=255)
                                                                    ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt   >  2:1736012/1‑76 (MQ=255)
                                                                     tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc  >  2:1939455/1‑76 (MQ=255)
                                                                     tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAATGCCTCAGAGTACACACCACGTtc  >  1:20555/1‑76 (MQ=25)
                                                                     tATGGGTTTTACCCCCGCGTGGTGGCTACATCTTCGATTCGTAAGCGAGAANTGCCTTAGACTCCATACCACTTtc  <  1:176761/76‑1 (MQ=37)
                                                                     tATGGGTTTTACCCCCGCGTAGTGCCTACTTCTTCGATTCGTAAGCGAGAACTACCATAAAGTCCATACCACGTtc  >  1:96169/1‑76 (MQ=37)
                                                                          |                                                                      
CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTC  >  REL606/2774122‑2774266

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.