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)

GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524453‑4524601
                                                                          |                                                                          
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                                                      >  1:822199‑M2/1‑54 (MQ=255)
                     aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                      >  2:843455‑M2/1‑54 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                     >  2:739678‑M2/48‑76 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                     <  1:545334‑M2/29‑1 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                     >  1:306819‑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:442344‑M2/20‑76 (MQ=255)
                                                  ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                         >  1:509852‑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             >  1:869081‑M2/8‑76 (MQ=255)
                                                              ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc             >  2:721277‑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)
                                                                       gCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcac    >  2:395364‑M2/1‑4 (MQ=255)
                                                                       gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc    >  1:269344/1‑76 (MQ=255)
                                                                        cAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacg   >  1:770420‑M2/1‑3 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt  >  1:607062/1‑76 (MQ=255)
                                                                          |                                                                          
GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524453‑4524601

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