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)28 (0.860) 26/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)32 (0.990) 28/144 0.2 96.7% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.030)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524452‑4524601
                                                                           |                                                                          
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                            <  2:375697‑M2/6‑1 (MQ=255)
catcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCaa                                                                            <  1:375042‑M2/6‑1 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacgcttaagttagcgcttatgggCCGCAAc                                                                           >  2:517855‑M2/70‑76 (MQ=255)
 atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                           >  1:53218‑M2/70‑76 (MQ=255)
  ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                          >  2:418067‑M2/1‑74 (MQ=255)
  ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                          <  1:487127‑M2/76‑3 (MQ=255)
   tCAGAAACAACCCTGAATAACGGAACCAACGCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                         <  1:1063841‑M2/76‑4 (MQ=255)
   cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                         <  1:439128‑M2/9‑1 (MQ=255)
    gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                        >  2:69006‑M2/67‑76 (MQ=255)
    cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccca                                                                        >  2:1036353‑M2/1‑72 (MQ=255)
     aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                       <  1:380417‑M2/76‑6 (MQ=255)
      gAAACTACCCTGAATAACTGTACCTATTCCATTCCGGTCGAGGCGCGTTATTTGAAAACCGGTGCCTCAAcccata                                                                      <  1:109140/76‑4 (MQ=255)
      gAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccata                                                                      >  1:476494‑M2/1‑70 (MQ=255)
         aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                   >  1:1030506‑M2/62‑76 (MQ=255)
          acgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgc                                                                  >  1:462978‑M2/61‑76 (MQ=255)
           cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                 >  1:538111‑M2/60‑76 (MQ=255)
           aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgc                                                                 <  1:557234‑M2/76‑12 (MQ=255)
            gaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                <  2:503273‑M2/18‑1 (MQ=255)
            gaggcaccgaggcgtcgcattctgcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgctg                                                                <  2:55014‑M2/18‑1 (MQ=255)
            acCCTGAATAACGGAATCAATACCATTCCGTTCTAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                >  2:727289‑M2/1‑64 (MQ=255)
            acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                                >  1:408673‑M2/1‑64 (MQ=255)
              ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcactaa                                                              >  2:376888‑M2/1‑62 (MQ=255)
              ccTGAATAACGGAACCAATACAATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                              <  1:462595‑M2/76‑15 (MQ=255)
               cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaac                                                             <  2:802700‑M2/76‑16 (MQ=255)
                naccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                            <  1:924464‑M2/22‑1 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaattt                                                           <  1:510603‑M2/76‑18 (MQ=255)
                 gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                           >  2:4775‑M2/1‑59 (MQ=255)
                  aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                          <  2:861715‑M2/76‑19 (MQ=255)
                   aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                         <  2:734588‑M2/76‑20 (MQ=255)
                    gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                        <  1:273724‑M2/26‑1 (MQ=255)
                     aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                       >  1:612679‑M2/1‑55 (MQ=255)
                      aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaaggg                                                      <  1:472731‑M2/76‑23 (MQ=255)
                          aaCCAATACCAGTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                  >  2:593554‑M2/1‑50 (MQ=255)
                                tcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCAGATGCGACCTTc                                            >  1:799552‑M2/39‑76 (MQ=255)
                                tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                            >  2:214741‑M2/1‑44 (MQ=255)
                                 aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                           >  2:1074505‑M2/1‑43 (MQ=255)
                                    aTTCCGTTCCAGGCGCGTTATTTTTCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                        <  2:194052‑M2/76‑37 (MQ=255)
                                       ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                     >  2:316169‑M2/1‑37 (MQ=255)
                                         ggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAg                                   >  2:204732‑M2/30‑76 (MQ=255)
                                         gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgat                                   <  2:64321‑M2/76‑42 (MQ=255)
                                          gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                  >  2:879794‑M2/29‑76 (MQ=255)
                                           tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                 <  1:638995‑M2/76‑44 (MQ=255)
                                            tcaacccttaagttagcncttatgggCCGCAACCCCGGGTGCTGCTAATGCNCATGAGACCATCAAGGTTCAGTAt                                >  2:896285‑M2/27‑76 (MQ=255)
                                            tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                                <  2:614444‑M2/50‑1 (MQ=255)
                                             caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                               <  2:1059667‑M2/51‑1 (MQ=255)
                                             caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                               >  2:5531‑M2/26‑76 (MQ=255)
                                                   ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                         <  2:565492‑M2/57‑1 (MQ=255)
                                                   cgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaac                         >  2:476162‑M2/1‑25 (MQ=255)
                                                    taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                        <  1:315148‑M2/58‑1 (MQ=255)
                                                     aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                       <  2:344267‑M2/59‑1 (MQ=255)
                                                      tATTTTGCAACCGGGGCCGCAAcctataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                      >  2:802499‑M2/1‑22 (MQ=255)
                                                      agttagcgcttatgggCCGCCACCCCGGGTGCTGCTAATGCGGATGCGACCGTCAAGGTTCAGTATCAATAACCTa                      <  1:815178‑M2/60‑1 (MQ=255)
                                                      agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCCAGGTTCAGCATCAATAACCTa                      >  2:567265‑M2/17‑76 (MQ=255)
                                                        ttttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgcc                    <  1:601354‑M2/76‑57 (MQ=255)
                                                        ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGCTTCAGTATCAATAACCTAcc                    >  1:812175‑M2/15‑76 (MQ=255)
                                                        ttagcgcttatgggCCGCAACCCCGGGTGCTGCGAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                    <  1:190987‑M2/62‑1 (MQ=255)
                                                         tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATACGACCTTCAAGGTTCAGTATCAATAACCTAccc                   >  1:861069‑M2/14‑76 (MQ=255)
                                                               ttatgggCCGCAACCCCGGGTACTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGCTTc             >  2:624850‑M2/8‑76 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTtcn            <  1:898105‑M2/70‑2 (MQ=255)
                                                                tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCACTATCAATAACCTACCCAGGTTCa            <  1:834531‑M2/70‑1 (MQ=255)
                                                                tatgggCCGCAACCCAGGGTGCTGCTAATGCGGATGCGACCTTCAAGATTCAGTATCAATAACCTACCCAGGTTCa            >  2:495339‑M2/7‑76 (MQ=255)
                                                                 cGGGGCCGCAAgccagacgcgctaacttcacggttatggtatnangcctgatntgatttaacgtgccgatgaatta           >  1:1001882‑M2/1‑11 (MQ=255)
                                                                      ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg      <  1:673961/76‑1 (MQ=255)
                                                                       cGCAAcccataagcgctaacttaagggctgtggtattacgcctgatatgatttaacgtgccgatgaattactctca     >  1:650043‑M2/1‑5 (MQ=255)
                                                                        gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc    <  1:157429/76‑1 (MQ=255)
                                                                          aacccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacga  >  1:254353‑M2/1‑2 (MQ=255)
                                                                           |                                                                          
AGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCAT  >  REL606/4524452‑4524601

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