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

ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGG  >  REL606/4524448‑4524594
                                                                          |                                                                        
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                         >  2:543813‑M2/75‑76 (MQ=255)
 aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                        >  2:1190474‑M2/74‑76 (MQ=255)
 aggcatcgtcttaacgaggcacagaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                        >  1:30974‑M2/74‑76 (MQ=255)
   gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                      <  1:1303481‑M2/5‑1 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                     >  2:1280826/1‑76 (MQ=255)
      ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                   <  2:1447619‑M2/76‑3 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                <  2:639235‑M2/76‑6 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                <  2:479770‑M2/76‑6 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                               >  2:396802‑M2/1‑70 (MQ=255)
          cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                               >  1:1300577‑M2/65‑76 (MQ=255)
           aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                              <  2:1363317‑M2/76‑8 (MQ=255)
            aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                             <  1:1410669‑M2/76‑9 (MQ=255)
             acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                            <  1:1069945‑M2/76‑10 (MQ=255)
              caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                           <  1:1336237‑M2/76‑11 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                          <  1:1086684‑M2/76‑12 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                          <  1:999504‑M2/76‑12 (MQ=255)
                gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                         <  1:898889‑M2/18‑1 (MQ=255)
                 cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                        >  1:790309‑M2/1‑63 (MQ=255)
                 aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                        >  2:1181‑M2/58‑76 (MQ=255)
                   gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                      <  1:15592‑M2/21‑1 (MQ=255)
                   cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                      >  2:290812‑M2/1‑61 (MQ=255)
                    tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                     <  2:1326881‑M2/76‑17 (MQ=255)
                    tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                     <  1:394451‑M2/76‑17 (MQ=255)
                     gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                    <  1:1132668‑M2/76‑18 (MQ=255)
                     gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                    <  1:864387‑M2/76‑18 (MQ=255)
                     accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                    >  1:1092128‑M2/54‑76 (MQ=255)
                      ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                   <  1:940479‑M2/24‑1 (MQ=255)
                         aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                <  1:1561505‑M2/27‑1 (MQ=255)
                         aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                <  1:374491‑M2/27‑1 (MQ=255)
                          ggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa                                               >  2:780595‑M2/49‑76 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                              <  1:1247957‑M2/29‑1 (MQ=255)
                            cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                             <  2:1047863‑M2/30‑1 (MQ=255)
                             gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                            <  1:640840‑M2/31‑1 (MQ=255)
                              tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                           <  1:95275‑M2/32‑1 (MQ=255)
                                gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                         <  2:1379515‑M2/34‑1 (MQ=255)
                                 cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                        >  1:1217000‑M2/42‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                       >  1:1141859‑M2/41‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                       >  1:1141865‑M2/41‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                       >  1:64429‑M2/41‑76 (MQ=255)
                                  aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                       <  1:1457310‑M2/76‑31 (MQ=255)
                                   ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                      >  2:273787‑M2/40‑76 (MQ=255)
                                    tACCATTCCGTTCCAGGCGCGTTATTTTACAACCGGGGCCGCAAcccataagcgctaactttaggggtgtggcatt                                     <  1:457697‑M2/76‑33 (MQ=255)
                                     cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                    <  2:1530532‑M2/39‑1 (MQ=255)
                                       tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                  >  1:1286782‑M2/36‑76 (MQ=255)
                                       cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                  >  2:1328905‑M2/1‑41 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                 >  2:404788‑M2/35‑76 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                 >  1:713766‑M2/35‑76 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                 >  1:1087544‑M2/35‑76 (MQ=255)
                                         ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                <  1:1409418‑M2/76‑38 (MQ=255)
                                         agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                >  2:132411‑M2/34‑76 (MQ=255)
                                         agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                >  2:354806‑M2/34‑76 (MQ=255)
                                          tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                               <  1:1325805‑M2/76‑39 (MQ=255)
                                          gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                               <  2:1307558‑M2/44‑1 (MQ=255)
                                            tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                             >  1:616223‑M2/31‑76 (MQ=255)
                                            cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                             <  1:1148045‑M2/76‑41 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                           >  1:117067‑M2/29‑76 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                           <  1:749656‑M2/48‑1 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                          <  1:580909‑M2/49‑1 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                          >  2:228124‑M2/28‑76 (MQ=255)
                                               tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                          >  2:746134‑M2/1‑33 (MQ=255)
                                               tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                          >  2:859060‑M2/1‑33 (MQ=255)
                                                tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                         <  2:67860‑M2/50‑1 (MQ=255)
                                                 caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                        >  1:1360066‑M2/26‑76 (MQ=255)
                                                  aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                       >  2:1412279‑M2/1‑30 (MQ=255)
                                                   acccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                      <  2:1484286‑M2/53‑1 (MQ=255)
                                                   acccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                      >  1:404735‑M2/24‑76 (MQ=255)
                                                    cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                     >  1:168077‑M2/23‑76 (MQ=255)
                                                    cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                     >  1:168090‑M2/23‑76 (MQ=255)
                                                       ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                  >  2:611637‑M2/20‑76 (MQ=255)
                                                       ctttatTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                  <  2:1506091‑M2/74‑52 (MQ=255)
                                                         ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                <  1:791464‑M2/76‑54 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa               >  2:925882‑M2/17‑76 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa               <  2:1080217‑M2/60‑1 (MQ=255)
                                                             tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg            <  1:668802‑M2/76‑58 (MQ=255)
                                                             tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg            >  2:668802‑M2/1‑19 (MQ=255)
                                                              ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga           >  2:578009‑M2/1‑18 (MQ=255)
                                                              ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga           >  2:830941‑M2/1‑18 (MQ=255)
                                                               gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg          <  1:355286‑M2/65‑1 (MQ=255)
                                                                gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg         >  2:1168500‑M2/1‑16 (MQ=255)
                                                                gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg         <  1:423895‑M2/76‑61 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt       <  2:1569413‑M2/68‑1 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt       <  2:10450‑M2/68‑1 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc      >  1:268809‑M2/8‑76 (MQ=255)
                                                                    tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa     <  1:1521471‑M2/70‑1 (MQ=255)
                                                                     cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta    >  2:525312‑M2/1‑11 (MQ=255)
                                                                     atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg    <  1:321066‑M2/71‑1 (MQ=255)
                                                                      tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg   <  1:1077496‑M2/72‑1 (MQ=255)
                                                                      tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg   >  1:469554‑M2/5‑76 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact  <  2:701060‑M2/76‑68 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact  >  2:806761‑M2/1‑9 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact  <  2:1358937‑M2/76‑68 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcccgatatgatttaacgtgccgatgaattact  <  2:679663‑M2/76‑68 (MQ=255)
                                                                          |                                                                        
ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGG  >  REL606/4524448‑4524594

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