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)25 (0.990) 21/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)25 (0.990) 23/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)

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGG  >  REL606/4524447‑4524593
                                                                           |                                                                       
cgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaaccctcaagttagcgcttatgggc                                                                         <  1:294126‑M2/1‑1 (MQ=255)
  tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                       >  1:91534/1‑76 (MQ=255)
     catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                    <  2:247255‑M2/6‑1 (MQ=255)
      atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                   <  1:254731‑M2/7‑1 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                              >  2:559767‑M2/65‑76 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                            <  2:425685‑M2/76‑9 (MQ=255)
              aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                           >  2:630027‑M2/62‑76 (MQ=255)
               caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                          <  1:682879‑M2/76‑11 (MQ=255)
                    gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                     <  2:430299‑M2/21‑1 (MQ=255)
                     caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                    >  2:548091‑M2/55‑76 (MQ=255)
                      gAATAACGGAACCAATACCATGACGTTNCANGCNCGATATTNTGCAACAGGGGCCGCAAcccataagcgctaactt                                                   >  1:734229‑M2/1‑59 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                   >  1:221399‑M2/54‑76 (MQ=255)
                      accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                   >  2:412850‑M2/54‑76 (MQ=255)
                       ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                  >  1:471000‑M2/53‑76 (MQ=255)
                       ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                  <  1:462170‑M2/24‑1 (MQ=255)
                        aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                 >  1:369936‑M2/1‑57 (MQ=255)
                            cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                             >  2:580995‑M2/1‑53 (MQ=255)
                             cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                            <  2:226910‑M2/30‑1 (MQ=255)
                              gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                           <  2:808222‑M2/76‑26 (MQ=255)
                                 gcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                        >  1:496104‑M2/43‑76 (MQ=255)
                                 ccAATCCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgg                                        >  2:648923‑M2/1‑48 (MQ=255)
                                  cattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc                                       >  1:250273‑M2/42‑76 (MQ=255)
                                  cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                       <  2:759517‑M2/76‑30 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                      <  1:463257‑M2/76‑31 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                      >  1:184464‑M2/1‑46 (MQ=255)
                                    aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                     <  2:97663‑M2/76‑32 (MQ=255)
                                    aTACCATTCCGTTCCAGGCGCGTTATTNTGCAACCGAGGACGCAAcccataagcgctaacttaagggttgtggtat                                     <  1:740658‑M2/76‑32 (MQ=255)
                                    aTACCATTCCGTGCCCGGCGCGTTATTTTGCAAGCGGGGCCGCAAcccataagcgctaactnaagggttgtggtat                                     <  1:719889‑M2/76‑32 (MQ=255)
                                     tcttcagttgattcaacctttaagttagcgcttgtgggCNNCAACCCCNGGTGCTGCTAATGCGGATGCGACCTTc                                    <  1:828007‑M2/38‑1 (MQ=255)
                                     tACCATTCCGTTCCATGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttcagggtcgtggtctt                                    <  2:70814‑M2/76‑33 (MQ=255)
                                       ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                  >  1:232874‑M2/37‑76 (MQ=255)
                                         canatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGATCTTCTNgg                                <  2:691774‑M2/42‑1 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                             >  1:466634‑M2/32‑76 (MQ=255)
                                             cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctga                            <  1:800441‑M2/76‑41 (MQ=255)
                                                  cAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatga                       <  2:114346‑M2/76‑46 (MQ=255)
                                                   aacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCa                      >  2:708100‑M2/25‑76 (MQ=255)
                                                     gcgcGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattt                    >  1:537776‑M2/1‑28 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                  >  2:447562‑M2/1‑26 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                  <  1:532888‑M2/56‑1 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                  >  2:635234‑M2/21‑76 (MQ=255)
                                                       cttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATaa                  <  2:76857‑M2/56‑1 (MQ=255)
                                                         gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                >  1:623590‑M2/1‑24 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataaacgctaacttaagggtngaggtattacgcctgntatgatctaacgcgccgatg        <  2:833321‑M2/76‑61 (MQ=255)
                                                                  cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga       <  1:379718‑M2/76‑62 (MQ=255)
                                                                   cttatgcgTCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt      >  1:782781‑M2/9‑76 (MQ=255)
                                                                     tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa    <  2:479626‑M2/70‑1 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatt    >  2:247828‑M2/1‑12 (MQ=255)
                                                                     ccGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatataatttaacgtgccgatgaatt    >  2:516747‑M2/1‑12 (MQ=255)
                                                                      atgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg   >  2:46045‑M2/6‑76 (MQ=255)
                                                                       ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac  >  1:262398‑M2/1‑10 (MQ=255)
                                                                           |                                                                       
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGG  >  REL606/4524447‑4524593

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.