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)100 (1.910) 65/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)95 (1.820) 68/144 0.0 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 202 total aligned reads displayed.

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                             >  2:467025‑M2/71‑76 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  1:394194‑M2/1‑73 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                          >  1:613136‑M2/1‑73 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                         <  2:790857‑M2/10‑1 (MQ=255)
     tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        >  2:1285279‑M2/66‑76 (MQ=255)
     tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        <  2:1701794‑M2/11‑1 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  2:1309317‑M2/1‑70 (MQ=255)
      cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                       <  2:1582679‑M2/12‑1 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  2:1083178‑M2/1‑69 (MQ=255)
        naacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     >  1:167273‑M2/63‑76 (MQ=255)
        aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataan                                                                     >  1:56617‑M2/1‑68 (MQ=255)
         acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    >  1:995879‑M2/1‑67 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  2:1298261‑M2/62‑76 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    <  2:1178747‑M2/15‑1 (MQ=255)
          caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                   >  1:810783‑M2/1‑66 (MQ=255)
          caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagag                                                                   <  2:415908‑M2/76‑11 (MQ=255)
          acgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgc                                                                   >  2:591192‑M2/61‑76 (MQ=255)
          acgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgc                                                                   >  1:1645638‑M2/61‑76 (MQ=255)
           cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCAGggcgcg                                                                  >  2:135033‑M2/60‑75 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCCAcccataagcgc                                                                  <  1:709462‑M2/76‑12 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                  <  2:757860‑M2/76‑12 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                  <  1:967119‑M2/76‑12 (MQ=255)
             cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  1:683087‑M2/1‑63 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             <  1:1477387‑M2/76‑17 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             >  1:251553‑M2/1‑60 (MQ=255)
                caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                             >  2:887678‑M2/55‑76 (MQ=255)
                 naatAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  1:384158‑M2/75‑18 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  2:223366‑M2/76‑19 (MQ=255)
                   cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                          <  1:640043‑M2/25‑1 (MQ=255)
                   aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          >  1:1347921‑M2/1‑57 (MQ=255)
                       gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:1597350‑M2/29‑1 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  1:1079975‑M2/76‑25 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTACTGCGGATGc                                                    >  2:698326‑M2/46‑76 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    >  1:1033585‑M2/46‑76 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    >  1:888759‑M2/46‑76 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  1:1392379‑M2/31‑1 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  1:176238‑M2/31‑1 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   >  2:1498067‑M2/45‑76 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   <  2:215856‑M2/32‑1 (MQ=255)
                          aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                   <  1:311508‑M2/76‑27 (MQ=255)
                           cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  1:9396‑M2/44‑76 (MQ=255)
                            gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                 >  2:1514651‑M2/43‑76 (MQ=255)
                            ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                                 >  1:426555‑M2/1‑48 (MQ=255)
                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               >  2:1639786‑M2/1‑46 (MQ=255)
                               ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                              <  1:1267101‑M2/37‑1 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                              >  2:103597‑M2/1‑45 (MQ=255)
                                tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                             <  2:1142871‑M2/38‑1 (MQ=255)
                                tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  2:892416‑M2/1‑44 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                            <  1:1492847‑M2/39‑1 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            >  2:815235‑M2/1‑43 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  1:465986‑M2/76‑34 (MQ=255)
                                  ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                           <  2:672670‑M2/40‑1 (MQ=255)
                                  ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  2:652725‑M2/1‑42 (MQ=255)
                                  ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:1141326‑M2/1‑42 (MQ=255)
                                   tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                          >  1:1445818‑M2/36‑76 (MQ=255)
                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCCTCAAgg                                         >  2:435489‑M2/35‑76 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  2:248893‑M2/76‑37 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         >  1:1209220‑M2/1‑40 (MQ=255)
                                     agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                        <  2:370625‑M2/43‑1 (MQ=255)
                                      tCCGTCCCAGGCGCGTTAGTGTGCAACCAGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                       >  2:903454‑M2/1‑38 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  2:745244‑M2/44‑1 (MQ=255)
                                        tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                     <  2:276024‑M2/46‑1 (MQ=255)
                                          ttCCAGGCGCGTTATTTTGCAGCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  1:739827‑M2/1‑34 (MQ=255)
                                          ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   >  2:1746366‑M2/1‑34 (MQ=255)
                                          ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                   <  1:1722349‑M2/76‑43 (MQ=255)
                                          gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   <  1:340222‑M2/48‑1 (MQ=255)
                                          gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                   <  1:351300‑M2/48‑1 (MQ=255)
                                            tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                 <  2:1052543‑M2/50‑1 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  1:1582198‑M2/55‑1 (MQ=255)
                                                  gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                           >  1:530389‑M2/1‑26 (MQ=255)
                                                   ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                          >  1:1383644‑M2/20‑76 (MQ=255)
                                                   ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                          >  1:1689814‑M2/20‑76 (MQ=255)
                                                    gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                         <  1:155618‑M2/76‑53 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  2:928236‑M2/18‑76 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        <  1:433911‑M2/59‑1 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       <  2:1087825‑M2/76‑55 (MQ=255)
                                                       aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                      <  1:702894‑M2/76‑56 (MQ=255)
                                                         tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg                    >  1:757860‑M2/1‑19 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                    >  1:850334‑M2/14‑76 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                    <  1:52677‑M2/63‑1 (MQ=255)
                                                            gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  1:1601371‑M2/1‑16 (MQ=255)
                                                             cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga                <  1:663231‑M2/76‑62 (MQ=255)
                                                              cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               <  2:324636‑M2/68‑1 (MQ=255)
                                                              cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               >  1:1154761‑M2/9‑76 (MQ=255)
                                                               ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc              <  2:762898‑M2/69‑1 (MQ=255)
                                                               aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat              <  1:1573351‑M2/76‑64 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  2:352373‑M2/7‑76 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:1245018‑M2/70‑1 (MQ=255)
                                                                 cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta            >  1:476745‑M2/1‑11 (MQ=255)
                                                                 atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            <  1:541580‑M2/71‑1 (MQ=255)
                                                                 atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            >  1:1277327‑M2/6‑76 (MQ=255)
                                                                   gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTCCAggg          >  1:204207‑M2/4‑76 (MQ=255)
                                                                    ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc         <  1:1526189‑M2/76‑69 (MQ=255)
                                                                    ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc         >  2:923040‑M2/1‑8 (MQ=255)
                                                                    ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         >  2:1020371‑M2/3‑76 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  1:1116725/76‑1 (MQ=255)
                                                                         cAAcccataagcgataacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  2:736874‑M2/76‑74 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga   >  1:228816‑M2/1‑2 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga   <  2:358639‑M2/76‑75 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  >  1:1231148/1‑76 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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