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)34 (0.760) 32/144 0.4 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)45 (1.010) 36/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
                                                                          |                                                                           
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            <  2:389241‑M2/2‑1 (MQ=255)
gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                            >  1:1366433‑M2/75‑76 (MQ=255)
 tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           >  1:1129643/1‑76 (MQ=255)
 tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           <  1:291006/76‑1 (MQ=255)
  ggcatcgtcttgacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGc                                                                          <  2:1132261‑M2/4‑1 (MQ=255)
  ggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGc                                                                          <  1:683877‑M2/4‑1 (MQ=255)
  ggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGc                                                                          <  1:272895‑M2/4‑1 (MQ=255)
   gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         <  1:465373‑M2/5‑1 (MQ=255)
   gcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCa                                                                         <  1:282810‑M2/5‑1 (MQ=255)
   cgcttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCa                                                                         >  2:19833/4‑76 (MQ=255)
    aGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                                                        <  2:1397192/76‑1 (MQ=255)
     atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                       >  1:720870‑M2/70‑76 (MQ=255)
      ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                      >  2:224009‑M2/1‑74 (MQ=255)
      tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                      <  1:817469‑M2/8‑1 (MQ=255)
       tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                     <  2:333363‑M2/76‑4 (MQ=255)
       cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                     >  1:877305‑M2/68‑76 (MQ=255)
          gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  >  1:1201247‑M2/1‑70 (MQ=255)
           aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                 >  2:805202‑M2/1‑69 (MQ=255)
            aacaacCCTGAATAACGGTACCAATTCCATTCCGTACCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                >  2:118037‑M2/1‑68 (MQ=255)
             aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                               >  1:16874‑M2/62‑76 (MQ=255)
             aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcctatgggCCGCAACCCCGGGtg                                                               >  2:1198812‑M2/62‑76 (MQ=255)
                gaggcaccgaggcttcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            <  2:276175‑M2/18‑1 (MQ=255)
                gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                            >  1:1340296‑M2/59‑76 (MQ=255)
                 cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           >  2:807201‑M2/1‑63 (MQ=255)
                  ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                          >  1:40956‑M2/57‑76 (MQ=255)
                  ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                          >  1:1159557‑M2/1‑62 (MQ=255)
                   cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                         <  1:1034132‑M2/76‑16 (MQ=255)
                    caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                        <  2:147392‑M2/22‑1 (MQ=255)
                     gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                       <  1:626982‑M2/76‑18 (MQ=255)
                     accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                       <  2:340243‑M2/23‑1 (MQ=255)
                       cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     >  1:765144‑M2/52‑76 (MQ=255)
                       cgaggcgtcgcattcgtcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     <  2:1061051‑M2/25‑1 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     >  1:529813‑M2/1‑57 (MQ=255)
                       aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  1:736530‑M2/76‑20 (MQ=255)
                         aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                   <  2:256569‑M2/27‑1 (MQ=255)
                          ggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa                                                  <  2:600544‑M2/28‑1 (MQ=255)
                           cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggg                                                 <  1:617157‑M2/76‑24 (MQ=255)
                            cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                >  2:1081149‑M2/47‑76 (MQ=255)
                              tcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              >  1:300955‑M2/45‑76 (MQ=255)
                              tcgcattcttcagatggctcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCGGCTAATGCGGATGCg                                              >  2:1061605‑M2/45‑76 (MQ=255)
                               cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              >  2:472291‑M2/44‑75 (MQ=255)
                               cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg                                              <  1:472291‑M2/32‑1 (MQ=255)
                                gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                            >  2:282869‑M2/43‑76 (MQ=255)
                                  aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                          <  1:1256615‑M2/76‑31 (MQ=255)
                                   tacttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt                                         <  2:1280416‑M2/37‑1 (MQ=255)
                                   aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                         <  1:1186616‑M2/76‑32 (MQ=255)
                                    tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                        <  2:1072921‑M2/38‑1 (MQ=255)
                                    tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  2:269813‑M2/1‑44 (MQ=255)
                                        cagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                    <  1:454040‑M2/42‑1 (MQ=255)
                                           ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 <  1:768010‑M2/76‑40 (MQ=255)
                                           atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                 >  2:815477‑M2/32‑76 (MQ=255)
                                             ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                               >  2:1424940‑M2/30‑76 (MQ=255)
                                              gttgaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAATGTTCAGt                              <  1:955589‑M2/48‑1 (MQ=255)
                                                tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            >  1:963529‑M2/27‑76 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        <  1:415207‑M2/76‑49 (MQ=255)
                                                    gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                        >  2:893199‑M2/1‑28 (MQ=255)
                                                    cccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                        >  2:309457‑M2/23‑76 (MQ=255)
                                                     cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                       >  1:1195341‑M2/1‑27 (MQ=255)
                                                      gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                      >  2:1466638‑M2/1‑26 (MQ=255)
                                                          tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  >  1:1264157‑M2/1‑22 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:901457‑M2/60‑1 (MQ=255)
                                                          agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:1474251‑M2/60‑1 (MQ=255)
                                                           gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTGTCAATAACCTAc                 <  1:1471739‑M2/61‑1 (MQ=255)
                                                           gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 <  2:845110‑M2/61‑1 (MQ=255)
                                                           gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                 <  2:502610‑M2/61‑1 (MQ=255)
                                                           aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                 <  1:328588‑M2/76‑56 (MQ=255)
                                                              ttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccga              >  2:1148667‑M2/1‑18 (MQ=255)
                                                              agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa              >  2:956778‑M2/13‑76 (MQ=255)
                                                               tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat             <  1:233202‑M2/76‑60 (MQ=255)
                                                               gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             >  2:222911‑M2/12‑76 (MQ=255)
                                                                gCAACCGGGGCCGCAAcccataagcgctaacttaagggctgtggtattacgcctgatatgatttaacgtgccgatg            <  2:766381‑M2/76‑61 (MQ=255)
                                                                 cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:729981‑M2/1‑15 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  1:627087‑M2/9‑76 (MQ=255)
                                                                  cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  1:1195211‑M2/9‑76 (MQ=255)
                                                                  aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa          <  2:822980‑M2/76‑63 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  2:561780‑M2/1‑12 (MQ=255)
                                                                    ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt        >  1:496781‑M2/1‑12 (MQ=255)
                                                                      tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg      <  2:108618‑M2/72‑1 (MQ=255)
                                                                       gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact     <  1:490570‑M2/76‑68 (MQ=255)
                                                                       gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg     <  1:1350403‑M2/73‑1 (MQ=255)
                                                                       cggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg     >  1:790731‑M2/4‑76 (MQ=255)
                                                                        ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    >  1:1263421‑M2/1‑8 (MQ=255)
                                                                        ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa    <  1:622315‑M2/74‑1 (MQ=255)
                                                                         gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc   <  2:519469‑M2/75‑1 (MQ=255)
                                                                          ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg  <  2:367659/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 < 28 ≤ ATCG/ATCG < 33 ≤ 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.