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)55 (1.150) 42/142 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)53 (1.110) 46/142 0.1 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCT  >  REL606/2774126‑2774269
                                                                      |                                                                         
gCATCGTCTTAGCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                       <  2:803093/75‑1 (MQ=25)
   tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                    >  1:1149176/1‑75 (MQ=40)
   tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                    >  2:1030057/1‑75 (MQ=40)
    cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccc                                                                   >  1:429548/1‑75 (MQ=255)
       cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcg                                                                >  2:319242/1‑75 (MQ=255)
           aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                            <  2:638997/75‑1 (MQ=255)
             gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                          <  1:70547/75‑1 (MQ=255)
              aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                         >  2:1414650/1‑75 (MQ=255)
               ggCACCGAGGCGTCGCATTCTTCAGATGGTTAAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                        >  2:32452/1‑75 (MQ=255)
                 cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                      <  1:698989/75‑1 (MQ=255)
                   ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAc                                                    >  1:654331/1‑75 (MQ=255)
                   ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTAc                                                    >  1:901168/1‑75 (MQ=255)
                             gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTCATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                          >  1:520075/1‑75 (MQ=255)
                              cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTc                                         >  1:795235/1‑75 (MQ=255)
                              cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTTTACCCCCGCGTGGTGGCTACTTCTTCGATTc                                         >  2:373279/1‑75 (MQ=255)
                                ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                       >  2:277234/1‑75 (MQ=255)
                                 tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                      >  1:473421/1‑75 (MQ=255)
                                      aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                                 >  2:1087923/1‑75 (MQ=255)
                                      aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                                 >  2:244418/1‑75 (MQ=255)
                                      aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                                 <  1:439893/75‑1 (MQ=255)
                                      aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                                 <  2:673366/75‑1 (MQ=255)
                                         tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                              >  2:132705/1‑75 (MQ=255)
                                          ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                             <  2:90000/75‑1 (MQ=255)
                                          ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                             >  2:1171830/1‑75 (MQ=255)
                                            ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGCTTCGTAAGCGAGAACTg                           >  1:661778/1‑75 (MQ=255)
                                               aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                        <  2:254051/75‑1 (MQ=255)
                                                aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                       <  2:266758/75‑1 (MQ=255)
                                                  ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                     <  2:377811/75‑1 (MQ=255)
                                                  ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                     <  1:1007888/75‑1 (MQ=255)
                                                  ccTTAAGTGAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                     <  2:1264689/75‑1 (MQ=255)
                                                    ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                   <  1:1438844/75‑1 (MQ=255)
                                                     tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                  >  1:1208238/1‑75 (MQ=255)
                                                         ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt              <  2:1171362/75‑1 (MQ=255)
                                                          tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa             >  1:681963/1‑75 (MQ=255)
                                                           aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc            >  2:376191/1‑75 (MQ=255)
                                                           aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc            <  1:1355857/75‑1 (MQ=255)
                                                            gcgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAcc           >  1:1607891/1‑75 (MQ=255)
                                                             cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCa          <  1:1392221/75‑1 (MQ=255)
                                                             cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCa          <  2:919658/75‑1 (MQ=255)
                                                              gcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc         >  2:196387/1‑75 (MQ=255)
                                                              gcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc         >  1:552808/1‑75 (MQ=255)
                                                                ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGt       <  2:141968/75‑1 (MQ=255)
                                                                tactGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGt       >  2:153926/4‑75 (MQ=255)
                                                                 tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt      >  1:311006/1‑75 (MQ=255)
                                                                 tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt      >  2:823642/1‑75 (MQ=255)
                                                                  aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc     <  1:925712/75‑1 (MQ=255)
                                                                     ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct  >  1:100737/1‑75 (MQ=255)
                                                                     ggTTTTACCCCCGCGTGGTGGCTACTTCTTCCATTCGTAAGCGAGAACTGCCTTACACTCCATACCACGTTctct  >  2:823636/1‑75 (MQ=255)
                                                                      |                                                                         
GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCT  >  REL606/2774126‑2774269

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