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)57 (1.310) 42/144 0.1 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.760) 30/144 0.5 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524452‑4524601
                                                                           |                                                                          
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                            <  1:391533‑M2/6‑1 (MQ=255)
   cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                         >  1:678495‑M2/68‑76 (MQ=255)
    cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                        >  1:445919‑M2/1‑72 (MQ=255)
     aGAAACAACGCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                       <  2:1229953‑M2/76‑6 (MQ=255)
     aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                       <  1:371611‑M2/76‑6 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                      >  1:899548‑M2/1‑70 (MQ=255)
      cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                      >  2:1251409‑M2/65‑76 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                     >  1:766527‑M2/1‑69 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                     >  1:53032‑M2/1‑69 (MQ=255)
       aaaCAACCCTGAAAAACGGAACCAATACCATTCNGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                     <  2:1214404‑M2/76‑8 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                    >  1:495761‑M2/63‑76 (MQ=255)
        aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                    >  2:366398‑M2/1‑68 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                 <  1:65344‑M2/76‑12 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                 <  1:142317‑M2/76‑12 (MQ=255)
             cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                               >  2:112356‑M2/1‑63 (MQ=255)
              ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                              >  2:973815‑M2/57‑76 (MQ=255)
                  ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                          <  2:725393‑M2/24‑1 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                          <  1:1107581‑M2/76‑19 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                          >  2:297011‑M2/1‑58 (MQ=255)
                   aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                         >  1:218779‑M2/1‑57 (MQ=255)
                       gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                     <  1:266944‑M2/29‑1 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                    <  1:404452‑M2/76‑25 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                    <  1:854140‑M2/76‑25 (MQ=255)
                        cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                    <  2:343014‑M2/30‑1 (MQ=255)
                        cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                    >  2:805885‑M2/47‑76 (MQ=255)
                          aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                               <  1:877823‑M2/63‑14 (MQ=255)
                          aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                               >  2:877823‑M2/1‑50 (MQ=255)
                           aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                 <  1:752588‑M2/76‑28 (MQ=255)
                           aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                 <  1:1105715‑M2/76‑28 (MQ=255)
                             cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                               >  2:1163682‑M2/1‑47 (MQ=255)
                              attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                              >  2:1385457‑M2/41‑76 (MQ=255)
                               ttcttcagatggttcaactcttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCCGCCtt                                             >  2:1381255‑M2/40‑76 (MQ=255)
                               ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                             <  2:1098659‑M2/37‑1 (MQ=255)
                               ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                             <  2:1010498‑M2/37‑1 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                             >  1:339391‑M2/1‑45 (MQ=255)
                               aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtan                                             >  1:1226148‑M2/1‑45 (MQ=255)
                                tACCATTCCGTTCCAGGCGTGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                            >  1:563636‑M2/1‑44 (MQ=255)
                                tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctnacttaagggttgtggtatt                                            >  1:1347374‑M2/1‑44 (MQ=255)
                                tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                            >  2:670099‑M2/1‑44 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                           <  2:1063749‑M2/39‑1 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                           <  1:429379‑M2/39‑1 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                           >  1:734576‑M2/38‑76 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                           >  1:1404907‑M2/1‑43 (MQ=255)
                                   cATTCCGTTCCAGNCGCGTTATTTTGCANCCNGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                         >  2:1302711‑M2/1‑41 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                         <  2:172331‑M2/76‑36 (MQ=255)
                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                        >  1:1119695‑M2/35‑76 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                        <  2:147896‑M2/76‑37 (MQ=255)
                                        tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                    <  1:1159025‑M2/46‑1 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   <  2:671550‑M2/76‑42 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   >  2:961975‑M2/1‑35 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   <  1:557540‑M2/76‑42 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   <  1:1318706‑M2/76‑42 (MQ=255)
                                         gTTCCAGGCGCGTAATTTTGCAAGCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   >  1:1236416‑M2/1‑35 (MQ=255)
                                          ttCCAGGCGTGTTGTTTTGCAACCGGGGCCGCAAccgataagcgctaacttaagggttgtggtattacgcctgata                                  >  1:372364‑M2/1‑34 (MQ=255)
                                          ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                  <  1:562691‑M2/76‑43 (MQ=255)
                                          gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                  >  1:974219‑M2/29‑76 (MQ=255)
                                           ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                 >  1:532694‑M2/28‑76 (MQ=255)
                                           tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                 >  1:1219536‑M2/1‑33 (MQ=255)
                                            tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                >  2:298842‑M2/27‑76 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                >  1:267962‑M2/1‑32 (MQ=255)
                                             caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                               <  1:61890‑M2/51‑1 (MQ=255)
                                             cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                               <  1:926123‑M2/76‑46 (MQ=255)
                                              aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                              <  2:1158362‑M2/76‑47 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                            >  1:681706‑M2/1‑28 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                            >  2:1258510‑M2/1‑28 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                            >  1:942693‑M2/1‑28 (MQ=255)
                                                gcgcGTTATTTTGCAACCGGGGCCGCAAccaataagngctaacgaaatggttgtggtattacgcctgatatnattt                            <  1:173998‑M2/76‑49 (MQ=255)
                                                cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                            <  1:316528‑M2/54‑1 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                          >  2:467310‑M2/21‑76 (MQ=255)
                                                  cttaagttagcgctaatgggCGGCAACCCCGGGTGCTGCTAATGCGGATGCGAACTTCTAGGTTCAGTATCAATaa                          <  2:224119‑M2/56‑1 (MQ=255)
                                                   ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                         >  2:620834‑M2/20‑76 (MQ=255)
                                                    gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                        <  2:242879‑M2/76‑53 (MQ=255)
                                                     ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                       >  2:1131380‑M2/1‑23 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                       <  1:345017‑M2/59‑1 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                      >  2:1023910‑M2/1‑22 (MQ=255)
                                                      agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                      <  2:1119487‑M2/60‑1 (MQ=255)
                                                         tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg                   >  2:274298‑M2/1‑19 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                 <  2:13840‑M2/65‑1 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGGTCAGTATCAATAACCTACCCAg                 <  1:5443‑M2/65‑1 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                <  1:634553‑M2/66‑1 (MQ=255)
                                                              aaCCGGGACCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa              <  1:799572‑M2/76‑63 (MQ=255)
                                                               ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATACACTACCCAGGTTc             >  1:1131546‑M2/8‑76 (MQ=255)
                                                                 cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta           <  2:757830‑M2/76‑66 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac          <  1:1392741‑M2/76‑67 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac          <  2:1287323‑M2/76‑67 (MQ=255)
                                                                    ggCCGCAANCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa        <  1:1271958‑M2/74‑1 (MQ=255)
                                                                      ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc      <  2:1289175‑M2/76‑71 (MQ=255)
                                                                       cGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctca     >  1:550568‑M2/1‑5 (MQ=255)
                                                                       cGCAACCCCGGGTGCTGCTAATGCGGATGCGTCCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt     <  2:171312/76‑1 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc    >  1:686432/1‑76 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc    <  2:164094/76‑1 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg   <  2:853647‑M2/76‑74 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga  <  1:1280474‑M2/76‑75 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt  >  2:482761/1‑76 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt  >  2:348345/1‑76 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt  >  1:1215390/1‑76 (MQ=255)
                                                                           |                                                                          
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524452‑4524601

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.