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)230 (1.200) 105/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)207 (1.080) 97/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 448 total aligned reads displayed.

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596
                                                                           |                                                                          
cgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggc                                                                            <  1:4972893‑M2/1‑1 (MQ=255)
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                            >  2:5830774/1‑76 (MQ=255)
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                            <  2:5390039/76‑1 (MQ=255)
 gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                           >  2:5055019‑M2/75‑76 (MQ=255)
 aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                           >  2:3407368/1‑76 (MQ=255)
 aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                           >  1:3227733/1‑76 (MQ=255)
  aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                          <  1:1274250‑M2/3‑1 (MQ=255)
    tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                        <  1:2028927/76‑1 (MQ=255)
    tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                        >  1:3720568/1‑76 (MQ=255)
    gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                        <  1:147559‑M2/5‑1 (MQ=255)
     aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                       >  1:2105497/1‑76 (MQ=255)
      gTTCAGAAACAACCCTGAATAACGGAACCAATCCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                      >  2:532943‑M2/1‑75 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                      >  1:4881347‑M2/70‑76 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     >  2:5741209‑M2/1‑74 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     <  1:6241335‑M2/76‑3 (MQ=255)
         gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                   >  1:1983232‑M2/67‑76 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                  >  2:705172‑M2/1‑71 (MQ=255)
           gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                 >  2:602619‑M2/1‑70 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                 <  1:2624242‑M2/12‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                               <  1:4736188‑M2/14‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                               <  2:2288076‑M2/14‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                               <  2:3804525‑M2/14‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                               <  2:5161526‑M2/14‑1 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                               >  2:2567654‑M2/1‑68 (MQ=255)
                cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                            >  2:4169610‑M2/60‑76 (MQ=255)
                  cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                          >  2:6107920‑M2/1‑63 (MQ=255)
                   ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                         <  2:2993568‑M2/20‑1 (MQ=255)
                   ccTGAAAAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                         >  1:4270662‑M2/1‑62 (MQ=255)
                     caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                       >  1:836663‑M2/55‑76 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                      <  1:5949655‑M2/76‑18 (MQ=255)
                       ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                     >  2:2773897‑M2/53‑76 (MQ=255)
                        cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                    <  2:4731500‑M2/25‑1 (MQ=255)
                           aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                 >  1:4616446‑M2/1‑54 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                <  2:1995678‑M2/29‑1 (MQ=255)
                             ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                               <  1:3427467‑M2/76‑25 (MQ=255)
                             cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                               >  1:2508828‑M2/47‑76 (MQ=255)
                             cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                               <  2:131496‑M2/30‑1 (MQ=255)
                               tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                             <  1:238798‑M2/32‑1 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             <  1:1728118‑M2/76‑27 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             <  1:285746‑M2/76‑27 (MQ=255)
                                aCCAATCCCATTCCGTTCCAGGCGCGTTATTTTACAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                            >  2:1074610‑M2/1‑49 (MQ=255)
                                 ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                           >  1:1055795‑M2/1‑48 (MQ=255)
                                  cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                          >  1:3008080‑M2/1‑47 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                         >  2:4280659‑M2/41‑76 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         <  1:6145150‑M2/76‑31 (MQ=255)
                                   aaGACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         >  2:149998‑M2/1‑46 (MQ=255)
                                    aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                        >  1:4046850‑M2/1‑45 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccaataagcgctaacttaagggttgtggtatt                                       <  2:2630375‑M2/76‑33 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCATTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                       >  1:5710860‑M2/1‑44 (MQ=255)
                                      cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                      <  2:2907706‑M2/39‑1 (MQ=255)
                                      cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                      <  2:773965‑M2/39‑1 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                      >  1:5467123‑M2/1‑43 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                     >  1:1443540‑M2/1‑42 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggcattacg                                    <  2:4929342‑M2/76‑36 (MQ=255)
                                         cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                   >  1:2265443‑M2/35‑76 (MQ=255)
                                          ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                  <  1:943900‑M2/76‑38 (MQ=255)
                                          agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                  >  2:3617968‑M2/34‑76 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                 >  1:246575‑M2/1‑38 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                 >  2:2039111‑M2/1‑38 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                 <  2:2716686‑M2/76‑39 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                >  2:6360791‑M2/32‑76 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                >  1:157090‑M2/32‑76 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                <  1:3118047‑M2/45‑1 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                >  2:3196802‑M2/32‑76 (MQ=255)
                                             tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                               >  2:1179197‑M2/31‑76 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            >  2:5682452‑M2/1‑33 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                          <  1:4403498‑M2/51‑1 (MQ=255)
                                                  cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                          >  1:5793330‑M2/1‑31 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                         <  1:4238521‑M2/52‑1 (MQ=255)
                                                     cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                       >  1:1113209‑M2/23‑76 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      >  1:3850757‑M2/1‑27 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                     >  1:3300842‑M2/1‑26 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                     >  2:1411324‑M2/1‑26 (MQ=255)
                                                        ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                    >  2:3460596‑M2/20‑76 (MQ=255)
                                                          aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                  >  1:3560137‑M2/18‑76 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                 >  2:1890347‑M2/1‑22 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                 <  1:3054074‑M2/76‑55 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                 <  1:4525306‑M2/76‑55 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 <  2:3132980‑M2/60‑1 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 <  1:4812794‑M2/60‑1 (MQ=255)
                                                            aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                <  1:1014511‑M2/76‑56 (MQ=255)
                                                              tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc              >  1:2757612‑M2/14‑76 (MQ=255)
                                                               ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga             >  1:6606856‑M2/1‑18 (MQ=255)
                                                               agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa             >  2:5561729‑M2/13‑76 (MQ=255)
                                                                 gCAGCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg           >  2:2222466‑M2/1‑16 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          >  2:5234441‑M2/1‑15 (MQ=255)
                                                                   aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa         <  2:1344431‑M2/76‑63 (MQ=255)
                                                                   aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa         <  2:3707663‑M2/76‑63 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc        <  2:4251407‑M2/69‑1 (MQ=255)
                                                                    aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat        >  1:3954686‑M2/1‑13 (MQ=255)
                                                                    aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat        >  2:1816937‑M2/1‑13 (MQ=255)
                                                                    aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat        >  1:6380943‑M2/1‑13 (MQ=255)
                                                                     tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCCGGTTCa       >  2:2175139‑M2/7‑76 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt       >  2:2430925‑M2/1‑12 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt       <  2:4064645‑M2/76‑65 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGGCCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     >  2:4582372‑M2/5‑76 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     <  1:5119988‑M2/72‑1 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc   <  1:6057253‑M2/76‑69 (MQ=255)
                                                                         ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa   <  2:6501024‑M2/74‑1 (MQ=255)
                                                                          gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct  >  1:359583‑M2/1‑7 (MQ=255)
                                                                           |                                                                          
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596

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