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)47 (1.150) 36/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)32 (0.780) 26/144 0.7 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
                                                                           |                                                                           
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                             >  2:1106379/1‑76 (MQ=255)
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                             <  1:840370/76‑1 (MQ=255)
   ttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                                                          >  1:321489/1‑76 (MQ=255)
     catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                        >  1:1261159‑M2/71‑76 (MQ=255)
     aGTTCAGAAACAACCCTTAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        <  1:835271/76‑1 (MQ=255)
      gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                       >  1:1066556‑M2/1‑75 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcggatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       <  2:1253017‑M2/7‑1 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgtttatcggCCGCAAc                                                                       >  2:81512‑M2/70‑76 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       >  2:899950‑M2/70‑76 (MQ=255)
         cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  1:708653‑M2/1‑72 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   >  2:725882‑M2/1‑71 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  >  2:393043‑M2/65‑76 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 <  1:999326‑M2/76‑8 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 <  1:610288‑M2/76‑8 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 <  1:311146‑M2/76‑8 (MQ=255)
             taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                <  2:691654‑M2/14‑1 (MQ=255)
              aacgaggcaccgagncgtcgcattcttcagatggttcaacccttaaggtagcgcttatgggCCGCACCTCNGGGtg                                                               <  2:1340299‑M2/15‑1 (MQ=255)
                aacCCTGTATTACGGAACGAATACCATGCCGTTTCAGGCNNGTTATTTNGCAACCGGGGCCGCAAcccataagcgc                                                             >  1:1318357‑M2/1‑65 (MQ=255)
                  aggcaccgaggcgtcgcattcttcagatggttcaacacttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                           >  2:1065080‑M2/58‑76 (MQ=255)
                   ccTGAATAACGGNACNAATACCATTCCGTTCCAGNCGCGTTTTTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                          >  1:1090262‑M2/1‑62 (MQ=255)
                    gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         <  2:471468‑M2/21‑1 (MQ=255)
                    gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         <  2:865408‑M2/21‑1 (MQ=255)
                      gAATAACGGAACCAATANCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGCCCGCAAcccacaaacgctaactt                                                       <  1:168268‑M2/76‑18 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       <  1:845268‑M2/76‑18 (MQ=255)
                      accgaggcgtcgcattcttcagttggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       >  2:471946‑M2/54‑76 (MQ=255)
                         gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    >  2:422858‑M2/51‑76 (MQ=255)
                           ggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa                                                  <  2:1185714‑M2/28‑1 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:1278554‑M2/29‑1 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:450479‑M2/29‑1 (MQ=255)
                            gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:813726‑M2/29‑1 (MQ=255)
                            cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                 >  1:992612‑M2/1‑53 (MQ=255)
                             ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                >  2:557076‑M2/1‑52 (MQ=255)
                             cgtcgcattcttcaggtggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                <  1:1109026‑M2/30‑1 (MQ=255)
                              gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                               >  2:848926‑M2/1‑51 (MQ=255)
                              gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                               <  2:827528‑M2/76‑26 (MQ=255)
                               aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                              <  1:1004857‑M2/76‑27 (MQ=255)
                                  cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGNATGCGAcc                                           <  1:1289746‑M2/35‑1 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaagttaagggttgtggta                                          <  1:1110625‑M2/76‑31 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          <  1:899929‑M2/76‑31 (MQ=255)
                                    ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  1:157806‑M2/37‑1 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCCAcccataagcgctaacttaagggttgaggtatt                                        <  1:544154‑M2/76‑33 (MQ=255)
                                     tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        <  1:962716‑M2/76‑33 (MQ=255)
                                      cttcagatggttcaacccttaagttagcgcttttgggCCGCAACCCCGGGTGCTGCTAATGCGGGTGCGACCTTCa                                       <  1:79372‑M2/39‑1 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  1:397240‑M2/1‑42 (MQ=255)
                                       ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  1:734973‑M2/1‑42 (MQ=255)
                                        tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                     <  1:562237‑M2/41‑1 (MQ=255)
                                        cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                     >  1:250913‑M2/1‑41 (MQ=255)
                                        cATACCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                     >  2:1037359‑M2/1‑41 (MQ=255)
                                         cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                    <  1:1275587‑M2/42‑1 (MQ=255)
                                          ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                   >  1:659267‑M2/1‑39 (MQ=255)
                                          tgatggttcaaccattaagttagcgcttatgggCCGCATCCCCGGGTGCTGCTAATGCGGGCGCTACCTTCACGGt                                   >  2:17078‑M2/34‑76 (MQ=255)
                                            ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 >  1:555847‑M2/1‑37 (MQ=255)
                                              gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtactacgcctgat                               <  1:1063205‑M2/76‑42 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  2:877880‑M2/1‑34 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            >  2:976692‑M2/1‑32 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccaaaaacgctaactnaagggttgtggtattacgcctgatatg                            >  1:1169743‑M2/1‑32 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           <  2:1197086‑M2/51‑1 (MQ=255)
                                                  caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           <  1:727125‑M2/51‑1 (MQ=255)
                                                  cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                           <  1:698938‑M2/76‑46 (MQ=255)
                                                   aGGCGCGTTATTTTGAAACCGGGGCCGCAAcccctaagcgctaacttaagggttgtggtatnacgcctgatatgat                          <  1:1144186‑M2/76‑47 (MQ=255)
                                                     cccttaagctagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGTGGATGCGACTTTCATGGTTCAGTATCAAt                        >  2:20974‑M2/23‑76 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGCGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       <  1:1136334‑M2/76‑50 (MQ=255)
                                                      ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  2:696308‑M2/55‑1 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:753230‑M2/1‑26 (MQ=255)
                                                        cgTTATTTTGNAGGCGGGGCCNCAAcccataagctctaacttaagggttgtggtattacgcctgatatgatttaac                     <  2:1324086‑M2/76‑52 (MQ=255)
                                                          ttATTTTGCAACGGGGGCCGCAAgccataagcgctaacttaagggttgtggtattacgnctgatatcatttaacgt                   >  1:166423‑M2/1‑23 (MQ=255)
                                                          aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                   <  2:932776‑M2/59‑1 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  >  1:417687‑M2/17‑76 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  2:679142‑M2/60‑1 (MQ=255)
                                                              tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg               >  1:642846‑M2/1‑19 (MQ=255)
                                                                 gCGACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  2:876490‑M2/1‑16 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg            >  2:1269068‑M2/1‑16 (MQ=255)
                                                                   cttatgagCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  2:741443‑M2/9‑76 (MQ=255)
                                                                    aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat         >  2:475787‑M2/1‑13 (MQ=255)
                                                                     tatgggCCGCCACCCCGAGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATATCCTACCCAGGTTCa        >  1:112438‑M2/7‑76 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcngatgaatt        >  1:1283467‑M2/1‑12 (MQ=255)
                                                                       ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac      <  1:726534‑M2/76‑67 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    >  2:924440‑M2/1‑8 (MQ=255)
                                                                           ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  <  1:657191‑M2/76‑71 (MQ=255)
                                                                           ccGCAAcccataagcgctaacttaagggttgtggaatnacgcctgatatgatttaacgtgccgatgaattactctc  <  1:1303164‑M2/76‑71 (MQ=255)
                                                                           |                                                                           
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524447‑4524597

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