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)302 (1.350) 117/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)238 (1.060) 100/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 572 total aligned reads displayed.

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602
                                                                           |                                                                           
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                             <  2:1380576‑M2/6‑1 (MQ=255)
aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                             >  1:4259514/1‑76 (MQ=255)
 gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                            >  2:6604215‑M2/1‑75 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                             >  2:4230762‑M2/70‑75 (MQ=255)
    cAGAAACAGCCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         >  2:6738046‑M2/1‑72 (MQ=255)
    cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                         >  2:4117051‑M2/1‑72 (MQ=255)
     tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                        <  1:7147827‑M2/11‑1 (MQ=255)
     aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  1:1611376‑M2/1‑71 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  2:6257813‑M2/1‑70 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  1:3586483‑M2/1‑70 (MQ=255)
      aaaacAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                       >  1:2044501‑M2/2‑70 (MQ=255)
       ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                      <  2:6791821‑M2/13‑1 (MQ=255)
       aaaCACCCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  1:3233326‑M2/1‑69 (MQ=255)
       aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  1:5064069‑M2/1‑69 (MQ=255)
        taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                     <  1:5286225‑M2/14‑1 (MQ=255)
         acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                    <  1:12831‑M2/76‑10 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    >  1:3822974‑M2/62‑76 (MQ=255)
          caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                   <  2:1468271‑M2/76‑11 (MQ=255)
          caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                   >  1:7075109‑M2/1‑66 (MQ=255)
            acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                 <  1:2378140‑M2/76‑13 (MQ=255)
             cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  1:5716967‑M2/1‑63 (MQ=255)
               gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                              <  1:1183875‑M2/21‑1 (MQ=255)
                tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             <  1:962799‑M2/76‑17 (MQ=255)
                cacggaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                             <  2:6332097‑M2/22‑1 (MQ=255)
                  ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                           <  2:2389020‑M2/24‑1 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  1:5501014‑M2/76‑19 (MQ=255)
                    tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                         <  2:514435‑M2/76‑21 (MQ=255)
                       gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:862687‑M2/29‑1 (MQ=255)
                       gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:7492060‑M2/29‑1 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  1:188397‑M2/76‑25 (MQ=255)
                        ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  2:6977646‑M2/76‑25 (MQ=255)
                        ggAACCAATACCATTCCGTGCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  1:2795758‑M2/1‑52 (MQ=255)
                        cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                     <  2:2985985‑M2/30‑1 (MQ=255)
                         gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    >  2:459893‑M2/46‑76 (MQ=255)
                          tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                                   <  1:2990395‑M2/32‑1 (MQ=255)
                          aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                   >  2:766751‑M2/1‑50 (MQ=255)
                          aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                   >  2:7077583‑M2/1‑50 (MQ=255)
                           cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                  >  2:2123601‑M2/44‑76 (MQ=255)
                           aCCAATACCTTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                  >  1:3969601‑M2/1‑49 (MQ=255)
                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                               <  1:3146915‑M2/76‑31 (MQ=255)
                              aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcacataagcgctaacttaagggttgtggaa                                               <  1:5788778‑M2/76‑31 (MQ=255)
                               ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                              <  2:7204509‑M2/37‑1 (MQ=255)
                                 cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                            <  2:574564‑M2/39‑1 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  1:6600740‑M2/76‑34 (MQ=255)
                                  ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:371542‑M2/1‑42 (MQ=255)
                                   tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                          <  1:1802208‑M2/41‑1 (MQ=255)
                                   cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          <  1:1065710‑M2/76‑36 (MQ=255)
                                    cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                         <  2:6020686‑M2/42‑1 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  1:4761103‑M2/76‑37 (MQ=255)
                                     agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                        >  1:6688314‑M2/34‑76 (MQ=255)
                                      gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  2:1149323‑M2/44‑1 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:4006289‑M2/1‑37 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  1:1382842‑M2/1‑37 (MQ=255)
                                       atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                      >  2:4986933‑M2/32‑76 (MQ=255)
                                          gtcCAGGCGCGTTATTTTGCAACCGGGGGCGCAAcccctaagcgctaacttaaaggttttggtattacgcctgata                                   >  2:5566126‑M2/2‑34 (MQ=255)
                                           ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                  >  2:31458‑M2/28‑76 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  1:5659349‑M2/1‑32 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  1:6340878‑M2/1‑32 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  2:2798973‑M2/1‑32 (MQ=255)
                                            ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 >  1:2915170‑M2/1‑32 (MQ=255)
                                              aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCCAGGGTCAGTATCa                               >  2:186372‑M2/25‑76 (MQ=255)
                                               acccttaagttagcgcttatgggCCGCATCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                              <  2:2020693‑M2/53‑1 (MQ=255)
                                                cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                             >  1:7273711‑M2/23‑76 (MQ=255)
                                                 ggcgtTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  2:5704363‑M2/2‑27 (MQ=255)
                                                 ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:1805768‑M2/55‑1 (MQ=255)
                                                  cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                           <  2:5326298‑M2/56‑1 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        <  2:4331914‑M2/59‑1 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:2639656‑M2/1‑22 (MQ=255)
                                                        ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                     >  1:6092600‑M2/15‑76 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                    <  2:365539‑M2/63‑1 (MQ=255)
                                                          agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                   <  1:4095900‑M2/64‑1 (MQ=255)
                                                           tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat                  <  1:2469579‑M2/76‑60 (MQ=255)
                                                           gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  1:2597822‑M2/65‑1 (MQ=255)
                                                            gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                 >  1:6207878‑M2/1‑16 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 >  1:4829522‑M2/11‑76 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 <  1:6759521‑M2/66‑1 (MQ=255)
                                                            cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                 <  2:5414631‑M2/66‑1 (MQ=255)
                                                               aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat              >  2:1317931‑M2/1‑13 (MQ=255)
                                                               aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat              >  1:2238172‑M2/1‑13 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  2:2515028‑M2/70‑1 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             >  1:418039‑M2/7‑76 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa             <  1:7671176‑M2/70‑1 (MQ=255)
                                                                ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             >  1:7408650‑M2/1‑12 (MQ=255)
                                                                ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt             >  2:3462854‑M2/1‑12 (MQ=255)
                                                                 atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg            >  2:4711803‑M2/6‑76 (MQ=255)
                                                                  tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg           >  2:4494263‑M2/5‑76 (MQ=255)
                                                                  tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg           <  1:4932985‑M2/72‑1 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           >  1:3489984‑M2/1‑10 (MQ=255)
                                                                  ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           <  1:6454976‑M2/76‑67 (MQ=255)
                                                                   gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact          <  1:4475589‑M2/76‑68 (MQ=255)
                                                                      ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc       <  2:1573296‑M2/76‑71 (MQ=255)
                                                                      ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg       <  1:3354277/76‑1 (MQ=255)
                                                                       cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt      >  2:5931552/1‑76 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  1:7004395/76‑1 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  2:4878692‑M2/76‑74 (MQ=255)
                                                                         cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg    <  2:7369282‑M2/76‑74 (MQ=255)
                                                                          aaCCCCGGGTGCTGCTAATGCGGATGCTACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt   >  1:5199299/1‑76 (MQ=255)
                                                                           acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  2:7471846‑M2/1‑1 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGTCCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  2:3856706/76‑1 (MQ=255)
                                                                           aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAtt  <  1:154301/76‑1 (MQ=255)
                                                                           |                                                                           
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524452‑4524602

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