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)39 (1.100) 35/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 (1.050) 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)

ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597
                                                                          |                                                                           
gaggcatcgtcttaacgaggcaccgaggcgttgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:953946‑M2/75‑76 (MQ=255)
    catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                        >  2:662453‑M2/71‑76 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  1:421504/1‑76 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  1:281419/1‑76 (MQ=255)
     atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       <  2:643502‑M2/7‑1 (MQ=255)
       tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                     <  1:810674‑M2/76‑4 (MQ=255)
         tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                   >  2:777053‑M2/66‑76 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  2:298935‑M2/1‑70 (MQ=255)
          cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  >  2:341747‑M2/65‑76 (MQ=255)
           ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                 <  2:128878‑M2/13‑1 (MQ=255)
           aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  2:653079‑M2/1‑69 (MQ=255)
            taacgaggcaccgagtcttcgcattcttcacatggctcaaccntnaagttggngcttatgggCCGCAACCCCGGGt                                                                <  1:1101989‑M2/14‑1 (MQ=255)
            taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgtatatgggCCGCAACCCcgggc                                                                >  2:1006136‑M2/63‑75 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTCGCAACCGGGGCCGCAAcccataagcgc                                                             <  1:714981‑M2/76‑12 (MQ=255)
                acCCTGAATAACGGAACCAATTCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                            >  1:671691‑M2/1‑64 (MQ=255)
                 cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           <  2:267020‑M2/76‑14 (MQ=255)
                    tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                        >  2:1147322‑M2/1‑60 (MQ=255)
                    tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                        <  1:701284‑M2/76‑17 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  2:1128168‑M2/1‑58 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  2:839730‑M2/76‑19 (MQ=255)
                       cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     <  1:42671‑M2/25‑1 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  2:899123‑M2/76‑20 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGACGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:1162384‑M2/76‑20 (MQ=255)
                        tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                    <  1:946099‑M2/76‑21 (MQ=255)
                        gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    >  1:825787‑M2/51‑76 (MQ=255)
                         aaCGGAACCAATACCATTCCGTTCCACGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaaattaagg                                                   <  2:657841‑M2/76‑22 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:88607‑M2/29‑1 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:691648‑M2/29‑1 (MQ=255)
                            ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                >  1:497487‑M2/1‑52 (MQ=255)
                              tcgcattcttcagatggttcaccccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              <  1:740388‑M2/32‑1 (MQ=255)
                              tcgcattcttcagatggttcaacccttaagttagcgcttatgagCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              >  1:947460‑M2/45‑76 (MQ=255)
                               cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                             >  1:545401‑M2/44‑76 (MQ=255)
                                  attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                          >  2:897328‑M2/41‑76 (MQ=255)
                                   ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  1:1142011‑M2/37‑1 (MQ=255)
                                   aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                         >  1:806686‑M2/1‑45 (MQ=255)
                                   aTACCATTCCGTTCCAGGCGCGTTACTTTGCAACNGGGGCCGCAAcccatangggctaacttaagggtngtggtat                                         <  2:1050475‑M2/76‑32 (MQ=255)
                                    tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        <  1:20442‑M2/76‑33 (MQ=255)
                                     ctgcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                       <  1:481373‑M2/39‑1 (MQ=255)
                                     aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataggcgctaacttaagggttgtggtatta                                       >  1:1050476‑M2/1‑43 (MQ=255)
                                      ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                      <  2:160308‑M2/40‑1 (MQ=255)
                                       tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     <  2:377720‑M2/41‑1 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                    >  1:256591‑M2/35‑76 (MQ=255)
                                          tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                  >  1:144945‑M2/1‑38 (MQ=255)
                                           ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccttaagcgctaacttaagggttgtggtgtnacgccgg                                 >  1:1047945‑M2/1‑37 (MQ=255)
                                            tggttcaacccttaagttagcgnttttgggCCGNAACCCCGGNTGCTGCTAATACGGATGCTACCNTCAAGGTTCa                                <  2:1037039‑M2/46‑1 (MQ=255)
                                            tggttcaacccttaagttagcgnttgtgggCCGNAACCCCGGNTGCTGCTAATACGGATGCGACCNTCAAGGTTCa                                <  2:1037035‑M2/46‑1 (MQ=255)
                                            cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                <  2:697945‑M2/76‑41 (MQ=255)
                                             gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                               <  1:219905‑M2/76‑42 (MQ=255)
                                              gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              <  2:681440‑M2/48‑1 (MQ=255)
                                                tcaacccttaagttagcgcttatgngCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            >  2:206198‑M2/27‑76 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  1:842497‑M2/76‑45 (MQ=255)
                                                 caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           >  2:1141060‑M2/26‑76 (MQ=255)
                                                 caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           >  2:565511‑M2/26‑76 (MQ=255)
                                                  aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                          >  1:573255‑M2/25‑76 (MQ=255)
                                                    ncgcatTATTTTGCAACCGGGNCCGGAAnccataggcgctaacttaagggttgtggtattacgcctgatatgattt                        >  1:1007192‑M2/2‑28 (MQ=255)
                                                    cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATcaan                        >  1:1043929‑M2/23‑75 (MQ=255)
                                                     ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       >  1:220916‑M2/22‑76 (MQ=255)
                                                       ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     <  1:800703‑M2/57‑1 (MQ=255)
                                                       ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     >  2:741104‑M2/20‑76 (MQ=255)
                                                       cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     <  1:1080167‑M2/76‑52 (MQ=255)
                                                       cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactcaagggntgtggtnttacgcctgatatgatttaan                     <  1:1000785‑M2/76‑52 (MQ=255)
                                                        taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                    >  1:783738‑M2/19‑76 (MQ=255)
                                                         ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                   >  1:482184‑M2/1‑23 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  1:303793‑M2/76‑55 (MQ=255)
                                                                gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  2:777444‑M2/1‑16 (MQ=255)
                                                                 gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt           >  1:885883‑M2/10‑76 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:440073‑M2/1‑15 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           <  2:238663‑M2/76‑62 (MQ=255)
                                                                  aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          >  1:1083913‑M2/1‑14 (MQ=255)
                                                                   aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat         <  1:701210‑M2/76‑64 (MQ=255)
                                                                    tatgggCCGCAACCCCGGGTGCTGCTAATTCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        <  1:789747‑M2/70‑1 (MQ=255)
                                                                    tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAATGTTCAGTATCAATAACCTACCCAGGTTCa        <  1:1163308‑M2/70‑1 (MQ=255)
                                                                    tatgggACGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        >  1:1116680‑M2/7‑76 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        <  1:744368‑M2/76‑65 (MQ=255)
                                                                    ccGGGACCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  2:939637‑M2/1‑12 (MQ=255)
                                                                      ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac      >  2:21734‑M2/1‑10 (MQ=255)
                                                                        ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    <  2:599352‑M2/74‑1 (MQ=255)
                                                                          ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  <  2:545401/76‑1 (MQ=255)
                                                                          |                                                                           
ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597

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