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)46 (0.900) 38/144 0.3 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)33 (0.640) 29/144 0.9 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597
                                                                         |                                                                           
tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                           <  1:997040/76‑1 (MQ=255)
 ttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGc                                                                          <  1:303050/76‑1 (MQ=255)
    gTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                                                       >  2:313864‑M2/1‑75 (MQ=255)
     tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                      >  1:901768‑M2/69‑76 (MQ=255)
       gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                    <  1:1529185‑M2/10‑1 (MQ=255)
       cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    <  1:96119‑M2/76‑5 (MQ=255)
       cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                    >  2:196077‑M2/1‑72 (MQ=255)
        tcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCg                                                                   <  1:1230343‑M2/11‑1 (MQ=255)
         gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  1:511535‑M2/76‑7 (MQ=255)
         gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                  <  2:891440‑M2/76‑7 (MQ=255)
           taacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGt                                                                <  1:304235‑M2/14‑1 (MQ=255)
           aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                >  1:1549856‑M2/1‑68 (MQ=255)
              cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                             >  1:1414841‑M2/60‑76 (MQ=255)
              aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                             <  1:1009358‑M2/76‑12 (MQ=255)
               acCCGGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                            >  2:367977‑M2/1‑64 (MQ=255)
                cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                           >  1:1295884‑M2/1‑63 (MQ=255)
                 ggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTgctgct                                                          <  1:1653344‑M2/20‑1 (MQ=255)
                  gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                         <  2:616698‑M2/21‑1 (MQ=255)
                  cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                         <  2:1042611‑M2/76‑16 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  2:1146159‑M2/1‑58 (MQ=255)
                     aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                      >  2:687455‑M2/1‑58 (MQ=255)
                      cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                     >  2:1283256‑M2/52‑76 (MQ=255)
                      aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                     <  2:650177‑M2/76‑20 (MQ=255)
                       gaggcgtcgcattcttcagatggttcaacccttaagtagacgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                    <  1:267867‑M2/26‑1 (MQ=255)
                        aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcacataagcgctaacttaagg                                                   <  2:1464058‑M2/76‑22 (MQ=255)
                         aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                  <  1:1482422‑M2/76‑23 (MQ=255)
                          gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                 <  2:636825‑M2/29‑1 (MQ=255)
                           ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                <  1:1532875‑M2/76‑25 (MQ=255)
                           ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                >  2:427213‑M2/1‑52 (MQ=255)
                            gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                               <  1:101852‑M2/31‑1 (MQ=255)
                              cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGcgg                                             >  2:1468387‑M2/44‑75 (MQ=255)
                              aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                             >  2:1658856‑M2/1‑49 (MQ=255)
                                   tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                        >  1:333108‑M2/1‑44 (MQ=255)
                                    aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                       >  2:728382‑M2/1‑43 (MQ=255)
                                     ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                      >  2:1192196‑M2/1‑42 (MQ=255)
                                        ttCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcc                                   >  1:369807‑M2/1‑39 (MQ=255)
                                         gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                  <  1:286170‑M2/44‑1 (MQ=255)
                                          ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                 <  1:223648‑M2/76‑40 (MQ=255)
                                           cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                                <  1:1196695‑M2/76‑41 (MQ=255)
                                             ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              <  1:1297053‑M2/76‑43 (MQ=255)
                                             ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  2:101807‑M2/1‑34 (MQ=255)
                                             ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  2:1044528‑M2/1‑34 (MQ=255)
                                             ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                              >  2:176209‑M2/1‑34 (MQ=255)
                                             gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                              >  2:698408‑M2/29‑76 (MQ=255)
                                              ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                             <  2:927256‑M2/49‑1 (MQ=255)
                                               tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                            <  2:874372‑M2/50‑1 (MQ=255)
                                               ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                            <  1:1480980‑M2/76‑45 (MQ=255)
                                                  ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatt                         <  2:1532364‑M2/76‑48 (MQ=255)
                                                    ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                       <  2:311742‑M2/55‑1 (MQ=255)
                                                     cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                      <  2:140911‑M2/56‑1 (MQ=255)
                                                      ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                     <  1:1123877‑M2/57‑1 (MQ=255)
                                                      cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                     >  1:825833‑M2/1‑25 (MQ=255)
                                                        aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAaccn                   <  1:508661‑M2/59‑2 (MQ=255)
                                                         tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  >  2:356819‑M2/1‑22 (MQ=255)
                                                         tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  >  2:498558‑M2/1‑22 (MQ=255)
                                                         tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                  <  1:51092‑M2/76‑55 (MQ=255)
                                                         agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:350931‑M2/60‑1 (MQ=255)
                                                         agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                  <  1:550243‑M2/60‑1 (MQ=255)
                                                          aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                 >  2:734396‑M2/1‑21 (MQ=255)
                                                           ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                >  2:978498‑M2/15‑76 (MQ=255)
                                                           ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                <  1:6780‑M2/62‑1 (MQ=255)
                                                              gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  1:1588658‑M2/65‑1 (MQ=255)
                                                              gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg             <  2:446780‑M2/65‑1 (MQ=255)
                                                               cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCTAgg            <  1:388350‑M2/66‑1 (MQ=255)
                                                               cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg            >  1:577258‑M2/11‑76 (MQ=255)
                                                               cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg            <  2:1530474‑M2/66‑1 (MQ=255)
                                                                cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           <  1:1461478‑M2/76‑62 (MQ=255)
                                                                cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           >  2:1598758‑M2/1‑15 (MQ=255)
                                                                cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga           <  2:402699‑M2/76‑62 (MQ=255)
                                                                 cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          <  2:1638381‑M2/68‑1 (MQ=255)
                                                                 cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt          >  2:1667645‑M2/9‑76 (MQ=255)
                                                                   tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa        >  2:1426832‑M2/7‑76 (MQ=255)
                                                                    cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta       >  1:981173‑M2/1‑11 (MQ=255)
                                                                    atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg       <  1:205203‑M2/71‑1 (MQ=255)
                                                                       ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc    <  2:1105619‑M2/76‑69 (MQ=255)
                                                                        gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct   >  2:1439437‑M2/1‑7 (MQ=255)
                                                                        gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAc   >  1:1178725‑M2/2‑76 (MQ=255)
                                                                         ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  2:884821‑M2/1‑6 (MQ=255)
                                                                         ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc  >  2:948195‑M2/1‑6 (MQ=255)
                                                                         |                                                                           
TTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACG  >  REL606/4524449‑4524597

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