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)79 (1.310) 58/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)62 (1.030) 46/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)

Only 100 of 146 total aligned reads displayed.

TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597
                                                                         |                                                                           
tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           >  2:1948865/1‑76 (MQ=255)
tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           <  2:195044/76‑1 (MQ=255)
gggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                           >  2:135359‑M2/74‑76 (MQ=255)
  gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         <  1:1228502‑M2/5‑1 (MQ=255)
  gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         >  1:380591‑M2/72‑76 (MQ=255)
   catcgtcttaacgaggcaccgaggcgtcgcattcttcggatgcttcaacccttaagttagcgcttatgggACGCaa                                                                        >  2:929820‑M2/71‑76 (MQ=255)
   aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGcga                                                                        >  1:2086/1‑74 (MQ=255)
    gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                       <  1:1673478‑M2/76‑2 (MQ=255)
    gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                       <  1:290406‑M2/76‑2 (MQ=255)
    atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       <  2:727060‑M2/7‑1 (MQ=255)
       gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    <  1:1859097‑M2/10‑1 (MQ=255)
       gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    >  2:382131‑M2/67‑76 (MQ=255)
       cAGAAACAACCCTGAATAACGGAACCAATATCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    <  2:1804951‑M2/76‑5 (MQ=255)
       cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  2:1724230‑M2/1‑72 (MQ=255)
        tcntaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaacccag                                                                   >  2:1404961‑M2/66‑74 (MQ=255)
        aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccac                                                                   <  1:928094‑M2/76‑6 (MQ=255)
         gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  1:870294‑M2/1‑70 (MQ=255)
         cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  >  2:895192‑M2/65‑76 (MQ=255)
          aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  2:1731668‑M2/1‑69 (MQ=255)
          aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  2:162155‑M2/1‑69 (MQ=255)
            acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                               >  1:152870‑M2/1‑67 (MQ=255)
            aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               <  2:530772‑M2/15‑1 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:690303‑M2/1‑66 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:1562223‑M2/1‑66 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:1616267‑M2/1‑66 (MQ=255)
             caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                              >  1:1036489‑M2/1‑66 (MQ=255)
              aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             >  1:193153‑M2/1‑65 (MQ=255)
               acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                            <  1:437944‑M2/76‑13 (MQ=255)
                cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           >  1:907780‑M2/1‑63 (MQ=255)
                cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           >  2:1218035‑M2/1‑63 (MQ=255)
                   caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        <  2:324445‑M2/22‑1 (MQ=255)
                   caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        <  1:1498513‑M2/22‑1 (MQ=255)
                   caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        >  2:900956‑M2/55‑76 (MQ=255)
                    taatAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       <  2:1522659‑M2/75‑18 (MQ=255)
                    accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  2:149850‑M2/23‑1 (MQ=255)
                    accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  2:1650116‑M2/23‑1 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  2:1454075‑M2/1‑58 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  2:1539225‑M2/1‑58 (MQ=255)
                      cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     >  1:626477‑M2/52‑76 (MQ=255)
                      aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     >  1:680660‑M2/1‑57 (MQ=255)
                       gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    <  2:1839015‑M2/26‑1 (MQ=255)
                         aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                  <  2:1228209‑M2/76‑23 (MQ=255)
                          gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:1441507‑M2/29‑1 (MQ=255)
                            gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                               >  1:1159095‑M2/46‑76 (MQ=255)
                             tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              <  1:448716‑M2/32‑1 (MQ=255)
                              cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             <  2:444207‑M2/33‑1 (MQ=255)
                              aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                             <  1:539296‑M2/76‑28 (MQ=255)
                                 attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          <  1:1621662‑M2/36‑1 (MQ=255)
                                   tgccatTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  1:1888875‑M2/3‑44 (MQ=255)
                                    cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  2:1614439‑M2/39‑1 (MQ=255)
                                    cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  2:946479‑M2/39‑1 (MQ=255)
                                    aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  1:1522660‑M2/1‑43 (MQ=255)
                                    aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  1:138333‑M2/1‑43 (MQ=255)
                                     ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                      <  2:528958‑M2/40‑1 (MQ=255)
                                      tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     <  2:1649705‑M2/41‑1 (MQ=255)
                                         tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                  <  2:789794‑M2/76‑39 (MQ=255)
                                           cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                >  2:40644‑M2/1‑36 (MQ=255)
                                            ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                               <  2:1835580‑M2/47‑1 (MQ=255)
                                            ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                               <  1:1186690‑M2/47‑1 (MQ=255)
                                             ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  1:139303‑M2/76‑43 (MQ=255)
                                             gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              <  1:888154‑M2/48‑1 (MQ=255)
                                             gttcaaccctgaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              <  1:888152‑M2/48‑1 (MQ=255)
                                              ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             >  1:992521‑M2/28‑76 (MQ=255)
                                              tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                             <  1:1441801‑M2/76‑44 (MQ=255)
                                               tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            >  1:152573‑M2/27‑76 (MQ=255)
                                               ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            >  2:63936‑M2/1‑32 (MQ=255)
                                                 aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          >  2:543731‑M2/25‑76 (MQ=255)
                                                   gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  1:1075151‑M2/1‑28 (MQ=255)
                                                   gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  2:1809570‑M2/1‑28 (MQ=255)
                                                   gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  1:1305684‑M2/1‑28 (MQ=255)
                                                    cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:716665‑M2/1‑27 (MQ=255)
                                                     cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      <  1:1731240‑M2/56‑1 (MQ=255)
                                                      ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     >  2:1145448‑M2/20‑76 (MQ=255)
                                                      cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     <  1:1115384‑M2/76‑52 (MQ=255)
                                                        aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                   <  2:1108976‑M2/59‑1 (MQ=255)
                                                         tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  1:835369‑M2/76‑55 (MQ=255)
                                                           ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                >  1:246450‑M2/1‑20 (MQ=255)
                                                           ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                <  2:1757456‑M2/62‑1 (MQ=255)
                                                            tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTTccc               <  2:794217‑M2/63‑1 (MQ=255)
                                                             ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga              <  1:1846721‑M2/76‑59 (MQ=255)
                                                               gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  2:1917688‑M2/1‑16 (MQ=255)
                                                                 cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  1:1409577‑M2/9‑76 (MQ=255)
                                                                 aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          >  2:6342‑M2/1‑14 (MQ=255)
                                                                  ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:1697329‑M2/8‑76 (MQ=255)
                                                                  aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat         >  1:10883‑M2/1‑13 (MQ=255)
                                                                   ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:105110‑M2/1‑12 (MQ=255)
                                                                   ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:1453600‑M2/1‑12 (MQ=255)
                                                                    cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:386126‑M2/1‑11 (MQ=255)
                                                                    atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGTCCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg       <  1:1606499‑M2/71‑1 (MQ=255)
                                                                     tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1841072‑M2/5‑76 (MQ=255)
                                                                     tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1778220‑M2/5‑76 (MQ=255)
                                                                     tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1000315‑M2/5‑76 (MQ=255)
                                                                     tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      >  2:1035153‑M2/5‑76 (MQ=255)
                                                                      gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg     <  1:1323575‑M2/73‑1 (MQ=255)
                                                                       ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    <  2:434027‑M2/76‑69 (MQ=255)
                                                                       ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    <  2:1839835‑M2/76‑69 (MQ=255)
                                                                        gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct   <  1:1834496‑M2/76‑70 (MQ=255)
                                                                         ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  <  1:1285886‑M2/76‑71 (MQ=255)
                                                                         ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  <  1:1866981/76‑1 (MQ=255)
                                                                         ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  >  1:446509/1‑76 (MQ=255)
                                                                         |                                                                           
TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597

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