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)46 (0.900) 38/144 0.3 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)33 (0.640) 29/144 0.9 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602
                                                                          |                                                                           
gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                            >  2:313864‑M2/1‑75 (MQ=255)
 tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                           >  1:901768‑M2/69‑76 (MQ=255)
   gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         <  1:1529185‑M2/10‑1 (MQ=255)
   cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         <  1:96119‑M2/76‑5 (MQ=255)
   cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         >  2:196077‑M2/1‑72 (MQ=255)
    tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        <  1:1230343‑M2/11‑1 (MQ=255)
     gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       <  1:511535‑M2/76‑7 (MQ=255)
     gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       <  2:891440‑M2/76‑7 (MQ=255)
       taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     <  1:304235‑M2/14‑1 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     >  1:1549856‑M2/1‑68 (MQ=255)
          cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  1:1414841‑M2/60‑76 (MQ=255)
          aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                  <  1:1009358‑M2/76‑12 (MQ=255)
           acCCGGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 >  2:367977‑M2/1‑64 (MQ=255)
            cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  1:1295884‑M2/1‑63 (MQ=255)
             ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                               <  1:1653344‑M2/20‑1 (MQ=255)
              gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                              <  2:616698‑M2/21‑1 (MQ=255)
              cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                              <  2:1042611‑M2/76‑16 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:1146159‑M2/1‑58 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:687455‑M2/1‑58 (MQ=255)
                  cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                          >  2:1283256‑M2/52‑76 (MQ=255)
                  aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          <  2:650177‑M2/76‑20 (MQ=255)
                   gaggcgtcgcattcttcagatggttcaacccttaagtagacgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         <  1:267867‑M2/26‑1 (MQ=255)
                    aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcacataagcgctaacttaagg                                                        <  2:1464058‑M2/76‑22 (MQ=255)
                     aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                       <  1:1482422‑M2/76‑23 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  2:636825‑M2/29‑1 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  1:1532875‑M2/76‑25 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  2:427213‑M2/1‑52 (MQ=255)
                        gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  1:101852‑M2/31‑1 (MQ=255)
                          cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGcgg                                                  >  2:1468387‑M2/44‑75 (MQ=255)
                          aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                  >  2:1658856‑M2/1‑49 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  1:333108‑M2/1‑44 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:728382‑M2/1‑43 (MQ=255)
                                 ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  2:1192196‑M2/1‑42 (MQ=255)
                                    ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                        >  1:369807‑M2/1‑39 (MQ=255)
                                     gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  1:286170‑M2/44‑1 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      <  1:223648‑M2/76‑40 (MQ=255)
                                       cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                     <  1:1196695‑M2/76‑41 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  2:1044528‑M2/1‑34 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  2:176209‑M2/1‑34 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  2:101807‑M2/1‑34 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   <  1:1297053‑M2/76‑43 (MQ=255)
                                         gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   >  2:698408‑M2/29‑76 (MQ=255)
                                          ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                  <  2:927256‑M2/49‑1 (MQ=255)
                                           tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 <  2:874372‑M2/50‑1 (MQ=255)
                                           ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 <  1:1480980‑M2/76‑45 (MQ=255)
                                              ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatt                              <  2:1532364‑M2/76‑48 (MQ=255)
                                                ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:311742‑M2/55‑1 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  2:140911‑M2/56‑1 (MQ=255)
                                                  ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                          <  1:1123877‑M2/57‑1 (MQ=255)
                                                  cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                          >  1:825833‑M2/1‑25 (MQ=255)
                                                    aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAaccn                        <  1:508661‑M2/59‑2 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       <  1:51092‑M2/76‑55 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:498558‑M2/1‑22 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:356819‑M2/1‑22 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:550243‑M2/60‑1 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       <  1:350931‑M2/60‑1 (MQ=255)
                                                      aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                      >  2:734396‑M2/1‑21 (MQ=255)
                                                       ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                     <  1:6780‑M2/62‑1 (MQ=255)
                                                       ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                     >  2:978498‑M2/15‑76 (MQ=255)
                                                          gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  2:446780‑M2/65‑1 (MQ=255)
                                                          gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  1:1588658‑M2/65‑1 (MQ=255)
                                                           cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCTAgg                 <  1:388350‑M2/66‑1 (MQ=255)
                                                           cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 <  2:1530474‑M2/66‑1 (MQ=255)
                                                           cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 >  1:577258‑M2/11‑76 (MQ=255)
                                                            cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga                <  2:402699‑M2/76‑62 (MQ=255)
                                                            cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga                >  2:1598758‑M2/1‑15 (MQ=255)
                                                            cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga                <  1:1461478‑M2/76‑62 (MQ=255)
                                                             cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               >  2:1667645‑M2/9‑76 (MQ=255)
                                                             cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               <  2:1638381‑M2/68‑1 (MQ=255)
                                                               tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  2:1426832‑M2/7‑76 (MQ=255)
                                                                cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            >  1:981173‑M2/1‑11 (MQ=255)
                                                                atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            <  1:205203‑M2/71‑1 (MQ=255)
                                                                   ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc         <  2:1105619‑M2/76‑69 (MQ=255)
                                                                    gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct        >  2:1439437‑M2/1‑7 (MQ=255)
                                                                    gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc        >  1:1178725‑M2/2‑76 (MQ=255)
                                                                     ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:948195‑M2/1‑6 (MQ=255)
                                                                     ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       >  2:884821‑M2/1‑6 (MQ=255)
                                                                      cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca      >  1:1185774‑M2/1‑5 (MQ=255)
                                                                      cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca      >  2:1299460‑M2/1‑5 (MQ=255)
                                                                      cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt      >  1:1580231/1‑76 (MQ=255)
                                                                       gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac     <  1:117983‑M2/76‑73 (MQ=255)
                                                                       gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac     >  1:403941‑M2/1‑4 (MQ=255)
                                                                       gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  1:1217280/76‑1 (MQ=255)
                                                                        cAACCCCGGGTGCTGCTAATGCGGATGTGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa    >  2:1556857/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  1:840482/1‑76 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  2:986368/1‑76 (MQ=255)
                                                                          aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  1:1314023/76‑1 (MQ=255)
                                                                          |                                                                           
GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602

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