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)40 (1.060) 36/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)37 (0.980) 31/144 0.3 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596
                                                                           |                                                                          
cgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggc                                                                            >  1:538026‑M2/76‑76 (MQ=255)
 gaggcatcgtcttaacgcggcaccgaggcgtcgcattcttcagatggttcaacccctaagttagcgcttatgggcc                                                                           >  1:1253828‑M2/75‑76 (MQ=255)
 gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                           >  2:1149682‑M2/75‑76 (MQ=255)
 gaggcatcgtcttaacgaggcaccgaagcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                           <  1:1205836‑M2/2‑1 (MQ=255)
 aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                           >  1:819837/1‑76 (MQ=255)
    tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                        >  1:675171/1‑76 (MQ=255)
    gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                        >  2:155306‑M2/72‑76 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     <  1:407902‑M2/76‑3 (MQ=255)
       tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                     <  1:134642‑M2/8‑1 (MQ=255)
        tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                    >  1:646269‑M2/1‑73 (MQ=255)
        cgtctcaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                    <  1:970249‑M2/9‑1 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                  <  2:328688‑M2/76‑6 (MQ=255)
           gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                 <  2:572936‑M2/76‑7 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                 <  1:812899‑M2/12‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcagcccttaagntagcgcttatgggCCGCAACCCCGGGt                                                               <  1:1115277‑M2/14‑1 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                               <  2:1048214‑M2/14‑1 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                              <  2:133903‑M2/76‑10 (MQ=255)
               caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                             >  1:527081‑M2/1‑66 (MQ=255)
               caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGGTATTTTGCAACCGGGGCCGCAAcccataagcg                                                             <  1:569992‑M2/76‑11 (MQ=255)
                  cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                          >  2:320381‑M2/1‑63 (MQ=255)
                  aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                          >  2:367094‑M2/58‑76 (MQ=255)
                   ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                         <  1:971007‑M2/20‑1 (MQ=255)
                   ggcaccgaggcgtcgcattcttcacatgnttcaacccttcagttagcgctnatcggCCGCAACCCCGGGTgctgct                                                         >  1:179839‑M2/57‑76 (MQ=255)
                   ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                         >  2:313059‑M2/1‑62 (MQ=255)
                     tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                       <  1:361444‑M2/76‑17 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                      <  1:1024156‑M2/76‑18 (MQ=255)
                        cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                    >  2:910153‑M2/52‑76 (MQ=255)
                         gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                   >  2:899464‑M2/51‑76 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  1:266120‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  >  1:679609‑M2/1‑55 (MQ=255)
                            cGGAACCNCTACCATTCCGTTCCAGGCGCGATATTTTGCAACCGGGGCCGCAAccaacaagagctaacttaagaga                                                <  2:1019068‑M2/76‑24 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                             <  2:977034‑M2/76‑27 (MQ=255)
                                  cAATACCATTCCGTTCTAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgatagctgacgggttgtcgt                                          <  1:14910‑M2/76‑30 (MQ=255)
                                  cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                          >  1:216394‑M2/1‑47 (MQ=255)
                                    ttcatcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                        <  2:128518‑M2/37‑1 (MQ=255)
                                     tcttcagatggttcaacccttaagttagggcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                       <  2:774568‑M2/38‑1 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                       >  1:681607‑M2/1‑44 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                      >  2:1139721‑M2/1‑43 (MQ=255)
                                        tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATACGACCTTCAAg                                    <  1:37601‑M2/41‑1 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcacataagcgctaacttaagggttgtggtattacg                                    <  2:70739‑M2/76‑36 (MQ=255)
                                            ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                <  1:526315‑M2/76‑40 (MQ=255)
                                             tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                               >  2:276090‑M2/31‑76 (MQ=255)
                                              gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                              <  1:507235‑M2/76‑42 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccanaagcgctaacttaagggttgtggtattacgcctgata                             >  1:1234362‑M2/1‑34 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            >  2:942913‑M2/1‑33 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                          <  1:630013‑M2/51‑1 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                          <  1:578528‑M2/51‑1 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                          >  1:319586‑M2/26‑76 (MQ=255)
                                                  cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                          <  1:640030‑M2/76‑46 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                         <  1:785173‑M2/52‑1 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                         <  2:505940‑M2/52‑1 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                         <  2:449224‑M2/52‑1 (MQ=255)
                                                   aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                         <  1:259417‑M2/76‑47 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      >  2:711837‑M2/1‑27 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                     <  2:1053630‑M2/76‑51 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                     >  2:758349‑M2/21‑76 (MQ=255)
                                                        cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                    <  2:367376‑M2/76‑52 (MQ=255)
                                                         gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                   >  2:1078822‑M2/1‑24 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 >  2:326378‑M2/17‑76 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 <  2:925535‑M2/60‑1 (MQ=255)
                                                            gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                <  1:382570‑M2/61‑1 (MQ=255)
                                                            gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                >  1:975487‑M2/16‑76 (MQ=255)
                                                            aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                <  1:571651‑M2/76‑56 (MQ=255)
                                                               agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa             >  2:423139‑M2/13‑76 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat            <  2:35264‑M2/76‑60 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggaattacgcctgatatgatttaacgtgccgatg           <  2:937149‑M2/76‑61 (MQ=255)
                                                                 cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg           >  2:863074‑M2/11‑76 (MQ=255)
                                                                  cAGCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          >  2:477763‑M2/1‑15 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          <  1:718868‑M2/76‑62 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga          <  1:670875‑M2/76‑62 (MQ=255)
                                                                   cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt         >  1:905255‑M2/9‑76 (MQ=255)
                                                                    ttatgggCCGCGACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc        >  1:1164946‑M2/8‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc        >  2:837434‑M2/8‑76 (MQ=255)
                                                                      cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaangtgccgatgaatta      >  1:1731‑M2/1‑11 (MQ=255)
                                                                      atgggccGTAAcccataagcgctaacttaggggttgtggtattacgcctgatatgatttaacgtgccgatgaatta      >  1:41966‑M2/3‑11 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTNAgg     <  1:1141404‑M2/72‑1 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     >  2:1169406‑M2/5‑76 (MQ=255)
                                                                          gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct  >  1:114591‑M2/1‑7 (MQ=255)
                                                                           |                                                                          
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596

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