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)27 (1.300) 26/144 0.0 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)23 (1.110) 20/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)

TCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524455‑4524601
                                                                        |                                                                          
cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                         >  1:498573‑M2/68‑76 (MQ=255)
   gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                      >  2:360047‑M2/1‑70 (MQ=255)
    ttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCggg                                                                     <  1:524399‑M2/13‑1 (MQ=255)
    aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                     >  1:311111‑M2/1‑69 (MQ=255)
      acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                                   >  2:522621‑M2/1‑67 (MQ=255)
       caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcg                                                                  >  1:280345‑M2/1‑66 (MQ=255)
       caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTCGCAACCGGGGCCGCAAcccatcagct                                                                  <  1:433929‑M2/76‑11 (MQ=255)
            gcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTa                                                             <  1:286163‑M2/21‑1 (MQ=255)
            cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                             <  2:26368‑M2/76‑16 (MQ=255)
              gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                           >  1:409759‑M2/1‑59 (MQ=255)
              accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                           <  2:325982‑M2/23‑1 (MQ=255)
              accgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAAt                                                           <  2:227806‑M2/23‑1 (MQ=255)
                cggggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                         <  1:285882‑M2/25‑1 (MQ=255)
                  aggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCgg                                                       >  2:69010‑M2/50‑76 (MQ=255)
                  aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                       <  2:381200‑M2/76‑22 (MQ=255)
                     ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                    <  2:301886‑M2/76‑25 (MQ=255)
                     ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                    >  2:34090‑M2/1‑52 (MQ=255)
                     cgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                    <  2:595941‑M2/30‑1 (MQ=255)
                      gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                   <  1:43438‑M2/31‑1 (MQ=255)
                      gAACCAATACCTTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                   >  1:439230‑M2/1‑51 (MQ=255)
                       aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                  <  1:549532‑M2/76‑27 (MQ=255)
                        cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                                 >  1:654976‑M2/44‑76 (MQ=255)
                          cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                               >  1:142007‑M2/1‑47 (MQ=255)
                          cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggt                                               <  1:545118‑M2/76‑30 (MQ=255)
                           aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                              <  1:107100‑M2/76‑31 (MQ=255)
                            aTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtat                                             >  2:3609‑M2/1‑45 (MQ=255)
                             tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                            >  1:162226‑M2/1‑44 (MQ=255)
                                  agatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGt                                       >  2:610733‑M2/34‑76 (MQ=255)
                                   tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                      >  2:125174‑M2/1‑38 (MQ=255)
                                       ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcatgata                                  >  2:191954‑M2/1‑34 (MQ=255)
                                        tCCAGGCGCGTTATTTTTCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                 >  1:534494‑M2/1‑33 (MQ=255)
                                        tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                 >  1:249692‑M2/1‑33 (MQ=255)
                                         tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                <  1:541502‑M2/50‑1 (MQ=255)
                                         tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCAGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                >  2:109498‑M2/27‑76 (MQ=255)
                                         ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                >  1:482307‑M2/1‑32 (MQ=255)
                                          caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                               <  1:215232‑M2/51‑1 (MQ=255)
                                          caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGAGACCTTCAAGGTTCAGTATc                               >  1:357783‑M2/26‑76 (MQ=255)
                                            ggCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatt                             >  1:382691‑M2/1‑29 (MQ=255)
                                              ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                           <  1:676813‑M2/55‑1 (MQ=255)
                                                    gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                     <  1:588009‑M2/61‑1 (MQ=255)
                                                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa                  <  2:377506‑M2/64‑1 (MQ=255)
                                                       agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCAACCTTCAAGGTTCAGTATCAATAACCTACCCa                  <  2:405724‑M2/64‑1 (MQ=255)
                                                         cgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAgg                <  1:566982‑M2/66‑1 (MQ=255)
                                                          cAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatga               >  2:574757‑M2/1‑15 (MQ=255)
                                                           cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt              >  1:622249‑M2/9‑76 (MQ=255)
                                                           aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa              <  1:501427‑M2/76‑63 (MQ=255)
                                                            aCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaat             <  1:531216‑M2/76‑64 (MQ=255)
                                                               tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg          >  2:444082‑M2/5‑76 (MQ=255)
                                                                   ccGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctc      >  2:310035‑M2/1‑6 (MQ=255)
                                                                   ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg      <  2:96627/76‑1 (MQ=255)
                                                                       aaCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACTCAGGTTCAGGGACGTCAt  >  1:424958/1‑76 (MQ=255)
                                                                        |                                                                          
TCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524455‑4524601

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