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)41 (1.130) 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)39 (1.070) 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                                                                           >  2:737093‑M2/70‑76 (MQ=255)
atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                           >  1:954177‑M2/70‑76 (MQ=255)
 tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                          <  2:624027‑M2/8‑1 (MQ=255)
    tcttaacgaggcaccgaggcgtcgcattcttcagatggtcaaacccttaagttagcgcttatgggCCGCAACCCCg                                                                       >  2:105381‑M2/66‑76 (MQ=255)
    aGAAACAACCCTGAATAATGGCACTAATACCATTCCGTTCCAGGCGCGTCATTTTGCAACCGGGGCCGCAAcccat                                                                       >  1:1146489‑M2/1‑71 (MQ=255)
     gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                      >  1:132728‑M2/1‑70 (MQ=255)
       taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCcgggg                                                                    >  2:266534‑M2/63‑75 (MQ=255)
       taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                    <  1:960370‑M2/14‑1 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                    >  1:700393‑M2/1‑68 (MQ=255)
        aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                   >  1:728785‑M2/62‑76 (MQ=255)
         caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                  <  2:646670‑M2/76‑11 (MQ=255)
          cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                 >  1:659789‑M2/60‑76 (MQ=255)
           gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                <  2:257526‑M2/18‑1 (MQ=255)
            aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                               <  2:548385‑M2/19‑1 (MQ=255)
            aggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctgc                                                               >  1:631197‑M2/58‑76 (MQ=255)
             ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                              <  2:341523‑M2/20‑1 (MQ=255)
             ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                              >  1:647842‑M2/57‑76 (MQ=255)
             ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                              <  2:900231‑M2/20‑1 (MQ=255)
               tGAATAACGGAACCAATCCCATCCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                            <  2:185169‑M2/76‑17 (MQ=255)
               tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                            <  1:847826‑M2/76‑17 (MQ=255)
                 aaTAACGGATCCAATACCATTCCGTTCCAGGCGAGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                          <  2:1042464‑M2/76‑19 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                          <  1:424146‑M2/76‑19 (MQ=255)
                  cgaggcgtcgcatncttcagatggttcaacccttaagttagcgctcatgggCCGCAACCCCGGGTGCTGCTAATGc                                                         <  2:1101448‑M2/25‑1 (MQ=255)
                     ggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCTGCAACCCCGGGTGCTGCTAATGCGGa                                                      <  2:1062904‑M2/28‑1 (MQ=255)
                      cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                     >  2:1067484‑M2/1‑53 (MQ=255)
                      cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                     >  1:1047259‑M2/1‑53 (MQ=255)
                       ggAACCAATATCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                    <  1:1114158‑M2/76‑25 (MQ=255)
                        gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                   <  2:992776‑M2/31‑1 (MQ=255)
                        gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                   <  1:251365‑M2/76‑26 (MQ=255)
                        gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                   <  2:97123‑M2/76‑26 (MQ=255)
                           gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                >  1:408282‑M2/43‑76 (MQ=255)
                             aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                              >  2:226771‑M2/1‑46 (MQ=255)
                              aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                             <  2:989289‑M2/76‑32 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                            >  2:488888‑M2/1‑44 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtata                                            >  2:953085‑M2/1‑44 (MQ=255)
                                cttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTATTGCGGATGCGACCTTCa                                           >  1:135364‑M2/38‑76 (MQ=255)
                                aCCATTCCNTGCCATGCGCNTTATTTTGCAACCGGGGCCGCANcccataagcgctaacttaagggttgtggtatta                                           >  2:1011443‑M2/1‑43 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                           >  2:292911‑M2/1‑43 (MQ=255)
                                  tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                         >  2:646912‑M2/36‑76 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                         >  1:439674‑M2/1‑41 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                         <  1:720997‑M2/76‑36 (MQ=255)
                                     tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                      >  1:447962‑M2/1‑38 (MQ=255)
                                     gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                      <  2:264701‑M2/44‑1 (MQ=255)
                                      atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                     <  2:554904‑M2/45‑1 (MQ=255)
                                      atagttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                     <  2:331205‑M2/45‑1 (MQ=255)
                                       cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                    <  2:1204279‑M2/76‑41 (MQ=255)
                                        ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAGCTTCAAGGTTCAg                                   <  1:77993‑M2/47‑1 (MQ=255)
                                        ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                   <  2:227380‑M2/47‑1 (MQ=255)
                                        gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   >  1:979642‑M2/1‑35 (MQ=255)
                                        gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggggtgtggtattaagcccgat                                   <  1:475807‑M2/76‑42 (MQ=255)
                                         ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                                  >  1:538569‑M2/1‑34 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                 <  2:578204‑M2/76‑44 (MQ=255)
                                            caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                               >  1:1037721‑M2/26‑76 (MQ=255)
                                             aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                              >  1:357228‑M2/25‑76 (MQ=255)
                                               gcgcGTTATTTTGCNACCGGGGCCGCAAtccatangcgctgacntaagggttgtggtattacgcctgatatgattt                            >  2:1003112‑M2/1‑28 (MQ=255)
                                               gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                            >  2:310929‑M2/1‑28 (MQ=255)
                                               gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                            >  2:574614‑M2/1‑28 (MQ=255)
                                                ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                           <  1:731871‑M2/55‑1 (MQ=255)
                                                 cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                          <  2:1164153‑M2/56‑1 (MQ=255)
                                                   taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAataaac                        <  2:367310‑M2/58‑3 (MQ=255)
                                                    aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATNNGGATGCGACCNNCAAGGTTCAGTATCAATAACCt                       >  1:1083932‑M2/18‑76 (MQ=255)
                                                    aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                       <  1:887116‑M2/59‑1 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                      <  1:749324‑M2/60‑1 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                      <  1:429436‑M2/60‑1 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                      <  2:49229‑M2/60‑1 (MQ=255)
                                                      aTTTTGCAACGGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                     <  1:1111195‑M2/76‑56 (MQ=255)
                                                      aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                     <  1:297817‑M2/76‑56 (MQ=255)
                                                        tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaaagtgccg                   <  1:920782‑M2/76‑58 (MQ=255)
                                                        tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattaaacgtgccg                   <  1:36459‑M2/76‑58 (MQ=255)
                                                         agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                  <  1:1006040‑M2/64‑1 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                >  1:784587‑M2/1‑16 (MQ=255)
                                                           gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg                <  1:954751‑M2/76‑61 (MQ=255)
                                                           cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                >  1:773288‑M2/11‑76 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt               >  1:491684‑M2/10‑76 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt               <  1:661691‑M2/67‑1 (MQ=255)
                                                              aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat             <  2:942602‑M2/76‑64 (MQ=255)
                                                               ccGGGGCCGCAAcccataagccctaacttaagggttgtggtattacgcctgatatgatgtaacgtgccgatgaatt            >  2:178577‑M2/1‑12 (MQ=255)
                                                                cGGGGCCGCAAcccataagcgctaacttaagggttgtggtatgacgcctgatatgatttaacgtgctgatgaatta           >  2:1040113‑M2/1‑11 (MQ=255)
                                                                  gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgcattact         <  2:1176132‑M2/76‑68 (MQ=255)
                                                                   ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc        <  1:522460‑M2/76‑69 (MQ=255)
                                                                    gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc       <  2:887823‑M2/75‑1 (MQ=255)
                                                                         aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAt  <  1:176782/76‑1 (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.