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)43 (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)34 (0.870) 29/144 0.4 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
                                                                          |                                                                           
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  2:896805‑M2/75‑76 (MQ=255)
aTTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGcc                                                                            >  1:384715/1‑76 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        >  2:385989/1‑76 (MQ=255)
        gtcttaacgaggcaccgaggcgtctcattcttcagatgggtcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    <  1:123952‑M2/10‑1 (MQ=255)
        cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    <  1:353318‑M2/76‑5 (MQ=255)
        cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  2:332773‑M2/1‑72 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                   <  1:92147‑M2/76‑6 (MQ=255)
         aGAAACAACCCTGAATAACGGAACCAATAACATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCACcccat                                                                   <  1:453078‑M2/76‑6 (MQ=255)
          tatacaaccctGANTAACGGAACCTATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  1:168141‑M2/4‑70 (MQ=255)
          cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                  >  2:1013308‑M2/65‑76 (MQ=255)
           aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  1:796809‑M2/1‑69 (MQ=255)
            taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                <  2:845024‑M2/14‑1 (MQ=255)
            taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                <  1:520155‑M2/14‑1 (MQ=255)
            taacgaggcaccgaggcgtcacnntcttcagatgnngcaaccctaaagttcgcgcttatgggCCGCAACCCCGGGt                                                                <  1:1174982‑M2/14‑1 (MQ=255)
            aacaacCCTGAATAACGGAACCAATACCAGTCCGTTTCAGGCGCGTTATTTTGCAACNGGGGCCGCAAcccataag                                                                >  1:1237477‑M2/1‑68 (MQ=255)
               cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                             >  1:278946‑M2/60‑76 (MQ=255)
               aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             <  2:105430‑M2/76‑12 (MQ=255)
                 cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCTGAGCCGCAAcccataagcgcta                                                           >  1:2021‑M2/1‑63 (MQ=255)
                  ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCTCAACCCCGGGTgctgct                                                          <  2:446399‑M2/20‑1 (MQ=255)
                    caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        >  2:158009‑M2/55‑76 (MQ=255)
                     gAATAACGGATCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       >  1:233676‑M2/1‑59 (MQ=255)
                     gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       >  1:233657‑M2/1‑59 (MQ=255)
                      ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                      <  1:310973‑M2/24‑1 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      <  1:618506‑M2/76‑19 (MQ=255)
                      aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCCCACcccataagcgctaactta                                                      <  1:1225262‑M2/76‑19 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:370525‑M2/76‑20 (MQ=255)
                         aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                   <  2:622106‑M2/76‑22 (MQ=255)
                           gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  1:260521‑M2/29‑1 (MQ=255)
                           gcgtcgcantattcagatgnttcaacccttaagttagcgcttntgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:1090958‑M2/29‑1 (MQ=255)
                             gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                               >  2:598153‑M2/46‑76 (MQ=255)
                              tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              <  2:608333‑M2/32‑1 (MQ=255)
                               aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                             >  2:16233‑M2/1‑49 (MQ=255)
                                gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                            >  1:211094‑M2/43‑76 (MQ=255)
                                ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                            >  1:451477‑M2/1‑48 (MQ=255)
                                ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                            <  2:867517‑M2/76‑29 (MQ=255)
                                   ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  1:492620‑M2/37‑1 (MQ=255)
                                   ttcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         >  2:773886‑M2/40‑76 (MQ=255)
                                     aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  2:1204510‑M2/1‑43 (MQ=255)
                                      ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                      <  2:871038‑M2/40‑1 (MQ=255)
                                        aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                    <  1:1041561‑M2/76‑37 (MQ=255)
                                           atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 <  2:388654‑M2/45‑1 (MQ=255)
                                           atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 <  1:669199‑M2/45‑1 (MQ=255)
                                            tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                <  1:625553‑M2/46‑1 (MQ=255)
                                               ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             >  2:597738‑M2/28‑76 (MQ=255)
                                                ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            >  2:809590‑M2/1‑32 (MQ=255)
                                                 caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                           >  1:218465‑M2/26‑76 (MQ=255)
                                                    cccttaagntagcgcttatnggCCGCAACCCCGGGTGCTGCTNATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                        <  2:1092743‑M2/54‑1 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  2:537053‑M2/1‑27 (MQ=255)
                                                      cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      <  2:490417‑M2/56‑1 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  1:249722‑M2/76‑55 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGGGCTGCTAATGCGGATGCGACCTTCAAGNTTCAGTATCAATAACCTa                  <  1:1244413‑M2/60‑1 (MQ=255)
                                                           aTTTTGCAACCGGGGCCGCAAaccataagcgctaacttaagggttgtggtattgcgccagatatgatttaacgtgc                 >  2:198380‑M2/1‑21 (MQ=255)
                                                            ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                >  1:1250870‑M2/15‑76 (MQ=255)
                                                             tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg               >  1:273608‑M2/1‑19 (MQ=255)
                                                              ttGCAACCGGGGCNGCAAcccataagcgcttacttaagggttgtgatattacgcctgatatgatttaacgtgccga              >  2:1181275‑M2/1‑18 (MQ=255)
                                                              agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa              <  1:329453‑M2/64‑1 (MQ=255)
                                                                cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       >  1:82397‑M2/11‑65 (MQ=255)
                                                                cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  2:82397‑M2/55‑1 (MQ=255)
                                                                 gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt           >  1:210175‑M2/10‑76 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:421015‑M2/1‑15 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           <  1:228706‑M2/76‑62 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  1:819393‑M2/1‑15 (MQ=255)
                                                                  cttatgggCCGCANCCCCGGGTGCTGCTNATNGGGATGCGACCTTCAACGTTCAGTATCAATAACCTACCCAGGtt          <  2:1169969‑M2/68‑1 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          <  1:198041‑M2/68‑1 (MQ=255)
                                                                  aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          >  2:618506‑M2/1‑14 (MQ=255)
                                                                   ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc         >  1:233402‑M2/8‑76 (MQ=255)
                                                                   aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat         <  1:874067‑M2/76‑64 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtggcgatgaatt        <  1:46438‑M2/76‑65 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:1088026‑M2/1‑12 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:1153473‑M2/1‑12 (MQ=255)
                                                                     cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  1:1077280‑M2/1‑11 (MQ=255)
                                                                     cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:834564‑M2/1‑11 (MQ=255)
                                                                     cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  2:489685‑M2/1‑11 (MQ=255)
                                                                      ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac      >  2:924765‑M2/1‑10 (MQ=255)
                                                                         gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct   >  2:514883‑M2/1‑7 (MQ=255)
                                                                          ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  1:85602‑M2/1‑6 (MQ=255)
                                                                          |                                                                           
ATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524448‑4524597

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