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)36 (0.990) 31/144 0.2 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)42 (1.160) 32/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)

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524447‑4524597
                                                                           |                                                                           
cgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggc                                                                             >  1:298129‑M2/76‑76 (MQ=255)
  tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           >  1:941553/1‑76 (MQ=255)
  aggcntcgtcttaacgaggcacctaggcatcgcattcttnagatggttcaacccttaagttagcgcttatgggCCg                                                                           <  1:152821‑M2/3‑1 (MQ=255)
    tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                         >  2:610771/1‑76 (MQ=255)
    gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         <  1:865779‑M2/5‑1 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       >  1:353146‑M2/70‑76 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                      >  1:844515‑M2/1‑74 (MQ=255)
         cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  1:243578‑M2/1‑72 (MQ=255)
          tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                   >  2:1052542‑M2/66‑76 (MQ=255)
           gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  1:559482‑M2/76‑7 (MQ=255)
           gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  2:126328‑M2/1‑70 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  >  2:922124‑M2/65‑76 (MQ=255)
            ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                 >  2:120090‑M2/64‑76 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                <  2:635597‑M2/76‑9 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                               >  2:896512‑M2/1‑67 (MQ=255)
              aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               <  1:955307‑M2/15‑1 (MQ=255)
                aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcga                                                             <  1:1156736‑M2/76‑12 (MQ=255)
                aacCCTGAATAACGGAACCAATACCATNCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             <  1:1209059‑M2/76‑12 (MQ=255)
                  cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTNTTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           >  2:1020781‑M2/1‑63 (MQ=255)
                   ggcaccgaggcgtcgcattcttcagatgggtcaacctttaagttagcgcttatgggCNGCAACCCCGGGTgctgct                                                          <  1:1149387‑M2/20‑1 (MQ=255)
                   ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                          <  2:415415‑M2/76‑15 (MQ=255)
                    gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         <  2:987768‑M2/21‑1 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCACAAcccataagcgctaactt                                                       >  1:1827‑M2/1‑59 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  2:249045‑M2/23‑1 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  2:396742‑M2/23‑1 (MQ=255)
                        cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     <  1:628022‑M2/25‑1 (MQ=255)
                           aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                  >  2:843455‑M2/1‑54 (MQ=255)
                           aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                  >  1:822199‑M2/1‑54 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 >  1:306819‑M2/48‑76 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:545334‑M2/29‑1 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 >  2:739678‑M2/48‑76 (MQ=255)
                            cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                 <  2:204219‑M2/76‑24 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              >  2:458152‑M2/1‑50 (MQ=255)
                                cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             <  2:381926‑M2/33‑1 (MQ=255)
                                cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             >  1:1016920‑M2/44‑76 (MQ=255)
                                cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             >  1:249875‑M2/44‑76 (MQ=255)
                                  cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                           <  1:813861‑M2/76‑30 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  2:100949‑M2/41‑76 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  2:13515‑M2/41‑76 (MQ=255)
                                    ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  1:1172798‑M2/37‑1 (MQ=255)
                                    aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                         >  1:958557‑M2/1‑45 (MQ=255)
                                    aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                         <  2:572928‑M2/76‑32 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                     >  2:456248‑M2/1‑41 (MQ=255)
                                          ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                   <  1:687842‑M2/76‑38 (MQ=255)
                                            atggttcaacccttaagatagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCNGATGCGACCTTCAAGGTTc                                 <  1:1121918‑M2/45‑1 (MQ=255)
                                             cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcctataagcgctaacttaagggttgtggtattacgtctga                                >  1:72245‑M2/1‑36 (MQ=255)
                                              ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACNTTCAAGGTTCAg                               <  1:16065‑M2/47‑1 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  2:776168‑M2/76‑43 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  1:1183155‑M2/76‑43 (MQ=255)
                                                ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             <  1:999921‑M2/49‑1 (MQ=255)
                                                 tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTCCTGTAt                            <  1:70502‑M2/50‑1 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  1:623306‑M2/76‑45 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgcggtattacgcctgatatg                            <  1:377542‑M2/76‑45 (MQ=255)
                                                  cAGGCGCGTTATTTTGCACCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                           <  2:454265‑M2/76‑46 (MQ=255)
                                                  cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                           <  1:54409‑M2/76‑46 (MQ=255)
                                                   aacccttaagttagcgcttntgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACTTTCAAGGTTCAGTATCa                          >  1:1124384‑M2/25‑76 (MQ=255)
                                                    acccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                         <  1:813610‑M2/53‑1 (MQ=255)
                                                     gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  1:301848‑M2/1‑28 (MQ=255)
                                                      cgngTTCTTTTCCAACCGGGGCCGCTAcccataagcgctaacttaagggttgtggtataacgcctgatatgtttta                       >  2:1140763‑M2/1‑27 (MQ=255)
                                                       ctttagttagcgcctatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      <  2:717527‑M2/56‑1 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      >  2:268116‑M2/21‑76 (MQ=255)
                                                        ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     >  1:509852‑M2/20‑76 (MQ=255)
                                                        ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     >  1:442344‑M2/20‑76 (MQ=255)
                                                         taagttagagtttatgngCCGCAACCCCGGATNCTGCTAATGCGGATGCGACCTTCAAGGTTCAGNATCAACNAnc                    >  1:90426‑M2/19‑76 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  1:354609‑M2/76‑55 (MQ=255)
                                                             ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                <  2:57348‑M2/62‑1 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  1:949694‑M2/1‑16 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            <  2:689668‑M2/76‑61 (MQ=255)
                                                                  gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGATCTTCAAGGTTCAGTATCAATAACGTACCCAGGt           <  2:1038270‑M2/67‑1 (MQ=255)
                                                                  gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt           <  1:310541‑M2/67‑1 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaaggattgtggtattacgcctgatatgatttaacgngccgatga           <  1:1060046‑M2/76‑62 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  2:721277‑M2/8‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:869081‑M2/8‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:591239‑M2/8‑76 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         <  1:344821‑M2/69‑1 (MQ=255)
                                                                    ttatgggCCGCAACCCCGGGTGCTGCTAATGAGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         <  2:200548‑M2/69‑1 (MQ=255)
                                                                    ttatgggCCACAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         <  1:234000‑M2/69‑1 (MQ=255)
                                                                    aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat         >  2:636189‑M2/1‑13 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  2:325763‑M2/1‑12 (MQ=255)
                                                                          gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc   <  1:286396‑M2/75‑1 (MQ=255)
                                                                           ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  >  2:72312/1‑76 (MQ=255)
                                                                           |                                                                           
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524447‑4524597

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